E. Wieprecht

PACS Evolutionary Probe (PEP) - A Herschel Key Program

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role within the entire set of Herschel surveys, and the field selection that includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, and EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high-redshift galaxy populations, thus testing and superseding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

The first Herschel view of the mass-SFR link in high-z galaxies

Aims. We exploit deep observations of the GOODS-N field taken with PACS, the Photodetector Array Camera and Spectrometer, onboard of Herschel, as part of the PACS evolutionary probe guaranteed time (PEP), to study the link between star formation and stellar mass in galaxies to z ∼ 2. Methods. Starting from a stellar mass – selected sample of ∼4500 galaxies with mag4.5 μm < 23.0 (AB), we identify ∼350 objects with a PACS detection at 100 or 160 μ ma nd∼ 1500 with only Spitzer 24 μm counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the spectral energy distributions (SEDs). Results. Consistently with other Herschel results, we find that LIR based only on 24 μm data is overestimated by a median factor ∼ 1. 8a tz ∼ 2, whereas it is underestimated (with our approach) up to a factor ∼ 1. 6a t 0.5 < z < 1.0. We then exploit this calibration to correct LIR based on the MIPS/Spitzer fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the star formation rate (SFR), of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the SFR, the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z < 1.0 and reaching a slope of α = −0.50 +0.13 −0.16 at z ∼ 2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ∼15 for

Star formation in AGN hosts in GOODS-N

Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high-z AGN hosts. Combining PACS 100 and 160 μm observations of the GOODS-N field with 2 Ms Chandra data, we detect ∼20% of X-ray AGN individually at >3σ. The host far-infrared luminosity of AGN with L2−10 keV ≈ 10 43 erg s −1 increases with redshift by an order of magnitude from z = 0 to z ∼ 1. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for L2−10 keV 1. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample which is dominated by L2−10 keV < 10 44 erg s −1 AGN. In conjunction with properties of local and luminous high-z AGN, we interpret these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star formation is traced locally over a wide range of luminosities and also extends to luminous high-z AGN. This correlation reflects an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive galaxies and shows little dependence on AGN luminosity. The level of this secular, non-merger driven star formation increasingly dominates over the correlation at increasing redshift.

Dissecting the cosmic infra-red background with Herschel/PEP

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 μ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σ depths of ∼3.0 mJy at 100 μ ma nd∼5.7 mJy at 160 μm. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 μ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 ≤ 0.5 we isolate a class of luminous sources (LIR ∼ 10 11 L� ), whose SEDs resemble late-spiral galaxies, peaking at ∼130 μ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± 1.67 and 6.58± 1.62 [nW m −2 sr −1 ] at 100 and 160 μm, resolving ∼45% and ∼52% of the CIB, respectively. When stacking 24 μm sources, the inferred CIB lies within 1.1σ and 0.5σ from direct measurements in the two bands, and fractions increase to 50% and 75%. Most of this resolved CIB fraction was radiated at z ≤ 1.0, with 160 μm sources found at higher redshift than 100 μm ones.

The star-formation rates of 1.5 < z < 2.5 massive galaxies

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z ∼ 2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5 12.2 L� . The SFGs and AGNs tend to exhibit the same 24 μm excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

Far-infrared properties of submillimeter and optically faint radio galaxies

We use deep observations obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space Observatory to study the far-infrared (FIR) properties of submillimeter and optically faint radio galaxies (SMGs and OFRGs). From literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. This sample is cross-matched with our PACS 100 μm and 160 μm multi-wavelength catalogs based on sources-extraction using prior detections at 24 μm. About half of the galaxies in our sample are detected in at least the PACS 160 μm bandpass. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 μm (only the upper limits for the OFRGs) flux densities with a single modified (β = 1.5) black body function. The median dust temperature of our SMG sample is Tdust = 36 ± 8 K while for our OFRG sample it is Tdust = 47 ± 3 K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005, ApJ, 622, 772) found a dust temperature of Tdust = 36± 7K for a large sample of SMGs assuming the validity of the FIR/radio correlation (i.e., q = log 10 (LFIR[W]/L1.4 GHz[W Hz −1 ]/3.75 × 10 12 )). The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we

New Wavelength Determinations of Mid-Infrared Fine-structure Lines by Infrared Space Observatory Short Wavelength Spectrometer

We report accurate new wavelengths for 29 mid-infrared ionic fine-structure lines, based on observations with the Short Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO). Our results originate from observations of NGC 7027, NGC 6543, NGC 6302, the Circinus galaxy, Sgr A West, and W51 IRS 2. The obtained accuracies (λ/Δλ) range from 3 × 104 to 1 × 105, depending on instrumental mode and uncertainty in radial velocities.

The effect of environment on star forming galaxies at redshift - I. First insight from PACS

We use deep 70, 100 and 160 μm observations taken with PACS, the Photodetector Array Camera and Spectrometer on board of Herschel, as part of the PACS Evolutionary Probe (PEP) guaranteed time, to study the relation between star formation rate and environment at redshift ∼1 in the GOODS-S and GOODS-N fields. We use the SDSS spectroscopic catalog to build the local analog and study the evolution of the star formation activity dependence on the environment. At z ∼ 1 we observe a reversal of the relation between star formation rate and local density, confirming the results based on Spitzer 24 μm data. However, due to the high accuracy provided by PACS in measuring the star formation rate also for AGN hosts, we identify in this class of objects the cause for the reversal of the density–SFR relation. Indeed, AGN hosts favor high stellar masses, dense regions and high star formation rates. Without the AGN contribution the relation flattens consistently with respect to the local analog in the same range of star formation rates. As in the local universe, the specific star formation rate anti-correlates with the density. This is due to mass segregation both at high and low redshift. The contribution of AGN hosts does not affect this anti-correlation, since AGN hosts exhibit the same specific star formation rate as star forming galaxies at the same mass. The same global trends and AGN contribution is observed once the relations are studied per morphological type. We study the specific star formation rate vs. stellar mass relation in three density regimes. Our data provides an indication that at M/M� > 10 11 the mean specific star formation rate tends to be higher at higher density, while the opposite trend is observed in the local SDSS star forming sample.

A Search for Broad Infrared Recombination Lines in NGC 1068

We report infrared spectroscopy of the prototypical Seyfert 2 galaxy NGC 1068, aiming at detection of broad components of hydrogen recombination lines that originate in the obscured broad-line region. Using the Short Wavelength Spectrometer on board the Infrared Space Observatory, we have observed for the first time the regions of Bβ 2.626 μm and Pfα 7.460 μm, and present improved data for Bα 4.052 μm. No significant broad components are detected, implying an equivalent visual extinction to the broad-line region of at least 50 mag and an obscuring column density of at least 1023 cm-2. While consistent with a highly obscured broad-line region, as required by the classical unified scenario, these limits are not yet significant enough to discriminate strongly between different torus models or to constrain properties of the gas causing the very large X-ray obscuration. We discuss the systematic limitations of infrared broad-line region searches and suggest that Bα may often be the most favorable transition for future searches.

Herschel deep far-infrared counts through Abell 2218 cluster-lens

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 μm observations toward the cluster lens Abell 2218 to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 μm and 2 mJy at 160 μm, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 μm. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94−35 mJy at 100 μ ma nd 1.47−35 mJy at 160 μm), we retrieve a cosmic infrared background surface brightness of ∼8.0 and ∼9.9 nW m −2 sr −1 , in the respective bands. These values correspond to 55 ± 24% and 77 ± 31% of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62 ± 25% and 88 ± 32% CIB total fractions, resolved at 100 and 160 μm, disregarding the high uncertainties of DIRBE absolute values.