Issue
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
Article Number L30
Number of page(s) 7
Section Letters
DOI https://doi.org/10.1051/0004-6361/201014610
Published online 16 July 2010
A&A 518, L30 (2010)

Herschel: the first science highlights

LETTER TO THE EDITOR

Dissecting the cosmic infra-red background with Herschel/PEP[*],[*]

S. Berta1 - B. Magnelli1 - D. Lutz1 - B. Altieri2 - H. Aussel3 - P. Andreani4,5 - O. Bauer1 - A. Bongiovanni6,7 - A. Cava6,7 - J. Cepa6,7 - A. Cimatti8 - E. Daddi3 - H. Dominguez9 - D. Elbaz3 - H. Feuchtgruber1 - N. M. Förster Schreiber1 - R. Genzel1 - C. Gruppioni9 - R. Katterloher1 - G. Magdis3 - R. Maiolino10 - R. Nordon1 - A. M. Pérez García6,7 - A. Poglitsch1 - P. Popesso1 - F. Pozzi8 - L. Riguccini3 - G. Rodighiero11 - A. Saintonge1 - P. Santini10 - M. Sanchez-Portal2 - L. Shao1 - E. Sturm1 - L. J. Tacconi1 - I. Valtchanov2 - M. Wetzstein1 - E. Wieprecht1

1 - Max-Planck-Institut für Extraterrestrische Physik (MPE), Postfach 1312, 85741 Garching, Germany
2 - Herschel Science Centre, ESAC, Villanueva de la Cañada, 28691 Madrid, Spain
3 - Laboratoire AIM, CEA/DSM-CNRS-Université Paris Diderot, IRFU/Service d'Astrophysique, Bât. 709, CEA-Saclay, 91191 Gif-sur-Yvette Cedex, France
4 - ESO, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany
5 - INAF - Osservatorio Astronomico di Trieste, via Tiepolo 11, 34143 Trieste, Italy
6 - Instituto de Astrofísica de Canarias, 38205 La Laguna, Spain
7 - Departamento de Astrofísica, Universidad de La Laguna, Spain
8 - Dipartimento di Astronomia, Università di Bologna, via Ranzani 1, 40127 Bologna, Italy
9 - INAF - Osservatorio Astronomico di Bologna, via Ranzani 1, 40127 Bologna, Italy
10 - INAF - Osservatorio Astronomico di Roma, via di Frascati 33, 00040 Monte Porzio Catone, Italy
11 - Dipartimento di Astronomia, Università di Padova, Vicolo dell'Osservatorio 3, 35122 Padova, Italy

Received 31 March 2010 / Accepted 21 April 2010

Abstract
The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 $\mu $m) by exploiting Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS evolutionary probe (PEP) guaranteed-time survey. The GOODS-N data reach 3$\sigma $ depths of $\sim$3.0 mJy at 100 $\mu $m and $\sim$5.7 mJy at 160 $\mu $m. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 $\mu $m. Differential number counts extend from a few mJy up to 100-200 mJy, and are approximated as a double power law, with the break lying between 5 and 10 mJy. The available ancillary information allows us to split number counts into redshift bins. At $z\le0.5$ we isolate a class of luminous sources ( $L_{\rm IR}\sim10^{11}~L_\odot$), whose SEDs resemble late-spiral galaxies, peaking at $\sim$130 $\mu $m restframe and significantly colder than what is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface brightness of $6.36\pm1.67$ and $6.58\pm1.62$ [nW m-2 sr-1] at 100 and 160 $\mu $m, resolving $\sim$45% and $\sim$52% of the CIB, respectively. When stacking 24 $\mu $m sources, the inferred CIB lies within 1.1$\sigma $ and 0.5$\sigma $ from direct measurements in the two bands, and fractions increase to 50% and 75%. Most of this resolved CIB fraction was radiated at $z\le1.0$, with 160 $\mu $m sources found at higher redshift than 100 $\mu $m ones.

Key words: infrared: diffuse background - infrared: galaxies - cosmic background radiation - galaxies: statistics - galaxies: evolution

1 Introduction

The cosmic IR background (CIB, Puget et al. 1996; Hauser et al. 1998) accounts for roughly half of the total extragalactic background light (EBL, Hauser & Dwek 2001; Lagache et al. 2005), i.e., half of the energy radiated by all galaxies, at all cosmic epochs, at any wavelength (Dole et al. 2006). It is therefore a crucial constraint on modes and times of galaxy formation.

Deep cosmological surveys carried out with the ISO (see Genzel & Cesarsky 2000, for a summary) and Spitzer Space Observatory (Soifer et al. 2008, for a review) produced large samples of mid-IR sources and deep number counts (Papovich et al. 2004; Elbaz et al. 2002). These surveys led to mid-IR CIB lower limits within a factor of two from the upper constraints set by TeV cosmic opacity measurements (e.g. Franceschini et al. 2010).

However, at CIB peak wavelengths (100-200 $\mu $m), the nature of individual galaxies building up the EBL is poorly known. Past surveys produced limited samples of distant far-IR objects (e.g. Frayer et al. 2009), mainly due to the small apertures of the available instruments and the low sensitivity in the far-IR. In the 160 $\mu $m Spitzer/MIPS band, $\sim$7% of the CIB was resolved into individually detected objects (Dole et al. 2004), and it was only through stacking the 24 $\mu $m sources that most (60-70%) of the far-IR CIB could be recovered (Dole et al. 2006; Béthermin et al. 2010). Similarly, at longer wavelengths, only stacking of 24 $\mu $m sources on BLAST maps could account for the majority of the EBL at 250, 350, and 500 $\mu $m (Marsden et al. 2009).

With the favorable diffraction limit of the large Herschel 3.5 m mirror (Pilbratt et al. 2010), and the high sensitivity of its Photodetector Array Camera and Spectrometer (PACS, 70, 100, 160 $\mu $m; Poglitsch et al. 2010), confusion and blending of sources are much less of a limitation. We are now able to resolve a large fraction of the CIB at its peak into individual galaxies.

The PACS evolutionary probe (PEP) extragalactic survey samples 4 different tiers: from the wide and shallow COSMOS field, through medium size areas like the Lockman hole, all the way down to the 160 and 100 $\mu $m confusion limit in the pencil-beam, very deep observations in GOODS-N and GOODS-S, and even beyond by exploiting gravitational lensing in low-redshift galaxy clusters (e.g. Abell 2218; Altieri et al. 2010).

\begin{figure}
\par\includegraphics[width=7.6cm,clip]{14610f1a.eps}\hspace*{4mm}...
...phics[width=7.6cm,clip]{14610f1b.eps} \vspace*{-1.8mm}
\vspace{-2mm}\end{figure} Figure 1:

Number counts at 100 ( left) and 160 $\mu $m ( right), normalized to the Euclidean slope. Filled/open symbols belong to flux bins above/below the 80% completeness limit. Models belong to Lagache et al. (2004), Franceschini et al. (2010), Rowan-Robinson (2009), Le Borgne et al. (2009), Valiante et al. (2009), Lacey et al. (2010), Gruppioni & Pozzi (in prep.). Shaded areas represent ISO and Spitzer data (Héraudeau et al. 2004; Béthermin et al. 2010; Rodighiero & Franceschini 2004). Inset: Collection of PACS number counts, including this work (red, limited to relative errors $\le $15% for clarity sake), PEP Abell 2218 (grey, Altieri et al. 2010) and HerMES-PACS (light-blue, Aussel et al., in prep.). The solid lines in the insets mark the trends expected for a non-evolving population of galaxies.

Open with DEXTER

Here we exploit the science demonstration phase (SDP) observations of the GOODS-N field, complemented with the COSMOS and Lockman hole (LH) wider and shallower layers of the survey, to build far-IR galaxy number counts and derive the fraction of CIB resolved by PEP. We then take advantage of the extensive multi-wavelength coverage of the GOODS-N field, with the aim of identifying the CIB contributors and the epoch when the bulk of the Universe IR energy budget was emitted.

2 Number counts

The area used to derive number counts includes observations in the GOODS-N ($\sim$140 arcmin2), COSMOS ($\sim$2 deg2) and LH ($\sim$450 arcmin2) fields. Appendix A describes data reduction and catalog extraction.

PACS number counts were computed in each field within 0.1 dex flux bins, and are reported in Fig. 1 and Table 1, normalized to the Euclidean slope ( ${\rm d}N/{\rm d}S\propto S^{-2.5}$). Incompleteness and spurious sources have been dealt with using Monte Carlo simulations (see Appendix A). Filled symbols belong to the conservative 80% completeness limit. Thanks to the method described in Chary et al. (2004), we extend the analysis down to the faintest objects detected at the 3$\sigma $ level. Error bars include Poisson statistics, flux calibration uncertainties, and photometric errors. The latter have been propagated into number counts via 104 realizations of random Gaussian flux errors applied to each PACS source, using a dispersion equal to the local measured noise. The 95 $\mu $m ISO counts (Héraudeau et al. 2004; Rodighiero & Franceschini 2004) and the 160 $\mu $m Spitzer results by Béthermin et al. (2010, including GOODS/FIDEL, COSMOS and SWIRE fields) are consistent with PACS data, excepted the ISO faintest flux bin.

Small insets in Fig. 1 show results from PACS guaranteed-time surveys: this work, Abell 2218 lensed counts (Altieri et al. 2010), and HerMES (Aussel et al., in prep.). The three datasets complement each other and show overall agreement, the 160 $\mu $m HerMES counts excepted. This discrepancy can be ascribed to the sheet of galaxies known to bias the HerMES LH North field, on which their counts are based (Owen & Morrison 2009).

Differential, Euclidean-normalized counts display various well-known features: below $\sim$100 mJy at 100 $\mu $m and $\sim$200 mJy at 160 $\mu $m, the slope is super-Euclidean, reaching a peak in the counts between $\sim$5-10 and 30 mJy, while at fainter fluxes the slope quickly becomes sub-Euclidean. These attributes are generally interpreted as indicating evolving properties in the IR galaxy population (see Sect. 4): the trend expected for no-evolution is shown in Fig. 1 insets.

The non-normalized differential counts can be represented by a power law of the type ${\rm d}N/{\rm d}S\propto S^{\alpha }$, with two distinct slopes at fluxes fainter/brighter than $S_{\rm break}$. The results of a weighted least squares fit are given in Table 2 for the different PEP fields. The $S_{\rm break}$ flux is $\sim$5.0 mJy at 100 $\mu $m and $\sim$8.5 mJy at 160 $\mu $m. At the bright end, the number of sources in the small fields is low, and error bars are dominated by Poisson statistics. Consequently, the large uncertainties allow for nearly-Euclidean fits in GOODS-N and the LH. In the case of COSMOS, the slope is significantly steeper than Euclidean at 160 $\mu $m, while the 100 $\mu $m counts are still very flat. At the faint end, the 160 $\mu $m slope derived from our counts (GOODS-N only) is consistent with what is found by stacking 24 $\mu $m sources on Spitzer 160 $\mu $m maps (Béthermin et al. 2010, $\alpha=-1.61\pm0.21$).

Table 2:   Power-law fit to PACS differential number counts in the form ${\rm d}N/{\rm d}S\propto S^{\alpha }$.

3 Resolved CIB fraction

The added contribution of resolved galaxies provides a lower limit to the IR background and can be compared to direct measurements of the total CIB. After the discovery of the CIB (Puget et al. 1996; Hauser et al. 1998), numerous authors attempted to directly measure its surface brightness from COBE/DIRBE maps (e.g. Renault et al. 2001; Lagache et al. 2000; Dole et al. 2006; Wright 2004). Here we adopt the most recent revision of the Lagache et al. (2000) DIRBE measurements, provided in Dole et al. (2006): $14.4\pm6.3$, $12.0\pm6.9$, and $12.3\pm2.5$ [nW m-2 sr-1], at 100, 140, and 240 $\mu $m, respectively. Interpolating between these values, one obtains a value of $12.8\pm6.4$ [nW m-2 sr-1] at 160 $\mu $m. These direct measurements are still affected by large uncertainties, mainly due to difficulties in defining an absolute flux scale and in removing zodiacal light. An alternative estimate of the total CIB is obtained by integrating the power-law extrapolation of our number counts, between 0.01 and 1000 mJy. We adopt the slope derived on GOODS-N data for $S<S_{\rm break}$ and an average between COSMOS and Lockman Hole at brighter fluxes (see Table 2 and Sect. 2). The total CIB estimates from this extrapolation are $\sim$12 and $\sim$13 [nW m-2 sr-1]in the green (100 $\mu $m) and red (160 $\mu $m) bands, within the errors of the direct measurements.

The contribution to the CIB of individual GOODS-N sources detected above 3$\sigma $ is $\nu I_\nu = 4.46\pm0.52$ at 100 $\mu $m and $4.41\pm 0.62$ [nW m-2 sr-1] at 160 $\mu $m, i.e., ${\sim}31\pm4$% and ${\sim}35\pm5$% of the direct DIRBE measurements (Fig. 2). Taking the effects of completeness into account (see Sect. 2), we can compute the CIB fraction due to the general FIR-population, overriding source extraction losses. We limit the analysis to the 3$\sigma $ detection limit again. The GOODS-N number counts cover the ranges 3-45 mJy (at 100 $\mu $m) and 5.5-72 mJy (at 160 $\mu $m); the COSMOS field allows upper boundaries to be extended to 142 and 179 mJy. Integrating the counts over the whole flux range, we obtain $\nu I_\nu=6.36\pm1.67$ and $6.58\pm1.62$ [nW m-2 sr-1], within 1.3$\sigma $ and 1.0$\sigma $ from the reference values at 100 and 160 $\mu $m, respectively (Fig. 2). These correspond to ${\sim}45\pm12$% and ${\sim}52\pm13$% of the Dole et al. (2006) estimate. The derived uncertainty now also includes the effect of Poisson statistics and not only photometric errors. The bright end of number counts, covered only by COSMOS, gives a small contribution ( ${\sim}4{-}6$%) to the total CIB surface brightness.

The PACS detection rate of Spitzer 24 $\mu $m sources ( $S_{24}\ge20$ $\mu $Jy, Magnelli et al. 2009) is roughly 15%. It is possible to derive a deeper CIB estimate, through stacking on the PEP maps at the positions of all 24 $\mu $m objects, including those not detected in the FIR (e.g. Dole et al. 2006). Uncertainties on the stacked fluxes are computed via a simple bootstrap procedure. The CIB surface brightness produced by 24 $\mu $m sources is $7.39\pm0.48$ and $9.57\pm 0.71$ [nW m-2 sr-1] at 100 and 160 $\mu $m, consistent with Béthermin et al. (2010), and providing $\sim$51% and $\sim$75% of the total background, within 1.1$\sigma $ and 0.5$\sigma $ from Dole et al. (2006).

\begin{figure}
\par\includegraphics[width=7.7cm,clip]{14610f2.eps}\vspace*{-0.8mm}
\vspace{-8mm}\end{figure} Figure 2:

Contribution of GOODS-N galaxies to the CIB. Red empty triangles represent the sum of individually detected sources; black filled triangles come from number counts above the 3$\sigma $ limit, and include completeness correction; blue upside-down empty triangles belong to stacking of 24 $\mu $m sources. Solid histograms denote the contribution of different redshift bins to the CIB, accounting for completeness correction. Hatched histograms are stacking results. Literature data include: DIRBE measurements revised by Dole et al. (2006, filled circles, 1)tex2html_wrap_inline2023# errors#; FIRAS spectrum (solid lines above 200 $\mu $m, Lagache et al. 2000,1999); Fixsen et al. (1998) modified Black Body (shaded area).

Open with DEXTER

\begin{figure}
\par\includegraphics[height=7.3cm,clip]{14610f3a.eps}\hspace*{4mm}
\includegraphics[height=7.3cm,clip]{14610f3b.eps} %
\vspace{-3mm}\end{figure} Figure 3:

Slicing the PEP GOODS-N population into redshift bins. Each PACS source has been associated to an optical/mid-IR counterpart as described in Appendix A. See Fig. 1 and text for details on the models shown.

Open with DEXTER

4 Discussion

In the attempt to reproduce the observed ISO and Spitzer number counts, several authors built ``backward'' evolutionary models, including luminosity and/or density evolution, as well as different galaxy populations. In Fig. 1, we overlay recent models onto the observed PACS counts. We include in this collection the Lacey et al. (2010) $\Lambda$-CDM semi-analytical model (SAM), complemented with radiative transfer dust reprocessing. The most successful models are the Valiante et al. (2009), including luminosity-dependent distribution functions for the galaxy IR SEDs and their AGN contribution, and Rowan-Robinson (2009), employing analytic evolutionary functions without discontinuities and 4 galaxy populations.

Gruppioni and Lagache et al. (2004) overestimate the amplitude of the number counts peak in both bands, while Franceschini et al. (2010) and Le Borgne et al. (2009) reproduce the counts fairly well only in one channel (100 and 160 $\mu $m, respectively). It is worth recalling that to date most of these models have been fine-tuned to reproduce mid-IR and sub-mm statistics, while a big gap in wavelength was affecting far-IR predictions. Most include a luminosity evolution $\propto(1+z)^{3.0-3.5}$, but the redshift limit for this slope, the details of density evolution, or the adopted galaxy zoo vary significantly from author to author.

Besides spanning a much wider range of observational data (UV, optical, near-IR luminosity functions, galaxy sizes, metallicity, etc.), the SAM approach suffers here for a limited flexibility in the choice of parameters, and significantly overestimates the bright end of PACS counts. Moreover, it cannot reproduce the peak, especially in the green band.

Thanks to the rich ancillary dataset in GOODS-N, we split the far-IR number counts into redshift bins (Fig. 3). This elaboration offers a remarkable chance to set detailed constraints on the evolution of the galaxy populations adopted in current recipes. This view highlights some new features and the main problems of the models under discussion.

First, in the lowest redshift bin, $0.0<z\le0.5$, the differential number counts, normalized to the Euclidean slope, monotonically increase as a function of flux, resembling the trend expected for a non-evolving population of galaxies. The consistency of models with data weakens at longer wavelengths. At the bright end, above 10 mJy and 20 mJy at 160 $\mu $m, PEP catalogs include 43 and 18 objects in this redshift range. Based on template fitting (Rodighiero et al. 2010), half of this sample shows a late-spiral SED (Sd-Sdm in the Polletta et al. 2007, template library) with typical luminosities above $10^{11}~L_\odot$ and a peak wavelength of $\sim$130 $\mu $m. Nevertheless, at these luminosities, Chary & Elbaz (2001) models, which adopt the local L-T relation, predict SEDs peaking around 90 $\mu $m restframe, thus indicating that this class is indeed colder than typical local IR-luminous sources. This is reminiscent of the cold sources found by ISO FIRBACK and serendipity surveys (e.g. Stickel et al. 2007; Dennefeld et al. 2005).

Significant differences between models arise at higher redshift. The bin $1.0<z\le2.0$ shows the largest discrepancies: apparently, most of the overestimation of total number counts by some models originates here, suggesting that the adopted evolution is either too protracted in time or too steep.

We are now able for the first time to infer how much of the resolved CIB peak was emitted at different cosmic epochs. Figure 2 depicts the amount of extragalactic IR background originating from the four redshift bins taken into account (see also Fig. 3). At 100 $\mu $m, at the GOODS-N 3$\sigma $ limit of $\sim$3 mJy, roughly 80% of the resolved CIB is emitted by objects detected at $z\le1$, equally distributed below and above z=0.5. The remainder belongs to more distant galaxies, with only 5% locked into individual objects at z>2. Most of the resolved 160 $\mu $m CIB is produced at higher redshift, with the $0.5<z\le1.0$ redshift bin dominating ($\sim$40%) the budget.

At the PEP depth, the CIB resolved into individual sources is mainly due to luminous and ultraluminous IR galaxies (LIRGs $10^{11}\le L_{\rm IR}<10^{12}~L_\odot$, ULIRGs $L_{\rm IR}\ge10^{12}~L_\odot$). The relative fraction of the two classes varies as a function of redshift, likely due to selection effects, with LIRGs providing 80-90% of the CIB resolved at $0.5<z\le1$ and ULIRGs dominating at higher redshift (60-70% at $1.0<z\le2.0$ and 100% above).

Finally, Fig. 2 also reports the CIB contribution of 24 $\mu $m sources, as obtained through stacking (see Sect. 3). The relative CIB fractions emitted in the four redshift bins change significantly at 100 $\mu $m, with respect to those obtained by integrating the observed counts. The redshift distribution of the CIB is now peaking between $0.5<z\le1.0$ and the relative contribution of z>1 increases by more than 10%. On the other hand, at 160 $\mu $m only the relative fraction of the highest redshift bin varies by more than 5%, while the others remain unchanged within 1-2%. Also in the case of stacking, the 160 $\mu $m CIB relative redshift distribution is more populated at high redshift than the 100 $\mu $m one. As expected, this indicates that PACS galaxies with redder observed colors lie on average at higher redshift than bluer ones.

Acknowledgements
PACS has been developed by a consortium of institutes led by MPE (Germany) and including UVIE (Austria); KU Leuven, CSL, IMEC (Belgium); CEA, LAM (France); MPIA (Germany); INAF-IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); IAC (Spain). This development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI/INAF (Italy), and CICYT/MCYT (Spain).

References

  1. Altieri, B., et al. 2010, A&A, 518, L17 Barger, A. J., Cowie, L. L., & Wang, W. 2008, ApJ, 689, 687 [NASA ADS] [CrossRef] [Google Scholar]
  2. Béthermin, M., Dole, H., Beelen, A., & Aussel, H. 2010, A&A, 512, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Brammer, G. B., van Dokkum, P. G., & Coppi, P. 2008, ApJ, 686, 1503 [NASA ADS] [CrossRef] [Google Scholar]
  4. Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 [NASA ADS] [CrossRef] [Google Scholar]
  5. Chary, R., Casertano, S., Dickinson, M. E., et al. 2004, ApJS, 154, 80 [NASA ADS] [CrossRef] [Google Scholar]
  6. Dennefeld, M., Lagache, G., Mei, S., et al. 2005, A&A, 440, 5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, A&AS, 147, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Dole, H., Lagache, G., & Puget, J. 2003, ApJ, 585, 617 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dole, H., Le Floc'h, E., Pérez-González, P. G., et al. 2004, ApJS, 154, 87 [NASA ADS] [CrossRef] [Google Scholar]
  10. Dole, H., Lagache, G., Puget, J., et al. 2006, A&A, 451, 417 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  11. Elbaz, D., Cesarsky, C. J., Chanial, P., et al. 2002, A&A, 384, 848 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Fixsen, D., Dwek, E., Mather, J., Bennett, C., & Shafer, R. 1998, ApJ, 508, 123 [NASA ADS] [CrossRef] [Google Scholar]
  13. Franceschini, A., Rodighiero, G., Vaccari, M., et al. 2010, A&A, in press [arXiv:0906.4264] [Google Scholar]
  14. Frayer, D. T., Sanders, D. B., Surace, J. A., et al. 2009, AJ, 138, 1261 [NASA ADS] [CrossRef] [Google Scholar]
  15. Fruchter, A. S., & Hook, R. N. 2002, PASP, 114, 144 [NASA ADS] [CrossRef] [Google Scholar]
  16. Genzel, R., & Cesarsky, C. J. 2000, ARA&A, 38, 761 [NASA ADS] [CrossRef] [Google Scholar]
  17. Giavalisco, M., Ferguson, H. C., Koekemoer, A. M., et al. 2004, ApJ, 600, L93 [NASA ADS] [CrossRef] [Google Scholar]
  18. Grazian, A., Fontana, A., de Santis, C., et al. 2006, A&A, 449, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Hauser, M. G., & Dwek, E. 2001, ARA&A, 39, 249 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hauser, M. G., Arendt, R. G., Kelsall, T., et al. 1998, ApJ, 508, 25 [NASA ADS] [CrossRef] [Google Scholar]
  21. Héraudeau, P., Oliver, S., del Burgo, C., et al. 2004, MNRAS, 354, 924 [NASA ADS] [CrossRef] [Google Scholar]
  22. Lacey, C. G., Baugh, C. M., Frenk, C. S., et al. 2010, MNRAS, 405, 2 [NASA ADS] [Google Scholar]
  23. Lagache, G., Abergel, A., Boulanger, F., Désert, & Puget. 1999, A&A, 344, 322 [NASA ADS] [Google Scholar]
  24. Lagache, G., Haffner, L., Reynolds, R., & Tufte, S. 2000, A&A, 354, 247 [Google Scholar]
  25. Lagache, G., Dole, H., & Puget, J. 2003, MNRAS, 338, 555 [NASA ADS] [CrossRef] [Google Scholar]
  26. Lagache, G., Dole, H., Puget, J., et al. 2004, ApJS, 154, 112 [NASA ADS] [CrossRef] [Google Scholar]
  27. Lagache, G., Puget, J., & Dole, H. 2005, ARA&A, 43, 727 [NASA ADS] [CrossRef] [Google Scholar]
  28. Le Borgne, D., Elbaz, D., Ocvirk, P., & Pichon, C. 2009, A&A, 504, 727 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Magnelli, B., Elbaz, D., Chary, R. R., et al. 2009, A&A, 496, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Marsden, G., Ade, P. A. R., Bock, J. J., et al. 2009, ApJ, 707, 1729 [NASA ADS] [CrossRef] [Google Scholar]
  31. Ott, S., et al. 2010, ASP Conf Series, ADASS XIX, Mizumoto et al. eds, in press [Google Scholar]
  32. Owen, F. N., & Morrison, G. E. 2009, ApJS, 182, 625 [NASA ADS] [CrossRef] [Google Scholar]
  33. Papovich, C., Dole, H., Egami, E., et al. 2004, ApJS, 154, 70 [NASA ADS] [CrossRef] [Google Scholar]
  34. Pilbratt, G. L., et al. 2010, A&A, 518, L1 [CrossRef] [EDP Sciences] [Google Scholar]
  35. Poglitsch, A., et al. 2010, A&A, 518, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81 [NASA ADS] [CrossRef] [Google Scholar]
  37. Puget, J., Abergel, A., Bernard, J., et al. 1996, A&A, 308, L5 [NASA ADS] [Google Scholar]
  38. Renault, C., Barrau, A., Lagache, G., & Puget, J. 2001, A&A, 371, 771 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Rodighiero, G., & Franceschini, A. 2004, A&A, 419, L55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Rodighiero, G., et al. 2010, A&A, 518, L25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Rowan-Robinson, M. 2009, MNRAS, 394, 117 [NASA ADS] [CrossRef] [Google Scholar]
  42. Soifer, B. T., Helou, G., & Werner, M. 2008, ARA&A, 46, 201 [NASA ADS] [CrossRef] [Google Scholar]
  43. Starck, J., & Murtagh, F. 1998, PASP, 110, 193 [NASA ADS] [CrossRef] [Google Scholar]
  44. Stickel, M., Klaas, U., & Lemke, D. 2007, A&A, 466, 831 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Sutherland, W., & Saunders, W. 1992, MNRAS, 259, 413 [NASA ADS] [CrossRef] [Google Scholar]
  46. Valiante, E., Lutz, D., Sturm, E., Genzel, R., & Chapin, E. 2009, ApJ, 701, 1814 [NASA ADS] [CrossRef] [Google Scholar]
  47. Wright, E. L. 2004, New Astron. Rev., 48, 465 [NASA ADS] [CrossRef] [Google Scholar]

Online Material

Table 1:   PEP number counts, normalized to the Euclidean slope.

Appendix A: PACS data

PEP science demonstration data include the GOODS-N and Abell 2218 fields. In addition to these, part of the other PEP blank fields have already been scheduled and observed: Lockman Hole and COSMOS ($\sim$85% of the planned depth). Observations of all fields were carried out by adopting the intermediate speed (20 arcsec/s) scan-map mode. Table A.1 lists the total exposure times already observed for these fields[*].

Data have been processed through the standard PACS reduction pipeline, version 2.0.1328, within the HCSS environment[*] (Ott et al. 2010). Additionally, we employed custom procedures aimed at removing of interference patterns, tracking anomalies, re-centering positional offsets, and mapping.

Glitch removal is based on multi-resolution median transform, developed by Starck & Murtagh (1998) to detect faint sources in ISOCAM data. The signal due to real sources and glitches show different signatures in the pixel timeline. These features are recognized using a multi-scale transform, separating the various frequencies of the signal. Once the glitch components are identified, they are replaced by interpolated values in the pixel timeline.

PACS photometers exhibit a noise with a roughly f-0.5 spectrum at relevant frequencies. To remove the bulk of the noise we apply a ``running-box'' high-pass median filter to each pixel timeline, but mask the position of bright sources. The objects mask is produced iteratively during the reduction by detecting sources on the final map, and then we mask them in a double-pass mapping scheme. Testing shows that this masked filtering method modifies the fluxes of point-like source by less than 5%.

Imperfections, drifts, and errors in the pointing accuracy of the Herschel satellite were corrected by re-centering the data on a grid of known 24 $\mu $m sources populating the fields. Such objects were stacked for all scan-legs in a given direction in a given map repetition (or for a subset of those, in the very large COSMOS maps). The stacking result was then used to compute the average offset to be applied to this set of scan-legs, for a given direction, in a given map repetition. This procedure also implicitly corrects for small timing offsets between pointing information and data. Absolute, systematic astrometric offsets turned out to be as high as 5 arcsec, while relative corrections between individual submaps are approximately 1 arcsec.

Table A.1:   PEP fields: total exposure times, noise properties, flux levels for 80% completeness and 30% spurious fraction, statistics of blind catalogs, and results of maximum-likelihood match to the multi-wavelength ancillary catalogs (labeled ``multi'', GOODS-N only).

Map reconstruction is done via simple image co-addition, based on a simplified version of the ``drizzle'' method (Fruchter & Hook 2002). Given the high data redundancy in the GOODS-N field, the drop size is set to 1/8 of the input array pixel size. This corresponds to 1/5 and $\sim$1/4 of the output pixel size at 100 $\mu $m and 160 $\mu $m, respectively, thus reducing the correlated noise in the final map. Fields with lower redundancy were mapped by adopting a smaller drop size (1/4 of the input PACS array pixel size). Images produced from each observation were weighted according to the effective exposure of each pixel and co-added to produce the final maps. The final error map was computed as the standard deviation of the weighted mean. Owing to the nature of scan maps, correlations exist between nearby pixels, in particular along the scan direction. These correlations are close to uniform across the final map, thus we derived a mean correlation correction factor which was then accounted for in the errors on the extracted fluxes.

PACS catalogs were extracted following two different approaches, optimized for the different scientific aims of the PEP project. We performed a blind extraction using the Starfinder PSF-fitting code (Diolaiti et al. 2000) and a guided extraction using 24 $\mu $m priors, following the method described in Magnelli et al. (2009). The two methods provide similar results: fluxes extracted in the two cases are consistent with each other, the prior extraction leading to slightly deeper - although possibly biased - catalogs. The number counts presented in this paper were built on the blind catalog. Point spread function (PSF) profiles were extracted from the final science maps, and turn out to have an FWHM of $\sim$7.5 and $\sim$11 arcsec in the 100 $\mu $m and 160 $\mu $m bands, respectively. Aperture corrections were characterized on calibration observations of the Vesta asteroid. Absolute flux calibration is based on $\gamma$Dra, $\alpha$Tau, and $\alpha$CMa and makes use of the standard calibration file embedded in the PACS pipeline (version 2.0.1328). Typical absolute flux calibration errors are $\sim$10% and include uncertainties on instrumental characterization, PSFs, and reference stars analysis.

Noise in PACS maps was measured with random aperture extractions on residual images and compared to the observed S/N ratio for the detected sources. The rms values thus obtained include both instrumental noise and confusion noise due to undetected sources (i.e., below the 3$\sigma $ threshold). This measured 1$\sigma $ noise is 1.00 mJy at 100 $\mu $m and 1.90 mJy at 160 $\mu $m for GOODS-N. Table A.1 includes the noise properties of the SDP fields, as well as the Lockman Hole and COSMOS.

To quantify the reliability of extracted fluxes, the level of incompleteness and the fraction of spurious sources, Monte Carlo simulations were performed, creating 500 images and adding 20 artificial objects onto science maps each, for a total of 10 000 sources. Input and output fluxes are consistent with each other within a few percent. Completeness is defined as the fraction of sources that have been detected with a photometric accuracy of at least 50% (Papovich et al. 2004). Spurious sources are defined as those extracted above 3$\sigma $ with an input flux lower than 3 $\sigma(Image)$. The latter is consistent with the spurious fraction inferred by blindly extracting from inverted maps. The GOODS-N blind catalog reaches 80% completeness at $\sim$5.5 mJy and $\sim$11.0 mJy in the two bands, and a 30% fraction of spurious detections at $\sim$2.5 and $\sim$7.0 mJy, in green and red respectively (see Table A.1).

The GOODS-N field benefits from an extensive multi-wavelength coverage. Adopting the Grazian et al. (2006) approach, the PEP Team built a reliable multi-wavelength, PSF-matched database, including ACS bviz (Giavalisco et al. 2004), Flamingos JHK[*] and Spitzer IRAC data. Moreover, MIPS 24 $\mu $m (Magnelli et al. 2009) and Barger et al. (2008) deep U, $K{\rm s}$, and spectroscopic redshifts have been added. When no spectroscopic redshifts were available, photometric redshifts have been derived using the EAZY code (Brammer et al. 2008).

The 100 and 160 $\mu $m blind catalogs were finally linked to this multi-wavelength catalog through a three-band maximum likelihood procedure (Sutherland & Saunders 1992), starting from the longest wavelength available (160 $\mu $m, PACS) and progressively matching 100 $\mu $m (PACS) and 24 $\mu $m (Spitzer/MIPS) data. Table A.1 summarizes the main properties of blind catalogs, as well as the results of the maximum-likelihood match in GOODS-N.

Confusion is a major concern in deep wide-beam observations, like far-IR or sub-mm imaging of blank fields. PEP Herschel observations are not dispensed from confusion: the high density of detected sources hinders the extraction of fainter objects, in the so-called ``source density confusion criterion'' (SDC, Dole et al. 2003). Adopting the Lagache et al. (2003) definition of beam (i.e., $\Omega=1.14\times\theta_{\it FWHM}^2$), the source density in PACS blind catalogs is 40 beams/source at 100 $\mu $m and 18 beams/source at 160 $\mu $m. This indicates that PEP GOODS-N is already hitting the SDC limit in the PACS red band, estimated to be 16.7 beams/source by Dole et al. (2003), while the green channel is not affected. By extrapolating integral number counts, we estimate that this limit will be reached at $\sim$2.0 mJy at 100 $\mu $m and $\sim$4.7 mJy at 160 $\mu $m. Nevertheless, this limit is known to be rather conservative, and some authors have already shown that source extraction can be reliably carried out to levels as low as 10 beam/source under favorable conditions (e.g. at 24 $\mu $m, Magnelli et al. 2009).

Footnotes

...Herschel/PEP[*]
Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
...[*]
Table 1 and Appendix A are only available in electronic form at http://www.aanda.org
... fields[*]
See the PEP web page for information about the other fields in the survey: http://www.mpe.mpg.de/ir/Research/PEP/.
... environment[*]
HCSS is a joint development by the Herschel Science Ground Segment Consortium, consisting of ESA, the NASA Herschel Science Center, and the HIFI, PACS and SPIRE consortia.
...JHK[*]
Kindly reduced by Kyoungsoo Lee.

All Tables

Table 2:   Power-law fit to PACS differential number counts in the form ${\rm d}N/{\rm d}S\propto S^{\alpha }$.

Table 1:   PEP number counts, normalized to the Euclidean slope.

Table A.1:   PEP fields: total exposure times, noise properties, flux levels for 80% completeness and 30% spurious fraction, statistics of blind catalogs, and results of maximum-likelihood match to the multi-wavelength ancillary catalogs (labeled ``multi'', GOODS-N only).

All Figures

  \begin{figure}
\par\includegraphics[width=7.6cm,clip]{14610f1a.eps}\hspace*{4mm}...
...phics[width=7.6cm,clip]{14610f1b.eps} \vspace*{-1.8mm}
\vspace{-2mm}\end{figure} Figure 1:

Number counts at 100 ( left) and 160 $\mu $m ( right), normalized to the Euclidean slope. Filled/open symbols belong to flux bins above/below the 80% completeness limit. Models belong to Lagache et al. (2004), Franceschini et al. (2010), Rowan-Robinson (2009), Le Borgne et al. (2009), Valiante et al. (2009), Lacey et al. (2010), Gruppioni & Pozzi (in prep.). Shaded areas represent ISO and Spitzer data (Héraudeau et al. 2004; Béthermin et al. 2010; Rodighiero & Franceschini 2004). Inset: Collection of PACS number counts, including this work (red, limited to relative errors $\le $15% for clarity sake), PEP Abell 2218 (grey, Altieri et al. 2010) and HerMES-PACS (light-blue, Aussel et al., in prep.). The solid lines in the insets mark the trends expected for a non-evolving population of galaxies.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{14610f2.eps}\vspace*{-0.8mm}
\vspace{-8mm}\end{figure} Figure 2:

Contribution of GOODS-N galaxies to the CIB. Red empty triangles represent the sum of individually detected sources; black filled triangles come from number counts above the 3$\sigma $ limit, and include completeness correction; blue upside-down empty triangles belong to stacking of 24 $\mu $m sources. Solid histograms denote the contribution of different redshift bins to the CIB, accounting for completeness correction. Hatched histograms are stacking results. Literature data include: DIRBE measurements revised by Dole et al. (2006, filled circles, 1)tex2html_wrap_inline2023# errors#; FIRAS spectrum (solid lines above 200 $\mu $m, Lagache et al. 2000,1999); Fixsen et al. (1998) modified Black Body (shaded area).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[height=7.3cm,clip]{14610f3a.eps}\hspace*{4mm}
\includegraphics[height=7.3cm,clip]{14610f3b.eps} %
\vspace{-3mm}\end{figure} Figure 3:

Slicing the PEP GOODS-N population into redshift bins. Each PACS source has been associated to an optical/mid-IR counterpart as described in Appendix A. See Fig. 1 and text for details on the models shown.

Open with DEXTER
In the text


Copyright ESO 2010

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.