sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Valtchanov, I.

Normalized to: Valtchanov, I.

204 article(s) in total. 1408 co-authors, from 1 to 103 common article(s). Median position in authors list is 24,0.

[1]  [pdf] - 1971468
SCUBA-2 observations of candidate starbursting protoclusters selected by Planck and Herschel-SPIRE
Comments: 30 pages, 8 figures, 5 tables, accepted for publication in MNRAS
Submitted: 2019-09-19, last modified: 2019-09-27
We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $2<z<3$, is consistent with that of known $z>2$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources.
[2]  [pdf] - 1887932
Confirming Herschel candidate protoclusters from ALMA/VLA CO observations
Comments: Accepted by ApJ. 24 pages, 9 figures
Submitted: 2019-05-21
ALMA 870$\mu$m continuum imaging has uncovered a population of blends of multiple dusty star-forming galaxies (DSFGs) in sources originally detected with the Herschel Space Observatory. However, their pairwise separations are much smaller that what is found by ALMA follow-up of other single-dish surveys or expected from theoretical simulations. Using ALMA and VLA, we have targeted three of these systems to confirm whether the multiple 870$\mu$m continuum sources lie at the same redshift, successfully detecting $^{12}$CO($J = 3$-2) and $^{12}$CO($J = 1$-0) lines and being able to confirm that in the three cases all the multiple DSFGs are likely physically associated within the same structure. Therefore, we report the discovery of two new gas-rich dusty protocluster cores (HELAISS02, $z = 2.171 \pm 0.004$; HXMM20, $z = 2.602 \pm 0.002$). The third target is located in the well known COSMOS overdensity at $z = 2.51$ (named CL J1001+0220 in the literature), for which we do not find any new secure CO(1-0) detection, although some of its members show only tentative detections and require further confirmation. From the gas, dust, and stellar properties of the two new protocluster cores, we find very large molecular gas fractions yet low stellar masses, pushing the sources above the main sequence, while not enhancing their star formation efficiency. We suggest that the sources might be newly formed galaxies migrating to the main sequence. The properties of the three systems compared to each other and to field galaxies may suggest a different evolutionary stage between systems.
[3]  [pdf] - 1830648
astroquery: An Astronomical Web-Querying Package in Python
Comments: Accepted to AJ. The documentation is at and the repository is at
Submitted: 2019-01-14
astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at
[4]  [pdf] - 1830594
Quantifying the suppression of the (un)-obscured star formation in galaxy cluster cores at 0.2$\lesssim$$z$$\lesssim$0.9
Comments: 34 pages, 13 figures, 15 tables; accepted for publication in MNRAS
Submitted: 2018-12-20
We quantify the star formation (SF) in the inner cores ($\mathcal{R}$/$R_{200}$$\leq$0.3) of 24 massive galaxy clusters at 0.2$\lesssim$$z$$\lesssim$0.9 observed by the $Herschel$ Lensing Survey and the Cluster Lensing and Supernova survey with $Hubble$. These programmes, covering the rest-frame ultraviolet to far-infrared regimes, allow us to accurately characterize stellar mass-limited ($\mathcal{M}_{*}$$>$$10^{10}$ $M_{\odot}$) samples of star-forming cluster members (not)-detected in the mid- and/or far-infrared. We release the catalogues with the photometry, photometric redshifts, and physical properties of these samples. We also quantify the SF displayed by comparable field samples from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We find that in intermediate-$z$ cluster cores, the SF activity is suppressed with respect the field in terms of both the fraction ($\mathcal{F}$) of star-forming galaxies (SFG) and the rate at which they form stars ($\mathcal{SFR}$ and $s\mathcal{SFR} = \mathcal{SFR}/\mathcal{M}_{*}$). On average, the $\mathcal{F}$ of SFGs is a factor $\sim$$2$ smaller in cluster cores than in the field. Furthermore, SFGs present average $\mathcal{SFR}$ and $s\mathcal{SFR}$ typically $\sim$0.3 dex smaller in the clusters than in the field along the whole redshift range probed. Our results favour long time-scale quenching physical processes as the main driver of SF suppression in the inner cores of clusters since $z$$\sim$0.9, with shorter time-scale processes being very likely responsible for a fraction of the missing SFG population.
[5]  [pdf] - 1791576
ESASky : A Science-Driven Discovery Portal for Space-Based Astronomy Missions
Comments: Published in Astronomy and Computing; 8 pages, 5 figures. ESASky: <>
Submitted: 2018-11-26
In the era of "big data" and with the advent of web 2.0 technologies, ESASky ( aims at providing a modern and visual way to access astronomical science-ready data products and metadata. The main goal of the application is to simplify the interaction between the scientific community and the ever-growing amount of data collected over the past decades from the most important astronomy missions. The ESASky concept is to offer a complementary scientific application to the more-traditional table-oriented exploitation of astronomical data, by allowing a more natural and visual approach and enabling the exploration of astronomical objects across the entire electromagnetic spectrum. To fulfill this goal, ESASky provides a multiwavelength interface to a set of astronomy data from an increasing number of missions and surveys, with the intention of becoming the single-point of entry to perform visual analysis and cross-matching among different energy ranges. A lot of effort has been invested on the design of a user-friendly, responsive Graphical User Interface (GUI) by the definition and optimisation of algorithms running behind each visual feature offered. In this contribution, we describe in detail the design and solutions adopted for the technical challenges arising during the development. We present the data services and features implemented in the latest version of ESASky (v2.1), including a Mission Planning Tool to support current James Webb Space Telescope (JWST) planning, the possibility to search for observations of Solar System Objects (planets, comets and moons) taken by astronomy missions, the integration of the SAO/NASA Astrophysics Data System (ADS) publication system in the ESASky GUI, and the retrieval of metadata and data products available within a specific region of the sky.
[6]  [pdf] - 1771717
The XXL Survey XX: The 365 cluster catalogue
Comments: A&A, in press Shortened abstract Full catalogues available online, see paper
Submitted: 2018-10-09
In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogue in its present state, as well as of associated catalogues as super-clusters and fossil groups. We release a sample containing 365 clusters in total. We give the details of the follow-up observations and explain the procedure adopted to validate the cluster spectroscopic redshifts. Considering the whole XXL cluster sample, we have provided two types of selection, both complete in a particular sense: one based on flux-morphology criteria, and an alternative based on the [0.5-2] keV flux within one arcmin of the cluster centre. We have also provided X-ray temperature measurements for 80$\%$ of the clusters having a flux larger than 9$\times$10$^{-15}$$\rm \thinspace erg \, s^{-1} \, cm^{-2}$. Our cluster sample extends from z$\sim$0 to z$\sim$1.2, with one cluster at z$\sim$2. Clusters were identified through a mean number of six spectroscopically confirmed cluster members. Our updated luminosity function and luminosity-temperature relation are compatible with our previous determinations based on the 100 brightest clusters, but show smaller uncertainties. We also present an enlarged list of super-clusters and a sample of 18 possible fossil groups. This intermediate publication is the last before the final release of the complete XXL cluster catalogue when the ongoing C2 cluster spectroscopic follow-up is complete. It provides a unique inventory of medium-mass clusters over a 50~\dd\ area out to z$\sim$1.
[7]  [pdf] - 1763192
The XXL Survey: XXVII. The 3XLSS point source catalogue
Comments: 19 pages, 11 figures, part of A&A Special Issue "The XXL Survey: second series"
Submitted: 2018-10-09
We present the version of the point source catalogue of the XXL Survey that was used, in part, in the first series of XXL papers. In this paper we release, in our database in Milan and at CDS: (i) the X-ray source catalogue with 26056 objects in two areas of 25 deg2; (ii) the associated multiwavelength catalogues with candidate counterparts of the X-ray sources in the infrared, near-infrared, optical, and ultraviolet (plus spectroscopic redshift when available); and (iii) a catalogue of spectroscopic redshifts recently obtained in the southern XXL area. We also present the basic properties of the X-ray point sources and their counterparts. Other catalogues described in the second series of XXL papers will be released contextually, and will constitute the second XXL data release.
[8]  [pdf] - 1760646
The XXL Survey XXV. Cosmological analysis of the C1 cluster number counts
Comments: 11 pages, 5 figures, 5 tables, accepted as a letter by Astronomy & Astrophysics
Submitted: 2018-10-03
Context. We present an estimation of cosmological parameters with clusters of galaxies. Aims. We constrain the $\Omega_m$, $\sigma_8$, and $w$ parameters from a stand-alone sample of X-ray clusters detected in the 50 deg$^2$ XMM-XXL survey with a well-defined selection function. Methods. We analyse the redshift distribution of a sample comprising 178 high S/N clusters out to a redshift of unity. The cluster sample scaling relations are determined in a self-consistent manner. Results. In a lambda cold dark matter ($\Lambda$CDM) model, the cosmology favoured by the XXL clusters compares well with results derived from the Planck S-Z clusters for a totally different sample (mass/redshift range, selection biases, and scaling relations). However, with this preliminary sample and current mass calibration uncertainty, we find no inconsistency with the Planck CMB cosmology. If we relax the $w$ parameter, the Planck CMB uncertainties increase by a factor of $\sim$10 and become comparable with those from XXL clusters. Combining the two probes allows us to put constraints on $\Omega_m$=0.316$\pm$0.060, $\sigma_8$=0.814$\pm$0.054, and $w$=-1.02$\pm$0.20. Conclusions. This first self-consistent cosmological analysis of a sample of serendipitous XMM clusters already provides interesting insights into the constraining power of the XXL survey. Subsequent analysis will use a larger sample extending to lower confidence detections and include additional observable information, potentially improving posterior uncertainties by roughly a factor of 3.
[9]  [pdf] - 1748121
The Herschel-PACS North Ecliptic Pole Survey
Comments: Accepted for publication in PASJ, 11 pages, 8 figures
Submitted: 2018-09-11
A detailed analysis of Herschel-PACS observations at the North Ecliptic Pole is presented. High quality maps, covering an area of 0.44 square degrees, are produced and then used to derive potential candidate source lists. A rigorous quality control pipeline has been used to create final legacy catalogues in the PACS Green 100 micron and Red 160 micron bands, containing 1384 and 630 sources respectively. These catalogues reach to more than twice the depth of the current archival Herschel/PACS Point Source Catalogue, detecting 400 and 270 more sources in the short and long wavelength bands respectively. Galaxy source counts are constructed that extend down to flux densities of 6mJy and 19mJy (50% completeness) in the Green 100 micron and Red 160 micron bands respectively. These source counts are consistent with previously published PACS number counts in other fields across the sky. The source counts are then compared with a galaxy evolution model identifying a population of luminous infrared galaxies as responsible for the bulk of the galaxy evolution over the flux range (5-100mJy) spanned by the observed counts, contributing approximate fractions of 50% and 60% to the cosmic infrared background (CIRB) at 100 microns and 160 microns respectively.
[10]  [pdf] - 1789453
The XXL Survey: XVI. The clustering of X-ray selected galaxy clusters at z~0.3
Comments: 10 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics
Submitted: 2018-07-12
Galaxy clusters trace the highest density peaks in the large-scale structure of the Universe. Their clustering provides a powerful probe that can be exploited in combination with cluster mass measurements to strengthen the cosmological constraints provided by cluster number counts. We investigate the spatial properties of a homogeneous sample of X-ray selected galaxy clusters from the XXL survey, the largest programme carried out by the XMM-Newton satellite. The measurements are compared to $\Lambda$-cold dark matter predictions, and used in combination with self-calibrated mass scaling relations to constrain the effective bias of the sample, $b_{eff}$, and the matter density contrast, $\Omega_{\rm M}$. We measured the angle-averaged two-point correlation function of the XXL cluster sample. The analysed catalogue consists of $182$ X-ray selected clusters from the XXL second data release, with median redshift $\langle z \rangle=0.317$ and median mass $\langle M_{500} \rangle\simeq1.3\cdot10^{14} M_\odot$. A Markov chain Monte Carlo analysis is performed to extract cosmological constraints using a likelihood function constructed to be independent of the cluster selection function. Modelling the redshift-space clustering in the scale range $10<r\,[$\Mpch$]<40$, we obtain $\Omega_{\rm M}=0.27_{-0.04}^{+0.06}$ and $b_{eff}=2.73_{-0.20}^{+0.18}$. This is the first time the two-point correlation function of an X-ray selected cluster catalogue at such relatively high redshifts and low masses has been measured. The XXL cluster clustering appears fully consistent with standard cosmological predictions. The analysis presented in this work demonstrates the feasibility of a cosmological exploitation of the XXL cluster clustering, paving the way for a combined analysis of XXL cluster number counts and clustering.
[11]  [pdf] - 1686697
The Herschel-ATLAS Data Release 2 Paper II: Catalogues of far-infrared and submillimetre sources in the fields at the south and north Galactic Poles
Submitted: 2017-12-19, last modified: 2018-03-19
The {\it Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 660 deg$^2$ with the PACS and SPIRE cameras in five photometric bands: 100, 160, 250, 350 and 500\mic. This is the second of three papers describing the data release for the large fields at the south and north Galactic poles (NGP and SGP). In this paper we describe the catalogues of far-infrared and submillimetre sources for the NGP and SGP, which cover 177 deg$^2$ and 303 deg$^2$, respectively. The catalogues contain 153,367 sources for the NGP field and 193,527 sources for the SGP field detected at more than 4$\sigma$ significance in any of the 250, 350 or 500\mic\ bands. The source detection is based on the 250\mic\ map, and we present photometry in all five bands for each source, including aperture photometry for sources known to be extended. The rms positional accuracy for the faintest sources is about 2.4 arc seconds in both right ascension and declination. We present a statistical analysis of the catalogues and discuss the practical issues -- completeness, reliability, flux boosting, accuracy of positions, accuracy of flux measurements -- necessary to use the catalogues for astronomical projects.
[12]  [pdf] - 1611911
Correcting the extended-source calibration for the Herschel-SPIRE Fourier-Transform Spectrometer
Comments: 9 pages, 6 figures, 1 table, MNRAS in press
Submitted: 2017-08-30, last modified: 2018-01-02
We describe an update to the Herschel-SPIRE Fourier-Transform Spectrometer (FTS) calibration for extended sources, which incorporates a correction for the frequency-dependent far-field feedhorn efficiency, $\eta_\mathrm{FF}$. This significant correction affects all FTS extended-source calibrated spectra in sparse or mapping mode, regardless of the spectral resolution. Line fluxes and continuum levels are underestimated by factors of 1.3-2 in the Spectrometer Long-Wavelength band (SLW, 447-1018 GHz; 671-294 $\mu$m) and 1.4-1.5 in the Spectrometer Short-Wavelength band (SSW, 944-1568 GHz; 318-191 $\mu$m). The correction was implemented in the FTS pipeline version 14.1 and has also been described in the SPIRE Handbook since Feb 2017. Studies based on extended-source calibrated spectra produced prior to this pipeline version should be critically reconsidered using the current products available in the Herschel Science Archive. Once the extended-source calibrated spectra are corrected for $\eta_\mathrm{FF}$, the synthetic photometry and the broadband intensities from SPIRE photometer maps agree within 2-4% -- similar levels to the comparison of point-source calibrated spectra and photometry from point-source calibrated maps. The two calibration schemes for the FTS are now self-consistent: the conversion between the corrected extended-source and point-source calibrated spectra can be achieved with the beam solid angle and a gain correction that accounts for the diffraction loss.
[13]  [pdf] - 1622394
Candidate high-z proto-clusters among the Planck compact sources, as revealed by Herschel-SPIRE
Comments: Resubmitted to MNRAS after request for minor revision
Submitted: 2017-12-19
By determining the nature of all the Planck compact sources within 808.4 deg^2 of large Herschel surveys, we have identified 27 candidate proto-clusters of dusty star forming galaxies (DSFGs) that are at least 3{\sigma} overdense in either 250, 350 or 500 $\mu$mm sources. We find roughly half of all the Planck compact sources are resolved by Herschel into multiple discrete objects, with the other half remaining unresolved by Herschel. We find a significant difference between versions of the Planck catalogues, with earlier releases hosting a larger fraction of candidate proto-clusters and Galactic Cirrus than later releases, which we ascribe to a difference in the filters used in the creation of the three catalogues. We find a surface density of DSFG candidate proto-clusters of (3.3 $\pm$ 0.7) x 10^(-2) sources deg^(-2), in good agreement with previous similar studies. We find that a Planck colour selection of S_{857}/S_{545} < 2 works well to select candidate proto-clusters, but can miss proto-clusters at z < 2. The Herschel colours of individual candidate proto-cluster members indicate our candidate proto-clusters all likely all lie at z > 1. Our candidate proto-clusters are a factor of 5 times brighter at 353 GHz than expected from simulations, even in the most conservative estimates. Further observations are needed to confirm whether these candidate proto-clusters are physical clusters, multiple proto-clusters along the line of sight, or chance alignments of unassociated sources.
[14]  [pdf] - 1602882
ESASky v.2.0: all the skies in your browser
Comments: 4 pages, 2 figures, ADASS 2017 conference proceedings
Submitted: 2017-12-11
With the goal of simplifying the access to science data to scientists and citizens, ESA recently released ESASky (, a new open-science easy-to-use portal with the science-ready Astronomy data from ESA and other major data providers. In this presentation, we announced version 2.0 of the application, which includes access to all science-ready images, catalogues and spectra, a feature to help planning of future JWST observations, the possibility to search for data of all (targeted and serendipitously observed) Solar System Objects in Astronomy images, a first support to mobile devices and several other smaller usability features. We also discussed the future evolution of the portal and the lessons learnt from the 1+ year of operations from the point of view of access, visualization and manipulation of big datasets (all sky maps, also called HiPS) and large catalogues (like e.g. the Gaia DR1 catalogues or the Hubble Source Catalogue) and the design and validation principles for the development of friendly GUIs for thin layer web clients aimed at scientists.
[15]  [pdf] - 1789392
The XXL Survey: XXIII. The Mass Scale of XXL Clusters from Ensemble Spectroscopy
Comments: 14 pages, 9 figures, 4 tables, accepted for publication in A&A
Submitted: 2017-11-19
An X-ray survey with the XMM-Newton telescope, XMM-XXL, has identified hundreds of galaxy groups and clusters in two 25 deg$^2$ fields. Combining spectroscopic and X-ray observations in one field, we determine how the kinetic energy of galaxies scales with hot gas temperature and also, by imposing prior constraints on the relative energies of galaxies and dark matter, infer a power-law scaling of total mass with temperature. Our goals are: i) to determine parameters of the scaling between galaxy velocity dispersion and X-ray temperature, $T_{\rm 300kpc}$, for the halos hosting XXL-selected clusters, and; ii) to infer the log-mean scaling of total halo mass with temperature, $\langle \ln M_{200} \, | \, T, z \rangle$. We apply an ensemble velocity likelihood to a sample of $> 1500$ spectroscopic redshifts within $132$ spectroscopically confirmed clusters with redshifts $z < 0.6$ to model, $\langle \ln \sigma_{\rm gal}\,|\,T,z\rangle$, where $\sigma_{\rm gal}$ is the velocity dispersion of XXL cluster member galaxies and $T$ is a 300 kpc aperture temperature. To infer total halo mass we use a precise virial relation for massive halos calibrated by N-body simulations along with a single degree of freedom summarizing galaxy velocity bias with respect to dark matter. For the XXL-N cluster sample, we find $\sigma_{\rm gal} \propto T^{0.63\pm0.05}$, a slope significantly steeper than the self-similar expectation of $0.5$. Assuming scale-independent galaxy velocity bias, we infer a mean logarithmic mass at a given X-ray temperature and redshift, $\langle\ln (E(z) M_{200}/10^{14}\,{\rm M}_{\odot})|T,z\rangle=\pi+\alpha \ln(T/T_p )+\beta\ln (E(z)/E(z_p) )$ using pivot values ${\rm k}T_{p}=2.2\,{\rm keV}$ and $z_p=0.25$, with normalization $\pi=0.45\pm0.24$ and slope $\alpha=1.89\pm0.15$. We obtain only weak constraints on redshift evolution, $\beta=-1.29\pm1.14$.
[16]  [pdf] - 1789389
The XXL Survey: XXII. The XXL-North spectrophotometric sample and galaxy stellar mass function in X-ray detected groups and clusters
Comments: 20 pages, 18 figures, 6 tables, accepted for publication in A&A
Submitted: 2017-10-12
The fraction of galaxies bound in groups in the nearby Universe is high (50% at z~0). Systematic studies of galaxy properties in groups are important in order to improve our understanding of the evolution of galaxies and of the physical phenomena occurring within this environment. We have built a complete spectrophotometric sample of galaxies within X-ray detected, optically spectroscopically confirmed groups and clusters (G&C), covering a wide range of halo masses at z<= 0.6. In the context of the XXL survey, we analyse a sample of 164 G&C in the XXL-North region (XXL-N), at z <= 0.6, with a wide range of virial masses (1.24 x 10^13 <=M_500 M_sun <= 6.63 x 10^14) and X-ray luminosities ( 2.27 x 10^41 <= L^XXL_500 (erg/s)<= 2.15 x10^44). The G&C are X-ray selected and spectroscopically confirmed. We describe the membership assignment and the spectroscopic completeness analysis, and compute stellar masses. As a first scientific exploitation of the sample, we study the dependence of the galaxy stellar mass function (GSMF) on global environment. We present a spectrophotometric characterisation of the G&C and their galaxies. The final sample contains 132 G&C, 22111 field galaxies and 2225 G&C galaxies with r-band magnitude <20. Of the G&C, 95% have at least three spectroscopic members, and 70% at least ten. The shape of the GSMF seems not to depend on environment (field versus G&C) or X-ray luminosity ( used as a proxy for the virial mass of the system). These results are confirmed by the study of the correlation between mean stellar mass of G&C members and L^XXL_500.We release the spectrophotometric catalogue of galaxies with all the quantities computed in this work. As a first homogeneous census of galaxies within X-ray spectroscopically confirmed G&C at these redshifts, this sample will allow environmental studies of the evolution of galaxy properties.
[17]  [pdf] - 1583012
Stacked Average Far-Infrared Spectrum of Dusty Star-Forming Galaxies from the Herschel/SPIRE Fourier Transform Spectrometer
Comments: 30 pages, 19 figures
Submitted: 2017-05-01, last modified: 2017-09-13
We present stacked average far-infrared spectra of a sample of 197 dusty, star-forming galaxies (DSFGs) at $0.005 < z < 4$ using close to 90% of the SPIRE Fourier Transform Spectrometer (FTS) extragalactic data archive from the Herschel Space Observatory based on 3.5 years of science operations. These spectra explore an observed-frame $\rm 447\,GHz-1568\,GHz$ ($\rm 191\,\mu m-671\,\mu m$) frequency (wavelength) range allowing us to observe the main atomic and molecular lines emitted by gas in the interstellar medium. The sample is sub-divided into five redshift bins at $0.005 < z < 0.05$, $0.05 < z < 0.2$, $0.2 < z < 0.5$, $0.8 < z <2$, and $2 < z < 4$. To study the dependence of observed spectral lines on total infrared luminosity, the sources in a subset of the redshift bins are stacked in luminosity bins. These stacked spectra are used to determine the average properties of the interstellar medium and dense molecular gas properties of DSFGs, in particular, the fine-structure line ([CII] 158 $\mu$m and [OI] 63 $\mu$m) luminosity ratios, and the line to far-IR luminosity ratios are used to model the gas density and radiation field strength in the photodissociation regions (PDRs). For the low-redshift sample, we additionally present the average spectral line energy distributions (SLED) of CO and $\rm{H_2O}$ rotational transitions and also consider PDR conditions based on a combination of [CI] 370 $\mu$m and 609 $\mu$m and $\rm CO (7-6)$ lines. For the high-z ($0.8 < z < 4$) sample PDR models suggest a molecular gas distribution in the presence of a radiation field that is at least a factor of 10$^3$ larger than the Milky-Way and with a neutral gas density of roughly 10$^3$ to 10$^5$ cm$^{-3}$. The corresponding PDR models for the low-z sample suggest a UV radiation field and gas density comparable to those at high-z.
[18]  [pdf] - 1580346
The XXL survey: first results and future
Comments: Proceeding of the XMM Next Decade Workshop held at ESAC, 9-11 May 2016
Submitted: 2016-10-06, last modified: 2017-07-11
The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.
[19]  [pdf] - 1598002
The implications of the surprising existence of a large, massive CO disk in a distant protocluster
Comments: Accepted by A&A. 15 pages, 8 figures and 4 tables
Submitted: 2017-01-18, last modified: 2017-07-09
It is not yet known if the properties of molecular gas in distant protocluster galaxies are significantly affected by their environment as galaxies are in local clusters. Through a deep, 64 hours of effective on-source integration with the Australian Telescope Compact Array (ATCA), we discovered a massive, M_mol=2.0+-0.2x10^11 M_sun, extended, ~40 kpc, CO(1-0)-emitting disk in the protocluster surrounding the radio galaxy, MRC1138-262. The galaxy, at z_CO=2.1478, is a clumpy, massive disk galaxy, M_star~5x10^11 M_sun, which lies 250 kpc in projection from MRC1138-262 and is a known H-alpha emitter, HAE229. HAE229 has a molecular gas fraction of ~30%. The CO emission has a kinematic gradient along its major axis, centered on the highest surface brightness rest-frame optical emission, consistent with HAE229 being a rotating disk. Surprisingly, a significant fraction of the CO emission lies outside of the UV/optical emission. In spite of this, HAE229 follows the same relation between star-formation rate and molecular gas mass as normal field galaxies. HAE229 is the first CO(1-0) detection of an ordinary, star-forming galaxy in a protocluster. We compare a sample of cluster members at z>0.4 that are detected in low-order CO transitions with a similar sample of sources drawn from the field. We confirm findings that the CO-luminosity and FWHM are correlated in starbursts and show that this relation is valid for normal high-z galaxies as well as those in overdensities. We do not find a clear dichotomy in the integrated Schmidt-Kennicutt relation for protocluster and field galaxies. Not finding any environmental dependence in the "star-formation efficiency" or the molecular gas content, especially for such an extended CO disk, suggests that environmentally-specific processes such as ram pressure stripping are not operating efficiently in (proto)clusters. (abridged)
[20]  [pdf] - 1584162
SPIRE Point Source Catalog Explanatory Supplement
Submitted: 2017-06-01
The Spectral and Photometric Imaging Receiver (SPIRE) was launched as one of the scientific instruments on board of the space observatory Herschel. The SPIRE photometer opened up an entirely new window in the Submillimeter domain for large scale mapping, that up to then was very difficult to observe. There are already several catalogs that were produced by individual Herschel science projects. Yet, we estimate that the objects of only a fraction of these maps will ever be systematically extracted and published by the science teams that originally proposed the observations. The SPIRE instrument performed its standard photometric observations in an optically very stable configuration, only moving the telescope across the sky, with variations in its configuration parameters limited to scan speed and sampling rate. This and the scarcity of features in the data that require special processing steps made this dataset very attractive for producing an expert reduced catalog of point sources that is being described in this document. The Catalog was extracted from a total of 6878 unmodified SPIRE scan map observations. The photometry was obtained by a systematic and homogeneous source extraction procedure, followed by a rigorous quality check that emphasized reliability over completeness. Having to exclude regions affected by strong Galactic emission, that pushed the limits of the four source extraction methods that were used, this catalog is aimed primarily at the extragalactic community. The result can serve as a pathfinder for ALMA and other Submillimeter and Far-Infrared facilities. 1,693,718 sources are included in the final catalog, splitting into 950688, 524734, 218296 objects for the 250\mu m, 350\mu m, and 500\mu m bands, respectively. The catalog comes with well characterized environments, reliability, completeness, and accuracies, that single programs typically cannot provide.
[21]  [pdf] - 1542889
The interstellar medium in high-redshift submillimeter galaxies as probed by infrared spectroscopy
Comments: Accepted for publication in ApJ
Submitted: 2017-01-11
Submillimeter galaxies (SMGs) at $z\gtrsim1$ are luminous in the far-infrared and have star-formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at z~1-4. We present 70-160 micron photometry and spectroscopy of the [OIV]26 micron, [FeII]26 micron, [SIII]33 micron, [SiII]34 micron, [OIII]52 micron, [NIII]57 micron, and [OI]63 micron fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected and we instead focus on stacking these spectra with observations of an additional 32 SMGs from the \herschel\ archive -- representing a complete compilation of PACS spectroscopy of SMGs. We detect [OI]63 micron, [SiII]34 micron, and [NIII]57 micron at >3sigma in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are $(0.36\pm0.12)\times10^{-3}$, $(0.84\pm0.17)\times10^{-3}$, and $(0.27\pm0.10)\times10^{-3}$, respectively. Using the [OIII]52/[NIII]57 emission line ratio we show that SMGs have average gas-phase metallicities $\gtrsim Z_{\rm sun}$. By using PDR modelling and combining the new spectral measurements with integrated far-infrared fluxes and existing [CII]158 micron data we show that SMGs have average gas densities, n, of $\sim10^{1-3}{\rm cm^{-3}}$ and FUV field strengths, $G_0\sim10^{2.2-4.5}$ (in Habing units: $1.6\times10^{-3}{\rm erg~cm^{-2}~s^{-1}}$), consistent with both local ULIRGs and lower luminosity star-forming galaxies.
[22]  [pdf] - 1521401
Calibration of Herschel SPIRE FTS observations at different spectral resolutions
Comments: 13 pages, 17 figures, 4 tables. Accepted for publication in MNRAS
Submitted: 2016-10-10
The SPIRE Fourier Transform Spectrometer on board the Herschel Space Observatory had two standard spectral resolution modes for science observations: high resolution (HR) and low resolution (LR), which could also be performed in sequence (H+LR). A comparison of the HR and LR resolution spectra taken in this sequential mode, revealed a systematic discrepancy in the continuum level. Analysing the data at different stages during standard pipeline processing, demonstrates the telescope and instrument emission affect HR and H+LR observations in a systematically different way. The origin of this difference is found to lie in the variation of both the telescope and instrument response functions, while it is triggered by fast variation of the instrument temperatures. As it is not possible to trace the evolution of the response functions through auxiliary housekeeping parameters, the calibration cannot be corrected analytically. Therefore an empirical correction for LR spectra has been developed, which removes the systematic noise introduced by the variation of the response functions.
[23]  [pdf] - 1486886
HELP: XID+, The Probabilistic De-blender for Herschel SPIRE maps
Comments: 14 pages, 10 Figures
Submitted: 2016-06-18
We have developed a new prior-based source extraction tool, XID+, to carry out photometry in the Herschel SPIRE maps at the positions of known sources. XID+ is developed using a probabilistic Bayesian framework which provides a natural framework in which to include prior information, and uses the Bayesian inference tool, Stan, to obtain the full posterior probability distribution on flux estimates. In this paper, we discuss the details of XID+ and demonstrate the basic capabilities and performance by running it on simulated SPIRE maps resembling the COSMOS field, and comparing to the current prior-based source extraction tool DESPHOT. We show that not only does XID+ perform better on metrics such as flux accuracy and flux uncertainty accuracy, we illustrate how obtaining the posterior probability distribution can help overcome some of the issues inherent with maximum likelihood based source extraction routines. We run XID+ on the COSMOS SPIRE maps from HerMES, using a 24 $\mathrm{\mu m}$ catalogue as a prior and show the marginalised SPIRE colour-colour plot and marginalised contribution to the cosmic infrared background at the SPIRE wavelengths. XID+ is a core tool arising from the Herschel Extragalactic Legacy Project (HELP) and we discuss how additional work within HELP providing prior information on fluxes can and will be utilised. The software is available at We also provide the data product for COSMOS. We believe this is the first time that the full posterior probability of galaxy photometry has been provided as a data product.
[24]  [pdf] - 1436413
The Footprint Database and Web Services of the Herschel Space Observatory
Comments: Accepted for publication in Experimental Astronomy
Submitted: 2016-06-13
Data from the Herschel Space Observatory is freely available to the public but no uniformly processed catalogue of the observations has been published so far. To date, the Herschel Science Archive does not contain the exact sky coverage (footprint) of individual observations and supports search for measurements based on bounding circles only. Drawing on previous experience in implementing footprint databases, we built the Herschel Footprint Database and Web Services for the Herschel Space Observatory to provide efficient search capabilities for typical astronomical queries. The database was designed with the following main goals in mind: (a) provide a unified data model for meta-data of all instruments and observational modes, (b) quickly find observations covering a selected object and its neighbourhood, (c) quickly find every observation in a larger area of the sky, (d) allow for finding solar system objects crossing observation fields. As a first step, we developed a unified data model of observations of all three Herschel instruments for all pointing and instrument modes. Then, using telescope pointing information and observational meta-data, we compiled a database of footprints. As opposed to methods using pixellation of the sphere, we represent sky coverage in an exact geometric form allowing for precise area calculations. For easier handling of Herschel observation footprints with rather complex shapes, two algorithms were implemented to reduce the outline. Furthermore, a new visualisation tool to plot footprints with various spherical projections was developed. Indexing of the footprints using Hierarchical Triangular Mesh makes it possible to quickly find observations based on sky coverage, time and meta-data. The database is accessible via a web site ( and also as a set of REST web service functions.
[25]  [pdf] - 1409453
H-ATLAS: A Candidate High Redshift Cluster/Protocluster of Star-Forming Galaxies
Comments: 16 pages, 7 figures, MNRAS in press
Submitted: 2016-05-20
We investigate the region around the Planck-detected z=3.26 gravitationally lensed galaxy HATLAS J114637.9-001132 (hereinafter HATLAS12-00) using both archival Herschel data from the H-ATLAS survey and using submm data obtained with both LABOCA and SCUBA2. The lensed source is found to be surrounded by a strong overdensity of both Herschel-SPIRE sources and submm sources. We detect 17 bright (S_870 >~7 mJy) sources at >4sigma closer than 5 arcmin to the lensed object at 850/870 microns. Ten of these sources have good cross-identifications with objects detected by Herschel-SPIRE which have redder colours than other sources in the field, with 350 micron flux > 250 micron flux, suggesting that they lie at high redshift. Submillimeter Array (SMA) observations localise one of these companions to ~1 arcsecond, allowing unambiguous cross identification with a 3.6 and 4.5 micron Spitzer source. The optical/near-IR spectral energy distribution (SED) of this source is measured by further observations and found to be consistent with z>2, but incompatible with lower redshifts. We conclude that this system may be a galaxy cluster/protocluster or larger scale structure that contains a number of galaxies undergoing starbursts at the same time.
[26]  [pdf] - 1399596
A complete census of Herschel-detected infrared sources within the HST Frontier Fields
Comments: 21 pages (+10 pages of tables), 18 figures, 9 tables. Published in MNRAS. All imaging, catalogues and SEDs are available on Rainbow (
Submitted: 2015-08-03, last modified: 2016-05-02
We present a complete census of all Herschel-detected sources within the six massive lensing clusters of the HST Frontier Fields (HFF). We provide a robust legacy catalogue of 263 sources with Herschel fluxes, primarily based on imaging from the Herschel Lensing Survey (HLS) and PEP/HerMES Key Programmes. We optimally combine Herschel, Spitzer and WISE infrared (IR) photometry with data from HST, VLA and ground-based observatories, identifying counterparts to gain source redshifts. For each Herschel-detected source we also present magnification factor (mu), intrinsic IR luminosity and characteristic dust temperature, providing a comprehensive view of dust-obscured star formation within the HFF. We demonstrate the utility of our catalogues through an exploratory overview of the magnified population, including more than 20 background sub-LIRGs unreachable by Herschel without the assistance gravitational lensing.
[27]  [pdf] - 1370441
The data processing pipeline for the Herschel SPIRE Fourier Transform Spectrometer
Comments: 14 pages, 13 figures, MNRAS in press
Submitted: 2016-03-01
We present the data processing pipeline to generate calibrated data products from the Spectral and Photometric Imaging Receiver (SPIRE) imaging Fourier Transform Spectrometer on the Herschel Space Observatory. The pipeline processes telemetry from SPIRE observations and produces calibrated spectra for all resolution modes. The spectrometer pipeline shares some elements with the SPIRE photometer pipeline, including the conversion of telemetry packets into data timelines and calculation of bolometer voltages. We present the following fundamental processing steps unique to the spectrometer: temporal and spatial interpolation of the scan mechanism and detector data to create interferograms; Fourier transformation; apodization; and creation of a data cube. We also describe the corrections for various instrumental effects including first- and second-level glitch identification and removal, correction of the effects due to emission from the Herschel telescope and from within the spectrometer instrument, interferogram baseline correction, temporal and spatial phase correction, non-linear response of the bolometers, and variation of instrument performance across the focal plane arrays. Astronomical calibration is based on combinations of observations of standard astronomical sources and regions of space known to contain minimal emission.
[28]  [pdf] - 1426198
The XXL Survey: I. Scientific motivations - XMM-Newton observing plan - Follow-up observations and simulation programme
Comments: 17 pages, accepted in A&A
Submitted: 2015-12-14
We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.
[29]  [pdf] - 1426197
The XXL Survey. II. The bright cluster sample: catalogue and luminosity function
Comments: Accepted by A&A, 25 pages, 17 figures, 7 tables. The XXL-100-GC Master catalogue is available from the XXL Milan database ( and through the CDS (
Submitted: 2015-12-14
Context. The XXL Survey is the largest survey carried out by the XMM-Newton satellite and covers a total area of 50 square degrees distributed over two fields. It primarily aims at investigating the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Aims. This article presents the XXL bright cluster sample, a subsample of 100 galaxy clusters selected from the full XXL catalogue by setting a lower limit of $3\times 10^{-14}\,\mathrm{erg \,s^{-1}cm^{-2}}$ on the source flux within a 1$^{\prime}$ aperture. Methods. The selection function was estimated using a mixture of Monte Carlo simulations and analytical recipes that closely reproduce the source selection process. An extensive spectroscopic follow-up provided redshifts for 97 of the 100 clusters. We derived accurate X-ray parameters for all the sources. Scaling relations were self-consistently derived from the same sample in other publications of the series. On this basis, we study the number density, luminosity function, and spatial distribution of the sample. Results. The bright cluster sample consists of systems with masses between $M_{500}=7\times 10^{13}$ and $3\times 10^{14} M_\odot$, mostly located between $z=0.1$ and 0.5. The observed sky density of clusters is slightly below the predictions from the WMAP9 model, and significantly below the predictions from the Planck 2015 cosmology. In general, within the current uncertainties of the cluster mass calibration, models with higher values of $\sigma_8$ and/or $\Omega_m$ appear more difficult to accommodate. We provide tight constraints on the cluster differential luminosity function and find no hint of evolution out to $z\sim1$. We also find strong evidence for the presence of large-scale structures in the XXL bright cluster sample and identify five new superclusters.
[30]  [pdf] - 1426200
The XXL Survey VII: A supercluster of galaxies at z=0.43
Comments: A&A, accepted; special XXL issue
Submitted: 2015-12-14
The XXL Survey is the largest homogeneous and contiguous survey carried out with XMM-Newton. Covering an area of 50 square degrees distributed over two fields, it primarily investigates the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Given its depth and sky coverage, XXL is particularly suited to systematically unveiling the clustering of X-ray clusters and to identifying superstructures in a homogeneous X-ray sample down to the typical mass scale of a local massive cluster. A friends-of-friends algorithm in three-dimensional physical space was run to identify large-scale structures. In this paper we report the discovery of the highest redshift supercluster of galaxies found in the XXL Survey. We describe the X-ray properties of the clusters members of the structure and the optical follow-up. The newly discovered supercluster is composed of six clusters of galaxies at a median redshift z around 0.43 and distributed across approximately 30 by 15 arc minutes (10 by 5 Mpc on sky) on the sky. This structure is very compact with all the clusters residing in one XMM pointing; for this reason this is the first supercluster discovered with the XXL Survey. Spectroscopic follow-up with WHT (William Herschel Telescope) and NTT (New Technology Telescope) confirmed a median redshift of z = 0.43. An estimate of the X-ray mass and luminosity of this supercluster and of its total gas mass put XLSSC-e at the average mass range of superclusters; its appearance, with two members of equal size, is quite unusual with respect to other superclusters and provides a unique view of the formation process of a massive structure.
[31]  [pdf] - 1347655
The HerMES sub-millimetre local and low-redshift luminosity functions
Comments: Accepted for publication by MNRAS on 17 Nov 2015
Submitted: 2015-11-19
We used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.
[32]  [pdf] - 1331107
Starbursting Brightest Cluster Galaxy: a Herschel view of the massive cluster MACS J1931.8-2634
Comments: 6 pages, 3 figures, MNRAS accepted
Submitted: 2015-11-10
We investigate the dust-obscured star formation properties of the massive, X-ray selected galaxy cluster MACS J1931.8-2634 at $z$=0.352. Using far-infrared (FIR) imaging in the range 100-500$\mu$m obtained with the \textit{Herschel} telescope, we extract 31 sources (2$\sigma$) within $r\sim$1 Mpc from the brightest cluster galaxy (BCG). Among these sources we identify six cluster members for which we perform an analysis of their spectral energy distributions (SEDs). We measure total infrared luminosity (L$_{IR}$), star formation rate (SFR) and dust temperature. The BCG, with L$_{IR}$=1.4$\times$10$^{12}$L$_\odot$ is an Ultra Luminous Infrared Galaxy and hosts a type II AGN. We decompose its FIR SED into AGN and starburst components and find equal contributions from AGN and starburst. We also recompute the SFR of the BCG finding SFR=150$\pm$15 M$_\odot$yr$^{-1}$. We search for an isobaric cooling flow in the cool core using {\sl Chandra} X-ray data, and find no evidence for gas colder than 1.8 keV in the inner 30 kpc, for an upper limit to the istantaneous mass-deposition rate of 58 M$_\odot$yr$^{-1}$ at 95 % c.l. This value is $3\times$ lower than the SFR in the BCG, suggesting that the on-going SF episode lasts longer than the ICM cooling events.
[33]  [pdf] - 1530353
The Herschel Point Source Catalogue
Comments: 2 pages, 1 figure
Submitted: 2015-10-27
The Herschel Space Observatory was the fourth cornerstone mission in the European Space Agency (ESA) science programme with excellent broad band imaging capabilities in the sub-mm and far-infrared part of the spectrum. Although the spacecraft finished its observations in 2013, it left a large legacy dataset that is far from having been fully scrutinised and still has a large potential for new scientific discoveries. This is specifically true for the photometric observations of the PACS and SPIRE instruments. Some source catalogues have already been produced by individual observing programs, but there are many observations that risk to remain unexplored. To maximise the science return of the SPIRE and PACS data sets, we are in the process of building the Herschel Point Source Catalogue (HPSC) from all primary and parallel mode observations. Our homogeneous source extraction enables a systematic and unbiased comparison of sensitivity across the different Herschel fields that single programs will generally not be able to provide. The catalogue will be made available online through archives like the Herschel Science Archive (HSA), the Infrared Science Archive (IRSA), and the Strasbourg Astronomical Data Center (CDS).
[34]  [pdf] - 1265714
HerMES: Current Cosmic Infrared Background Estimates Can be Explained by Known Galaxies and their Faint Companions at z < 4
Comments: Accepted to ApJL. 6 Pages, 3 figures
Submitted: 2015-05-22, last modified: 2015-08-08
We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 {\mu}m from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 {\mu}m at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 {\mu}m, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.
[35]  [pdf] - 1249120
New XMM-Newton observation of the Phoenix cluster: properties of the cool core
Comments: A&A in press, typos corrected, revised text according to published version
Submitted: 2015-05-26, last modified: 2015-07-01
(Abridged) We present a spectral analysis of a deep (220 ks) XMM-Newton observation of the Phoenix cluster (SPT-CL J2344-4243), which we also combine with Chandra archival ACIS-I data. We extract CCD and RGS X-ray spectra from the core region to search for the signature of cold gas, and constrain the mass deposition rate in the cooling flow which is thought to be responsible of the massive star formation episode observed in the BCG. We find an average mass deposition rate of $\dot M = 620 (-190 +200)_{stat} (-50 +150)_{syst} M_\odot$/yr in the temperature range 0.3-3.0 keV from MOS data. A temperature-resolved analysis shows that a significant amount of gas is deposited only above 1.8 keV, while upper limits of the order of hundreds of $M_\odot$/yr can be put in the 0.3-1.8 keV temperature range. From pn data we obtain $\dot M = 210 (-80 +85)_{stat} ( -35 +60)_{syst} M_\odot$/yr, and the upper limits from the temperature-resolved analysis are typically a factor of 3 lower than MOS data. In the RGS spectrum, no line emission from ionization states below Fe XXIII is seen above $12 \AA$, and the amount of gas cooling below $\sim 3$ keV has a best-fit value $\dot M = 122_{-122}^{+343}$ $M_{\odot}$/yr. In addition, our analysis of the FIR SED of the BCG based on Herschel data provides $SFR = (530 \pm 50) M_\odot$/yr, significantly lower than previous estimates by a factor 1.5. Current data are able to firmly identify substantial amount of cooling gas only above 1.8 keV in the core of the Phoenix cluster. While MOS data analysis is consistent with values as high as $\dot M \sim 1000$ within $1 \sigma$, pn data provide $\dot M < 500 M_\odot$ yr$^{-1}$ at $3\sigma$ c.l. at temperature below 1.8 keV. At present, this discrepancy cannot be explained on the basis of known calibration uncertainties or other sources of statistical noise.
[36]  [pdf] - 1277049
Planck's Dusty GEMS: Gravitationally lensed high-redshift galaxies discovered with the Planck survey
Comments: A&A accepted
Submitted: 2015-06-05
We present an analysis of 11 bright far-IR/submm sources discovered through a combination of the Planck survey and follow-up Herschel-SPIRE imaging. Each source has a redshift z=2.2-3.6 obtained through a blind redshift search with EMIR at the IRAM 30-m telescope. Interferometry obtained at IRAM and the SMA, and optical/near-infrared imaging obtained at the CFHT and the VLT reveal morphologies consistent with strongly gravitationally lensed sources. Additional photometry was obtained with JCMT/SCUBA-2 and IRAM/GISMO at 850 um and 2 mm, respectively. All objects are bright, isolated point sources in the 18 arcsec beam of SPIRE at 250 um, with spectral energy distributions peaking either near the 350 um or the 500 um bands of SPIRE, and with apparent far-infrared luminosities of up to 3x10^14 L_sun. Their morphologies and sizes, CO line widths and luminosities, dust temperatures, and far-infrared luminosities provide additional empirical evidence that these are strongly gravitationally lensed high-redshift galaxies. We discuss their dust masses and temperatures, and use additional WISE 22-um photometry and template fitting to rule out a significant contribution of AGN heating to the total infrared luminosity. Six sources are detected in FIRST at 1.4 GHz. Four have flux densities brighter than expected from the local far-infrared-radio correlation, but in the range previously found for high-z submm galaxies, one has a deficit of FIR emission, and 6 are consistent with the local correlation. The global dust-to-gas ratios and star-formation efficiencies of our sources are predominantly in the range expected from massive, metal-rich, intense, high-redshift starbursts. An extensive multi-wavelength follow-up programme is being carried out to further characterize these sources and the intense star-formation within them.
[37]  [pdf] - 1288345
Planck intermediate results. XXVII. High-redshift infrared galaxy overdensity candidates and lensed sources discovered by Planck and confirmed by Herschel-SPIRE
Planck Collaboration; Aghanim, N.; Altieri, B.; Arnaud, M.; Ashdown, M.; Aumont, J.; Baccigalupi, C.; Banday, A. J.; Barreiro, R. B.; Bartolo, N.; Battaner, E.; Beelen, A.; Benabed, K.; Benoit-Lévy, A.; Bernard, J. -P.; Bersanelli, M.; Bethermin, M.; Bielewicz, P.; Bonavera, L.; Bond, J. R.; Borrill, J.; Bouchet, F. R.; Boulanger, F.; Burigana, C.; Calabrese, E.; Canameras, R.; Cardoso, J. -F.; Catalano, A.; Chamballu, A.; Chary, R. -R.; Chiang, H. C.; Christensen, P. R.; Clements, D. L.; Colombi, S.; Couchot, F.; Crill, B. P.; Curto, A.; Danese, L.; Dassas, K.; Davies, R. D.; Davis, R. J.; de Bernardis, P.; de Rosa, A.; de Zotti, G.; Delabrouille, J.; Diego, J. M.; Dole, H.; Donzelli, S.; Doré, O.; Douspis, M.; Ducout, A.; Dupac, X.; Efstathiou, G.; Elsner, F.; Enßlin, T. A.; Falgarone, E.; Flores-Cacho, I.; Forni, O.; Frailis, M.; Fraisse, A. A.; Franceschi, E.; Frejsel, A.; Frye, B.; Galeotta, S.; Galli, S.; Ganga, K.; Giard, M.; Gjerløw, E.; González-Nuevo, J.; Górski, K. M.; Gregorio, A.; Gruppuso, A.; Guéry, D.; Hansen, F. K.; Hanson, D.; Harrison, D. L.; Helou, G.; Hernández-Monteagudo, C.; Hildebrandt, S. R.; Hivon, E.; Hobson, M.; Holmes, W. A.; Hovest, W.; Huffenberger, K. M.; Hurier, G.; Jaffe, A. H.; Jaffe, T. R.; Keihänen, E.; Keskitalo, R.; Kisner, T. S.; Kneissl, R.; Knoche, J.; Kunz, M.; Kurki-Suonio, H.; Lagache, G.; Lamarre, J. -M.; Lasenby, A.; Lattanzi, M.; Lawrence, C. R.; Floc'h, E. Le; Leonardi, R.; Levrier, F.; Liguori, M.; Lilje, P. B.; Linden-Vørnle, M.; López-Caniego, M.; Lubin, P. M.; Macías-Pérez, J. F.; MacKenzie, T.; Maffei, B.; Mandolesi, N.; Maris, M.; Martin, P. G.; Martinache, C.; Martínez-González, E.; Masi, S.; Matarrese, S.; Mazzotta, P.; Melchiorri, A.; Mennella, A.; Migliaccio, M.; Moneti, A.; Montier, L.; Morgante, G.; Mortlock, D.; Munshi, D.; Murphy, J. A.; Natoli, P.; Negrello, M.; Nesvadba, N. P. H.; Novikov, D.; Novikov, I.; Omont, A.; Pagano, L.; Pajot, F.; Pasian, F.; Perdereau, O.; Perotto, L.; Perrotta, F.; Pettorino, V.; Piacentini, F.; Piat, M.; Plaszczynski, S.; Pointecouteau, E.; Polenta, G.; Popa, L.; Pratt, G. W.; Prunet, S.; Puget, J. -L.; Rachen, J. P.; Reach, W. T.; Reinecke, M.; Remazeilles, M.; Renault, C.; Ristorcelli, I.; Rocha, G.; Roudier, G.; Rusholme, B.; Sandri, M.; Santos, D.; Savini, G.; Scott, D.; Spencer, L. D.; Stolyarov, V.; Sunyaev, R.; Sutton, D.; Sygnet, J. -F.; Tauber, J. A.; Terenzi, L.; Toffolatti, L.; Tomasi, M.; Tristram, M.; Tucci, M.; Umana, G.; Valenziano, L.; Valiviita, J.; Valtchanov, I.; Van Tent, B.; Vieira, J. D.; Vielva, P.; Wade, L. A.; Wandelt, B. D.; Wehus, I. K.; Welikala, N.; Zacchei, A.; Zonca, A.
Comments: Updated author list, in press and available at A&A ; Full-res paper (3.6Mb) at A&A ; Check also web story at ESA and NASA
Submitted: 2015-03-30, last modified: 2015-04-08
[Abridged] We use the Planck all-sky submm and mm maps to search for rare sources distinguished by extreme brightness, a few hundreds of mJy, and their potential for being situated at high redshift. These "cold" Planck sources, selected using the High Frequency Instrument (HFI) directly from the maps and from the Planck Catalogue of Compact Sources (PCCS), all satisfy the criterion of having their rest-frame far-infrared peak redshifted to the frequency range 353 and 857 GHz. This colour-selection favours galaxies in the redshift range z=2-4, which we consider as cold peaks in the cosmic infrared background (CIB). We perform a dedicated Herschel-SPIRE follow-up of 234 such Planck targets, finding a significant excess of red 350 and 500um sources, in comparison to reference SPIRE fields. About 94% of the SPIRE sources in the Planck fields are consistent with being overdensities of galaxies peaking at 350um. About 3% are candidate lensed systems, all 12 of which have secure spectroscopic confirmations, placing them at redshifts z>2.2. The galaxy overdensities are detected with high significance, half of the sample showing statistical significance above 10sigma. The SPIRE photometric redshifts of galaxies in overdensities suggest a peak at z~2. Under the Td=35K assumption, we derive an infrared (IR) luminosity for each SPIRE source of about 4x10^12 Lsun, yielding star formation rates of typically 700 Msun.yr^-1. If the observed overdensities are actual gravitationally-bound structures, the total total star formation rates reaches 7x10^3 Msun.yr^-1. Taken together, these sources show the signatures of high-z (z>$) protoclusters of intensively star-forming galaxies. All these observations confirm the uniqueness of our sample and demonstrate the ability of the all-sky Planck-HFI cold sources to select populations of cosmological and astrophysical interest for structure formation studies.
[38]  [pdf] - 1242234
GOODS-$Herschel$: identification of the individual galaxies responsible for the 80-290$\mu$m cosmic infrared background
Comments: 22 pages, 16 figures. Accepted for publication by Astronomy and Astrophysics
Submitted: 2015-03-19
We propose a new method of pushing $Herschel$ to its faintest detection limits using universal trends in the redshift evolution of the far infrared over 24$\mu$m colours in the well-sampled GOODS-North field. An extension to other fields with less multi-wavelength information is presented. This method is applied here to raise the contribution of individually detected $Herschel$ sources to the cosmic infrared background (CIRB) by a factor 5 close to its peak at 250$\mu$m and more than 3 in the 350$\mu$m and 500$\mu$m bands. We produce realistic mock $Herschel$ images of the deep PACS and SPIRE images of the GOODS-North field from the GOODS-$Herschel$ Key Program and use them to quantify the confusion noise at the position of individual sources, i.e., estimate a "local confusion noise". Two methods are used to identify sources with reliable photometric accuracy extracted using 24$\mu$m prior positions. The clean index (CI), previously defined but validated here with simulations, which measures the presence of bright 24$\mu$m neighbours and the photometric accuracy index (PAI) directly extracted from the mock $Herschel$ images. After correction for completeness, thanks to our mock $Herschel$ images, individually detected sources make up as much as 54% and 60% of the CIRB in the PACS bands down to 1.1 mJy at 100$\mu$m and 2.2 mJy at 160$\mu$m and 55, 33, and 13% of the CIRB in the SPIRE bands down to 2.5, 5, and 9 mJy at 250$\mu$m, 350$\mu$m, and 500$\mu$m, respectively. The latter depths improve the detection limits of $Herschel$ by factors of 5 at 250$\mu$m, and 3 at 350$\mu$m and 500$\mu$m as compared to the standard confusion limit. Interestingly, the dominant contributors to the CIRB in all $Herschel$ bands appear to be distant siblings of the Milky Way ($z$$\sim$0.96 for $\lambda$$<$300$\mu$m) with a stellar mass of $M_{\star}$$\sim$9$\times$10$^{10}$M$_{\odot}$.
[39]  [pdf] - 943876
A multi-wavelength exploration of the [CII]/IR ratio in H-ATLAS/GAMA galaxies out to z=0.2
Comments: 16 pages, 10 figures, 3 tables
Submitted: 2015-03-03
We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).
[40]  [pdf] - 1224510
Systematic characterisation of the Herschel SPIRE Fourier Transform Spectrometer
Comments: 41 pages, 37 figures, 32 tables. Accepted for publication in MNRAS
Submitted: 2015-02-19
A systematic programme of calibration observations was carried out to monitor the performance of the SPIRE FTS instrument on board the Herschel Space Observatory. Observations of planets (including the prime point-source calibrator, Uranus), asteroids, line sources, dark sky, and cross-calibration sources were made in order to monitor repeatability and sensitivity, and to improve FTS calibration. We present a complete analysis of the full set of calibration observations and use them to assess the performance of the FTS. Particular care is taken to understand and separate out the effect of pointing uncertainties, including the position of the internal beam steering mirror for sparse observations in the early part of the mission. The repeatability of spectral line centre positions is <5km/s, for lines with signal-to-noise ratios >40, corresponding to <0.5-2.0% of a resolution element. For spectral line flux, the repeatability is better than 6%, which improves to 1-2% for spectra corrected for pointing offsets. The continuum repeatability is 4.4% for the SLW band and 13.6% for the SSW band, which reduces to ~1% once the data have been corrected for pointing offsets. Observations of dark sky were used to assess the sensitivity and the systematic offset in the continuum, both of which were found to be consistent across the FTS detector arrays. The average point-source calibrated sensitivity for the centre detectors is 0.20 and 0.21 Jy [1 sigma; 1 hour], for SLW and SSW. The average continuum offset is 0.40 Jy for the SLW band and 0.28 Jy for the SSW band.
[41]  [pdf] - 1224456
The nuclear and extended infrared emission of the Seyfert galaxy NGC 2992 and the interacting system Arp 245
Comments: 20 pages, 12 figures, accepted by MNRAS
Submitted: 2015-02-16
We present subarcsecond resolution infrared (IR) imaging and mid-IR spectroscopic observations of the Seyfert 1.9 galaxy NGC 2992, obtained with the Gemini North Telescope and the Gran Telescopio CANARIAS (GTC). The N-band image reveals faint extended emission out to ~3 kpc, and the PAH features detected in the GTC/CanariCam 7.5-13 micron spectrum indicate that the bulk of this extended emission is dust heated by star formation. We also report arcsecond resolution MIR and far-IR imaging of the interacting system Arp 245, taken with the Spitzer Space Telescope and the Herschel Space Observatory. Using these data, we obtain nuclear fluxes using different methods and find that we can only recover the nuclear fluxes obtained from the subarcsecond data at 20-25 micron, where the AGN emission dominates. We fitted the nuclear IR spectral energy distribution of NGC 2992, including the GTC/CanariCam nuclear spectrum (~50 pc), with clumpy torus models. We then used the best-fitting torus model to decompose the Spitzer/IRS 5-30 spectrum (~630 pc) in AGN and starburst components, using different starburst templates. We find that, whereas at shorter mid-IR wavelengths the starburst component dominates (64% at 6 micron), the AGN component reaches 90% at 20 micron. We finally obtained dust masses, temperatures and star formation rates for the different components of the Arp 245 system and find similar values for NGC 2992 and NGC 2993. These measurements are within those reported for other interacting systems in the first stages of the interaction.
[42]  [pdf] - 1055782
GLACE survey: OSIRIS/GTC Tuneable Filter H$\alpha$ imaging of the rich galaxy cluster ZwCl 0024.0+1652 at z = 0.395. Part I -- Survey presentation, TF data reduction techniques and catalogue
Comments: 23 pages, 13 figures, accepted for publication in A&A. Astro-ph abstract abridged
Submitted: 2015-02-10
The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.
[43]  [pdf] - 1223676
The reversal of the SF-density relation in a massive, X-ray selected galaxy cluster at z=1.58: results from Herschel
Comments: Letter accepted for publication in MNRAS, ESA Press Release on 18 December 2014
Submitted: 2014-12-16
Dusty, star-forming galaxies have a critical role in the formation and evolution of massive galaxies in the Universe. Using deep far-infrared imaging in the range 100-500um obtained with the Herschel telescope, we investigate the dust-obscured star formation in the galaxy cluster XDCP J0044.0-2033 at z=1.58, the most massive cluster at z >1.5, with a measured mass M200= 4.7x10$^{14}$ Msun. We perform an analysis of the spectral energy distributions (SEDs) of 12 cluster members (5 spectroscopically confirmed) detected with >3$\sigma$ significance in the PACS maps, all ULIRGs. The individual star formation rates (SFRs) lie in the range 155-824 Ms/yr, with dust temperatures of 24$\pm$35 K. We measure a strikingly high amount of star formation (SF) in the cluster core, SFR (< 250 kpc) > 1875$\pm$158 Ms/yr, 4x higher than the amount of star formation in the cluster outskirts. This scenario is unprecedented in a galaxy cluster, showing for the first time a reversal of the SF-density relation at z~1.6 in a massive cluster.
[44]  [pdf] - 908619
The far-infrared/radio correlation and radio spectral index of galaxies in the SFR-M* plane up to z 2
Comments: Accepted for publication in A&A; 18 pages, 10 figures
Submitted: 2014-10-27
[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.
[45]  [pdf] - 882044
An excess of dusty starbursts related to the Spiderweb galaxy
Comments: Accepted by A&A. ESO Press Release on 15 October 2014. 19 pages, 9 figures and 4 tables
Submitted: 2014-10-14
We present APEX LABOCA 870 micron observations of the field around the high-redshift radio galaxy MRC1138-262 at z=2.16. We detect 16 submillimeter galaxies in this ~140 square arcmin bolometer map with flux densities in the range 3-11 mJy. The raw number counts indicate a density of submillimeter galaxies (SMGs) that is up to four times that of blank field surveys. Based on an exquisite multiwavelength database, including VLA 1.4 GHz radio and infrared observations, we investigate whether these sources are members of the protocluster structure at z=2.2. Using Herschel PACS+SPIRE and Spitzer MIPS photometry, we derive reliable far-infrared photometric redshifts for all sources. Follow-up VLT ISAAC and SINFONI near-infrared spectra confirm that four of these SMGs have redshifts of z=2.2. We also present evidence that another SMG in this field, detected earlier at 850 micron, has a counterpart that exhibits Halpha and CO(1-0) emission at z=2.15. Including the radio galaxy and two SMGs with far-IR photometric redshifts at z=2.2, we conclude that at least eight submm sources are part of the protocluster at z=2.16 associated with the radio galaxy MRC1138-262. We measure a star formation rate density SFRD ~1500 Msun yr^-1 Mpc^-3, four magnitudes higher than the global SFRD of blank fields at this redshift. Strikingly, these eight sources are concentrated within a region of 2 Mpc (the typical size of clusters in the local universe) and are distributed within the filaments traced by the Halpha emitters at z=2.2. This concentration of massive, dusty starbursts is not centered on the submillimeter-bright radio galaxy which could support the infalling of these sources into the cluster center. Approximately half (6/11) of the SMGs that are covered by the Halpha imaging data are associated with Halpha emitters, demonstrating the potential of tracing SMG counterparts with this population (abridged).
[46]  [pdf] - 1217058
A Far-Infrared Spectroscopic Survey of Intermediate Redshift (Ultra) Luminous Infrared Galaxies
Comments: Accepted for publication in ApJ
Submitted: 2014-09-19
We present Herschel far-IR photometry and spectroscopy as well as ground based CO observations of an intermediate redshift (0.21 < z < 0.88) sample of Herschel-selected (ultra)-luminous infrared galaxies (L_IR > 10^11.5L_sun). With these measurements we trace the dust continuum, far-IR atomic line emission, in particular [CII]\,157.7microns, as well as the molecular gas of z~0.3 (U)LIRGs and perform a detailed investigation of the interstellar medium of the population. We find that the majority of Herschel-selected intermediate redshift (U)LIRGs have L_CII/L_FIR ratios that are a factor of about 10 higher than that of local ULIRGs and comparable to that of local normal and high-$z$ star forming galaxies. Using our sample to bridge local and high-z [CII] observations, we find that the majority of galaxies at all redshifts and all luminosities follow a L_CII-L_FIR relation with a slope of unity, from which local ULIRGs and high-z AGN dominated sources are clear outliers. We also confirm that the strong anti-correlation between the L_CII/L_FIR ratio and the far-IR color L_60/L_100 observed in the local Universe holds over a broad range of redshifts and luminosities, in the sense that warmer sources exhibit lower L_CII/L_FIR at any epoch. Intermediate redshift ULIRGs are also characterised by large molecular gas reservoirs and by lower star formation efficiencies compared to that of local ULIRGs. The high L_CII/L_FIR ratios, the moderate star formation efficiencies (L_LIR/L_CO or L_IR/M_gas) and the relatively low dust temperatures of our sample (which are also common characteristics of high-z star forming galaxies with ULIRG-like luminosities) indicate that the evolution of the physical properties of (U)LIRGs between the present day and z > 1 is already significant by z ~ 0.3.
[47]  [pdf] - 1253849
Detailed modelling of a large sample of Herschel sources in the Lockman Hole: identification of cold dust and of lensing candidates through their anomalous SEDs
Comments: 18 pages, 46 figures, accepted for publication MNRAS. arXiv admin note: text overlap with arXiv:1006.0120
Submitted: 2014-09-17
We have studied in detail a sample of 967 SPIRE sources with 5-sigma detections at 350 and 500 micron and associations with Spitzer-SWIRE 24 micron galaxies in the HerMES-Lockman survey area, fitting their mid- and far-infrared, and submillimetre, SEDs in an automatic search with a set of six infrared templates. For almost 300 galaxies we have modelled their SEDs individually to ensure the physicality of the fits. We confirm the need for the new cool and cold cirrus templates, and also of the young starburst template, introduced in earlier work. We also identify 109 lensing candidates via their anomalous SEDs and provide a set of colour-redshift constraints which allow lensing candidates to be identified from combined Herschel and Spitzer data. The picture that emerges of the submillimetre galaxy population is complex, comprising ultraluminous and hyperluminous starbursts, lower luminosity galaxies dominated by interstellar dust emission, lensed galaxies and galaxies with surprisingly cold (10-13K) dust. 11 % of 500$\mu$m selected sources are lensing candidates. 70 % of the unlensed sources are ultraluminous infrared galaxies and 26 % are hyperluminous. 34 % are dominated by optically thin interstellar dust ('cirrus') emission, but most of these are due to cooler dust than is characteristic of our Galaxy. At the highest infrared luminosities we see SEDs dominated by M82, Arp220 and young starburst types, in roughly equal proportions.
[48]  [pdf] - 862347
Flux Calibration of Broadband Far Infrared and Submillimetre Photometric Instruments: Theory and Application to Herschel-SPIRE
Comments: 14 pages, 16 figures, published in MNRAS (2013). Corrections to Eq.28 & Fig11 to match published version. Some numbers have very small adjustments (relative to published version) to agree with v12.3 of the SPIRE Calibration Tree. Further information about SPIRE available at calibration at
Submitted: 2013-06-07, last modified: 2014-08-11
Photometric instruments operating at far infrared to millimetre wavelengths often have broad spectral passbands (central wavelength/bandwidth ~ 3 or less), especially those operating in space. A broad passband can result in significant variation of the beam profile and aperture efficiency across the passband, effects which thus far have not generally been taken into account in the flux calibration of such instruments. With absolute calibration uncertainties associated with the brightness of primary calibration standards now in the region of 5% or less, variation of the beam properties across the passband can be a significant contributor to the overall calibration accuracy for extended emission. We present a calibration framework which takes such variations into account for both antenna-coupled and absorber-coupled focal plane architectures. The scheme covers point source and extended source cases, and also the intermediate case of a semi-extended source profile. We apply the new method to the Herschel-SPIRE space-borne photometer.
[49]  [pdf] - 1203557
Herschel reveals a molecular outflow in a z = 2.3 ULIRG
Comments: 7 pages, 4 figures; accepted for publication in MNRAS
Submitted: 2014-02-25, last modified: 2014-06-04
We report the results from a 19-h integration with the SPIRE Fourier Transform Spectrometer aboard the Herschel Space Observatory which has revealed the presence of a molecular outflow from the Cosmic Eyelash (SMM J2135-0102) via the detection of blueshifted OH absorption. Detections of several fine-structure emission lines indicate low-excitation HII regions contribute strongly to the [CII] luminosity in this z = 2.3 ULIRG. The OH feature suggests a maximum wind velocity of 700 km/s, which is lower than the expected escape velocity of the host dark matter halo, ~1000 km/s. A large fraction of the available molecular gas could thus be converted into stars via a burst protracted by the resulting gas fountain, until an AGN-driven outflow can eject the remaining gas.
[50]  [pdf] - 1209361
Star formation in the massive cluster merger Abell 2744
Comments: Accepted by MNRAS. 12 pages, 7 figures (high resolution versions of Figs. 1 & 2 are available in the published PDF)
Submitted: 2014-05-05
We present a comprehensive study of star-forming (SF) galaxies in the HST Frontier Field recent cluster merger A2744 (z=0.308). Wide-field, ultraviolet-infrared (UV-IR) imaging enables a direct constraint of the total star formation rate (SFR) for 53 cluster galaxies, with SFR{UV+IR}=343+/-10 Msun/yr. Within the central 4 arcmin (1.1 Mpc) radius, the integrated SFR is complete, yielding a total SFR{UV+IR}=201+/-9 Msun/yr. Focussing on obscured star formation, this core region exhibits a total SFR{IR}=138+/-8 Msun/yr, a mass-normalised SFR{IR} of Sigma{SFR}=11.2+/-0.7 Msun/yr per 10^14 Msun and a fraction of IR-detected SF galaxies f{SF}=0.080(+0.010,-0.037). Overall, the cluster population at z~0.3 exhibits significant intrinsic scatter in IR properties (total SFR{IR}, Tdust distribution) apparently unrelated to the dynamical state: A2744 is noticeably different to the merging Bullet cluster, but similar to several relaxed clusters. However, in A2744 we identify a trail of SF sources including jellyfish galaxies with substantial unobscured SF due to extreme stripping (SFR{UV}/SFR{IR} up to 3.3). The orientation of the trail, and of material stripped from constituent galaxies, indicates that the passing shock front of the cluster merger was the trigger. Constraints on star formation from both IR and UV are crucial for understanding galaxy evolution within the densest environments.
[51]  [pdf] - 1208191
Calibration of the Herschel SPIRE Fourier Transform Spectrometer
Comments: 20 pages, 18 figures, accepted for publication in MNRAS
Submitted: 2014-03-05
The Herschel SPIRE instrument consists of an imaging photometric camera and an imaging Fourier Transform Spectrometer (FTS), both operating over a frequency range of 450-1550 GHz. In this paper, we briefly review the FTS design, operation, and data reduction, and describe in detail the approach taken to relative calibration (removal of instrument signatures) and absolute calibration against standard astronomical sources. The calibration scheme assumes a spatially extended source and uses the Herschel telescope as primary calibrator. Conversion from extended to point-source calibration is carried out using observations of the planet Uranus. The model of the telescope emission is shown to be accurate to within 6% and repeatable to better than 0.06% and, by comparison with models of Mars and Neptune, the Uranus model is shown to be accurate to within 3%. Multiple observations of a number of point-like sources show that the repeatability of the calibration is better than 1%, if the effects of the satellite absolute pointing error (APE) are corrected. The satellite APE leads to a decrement in the derived flux, which can be up to ~10% (1 sigma) at the high-frequency end of the SPIRE range in the first part of the mission, and ~4% after Herschel operational day 1011. The lower frequency range of the SPIRE band is unaffected by this pointing error due to the larger beam size. Overall, for well-pointed, point-like sources, the absolute flux calibration is better than 6%, and for extended sources where mapping is required it is better than 7%.
[52]  [pdf] - 1179965
[CII] and 12CO(1-0) Emission Maps in HLSJ091828.6+514223: A Strongly Lensed Interacting System at z=5.24
Comments: 28 pages, 9 figures, accepted in ApJ
Submitted: 2013-10-15, last modified: 2014-01-15
We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.
[53]  [pdf] - 1202667
Herschel Observations of Far-Infrared Cooling Lines in intermediate Redshift (Ultra)-luminous Infrared Galaxies
Comments: Accepted for publication in ApJL, ApJ Volumne 781, L15
Submitted: 2014-01-10
We report the first results from a spectroscopic survey of the [CII] 158um line from a sample of intermediate redshift (0.2<z<0.8) (ultra)-luminous infrared galaxies, (U)LIRGs (LIR>10^11.5 Lsun), using the SPIRE-Fourier Transform Spectrometer (FTS) on board the Herschel Space Observatory. This is the first survey of [CII] emission, an important tracer of star-formation, at a redshift range where the star-formation rate density of the Universe increases rapidly. We detect strong [CII] 158um line emission from over 80% of the sample. We find that the [CII] line is luminous, in the range (0.8-4)x10^(-3) of the far-infrared continuum luminosity of our sources, and appears to arise from photodissociation regions on the surface of molecular clouds. The L[CII]/LIR ratio in our intermediate redshift (U)LIRGs is on average ~10 times larger than that of local ULIRGs. Furthermore, we find that the L[CII]/LIR and L[CII]/LCO(1-0) ratios in our sample are similar to those of local normal galaxies and high-z star-forming galaxies. ULIRGs at z~0.5 show many similarities to the properties of local normal and high-z star forming galaxies. Our findings strongly suggest that rapid evolution in the properties of the star forming regions of luminous infrared galaxies is likely to have occurred in the last 5 billion years.
[54]  [pdf] - 769842
Relative pointing offset analysis of calibration targets with repeated observations with Herschel-SPIRE Fourier-Transform Spectrometer
Comments: 17 pages, 19 figures, accepted for publication in Experimental Astronomy
Submitted: 2014-01-09
We present a method to derive the relative pointing offsets for SPIRE Fourier-Transform Spectrometer (FTS) solar system object (SSO) calibration targets, which were observed regularly throughout the Herschel mission. We construct ratios of the spectra for all observations of a given source with respect to a reference. The reference observation is selected iteratively to be the one with the highest observed continuum. Assuming that any pointing offset leads to an overall shift of the continuum level, then these ratios represent the relative flux loss due to mispointing. The mispointing effects are more pronounced for a smaller beam, so we consider only the FTS short wavelength array (SSW, 958-1546 GHz) to derive a pointing correction. We obtain the relative pointing offset by comparing the ratio to a grid of expected losses for a model source at different distances from the centre of the beam, under the assumption that the SSW FTS beam can be well approximated by a Gaussian. In order to avoid dependency on the point source flux conversion, which uses a particular observation of Uranus, we use extended source flux calibrated spectra to construct the ratios for the SSOs. In order to account for continuum variability, due to the changing distance from the Herschel telescope, the SSO ratios are normalised by the expected model ratios for the corresponding observing epoch. We confirm the accuracy of the derived pointing offset by comparing the results with a number of control observations, where the actual pointing of Herschel is known with good precision. Using the method we derived pointing offsets for repeated observations of Uranus (including observations centred on off-axis detectors), Neptune, Ceres and NGC7027. The results are used to validate and improve the point-source flux calibration of the FTS.
[55]  [pdf] - 1202650
Herschel SPIRE FTS Relative Spectral Response Calibration
Comments: 15 pages, 13 figures, accepted for publication in Experimental Astronomy
Submitted: 2014-01-09
Herschel/SPIRE Fourier transform spectrometer (FTS) observations contain emission from both the Herschel Telescope and the SPIRE Instrument itself, both of which are typically orders of magnitude greater than the emission from the astronomical source, and must be removed in order to recover the source spectrum. The effects of the Herschel Telescope and the SPIRE Instrument are removed during data reduction using relative spectral response calibration curves and emission models. We present the evolution of the methods used to derive the relative spectral response calibration curves for the SPIRE FTS. The relationship between the calibration curves and the ultimate sensitivity of calibrated SPIRE FTS data is discussed and the results from the derivation methods are compared. These comparisons show that the latest derivation methods result in calibration curves that impart a factor of between 2 and 100 less noise to the overall error budget, which results in calibrated spectra for individual observations whose noise is reduced by a factor of 2-3, with a gain in the overall spectral sensitivity of 23% and 21% for the two detector bands, respectively.
[56]  [pdf] - 1202649
Herschel SPIRE FTS Spectral Mapping Calibration
Comments: 12 pages, 9 figures, accepted for publication in Experimental Astronomy
Submitted: 2014-01-09
The Herschel SPIRE Fourier transform spectrometer (FTS) performs spectral imaging in the 447-1546 GHz band. It can observe in three spatial sampling modes: sparse mode, with a single pointing on sky, or intermediate or full modes with 1 and 1/2 beam spacing, respectively. In this paper, we investigate the uncertainty and repeatability for fully sampled FTS mapping observations. The repeatability is characterised using nine observations of the Orion Bar. Metrics are derived based on the ratio of the measured intensity in each observation compared to that in the combined spectral cube from all observations. The mean relative deviation is determined to be within 2%, and the pixel-by-pixel scatter is ~7%. The scatter increases towards the edges of the maps. The uncertainty in the frequency scale is also studied, and the spread in the line centre velocity across the maps is found to be ~15 km/s. Other causes of uncertainty are also discussed including the effect of pointing and the additive uncertainty in the continuum.
[57]  [pdf] - 769844
Herschel SPIRE FTS telescope model correction
Comments: 12 pages, 14 figures, accepted for publication in Experimental Astronomy
Submitted: 2014-01-09
Emission from the Herschel telescope is the dominant source of radiation for the majority of SPIRE Fourier transform spectrometer (FTS) observations, despite the exceptionally low emissivity of the primary and secondary mirrors. Accurate modelling and removal of the telescope contribution is, therefore, an important and challenging aspect of FTS calibration and data reduction pipeline. A dust-contaminated telescope model with time invariant mirror emissivity was adopted before the Herschel launch. However, measured FTS spectra show a clear evolution of the telescope contribution over the mission and strong need for a correction to the standard telescope model in order to reduce residual background (of up to 7 Jy) in the final data products. Systematic changes in observations of dark sky, taken over the course of the mission, provide a measure of the evolution between observed telescope emission and the telescope model. These dark sky observations have been used to derive a time dependent correction to the telescope emissivity that reduces the systematic error in the continuum of the final FTS spectra to ~0.35 Jy.
[58]  [pdf] - 769843
Herschel SPIRE Fourier Transform Spectrometer: Calibration of its Bright-source Mode
Comments: 15 pages, 16 figures, accepted for publication in Experimental Astronomy
Submitted: 2014-01-09
The Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) on board the ESA Herschel Space Observatory has two detector setting modes: (a) a nominal mode, which is optimized for observing moderately bright to faint astronomical targets, and (b) a bright-source mode recommended for sources significantly brighter than 500 Jy, within the SPIRE FTS bandwidth of 446.7-1544 GHz (or 194-671 microns in wavelength), which employs a reduced detector responsivity and out-of-phase analog signal amplifier/demodulator. We address in detail the calibration issues unique to the bright-source mode, describe the integration of the bright-mode data processing into the existing pipeline for the nominal mode, and show that the flux calibration accuracy of the bright-source mode is generally within 2% of that of the nominal mode, and that the bright-source mode is 3 to 4 times less sensitive than the nominal mode.
[59]  [pdf] - 1202077
Star formation in the cluster CLG0218.3-0510 at z=1.62 and its large-scale environment: the infrared perspective
Comments: accepted for publication in MNRAS. v2: minor corrections
Submitted: 2013-12-05, last modified: 2013-12-11
The galaxy cluster CLG0218.3-0510 at z=1.62 is one of the most distant galaxy clusters known, with a rich muti-wavelength data set that confirms a mature galaxy population already in place. Using very deep, wide area (20x20 Mpc) imaging by Spitzer/MIPS at 24um, in conjunction with Herschel 5-band imaging from 100-500um, we investigate the dust-obscured, star-formation properties in the cluster and its associated large scale environment. Our galaxy sample of 693 galaxies at z=1.62 detected at 24um (10 spectroscopic and 683 photo-z) includes both cluster galaxies (i.e. within r <1 Mpc projected clustercentric radius) and field galaxies, defined as the region beyond a radius of 3 Mpc. The star-formation rates (SFRs) derived from the measured infrared luminosity range from 18 to 2500 Ms/yr, with a median of 55 Ms/yr, over the entire radial range (10 Mpc). The cluster brightest FIR galaxy, taken as the centre of the galaxy system, is vigorously forming stars at a rate of 256$\pm$70 Ms/yr, and the total cluster SFR enclosed in a circle of 1 Mpc is 1161$\pm$96 Ms/yr. We estimate a dust extinction of about 3 magnitudes by comparing the SFRs derived from [OII] luminosity with the ones computed from the 24um fluxes. We find that the in-falling region (1-3 Mpc) is special: there is a significant decrement (3.5x) of passive relative to star-forming galaxies in this region, and the total SFR of the galaxies located in this region is lower (130 Ms/yr/Mpc2) than anywhere in the cluster or field, regardless of their stellar mass. In a complementary approach we compute the local galaxy density, Sigma5, and find no trend between SFR and Sigma5. However, we measure an excess of star-forming galaxies in the cluster relative to the field by a factor 1.7, that lends support to a reversal of the SF-density relation in CLG0218.
[60]  [pdf] - 1202018
Disentangling a group of lensed submm galaxies at z~2.9
Submitted: 2013-12-03
MS$\,$0451.6$-$0305 is a rich galaxy cluster whose strong lensing is particularly prominent at submm wavelengths. We combine new SCUBA-2 data with imaging from Herschel SPIRE and PACS and HST in order to try to understand the nature of the sources being lensed. In the region of the "giant submm arc," we uncover seven multiply imaged galaxies (up from the previously known three), of which six are found to be at a redshift of $z\sim2.9$, and possibly constitute an interacting system. Using a novel forward-modelling approach, we are able to simultaneously deblend and fit SEDs to the individual galaxies that contribute to the giant submm arc, constraining their dust temperatures, far infrared luminosities and star formation rates. The submm arc first identified by SCUBA can now be seen to be composed of at least five distinct sources, four of these within the galaxy group at $z\sim2.9$. The total unlensed luminosity for this galaxy group is $(3.1\pm0.3) \times 10^{12}\,\mathrm{L}_\odot$, which gives an unlensed star formation rate of $(450\pm50)$ M$_\odot$ yr$^{-1}$. From the properties of this system, we see no evidence of evolution towards lower temperatures in the dust temperature versus far-infrared luminosity relation for high redshift galaxies.
[61]  [pdf] - 1180673
The evolution of the dust temperatures of galaxies in the SFR$-M_{\ast}$ plane up to $z$$\,\thicksim\,$$2$
Comments: Accepted for publication in A&A; A simple IDL code implementing our results is available at ; This code predicts the dust temperature and the FIR/mm flux densities of a galaxy from its SFR, M* and z
Submitted: 2013-11-12
[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.
[62]  [pdf] - 1180301
HerMES: Candidate High-Redshift Galaxies Discovered with Herschel/SPIRE
Comments: Accepted for publication in ApJ
Submitted: 2013-10-28
We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.
[63]  [pdf] - 732631
The suppression of star formation by powerful active galactic nuclei
Comments: Published in Nature, May 2012. Includes supplementary information
Submitted: 2013-10-15
The old, red stars which 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 from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.
[64]  [pdf] - 1178901
An extended Herschel drop-out source in the center of AS1063, a 'normal' dusty galaxy at z=6.1 or SZ substructures?
Comments: accepted for publication in A&A
Submitted: 2013-08-30, last modified: 2013-10-03
In the course of our 870um APEX/LABOCA follow up of the Herschel Lensing Survey we have detected a source in AS1063 (RXC J2248.7-4431), that has no counterparts in any of the Herschel PACS/SPIRE bands, it is a Herschel 'drop-out' with S_870/S_500>0.5. The 870um emission is extended and centered on the brightest cluster galaxy suggesting either a multiply imaged background source or substructure in the Sunyaev-Zel'dovich (SZ) increment due to inhomogeneities in the hot cluster gas of this merging cluster. We discuss both interpretations with emphasis on the putative lensed source. Based on the observed properties and on our lens model we find that this source could be the first SMG with a moderate far infrared luminosity (L_FIR<10^12 L_sol) detected so far at z>4. In deep HST observations we identified a multiply imaged z~6 source and we measured its spectroscopic redshift z=6.107 with VLT/FORS. This source could be associated with the putative SMG but it is most likely offset spatially by 10-30kpc and they could be interacting galaxies. With a FIR luminosity in the range [5-15]x10^{11} L_sol corresponding to a star formation rate in the range [80-260]M_sol/yr, this SMG would be more representative than the extreme starbursts usually detected at z>4. With a total magnification of ~25 it would open a unique window to the 'normal' dusty galaxies at the end of the epoch of reionization.
[65]  [pdf] - 765760
Serendipitous detection of an overdensity of Herschel-SPIRE 250 micron sources south of MRC1138-26
Comments: 10 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2013-09-17
We report the serendipitous detection of a significant overdensity of Herschel-SPIRE 250 micron sources in the vicinity of MRC1138-26. We use an adaptive kernel density estimate to quantify the significance, including a comparison with other fields. The overdensity has a size of ~3.5-4' and stands out at ~5sigma with respect to the background estimate. No features with similar significance were found in four extragalactic control fields: GOODS-North, Lockman, COSMOS and UDS. The chance of having a similar overdensity in a field with the same number but randomly distributed sources is less than 2%. The clump is also visible as a low surface brightness feature in the Planck 857 GHz map. We detect 76 sources at 250 micron (with a signal-to-noise ratio greater than 3), in a region of 4' radius; 43 of those are above a flux density limit of 20 mJy. This is a factor of 3.6 in excess over the average in the four control fields, considering only the sources above 20 mJy. We also find an excess in the number counts of sources with 250 micron flux densities between 30 and 40 mJy, compared to deep extragalactic blank-field number counts. Assuming a fixed dust temperature (30 K) and emissivity (beta=1.5) a crude, blackbody-derived redshift distribution, zBB, of the detected sources is significantly different from the distributions in the control fields and exhibits a significant peak at zBB ~ 1.5, although the actual peak redshift is highly degenerate with the temperature. We tentatively suggest, based on zBB and the similar S250/S350 colours of the sources within the peak, that a significant fraction of the sources in the clump may be at a similar redshift. Since the overdensity lies ~7' south of the z=2.16 Spiderweb protocluster MRC1138-26, an intriguing possibility (that is presently unverifiable given the data in hand) is that it lies within the same large-scale structure.(abridged)
[66]  [pdf] - 1179028
Multi-wavelength landscape of the young galaxy cluster RXJ1257.2+4738 at z=0.866: I. The infrared view
Comments: 16 pages, 14 figures, accepted for publication in A&A
Submitted: 2013-09-05, last modified: 2013-09-09
We performed a thorough analysis of the star formation activity in the young massive galaxy cluster RXJ1257+4738 at z=0.866, with emphasis on the relationship between the local environment of the cluster galaxies and their star formation activity. We present an optical and IR study that benefited from the large amount of data available for this cluster, including new OSIRIS/GTC and Herschel imaging observations. Using a optical-to-NIR multi-wavelength catalogue, we measured photometric redshifts through a chi2 SED-fitting procedure. We implemented a reliable and carefully chosen cluster membership selection criterion including Monte Carlo simulations and derived a sample of 292 reliable cluster member galaxies for which we measured the following properties: optical colours, stellar masses, ages, ultraviolet luminosities and local densities. Using the MIPS 24um and Herschel data, we measured total IR luminosities and SFR. Of the sample of 292 cluster galaxies, 38 show FIR emission with an SFR between 0.5 and 45 Msun/yr. The spatial distribution of the FIR emitters within the cluster density map and the filament-like overdensities observed suggest that RXJ1257 is not virialised, but is in the process of assembly. The average star formation as a function of the cluster environment parametrised by the local density of galaxies does not show any clear trend. However, the fraction of SF galaxies unveils that the cluster intermediate-density regions is preferred for the SF activity to enhance, since we observe a significant increase of the FIR-emitter fraction in this environment. Focusing on the optically red SF galaxies, we can support the interpretation of this population as dusty red galaxies, since we observe an appreciable difference in their extinction compared with the blue population.
[67]  [pdf] - 1178788
H-ATLAS: estimating redshifts of Herschel sources from sub-mm fluxes
Comments: 13 pages, 11 figures
Submitted: 2013-08-26
Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}
[68]  [pdf] - 1173518
Far-infrared spectroscopy of a lensed starburst: a blind redshift from Herschel
Comments: 6 pages, 5 figures, accepted for publication in MNRAS as a Letter
Submitted: 2013-08-20
We report the redshift of HATLAS J132427.0+284452 (hereafter HATLAS J132427), a gravitationally lensed starburst galaxy, the first determined 'blind' by the Herschel Space Observatory. This is achieved via the detection of [C II] consistent with z = 1.68 in a far-infrared spectrum taken with the SPIRE Fourier Transform Spectrometer. We demonstrate that the [C II] redshift is secure via detections of CO J = 2 - 1 and 3 - 2 using the Combined Array for Research in Millimeter-wave Astronomy and the Institut de Radioastronomie Millimetrique's Plateau de Bure Interferometer. The intrinsic properties appear typical of high-redshift starbursts despite the high lensing-amplified fluxes, proving the ability of the FTS to probe this population with the aid of lensing. The blind detection of [C II] demonstrates the potential of the SAFARI imaging spectrometer, proposed for the much more sensitive SPICA mission, to determine redshifts of multiple dusty galaxies simultaneously without the benefit of lensing.
[69]  [pdf] - 700899
A new HST/Herschel deep field at the North Ecliptic Pole: preparing the way for JWST, SPICA and Euclid
Comments: White paper submitted to the HST / Herschel Deep Fields working group. CANDELS reference added
Submitted: 2012-09-17, last modified: 2013-07-31
We propose a co-ordinated multi-observatory survey at the North Ecliptic Pole. This field is the natural extragalactic deep field location for most space observatories (e.g. containing the deepest Planck, WISE and eROSITA data), is in the continuous viewing zones for e.g. Herschel, HST, JWST, and is a natural high-visibility field for the L2 halo orbit of SPICA with deep and wide-field legacy surveys already planned. The field is also a likely deep survey location for the forthcoming Euclid mission. It is already a multi-wavelength legacy field in its own right (e.g. AKARI, LOFAR, SCUBA-2): the outstanding and unparalleled continuous mid-IR photometric coverage in this field and nowhere else enables a wide range of galaxy evolution diagnostics unachievable in any other survey field, by spanning the wavelengths of redshifted PAH and silicate features and the peak energy output of AGN hot dust. We argue from the science needs of Euclid and JWST, and from the comparative multiwavelength depths, that the logical approach is (1) a deep (H-UDF) UV/optical tile in the NEP over ~10 square arcminutes, and (2) an overlapping wide-field UV/optical HST survey tier covering >100 square arcminutes, with co-ordinated submm SPIRE mapping up to or beyond the submm point source confusion limit over a wider area and PACS data over the shallower HST tier.
[70]  [pdf] - 1164763
The Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4
Comments: 32 pages, 20 figures, 9 tables. Published in MNRAS. Replaced Fig. 8 (small scaling bug in the previous version)
Submitted: 2013-02-21, last modified: 2013-07-25
We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.
[71]  [pdf] - 1165564
HerMES: A Deficit in the Surface Brightness of the Cosmic Infrared Background Due to Galaxy Cluster Gravitational Lensing
Comments: 5 pages, 3 figures, replacing with version matching published manuscript
Submitted: 2013-03-26, last modified: 2013-07-18
We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB.
[72]  [pdf] - 1171826
Lyman break and UV-selected galaxies at z ~ 1: II. PACS-100um/160um FIR detections
Comments: Accepted for publication in MNRAS
Submitted: 2013-06-05, last modified: 2013-07-17
We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.
[73]  [pdf] - 1172267
Observing Extended Sources with the \Herschel SPIRE Fourier Transform Spectrometer
Comments: Accepted for publication by A&A
Submitted: 2013-06-24
The Spectral and Photometric Imaging Receiver (SPIRE) on the European Space Agency's Herschel Space Observatory utilizes a pioneering design for its imaging spectrometer in the form of a Fourier Transform Spectrometer (FTS). The standard FTS data reduction and calibration schemes are aimed at objects with either a spatial extent much larger than the beam size or a source that can be approximated as a point source within the beam. However, when sources are of intermediate spatial extent, neither of these calibrations schemes is appropriate and both the spatial response of the instrument and the source's light profile must be taken into account and the coupling between them explicitly derived. To that end, we derive the necessary corrections using an observed spectrum of a fully extended source with the beam profile and the source's light profile taken into account. We apply the derived correction to several observations of planets and compare the corrected spectra with their spectral models to study the beam coupling efficiency of the instrument in the case of partially extended sources. We find that we can apply these correction factors for sources with angular sizes up to \theta_{D} ~ 17". We demonstrate how the angular size of an extended source can be estimated using the difference between the sub-spectra observed at the overlap bandwidth of the two frequency channels in the spectrometer, at 959<\nu<989 GHz. Using this technique on an observation of Saturn, we estimate a size of 17.2", which is 3% larger than its true size on the day of observation. Finally, we show the results of the correction applied on observations of a nearby galaxy, M82, and the compact core of a Galactic molecular cloud, Sgr B2.
[74]  [pdf] - 1164699
Herschel-ATLAS: a binary HyLIRG pinpointing a cluster of starbursting proto-ellipticals
Comments: 15 pages, 8 colour figures; accepted for publication in ApJ
Submitted: 2013-02-18, last modified: 2013-06-12
Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.
[75]  [pdf] - 1171835
Flux calibration of the Herschel-SPIRE photometer
Comments: 20 pages, 17 figures, accepted for publication in MNRAS
Submitted: 2013-06-05
We describe the procedure used to flux calibrate the three-band submillimetre photometer in the Spectral and Photometric Imaging REceiver (SPIRE) 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/V, we then use observations of Neptune in which the beam of each bolometer is mapped using Neptune observed in a very fine scan pattern. The total instrumental uncertainties in the flux calibration for the individual bolometers is ~0.5% for most bolometers, although a few bolometers have uncertainties of ~1-5% 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%. This is considerably less than the absolute calibration uncertainty associated with the model of Neptune, which is estimated at 4%.
[76]  [pdf] - 1150890
HerMES: Cosmic Infrared Background Anisotropies and the Clustering of Dusty Star-Forming Galaxies
Comments: 29 pages, 15 figures, 10 tables. Accepted for publication in the Astrophysical Journal. Tabulated data and other software available at
Submitted: 2012-08-24, last modified: 2013-05-29
We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (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_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. 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 towards 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_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.
[77]  [pdf] - 747035
The rapid assembly of an elliptical galaxy of 400 billion solar masses at a redshift of 2.3
Comments: Appearing in Nature online on May 22 and in print on May 30. Submitted here is the accepted version (including the Supplementary Information), see for the final version
Submitted: 2013-05-21
Stellar archeology shows that massive elliptical galaxies today formed rapidly about ten billion years ago with star formation rates above several hundreds solar masses per year (M_sun/yr). Their progenitors are likely the sub-millimeter-bright galaxies (SMGs) at redshifts (z) greater than 2. While SMGs' mean molecular gas mass of 5x10^10 M_sun can explain the formation of typical elliptical galaxies, it is inadequate to form ellipticals that already have stellar masses above 2x10^11 M_sun at z ~ 2. Here we report multi-wavelength high-resolution observations of a rare merger of two massive SMGs at z = 2.3. The system is currently forming stars at a tremendous rate of 2,000 M_sun/yr. With a star formation efficiency an order-of-magnitude greater than that of normal galaxies, it will quench the star formation by exhausting the gas reservoir in only ~200 million years. At a projected separation of 19 kiloparsecs, the two massive starbursts are about to merge and form a passive elliptical galaxy with a stellar mass of ~4x10^11 M_sun. Our observations show that gas-rich major galaxy mergers, concurrent with intense star formation, can form the most massive elliptical galaxies by z ~ 1.5.
[78]  [pdf] - 1166499
Dust-obscured star-formation in the outskirts of XMMU J2235.3-2557, a massive galaxy cluster at z=1.4
Comments: 14 pages, accepted for publication in MNRAS
Submitted: 2013-05-08
Star-formation in the galaxy populations of local massive clusters is reduced with respect to field galaxies, and tends to be suppressed in the core region. Indications of a reversal of the star-formation--density relation have been observed in a few z >1.4 clusters. Using deep imaging from 100-500um from PACS and SPIRE onboard Herschel, we investigate the infrared properties of spectroscopic and photo-z cluster members, and of Halpha emitters in XMMU J2235.3-2557, one of the most massive, distant, X-ray selected clusters known. Our analysis is based mostly on fitting of the galaxies spectral energy distribution in the rest-frame 8-1000um. We measure total IR luminosity, deriving star formation rates (SFRs) ranging from 89-463 Msun/yr for 13 galaxies individually detected by Herschel, all located beyond the core region (r >250 kpc). We perform a stacking analysis of nine star-forming members not detected by PACS, yielding a detection with SFR=48 Msun/yr. Using a color criterion based on a star-forming galaxy SED at the cluster redshift we select 41 PACS sources as candidate star-forming cluster members. We characterize a population of highly obscured SF galaxies in the outskirts of XMMU J2235.3-2557. We do not find evidence for a reversal of the SF-density relation in this massive, distant cluster.
[79]  [pdf] - 1166242
The Herschel-PEP survey: evidence for downsizing in the hosts of dusty star-forming systems
Comments: 12 pages, 8 figures, to appear in MNRAS
Submitted: 2013-04-26
By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)
[80]  [pdf] - 1165355
The deepest Herschel-PACS far-infrared survey: number counts and infrared luminosity functions from combined PEP/GOODS-H observations
Comments: Accepted for publication in A&A; 22 pages, 13 figures; V2: updated to match accepted version
Submitted: 2013-03-18, last modified: 2013-04-22
We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.
[81]  [pdf] - 1117175
HerMES: deep number counts at 250, 350, and 500 microns in the COSMOS and GOODS-N fields and the build-up of the cosmic infrared background
Comments: 25 pages, 17 figures, 15 tables, accepted by A&A, bug with Fig. 3 fixed
Submitted: 2012-03-08, last modified: 2013-04-17
ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. 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. 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 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. 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, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.
[82]  [pdf] - 653862
A Dust-Obscured Massive Maximum-Starburst Galaxy at a Redshift of 6.34
Comments: 35 pages, 19 figures, 5 tables, including Supplementary Information (version formatted by the authors). To appear in Nature. Under press embargo until 18:00 London time/13:00 US Eastern Time on 17 April 2013
Submitted: 2013-04-15
Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.
[83]  [pdf] - 1164740
The Herschel census of infrared SEDs through cosmic time
Comments: 27 pages, 23 figures, 2 tables, accepted for publication in MNRAS
Submitted: 2013-02-20
Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), 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[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=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 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised 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<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.
[84]  [pdf] - 622856
The Mean Star Formation Rate of X-ray selected Active Galaxies and its Evolution from z=2.5: Results from PEP-Herschel
Comments: 18 pages, 9 figures, 3 tables
Submitted: 2012-03-27, last modified: 2013-02-05
(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.
[85]  [pdf] - 1159146
Panchromatic spectral energy distributions of Herschel sources
Comments: Accepted for publication in A&A. Some figures are presented in low resolution. The new galaxy templates are available for download at the address
Submitted: 2013-01-18
(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is then fitted with a modified version of the MAGPHYS code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an AGN. The defined Gaussian grouping is also used to identify rare sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z~1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 micron detected sources in PEP GOODS fields. AGN appear to be distributed in the stellar mass (M*) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M*-SFR-z space.
[86]  [pdf] - 1158570
Distant galaxy clusters in the XMM Large Scale Structure survey
Comments: Accepted by MNRAS, 25 pages, 15 figures
Submitted: 2012-12-17
(Abridged) Distant galaxy clusters provide important tests of the growth of large scale structure in addition to highlighting the process of galaxy evolution in a consistently defined environment at large look back time. We present a sample of 22 distant (z>0.8) galaxy clusters and cluster candidates selected from the 9 deg2 footprint of the overlapping X-ray Multi Mirror (XMM) Large Scale Structure (LSS), CFHTLS Wide and Spitzer SWIRE surveys. Clusters are selected as extended X-ray sources with an accompanying overdensity of galaxies displaying optical to mid-infrared photometry consistent with z>0.8. Nine clusters have confirmed spectroscopic redshifts in the interval 0.8<z<1.2, four of which are presented here for the first time. A further 11 candidate clusters have between 8 and 10 band photometric redshifts in the interval 0.8<z<2.2, while the remaining two candidates do not have information in sufficient wavebands to generate a reliable photometric redshift. All of the candidate clusters reported in this paper are presented for the first time. Those confirmed and candidate clusters with available near infrared photometry display evidence for a red sequence galaxy population, determined either individually or via a stacking analysis, whose colour is consistent with the expectation of an old, coeval stellar population observed at the cluster redshift. We further note that the sample displays a large range of red fraction values indicating that the clusters may be at different stages of red sequence assembly. We compare the observed X-ray emission to the flux expected from a suite of model clusters and find that the sample displays an effective mass limit M200 ~ 1e14 Msolar with all clusters displaying masses consistent with M200 < 5e14 Msolar. This XMM distant cluster study represents a complete sample of X-ray selected z>0.8 clusters.
[87]  [pdf] - 1158476
GOODS-Herschel: Separating High Redshift active galactic Nuclei and star forming galaxies Using Infrared Color Diagnostics
Comments: Accepted for publication in ApJ. 13 pages, 8 figures
Submitted: 2012-12-12
We have compiled a large sample of 151 high redshift (z=0.5-4) galaxies selected at 24 microns (S24>100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-infrared spectrum into contributions from star formation and activity in the galactic nuclei. In addition, we have a wealth of photometric data from Spitzer IRAC/MIPS and Herschel PACS/SPIRE. We explore how effective different infrared color combinations are at separating our mid-IR spectroscopically determined active galactic nuclei from our star forming galaxies. We look in depth at existing IRAC color diagnostics, and we explore new color-color diagnostics combining mid-IR, far-IR, and near-IR photometry, since these combinations provide the most detail about the shape of a source's IR spectrum. An added benefit of using a color that combines far-IR and mid-IR photometry is that it is indicative of the power source driving the IR luminosity. For our data set, the optimal color selections are S250/S24 vs. S8.0/S3.6 and S100/S24 vs. S8.0/S3.6; both diagnostics have ~10% contamination rate in the regions occupied primarily by star forming galaxies and active galactic nuclei, respectively. Based on the low contamination rate, these two new IR color-color diagnostics are ideal for estimating both the mid-IR power source of a galaxy when spectroscopy is unavailable and the dominant power source contributing to the IR luminosity. In the absence of far-IR data, we present color diagnostics using the WISE mid-IR bands which can efficiently select out high z (z~2) star forming galaxies.
[88]  [pdf] - 1158441
Inferring the mass of sub-millimetre galaxies by exploiting their gravitational magnification of background galaxies
Comments: 9 pages, 3 figures, accepted by MNRAS
Submitted: 2012-12-11
Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star formation in the Universe. In order to understand the physical processes involved in this mode of star formation, it is essential to gain knowledge about the dark matter structures - most importantly their masses - that sub-millimetre galaxies live in. Here we use the magnification effect of gravitational lensing to determine the average mass and dust content of sub-millimetre galaxies with 250mu flux densities of S_250>15mJy selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions of background Lyman-break galaxies at z~3-5 selected using optical data from the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation signal at the 7-sigma level over a sky area of one square degree, with ~80% of this signal being due to magnification, whereas the remaining ~20% comes from dust extinction. Adopting some simple assumptions for the dark matter and dust profiles and the redshift distribution enables us to estimate the average mass of the halos hosting the sub-millimetre galaxies to be log(M_200/M_sun)=13.17+0.05-0.08(stat.) and their average dust mass fraction (at radii of >10kpc) to be M_dust/M_200~6x10^-5. This supports the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside in massive group-sized halos, and are a crucial phase in the assembly and evolution of structure in the Universe.
[89]  [pdf] - 582610
SAFIR: testing the coexistence of AGN and star formation activity and the nature of the dusty torus in the local universe
Comments: 6 pages, 3 figures. To appear in "Highlights of Spanish Astrophysics VII", Proceedings of the X Scientific Meeting held in Valencia, Spain on 2012, July 9-13
Submitted: 2012-10-28
We present the Seyfert and star formation Activity in the Far-InfraRed (SAFIR) project, a small (15.1h) Herschel guaranteed time proposal performing PACS and SPIRE imaging of a small sample of nearby Seyfert galaxies. This project is aimed at studying the physical nature of the nuclear IR emission by means of multi-component spectral energy distribution (SED) fitting and the star formation properties of AGN hosts, as traced by cold dust. We summarize the results achieved so far and outline the on-going work.
[90]  [pdf] - 1152228
A Redshift Survey of Herschel Far-Infrared Selected Starbursts and Implications for Obscured Star Formation
Comments: 45 pages (23 text), 23 figures; ApJ accepted
Submitted: 2012-10-17, last modified: 2012-10-25
We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]
[91]  [pdf] - 1152231
A Population of z>2 Far-Infrared Herschel-SPIRE selected Starbursts
Comments: 25 pages, 10 figures; ApJ accepted
Submitted: 2012-10-17, last modified: 2012-10-25
We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.
[92]  [pdf] - 1151568
GOODS-Herschel: Impact of Active Galactic Nuclei and Star Formation Activity on Infrared Spectral Energy Distributions at High Redshift
Comments: Accepted for publication in ApJ. 24 pages, 11 figures. High resolution version of the figures, and composite SEDs presented in the paper, are available at
Submitted: 2012-09-21
We explore the effects of active galactic nuclei (AGN) and star formation activity on the infrared (0.3-1000 microns) spectral energy distributions of luminous infrared galaxies from z = 0.5 to 4.0. We have compiled a large sample of 151 galaxies selected at 24 microns (S24 > 100 uJy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spectroscopy, allowing us to decompose the mid-IR spectrum into contributions from star formation and AGN activity. A significant portion (~25%) of our sample is dominated by an AGN in the mid-IR. Based on the mid-IR classification, we divide our full sample into four sub-samples: z~1 star-forming (SF) sources; z~2 SF sources; AGN with clear 9.7 micron silicate absorption; and AGN with featureless mid-IR spectra. From our large spectroscopic sample and wealth of multi-wavelength data, including deep Herschel imaging at 100, 160, 250, 350, and 500 microns, we use 95 galaxies with complete spectral coverage to create a composite spectral energy distribution (SED) for each sub-sample. We then fit a two-temperature component modified blackbody to the SEDs. We find that the IR SEDs have similar cold dust temperatures, regardless of the mid-IR power source, but display a marked difference in the warmer dust temperatures. We calculate the average effective temperature of the dust in each sub-sample and find a significant (~20 K) difference between the SF and AGN systems. We compare our composite SEDs to local templates and find that local templates do not accurately reproduce the mid-IR features and dust temperatures of our high redshift systems. High redshift IR luminous galaxies contain significantly more cool dust than their local counterparts. We find that a full suite of photometry spanning the IR peak is necessary to accurately account for the dominant dust temperature components in high redshift IR luminous galaxies.
[93]  [pdf] - 1151411
Herschel-ATLAS/GAMA: spatial clustering of low-redshift sub-mm galaxies
Comments: 10 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2012-09-14
We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.
[94]  [pdf] - 1151366
Evidence for a wide range of UV obscuration in z ~ 2 dusty galaxies from the GOODS-Herschel survey
Comments: 9 pages, 9 figures. Accepted by ApJ
Submitted: 2012-09-12
Dusty galaxies at z ~ 2 span a wide range of relative brightness between rest-frame mid-infrared (8um) and ultraviolet wavelengths. We attempt to determine the physical mechanism responsible for this diversity. Dust-obscured galaxies (DOGs), which have rest-frame mid-IR to UV flux density ratios > 1000, might be abnormally bright in the mid-IR, perhaps due to prominent AGN and/or PAH emission, or abnormally faint in the UV. We use far-infrared data from the GOODS-Herschel survey to show that most DOGs with 10^12 L_Sun < L_IR < 10^13 L_Sun are not abnormally bright in the mid-IR when compared to other dusty galaxies with similar IR (8--1000um) luminosities. We observe a relation between the median IR to UV luminosity ratios and the median UV continuum power-law indices for these galaxies, and we find that only 24% have specific star formation rates which indicate the dominance of compact star-forming regions. This circumstantial evidence supports the idea that the UV- and IR-emitting regions in these galaxies are spatially coincident, which implies a connection between the abnormal UV faintness of DOGs and dust obscuration. We conclude that the range in rest-frame mid-IR to UV flux density ratios spanned by dusty galaxies at z ~ 2 is due to differing amounts of UV obscuration. Of galaxies with these IR luminosities, DOGs are the most obscured. We attribute differences in UV obscuration to either: 1) differences in the degree of alignment between the spatial distributions of dust and massive stars, or 2) differences in the total dust content.
[95]  [pdf] - 1040951
Measurements of CO redshifts with Z-Spec for lensed submillimeter galaxies discovered in the H-ATLAS survey
Comments: 28 pages, 15 figures, accepted for publication in ApJ, October 1, 2012, Issue 757 - 2
Submitted: 2010-09-29, last modified: 2012-08-17
We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.
[96]  [pdf] - 1125199
GOODS-Herschel: Ultra-deep XMM-Newton observations reveal AGN/star-formation connection
Comments: Accepted for publication in A&A
Submitted: 2012-07-30
Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)
[97]  [pdf] - 1124700
Discovery of "Warm Dust" Galaxies in Clusters at z~0.3: Evidence for Stripping of Cool Dust in the Dense Environment?
Comments: 25 pages, 13 figures, accepted for publication in ApJ
Submitted: 2012-07-09
Using far-infrared imaging from the "Herschel Lensing Survey", we derive dust properties of spectroscopically-confirmed cluster member galaxies within two massive systems at z~0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30K). Several sub-LIRG (L_IR < 10^11 L_sun) Bullet Cluster members are much warmer (T_dust > 37K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR-luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (> 3sigma), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e. preferentially located in the cluster periphery, although the galaxy hosts tend towards lower stellar masses (M_* < 10^10 M_sun). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover a SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR-luminosity or dust mass than the currently observed normal star-forming population.
[98]  [pdf] - 1124029
Herschel observations of a z~2 stellar mass selected galaxy sample drawn from the GOODS NICMOS Survey
Comments: Accepted for publication in MNRAS, 17 pages, 17 figures; added/updated references
Submitted: 2012-06-11, last modified: 2012-07-09
We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.
[99]  [pdf] - 1124380
Rapid Coeval Black Hole and Host Galaxy Growth in MRC 1138-262: The Hungry Spider
Comments: 11 pages, 2 figures, accepted for publication in ApJ
Submitted: 2012-06-25
We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.
[100]  [pdf] - 1084267
The cosmological analysis of X-ray cluster surveys: II- Application of the CR-HR method to the XMM archive
Comments: 26 pages, 24 figures, 9 tables. Accepted for publication in MNRAS (minor changes with respect to submitted version). The corresponding galaxy cluster catalogue is available at
Submitted: 2011-09-20, last modified: 2012-05-10
We have processed 2774 high-galactic observations from the XMM archive (as of May 2010) and extracted a serendipitous catalogue of some 850 clusters of galaxies based on purely X-ray criteria, following the methodology developed for the XMM-LSS survey. Restricting the sample to the highest signal-to-noise objects (347 clusters), we perform a cosmological analysis using the X-ray information only. The analysis consists in the modelling of the observed colour-magnitude (CR-HR) diagram constructed from cluster instrumental count-rates measured in the [0.5-2], [1-2] and [0.5-1] keV bands. A MCMC procedure simultaneously fits the cosmological parameters, the evolution of the cluster scaling laws and the selection effects. Our results are consistent with the sigma_8 and Omega_m values obtained by WMAP-5 and point toward a negative evolution of the cluster scaling relations with respect to the self-similar expectation. We are further able to constrain the cluster fractional radius xc0=r_c/r500c, to xc0=0.24 +/- 0.04. This study stresses again the critical role of selection effects in deriving cluster scaling relations, even in the local universe. Finally, we show that CR-HR method applied to the eRosita all-sky survey - provided that cluster photometric redshifts are available - will enable the determination of the equation of state of the dark energy at the level of the DETF stage IV predictions; simultaneously, the evolution of the cluster scaling-relations will be unambiguously determined. The XMM CLuster Archive Super Survey (XCLASS) serendipitous cluster catalogue is available online at:
[101]  [pdf] - 1530086
Herschel-ATLAS: Planck sources in the Phase 1 fields
Comments: 14 pages, 12 figures, 2 tables. Submitted to A&A
Submitted: 2012-04-17
We present the results of a cross-correlation of the Planck Early Release Compact Source Catalog (ERCSC) with the catalog of Herschel-ATLAS sources detected in the Phase 1 fields, covering 134.55 deg2. There are 28 ERCSC sources detected by Planck at 857 GHz in this area. As many as 16 of them are probably high Galactic latitude cirrus; 10 additional sources can be clearly identified as bright, low-z galaxies; one further source is resolved by Herschel as two relatively bright sources; and the last is resolved into an unusual condensation of low-flux, probably high-redshift point sources, around a strongly lensed Herschel-ATLAS source at z = 3.26. Our results demonstrate that the higher sensitivity and higher angular resolution H-ATLAS maps provide essential information for the interpretation of candidate sources extracted from Planck sub-mm maps.
[102]  [pdf] - 1117742
Physical properties of Lyman-alpha emitters at $z\sim 0.3$ from UV-to-FIR measurements
Comments: Accepted for publication in ApJ
Submitted: 2012-04-03
The analysis of the physical properties of low-redshift Ly$\alpha$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$\alpha$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$\alpha$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$\alpha$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.
[103]  [pdf] - 1093147
Enhanced star formation rates in AGN hosts with respect to inactive galaxies from PEP-Herschel observations
Comments: Final version
Submitted: 2012-01-20, last modified: 2012-03-23
We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.
[104]  [pdf] - 1117226
The Herschel Multi-tiered Extragalactic Survey: HerMES
Comments: 23 pages, 17 figures, 9 Tables, MNRAS accepted
Submitted: 2012-03-12
The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.
[105]  [pdf] - 1116784
The starburst-AGN connection in the merger galaxy Mrk 938: an infrared and X-ray view
Comments: 12 pages, 6 figures, 4 tables, accepted for publication in MNRAS
Submitted: 2012-02-21
Mrk938 is a luminous infrared galaxy in the local Universe believed to be the remnant of a galaxy merger. It shows a Seyfert 2 nucleus and intense star formation according to optical spectroscopic observations. We have studied this galaxy using new Herschel far-IR imaging data in addition to archival X-ray, UV, optical, near-IR and mid-IR data. Mid- and far-IR data are crucial to characterise the starburst contribution, allowing us to shed new light on its nature and to study the coexistence of AGN and starburst activity in the local Universe. The decomposition of the mid-IR Spitzer spectrum shows that the AGN bolometric contribution to the mid-IR and total infrared luminosity is small (Lbol(AGN)/LIR~0.02), which agrees with previous estimations. We have characterised the physical nature of its strong infrared emission and constrained it to a relatively compact emitting region of <2kpc. It is in this obscured region where most of the current star formation activity is taking place as expected for LIRGs. We have used Herschel imaging data for the first time to constrain the cold dust emission with unprecedented accuracy. We have fitted the integrated far-IR spectral energy distribution and derived the properties of the dust, obtaining a dust mass of 3x10^7Msun. The far-IR is dominated by emission at 35K, consistent with dust heated by the on-going star formation activity.
[106]  [pdf] - 662540
PACS-Herschel FIR detections of Lyman-alpha emitters at 2.0<z<3.5
Comments: Accepted for publication in A&A
Submitted: 2012-02-05
In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$\alpha$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$\mu$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$\alpha$ emitting galaxies whose rest-frame Ly$\alpha$ equivalent widths (Ly$\alpha$ EW$_{rest-frame}$) are above 20\AA, the typical threshold in narrow-band searches. Ly$\alpha$ emitting galaxies with Ly$\alpha$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$\alpha$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$\alpha$ emission.
[107]  [pdf] - 1116379
A Herschel view of the far-infrared properties of submillimetre galaxies
Comments: Accepted for publication in Astronomy and Astrophysics (A&A); 37 pages, 17 figures
Submitted: 2012-02-03
We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]
[108]  [pdf] - 1092759
The Relation Between Cool Cluster Cores and Herschel-Detected Star Formation in Brightest Cluster Galaxies
Comments: 26 pages, 11 figures, accepted for publication in ApJ
Submitted: 2012-01-05
We present far-infrared (FIR) analysis of 68 Brightest Cluster Galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500um), we calculate the obscured star formation rate (SFR). 22(+6.2,-5.3)% of the BCGs are detected in the far-infrared, with SFR= 1-150 M_sun/yr. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Ha emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2x10^11 L_sun), only a small (<0.4 mag) reddening correction is required for SFR(Ha) to agree with SFR_FIR. The relatively low Ha extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate form normal stellar mass loss.
[109]  [pdf] - 1648465
FIR/submm spectroscopy with Herschel: first results from the VNGS and H-ATLAS surveys
Submitted: 2011-12-06
The FIR/submm window is one of the least-studied regions of the electromagnetic spectrum, yet this wavelength range is absolutely crucial for understanding the physical processes and properties of the ISM in galaxies. The advent of the Herschel Space Observatory has opened up the entire FIR/submm window for spectroscopic studies. We present the first FIR/submm spectroscopic results on both nearby and distant galaxies obtained in the frame of two Herschel key programs: the Very Nearby Galaxies Survey and the Herschel ATLAS.
[110]  [pdf] - 1077108
The impact of evolving infrared spectral energy distributions of galaxies on star formation rate estimates
Comments: 25 pages, 15 figures, accepted to ApJ; v3 include many minor modifications following the referee report, added a section comparing various z~2 IR template calibrations
Submitted: 2011-06-06, last modified: 2011-12-01
We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7<z<2.5 normal star forming galaxies around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main sequence galaxies tend to have a similar nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8--60 um as a function of log(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z~1 and z~2. Corresponding SED template calibrations are provided for use if no stellar masses are in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.
[111]  [pdf] - 1084056
The Herschel Multi-Tiered Extragalactic Survey: SPIRE-mm Photometric Redshifts
Comments: 17 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2011-09-13, last modified: 2011-11-18
We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.
[112]  [pdf] - 1084852
The evolution of the star formation activity per halo mass up to redshift ~ 1.6 as seen by Herschel
Comments: 8 pages, 3 figures, accepted for publication on A&A
Submitted: 2011-10-13
Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (\Sigma(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field \Sigma(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field \Sigma(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the \Sigma(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher \Sigma(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar \Sigma(SFR/M), while massive clusters have significantly lower \Sigma(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group \Sigma(SFR/M) lies above the field \Sigma(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.
[113]  [pdf] - 1084356
HerMES: point source catalogues from deep Herschel-SPIRE observations
Comments: 15 pages, 6 figures, accepted for publication in MNRAS. For associated data release, see
Submitted: 2011-09-23
We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.
[114]  [pdf] - 1076591
GOODS-Herschel: an infrared main sequence for star-forming galaxies
Comments: Astronomy and Astrophysics, vol.533, A119 (version with updated figures 2 and 3). [26 pages, 25 figures]
Submitted: 2011-05-12, last modified: 2011-09-13
We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.
[115]  [pdf] - 1083953
GOODS-Herschel: the impact of galaxy-galaxy interactions on the far-infrared properties of galaxies
Comments: 10 pages, 6 figures. To appear in A&A
Submitted: 2011-09-09
We study the impact of galaxy-galaxy interactions on the FIR properties of galaxies and its evolution at 0<z<1.2. Using the high-z galaxies in the fields of GOODS observed by Herschel in the framework of the GOODS-Herschel key program and the local IRAS or AKARI-selected galaxies in the field of SDSS DR7, we investigate the dependence of galaxy properties on the morphology of and the distance to the nearest neighbor galaxy. We find that the SFRs and the SSFRs of galaxies, on average, depend on the morphology of and the distance to the nearest neighbor galaxy in this redshift range. When a late-type galaxy has a close neighbor galaxy, the SFR and the SSFR increase as it approaches a late-type neighbor, which is supported by K-S and MC tests with a significance level of >99%. However, the SFR and the SSFR decrease or do not change much as it approaches an early-type neighbor. The bifurcations of SFRs and SSFRs depending on the neighbor's morphology seem to occur at R_n~0.5r_vir,nei, which is supported by K-S and MC tests with a significance level of >98%. For all redshift bins, the SSFRs of late-type galaxies interacting with late-type neighbors are increased by factors of about 1.8+/-0.7 and 4.0+/-1.2 compared to those of non-interacting galaxies when the pair separation is smaller than 0.5r_vir,nei and 0.1r_vir,nei, respectively. The Tdust of both local and high-z late-type galaxies strongly interacting with late-type neighbors (i.e. R_n < 0.1r_vir,nei) appears to be higher than that of non-interacting galaxies with a significance level of 96-99%. However, Tdust of local late-type galaxies strongly interacting with early-type neighbors seems to be lower than or similar to that of non-interacting galaxies. Our results suggest that galaxy-galaxy interactions and mergers have been strongly affecting the SFR and the dust properties of star-forming galaxies over at least 8 billion years.
[116]  [pdf] - 1077452
GOODS-Herschel: The far-infrared view of star formation in AGN host galaxies since z~3
Comments: 21 pages; 14 Figures; MNRAS accepted; some figures have been degraded, full resolution version at:
Submitted: 2011-06-21, last modified: 2011-09-07
We present a study of the infrared properties of X-ray selected, moderate luminosity (Lx=10^{42}-10^{44}ergs/s) active galactic nuclei (AGNs) up to z~3, to explore the links between star formation in galaxies and accretion onto their central black holes. We use 100um and 160um fluxes from GOODS-Herschel -the deepest survey yet undertaken by the Herschel telescope- and show that in >94 per cent of cases these fluxes are dominated by the host. We find no evidence of any correlation between the X-ray and infrared luminosities of moderate AGNs at any redshift, suggesting that star-formation is decoupled from nuclear (AGN) activity. The star formation rates of AGN hosts increase strongly with redshift; by a factor of 43 from z<0.1 to z=2-3 for AGNs with the same X-ray luminosities. This increase is consistent with the factor of 25-50 increase in the specific star formation rates (SSFRs) of normal, star-forming (main-sequence) galaxies. Indeed, the average SSFRs of AGN hosts are only marginally (20 per cent) lower than those of main-sequence galaxies, with this small deficit being due to a fraction of AGNs residing in quiescent (low-SSFR) galaxies. We estimate 79+/-10 per cent of moderate AGNs are hosted in main-sequence galaxies, 15+/-7 per cent in quiescent galaxies and <10 per cent in strongly starbursting galaxies. The fractions of all main sequence galaxies at z<2 experiencing a period of moderate nuclear activity is strongly dependent on galaxy stellar mass (Mstars); rising from a few per cent at Mstars~10^{10}Msun to >20 per cent at Mstars>10^{11}Msun. Our results indicate that it is galaxy stellar mass that is most important in dictating whether a galaxy hosts a moderate luminosity AGN. We argue that the majority of moderate nuclear activity is fuelled by internal mechanisms rather than violent mergers, suggesting that disk instabilities could be an important AGN feeding mechanism.
[117]  [pdf] - 1083410
FIR measurements of Ly-$\alpha$ emitters at z$\lesssim$1.0: dust attenuation from PACS-\emph{Herschel}
Comments: Accepted for publication in ApJL
Submitted: 2011-08-15
One remaining open question regarding the physical properties of Ly$\alpha$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$\mu$m, PACS-160$\mu$m and \emph{Spitzer} MIPS-24$\mu$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$\mu$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.0, with a median value of A$_{1200\textrm{\AA}}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$\mu$m and PACS, with $A_{1200\AA}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.
[118]  [pdf] - 1083185
GOODS-Herschel A population of 24micron dropout sources at z < 2
Comments: Accepted for publication at A&A
Submitted: 2011-08-03, last modified: 2011-08-05
Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.
[119]  [pdf] - 389420
Physical conditions of the interstellar medium of high-redshift, strongly lensed submillimetre galaxies from the Herschel-ATLAS
Comments: 13 pages, 6 figures, published in MNRAS
Submitted: 2011-05-19, last modified: 2011-07-22
We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.
[120]  [pdf] - 1077977
Gas and dust in a submillimeter galaxy at z = 4.24 from the Herschel ATLAS
Comments: To appear in ApJ
Submitted: 2011-07-14
We report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.
[121]  [pdf] - 1077928
Resolving the nuclear dust distribution of the Seyfert 2 galaxy NGC 3081
Comments: 6 pages, 4 figures. Accepted for publication in MNRAS letters
Submitted: 2011-07-12
We report far-infrared (FIR) imaging of the Seyfert 2 galaxy NGC 3081 in the range 70-500 micron, obtained with an unprecedented angular resolution, using the Herschel Space Observatory instruments PACS and SPIRE. The 11 kpc (~70 arcsec) diameter star-forming ring of the galaxy appears resolved up to 250 micron. We extracted infrared (1.6-500 micron) nuclear fluxes, that is active nucleus-dominated fluxes, and fitted them with clumpy torus models, which successfully reproduce the FIR emission with small torus sizes. Adding the FIR data to the near- and mid-infrared spectral energy distribution (SED) results in a torus radial extent of Ro=4(+2/-1) pc, as well as in a flat radial distribution of the clouds (i.e. the q parameter). At wavelengths beyond 200 micron, cold dust emission at T=28+/-1 K from the circumnuclear star-forming ring of 2.3 kpc (~15 arcsec) in diameter starts making a contribution to the nuclear emission. The dust in the outer parts of the galaxy is heated by the interstellar radiation field (19+/-3 K).
[122]  [pdf] - 1077587
On Star Formation Rates and Star Formation Histories of Galaxies out to z ~ 3
Comments: Accepted for publication in The Astrophysical Journal, 19 pages, 15 figures
Submitted: 2011-06-27
We compare multi-wavelength SFR indicators out to z~3 in GOODS-South. Our analysis uniquely combines U-to-8um photometry from FIREWORKS, MIPS 24um and PACS 70, 100, and 160um photometry from the PEP survey, and Ha spectroscopy from the SINS survey. We describe a set of conversions that lead to a continuity across SFR indicators. A luminosity-independent conversion from 24um to total infrared luminosity yields estimates of LIR that are in the median consistent with the LIR derived from PACS photometry, albeit with significant scatter. Dust correction methods perform well at low to intermediate levels of star formation. They fail to recover the total amount of star formation in systems with large SFR_IR/SFR_UV ratios, typically occuring at the highest SFRs (SFR_UV+IR \gtrsim 100 Msun/yr) and redshifts (z \gtrsim 2.5) probed. Finally, we confirm that Ha-based SFRs at 1.5<z<2.6 are consistent with SFR_SED and SFR_UV+IR provided extra attenuation towards HII regions is taken into account (Av,neb = Av,continuum / 0.44). With the cross-calibrated SFR indicators in hand, we perform a consistency check on the star formation histories inferred from SED modeling. We compare the observed SFR-M relations and mass functions at a range of redshifts to equivalents that are computed by evolving lower redshift galaxies backwards in time. We find evidence for underestimated stellar ages when no stringent constraints on formation epoch are applied. We demonstrate how resolved SED modeling, or alternatively deep UV data, may help to overcome this bias. The age bias is most severe for galaxies with young stellar populations, and reduces towards older systems. Finally, our analysis suggests that SFHs typically vary on timescales that are long (at least several 100 Myr) compared to the galaxies' dynamical time.
[123]  [pdf] - 1077352
PACS Evolutionary Probe (PEP) - A Herschel Key Program
Comments: 13 pages, 12 figures, accepted for publication in A&A
Submitted: 2011-06-16
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 in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, 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 superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.
[124]  [pdf] - 1077315
Building the cosmic infrared background brick by brick with Herschel/PEP
Comments: Proposed for acceptance on A&A
Submitted: 2011-06-15
The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).
[125]  [pdf] - 1076763
The PEP survey: clustering of infrared-selected galaxies and structure formation at z~2 in the GOODS South
Comments: 14 pages, 8 figures, MNRAS accepted
Submitted: 2011-05-20
ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.
[126]  [pdf] - 1053496
Discovery of a Multiply-Lensed Submillimeter Galaxy in Early HerMES Herschel/SPIRE Data
Comments: Accepted for publication in ApJL
Submitted: 2011-04-20
We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.
[127]  [pdf] - 1053500
Modeling of the HERMES J105751.1+573027 submillimeter source lensed by a dark matter dominated foreground group of galaxies
Comments: Submitted to ApJ
Submitted: 2011-04-20
We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.
[128]  [pdf] - 1042131
Spitzer Imaging of Herschel-ATLAS Gravitationally Lensed Submillimeter Sources
Comments: 6 pages, 4 figures, 1 table. Accepted for Publication in ApJ Letters. Minor corrections to address referee's comments
Submitted: 2010-11-22, last modified: 2011-01-13
We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.
[129]  [pdf] - 955849
Sub-millimetre galaxies reside in dark matter halos with masses greater than 3x10^11 solar masses
Comments: Accepted for publication as a Letter to Editor in Nature; This paper is under press embargo until Nature publication (early February). Supplementary material and electronic tables will be available from at the time of publication. For HerMES team information see and for data releases see
Submitted: 2011-01-05
The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. 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 a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. 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.
[130]  [pdf] - 275699
The GaLAxy Cluster Evolution Survey (GLACE): introduction and first results
Comments: 8 pages, 4 figures, submitted to the proceedings of the IX Scientific Meeting of the Spanish Astronomical Society (SEA)
Submitted: 2010-12-08
Aimed at understanding the evolution of galaxies in clusters, the GLACE survey is mapping a set of optical lines ([OII]3727, [OIII]5007, Hbeta and Halpha/[NII] when possible) in several galaxy clusters at redshift around 0.40, 0.63 and 0.86, using the Tuneable Filters (TF) of the OSIRIS instrument (Cepa et al. 2005) at the 10.4m GTC telescope. This study will address key questions about the physical processes acting upon the infalling galaxies during the course of hierarchical growth of clusters. GLACE is already ongoing: we present some preliminary results on our observations of the galaxy cluster Cl0024+1654 at z = 0.395; on the other hand, GLACE@0.86 has been approved as ESO/GTC large project to be started in 2011.
[131]  [pdf] - 1034896
GBT Zpectrometer CO(1-0) Observations of the Strongly-Lensed Submillimeter Galaxies from the Herschel ATLAS
Comments: v1 Submitted to ApJL, 2010 Sept. 11; v2 Accepted by ApJL, 2010 Nov. 19; 15 pages including tables and figures
Submitted: 2010-09-11, last modified: 2010-11-23
The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) has uncovered a population of strongly-lensed submillimeter galaxies (SMGs). The Zpectrometer instrument on the Green Bank Telescope (GBT) was used to measure the redshifts and constrain the masses of the cold molecular gas reservoirs for two candidate high-redshift lensed sources. We derive CO(1-0) redshifts of z=3.042+/-0.001 and z=2.625+/-0.001, and measure molecular gas masses of (1--3)x10^{10}Msun, corrected for lens amplification and assuming a conversion factor of alhpa=0.8 Msun(K km/s pc^2)^{-1}. We find typical L(IR)/L'(CO) ratios of 120+/-40 and 140+/-50 Lsun (K km/s pc^2)^{-1}, which are consistent with those found for local ULIRGs and other high-redshift SMGs. From analysis of published data, we find no evidence for enhanced L(IR)/L'(CO(1-0)) ratios for the SMG population in comparison to local ULIRGs. The GBT results highlight the power of using the CO lines to derive blind redshifts, which is challenging for the SMG population at optical wavelengths given their high obscuration.
[132]  [pdf] - 254499
The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies
Comments: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science VOL 330, (Nov 5th, 2010), doi:10.1126/science.1193420
Submitted: 2010-11-04
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 starforming 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.
[133]  [pdf] - 1041621
The XMM-LSS survey: optical assessment and properties of different X-ray selected cluster classes
Comments: Seven figures given in jpg format due to space limitation. The full sample of X-ray clusters with associated optical spectroscopic data is available via the L3SDB database (
Submitted: 2010-10-29
XMM and Chandra opened a new area for the study of clusters of galaxies. Not only for cluster physics but also, for the detection of faint and distant clusters that were inaccessible with previous missions. This article presents 66 spectroscopically confirmed clusters (0.05<z<1.5) within an area of 6 deg2 enclosed in the XMM-LSS survey. Almost two thirds have been confirmed with dedicated spectroscopy only and 10% have been confirmed with dedicated spectroscopy supplemented by literature redshifts. Sub-samples, or classes, of extended-sources are defined in a two-dimensional X-ray parameter space allowing for various degrees of completeness and contamination. We describe the procedure developed to assess the reality of these cluster candidates using the CFHTLS photometric data and spectroscopic information from our own follow-up campaigns. Most of these objects are low mass clusters, hence constituting a still poorly studied population. In a second step, we quantify correlations between the optical properties such as richness or velocity dispersion and the cluster X-ray luminosities. We examine the relation of the clusters to the cosmic web. Finally, we review peculiar structures in the surveyed area like very distant clusters and fossil groups.
[134]  [pdf] - 249436
HerMES : SPIRE detection of high redshift massive compact galaxies in GOODS-N field
Comments: 6 pages, 3 figure, 1 table, MNRAS Herschel issue, in V2 added references and corrected typos
Submitted: 2010-10-04, last modified: 2010-10-28
We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).
[135]  [pdf] - 1041584
The first release of data from the Herschel ATLAS: the SPIRE images
Comments: 8 pages, 4 figures, submitted to MNRAS
Submitted: 2010-10-27
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 sheer size of the data, 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 latter is estimated to 5.3, 6.4, and 6.7 mJy/beam (1-{\sigma}), at 250, 350, and 500 \mu{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 dataset is available from
[136]  [pdf] - 1040915
HerMES: Deep Galaxy Number Counts from a P(D) Fluctuation Analysis of SPIRE Science Demonstration Phase Observations
Comments: Accepted for publication in MNRAS
Submitted: 2010-09-28
Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-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 analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) 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 counts 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.
[137]  [pdf] - 1034824
The Herschel Multi-Tiered Extragalactic Survey: Source Extraction and Cross-Identifications in Confusion-Dominated SPIRE Images
Comments: Accepted for publication in MNRAS, 19 pages, 11 figures
Submitted: 2010-09-08
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 minimise 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 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration 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; Savage & Oliver et al. 2006) by Smith et al. 2010. 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 dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.
[138]  [pdf] - 1034774
The Deep SPIRE HerMES Survey: Spectral Energy Distributions and their Astrophysical Indications at High Redshift
Comments: 9 pages, 3 large tables, 2 tables in text, 5 figures Accepted for publication in MNRAS
Submitted: 2010-09-07, last modified: 2010-09-08
The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.
[139]  [pdf] - 1034730
Evolution of Dust Temperature of Galaxies through Cosmic Time as seen by Herschel
Hwang, H. S.; Elbaz, D.; Magdis, G. E.; Daddi, E.; Symeonidis, M.; Altieri, B.; Amblard, A.; Andreani, P.; Arumugam, V.; Auld, R.; Aussel, H.; Babbedge, T.; Berta, S.; Blain, A.; Bock, J.; Bongiovanni, A.; Boselli, A.; Buat, V.; Burgarella, D.; Castro-Rodriguez, N.; Cava, A.; Cepa, J.; Chanial, P.; Chapin, E.; Chary, R. -R.; Cimatti, A.; Clements, D. L.; Conley, A.; Conversi, L.; Cooray, A.; Dannerbauer, H.; Dickinson, M.; Dominguez, H.; Dowell, C. D.; Dunlop, J. S.; Dwek, E.; Eales, S.; Farrah, D.; Schreiber, N. Forster; Fox, M.; Franceschini, A.; Gear, W.; Genzel, R.; Glenn, J.; Griffin, M.; Gruppioni, C.; Halpern, M.; Hatziminaoglou, E.; Ibar, E.; Isaak, K.; Ivison, R. J.; Jeong, W. -S.; Lagache, G.; Borgne, D. Le; Floc'h, E. Le; Lee, H. M.; Lee, J. C.; Lee, M. G.; Levenson, L.; Lu, N.; Lutz, D.; Madden, S.; Maffei, B.; Magnelli, B.; Mainetti, G.; Maiolino, R.; Marchetti, L.; Mortier, A. M. J.; Nguyen, H. T.; Nordon, R.; O'Halloran, B.; Okumura, K.; Oliver, S. J.; Omont, A.; Page, M. J.; Panuzzo, P.; Papageorgiou, A.; Pearson, C. P.; Perez-Fournon, I.; Garcia, A. M. Perez; Poglitsch, A.; Pohlen, M.; Popesso, P.; Pozzi, F.; Rawlings, J. I.; Rigopoulou, D.; Riguccini, L.; Rizzo, D.; Rodighiero, G.; Roseboom, I. G.; Rowan-Robinson, M.; Saintonge, A.; Portal, M. Sanchez; Santini, P.; Sauvage, M.; Schulz, B.; Scott, Douglas; Seymour, N.; Shao, L.; Shupe, D. L.; Smith, A. J.; Stevens, J. A.; Sturm, E.; Tacconi, L.; Trichas, M.; Tugwell, K. E.; Vaccari, M.; Valtchanov, I.; Vieira, J. D.; Vigroux, L.; Wang, L.; Ward, R.; Wright, G.; Xu, C. K.; Zemcov, M.
Comments: 9 pages, 6 figures, To appear in MNRAS
Submitted: 2010-09-06
We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.
[140]  [pdf] - 711261
Herschel FIR counterparts of selected Ly-alpha emitters at z~2.2. Fast evolution since z~3 or missed obscured AGNs?
Comments: 11 pages (including Appendices), 6 figures. Accepted for publication in Astronomy & Astrophysics Letters (two references added)
Submitted: 2010-08-10, last modified: 2010-08-19
Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.
[141]  [pdf] - 1033834
A First Glimpse into the far-IR properties of high-z UV-selected Galaxies: Herschel/PACS observations of z~3 LBGs
Comments: Accepted for publication at ApJL
Submitted: 2010-07-22
We present first insights into the far-IR properties for a sample of IRAC and MIPS-24um detected Lyman Break Galaxies (LBGs) at z ~ 3, as derived from observations in the northern field of the Great Observatories Origins Survey (GOODS-N) carried out with the PACS instrument on board the Herschel Space Observatory. Although none of our galaxies are detected by Herschel, we employ a stacking technique to construct, for the first time, the average spectral energy distribution of infrared luminous LBGs from UV to radio wavelengths. We derive a median IR luminosity of L_{IR} = 1.6 x 10^12 Lo, placing the population in the class of ultra luminous infrared galaxies (ULIRGs). Complementing our study with existing multi-wavelength data, we put constraints on the dust temperature of the population and find that for their L_{IR}, MIPS-LBGs are warmer than submm-luminous galaxies while they fall in the locus of the L_{IR}-T_{d} relation of the local ULIRGs. This, along with estimates based on the average SED, explains the marginal detection of LBGs in current sub-mm surveys and suggests that these latter studies introduce a bias towards the detection of colder ULIRGs in the high-z universe, while missing high-z ULIRGS with warmer dust.
[142]  [pdf] - 205573
The Herschel-SPIRE Legacy Survey (HSLS): the scientific goals of a shallow and wide submillimeter imaging survey with SPIRE
Cooray, Asantha; Eales, Steve; Chapman, Scott; Clements, David L.; Dore, Olivier; Farrah, Duncan; Jarvis, Matt J.; Kaplinghat, Manoj; Negrello, Mattia; Melchiorri, Alessandro; Peiris, Hiranya; Pope, Alexandra; Santos, Mario G.; Serjeant, Stephen; Thompson, Mark; White, Glenn; Amblard, Alexandre; Banerji, Manda; Corasaniti, Pier-Stefano; Das, Sudeep; de_Bernardis, Francesco; de_Zotti, Gianfranco; Giannantonio, Tommaso; Gonzalez, Joaquin Gonzalez-Nuevo; Khostovan, Ali Ahmad; Mitchell-Wynne, Ketron; Serra, Paolo; Song, Yong-Seon; Vieira, Joaquin; Wang, Lingyu; Zemcov, Michael; Abdalla, Filipe; Afonso, Jose; Aghanim, Nabila; Andreani, Paola; Aretxaga, Itziar; Auld, Robbie; Baes, Maarten; Baker, Andrew; Barkats, Denis; Barreiro, R. Belen; Bartolo, Nicola; Barton, Elizabeth; Barway, Sudhanshu; Battistelli, Elia Stefano; Baugh, Carlton; Beelen, Alexander; Benabed, Karim; Blain, Andrew; Bland-Hawthorn, Joss; Bock, James~J.; Bond, J. Richard; Borrill, Julian; Borys, Colin; Boselli, Alessandro; Bouchet, Francois R.; Bridge, Carrie; Brighenti, Fabrizio; Buat, Veronique; Buote, David; Burgarella, Denis; Bussmann, Robert; Calabrese, Erminia; Cantalupo, Christopher; Carlberg, Raymond; Carvalho, Carla Sofia; Casey, Caitlin; Cava, Antonio; Cepa, Jordi; Chapin, Edward; Chary, Ranga Ram; Chen, Xuelei; Colafrancesco, Sergio; Cole, Shaun; Coles, Peter; Conley, Alexander; Conversi, Luca; Cooke, Jeff; Crawford, Steven; Cress, Catherine; da Cunha, Elisabete; Dalton, Gavin; Danese, Luigi; Dannerbauer, Helmut; Davies, Jonathan; de Bernardis, Paolo; de Putter, Roland; Devlin, Mark; Diego, Jose M.; Dole, Herve; Douspis, Marian; Dunkley, Joanna; Dunlop, James; Dunne, Loretta; Dunner, Rolando; Dye, Simon; Efstathiou, George; Egami, Eiichi; Fang, Taotao; Ferrero, Patrizia; Franceschini, Alberto; Frazer, Christopher C.; Frayer, David; Frenk, Carlos; Ganga, Ken; Gavazzi, Raphael; Glenn, Jason; Gong, Yan; Gonzalez-Solares, Eduardo; Griffin, Matt; Guo, Qi; Gurwell, Mark; Hajian, Amir; Halpern, Mark; Hanson, Duncan; Hardcastle, Martin; Hatziminaoglou, Evanthia; Heavens, Alan; Heinis, Sebastien; Herranz, Diego; Hilton, Matt; Ho, Shirley; Holwerda, Benne W.; Hopwood, Rosalind; Horner, Jonathan; Huffenberger, Kevin; Hughes, David H.; Hughes, John P.; Ibar, Edo; Ivison, Rob; Jackson, Neal; Jaffe, Andrew; Jenness, Timothy; Joncas, Gilles; Joudaki, Shahab; Kaviraj, Sugata; Kim, Sam; King, Lindsay; Kisner, Theodore; Knapen, Johan; Kniazev, Alexei; Komatsu, Eiichiro; Koopmans, Leon; Kuo, Chao-Lin; Lacey, Cedric; Lahav, Ofer; Lasenby, Anthony N.; Lawrence, Andy; Lee, Myung Gyoon; Leeuw, Lerothodi L.; Levenson, Louis R.; Lewis, Geraint; Loaring, Nicola; Lopez-Caniego, Marcos; Maddox, Steve; Marriage, Tobias; Marsden, Gaelen; Martinez-Gonzalez, Enrique; Masi, Silvia; Matarrese, Sabino; Mathews, William G.; Matsuura, Shuji; McMahon, Richard; Mellier, Yannick; Menanteau, Felipe; Michalowski, Michal~J.; Millea, Marius; Mobasher, Bahram; Mohanty, Subhanjoy; Montier, Ludovic; Moodley, Kavilan; Moriarty-Schieven, Gerald H.; Mortier, Angela; Munshi, Dipak; Murphy, Eric; Nandra, Kirpal; Natoli, Paolo; Nguyen, Hien; Oliver, Seb; Omont, Alain; Page, Lyman; Page, Mathew; Paladini, Roberta; Pandolfi, Stefania; Pascale, Enzo; Patanchon, Guillaume; Peacock, John; Pearson, Chris; Perez-Fournon, Ismael; Perez-Gonz, Pablo G.; Piacentini, Francesco; Pierpaoli, Elena; Pohlen, Michael; Pointecouteau, Etienne; Polenta, Gianluca; Rawlings, Jason; Reese, Erik~D.; Rigby, Emma; Rodighiero, Giulia; Romero-Colmenero, Encarni; Roseboom, Isaac; Rowan-Robinson, Michael; Sanchez-Portal, Miguel; Schmidt, Fabian; Schneider, Michael; Schulz, Bernhard; Scott, Douglas; Sedgwick, Chris; Sehgal, Neelima; Seymour, Nick; Sherwin, Blake~D.; Short, Jo; Shupe, David; Sievers, Jonathan; Skibba, Ramin; Smidt, Joseph; Smith, Anthony; Smith, Daniel J. B.; Smith, Matthew W. L.; Spergel, David; Staggs, Suzanne; Stevens, Jason; Switzer, Eric; Takagi, Toshinobu; Takeuchi, Tsutomu; Temi, Pasquale; Trichas, Markos; Trigilio, Corrado; Tugwell, Katherine; Umana, Grazia; Vacca, William; Vaccari, Mattia; Vaisanen, Petri; Valtchanov, Ivan; van der Heyden, Kurt; van der Werf, Paul P.; van_Kampen, Eelco; van_Waerbeke, Ludovic; Vegetti, Simona; Veneziani, Marcella; Verde, Licia; Verma, Aprajita; Vielva, Patricio; Viero, Marco P.; Vilaro, Baltasar Vila; Wardlow, Julie; Wilson, Grant; Wright, Edward L.; Xu, C. Kevin; Yun, Min S.
Comments: White paper supplement to the proposal submitted by the HSLS science team to ESA for Herschel open-time programs
Submitted: 2010-07-20, last modified: 2010-07-22
A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.
[143]  [pdf] - 1033625
Unveiling Far-Infrared Counterparts of Bright Submillimeter Galaxies Using PACS Imaging
Comments: 6 pages, 4 figures. ApJ Letters in press
Submitted: 2010-07-12
We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of >=3 at 100 micron down to 3.0 mJy, and/or at 160 micron down to 5.7 mJy. The fraction of SMGs seen at 160 micron is higher than that at 100 micron. About 50% of radio-identified SMGs are associated with PACS sources. We find a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMG sample) remain unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength may imply higher redshifts for them than for the average SMG population (e.g., z>3-4). The Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far-infrared spectral energy distributions have significantly different shapes than template libraries of local infrared galaxies.
[144]  [pdf] - 1033621
Measures of star formation rates from Infrared (Herschel) and UV (GALEX) emissions of galaxies in the HerMES fields
Comments: accepted for publication in MNRAS, Herschel special issue
Submitted: 2010-07-12
The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .
[145]  [pdf] - 1032842
The dust content of high-z submillimeter galaxies revealed by Herschel
Comments: Accepted for publication in Astronomy & Astrophysics Letters. One reference updated
Submitted: 2010-05-31, last modified: 2010-06-07
We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.
[146]  [pdf] - 1032794
In-flight calibration of the Herschel-SPIRE instrument
Comments: Accepted for publication in Astronomy&Astrophysics, Herschel First Results special issue
Submitted: 2010-05-27
SPIRE, the Spectral and Photometric Imaging Receiver, is the Herschel Space Observatory's submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 {\mu}m, and an imaging Fourier transform spectrometer (FTS) covering 194-671 {\mu}m (447-1550 GHz). In this paper we describe the initial approach taken to the absolute calibration of the SPIRE instrument using a combination of the emission from the Herschel telescope itself and the modelled continuum emission from solar system objects and other astronomical targets. We present the photometric, spectroscopic and spatial accuracy that is obtainable in data processed through the "standard" pipelines. The overall photometric accuracy at this stage of the mission is estimated as 15% for the photometer and between 15 and 50% for the spectrometer. However, there remain issues with the photometric accuracy of the spectra of low flux sources in the longest wavelength part of the SPIRE spectrometer band. The spectrometer wavelength accuracy is determined to be better than 1/10th of the line FWHM. The astrometric accuracy in SPIRE maps is found to be 2 arcsec when the latest calibration data are used. The photometric calibration of the SPIRE instrument is currently determined by a combination of uncertainties in the model spectra of the astronomical standards and the data processing methods employed for map and spectrum calibration. Improvements in processing techniques and a better understanding of the instrument performance will lead to the final calibration accuracy of SPIRE being determined only by uncertainties in the models of astronomical standards.
[147]  [pdf] - 1032799
The Herschel-SPIRE instrument and its in-flight performance
Griffin, M. J.; Abergel, A.; Abreu, A.; Ade, P. A. R.; André, P.; Augueres, J. -L.; Babbedge, T.; Bae, Y.; Baillie, T.; Baluteau, J. -P.; Barlow, M. J.; Bendo, G.; Benielli, D.; Bock, J. J.; Bonhomme, P.; Brisbin, D.; Brockley-Blatt, C.; Caldwell, M.; Cara, C.; Castro-Rodriguez, N.; Cerulli, R.; Chanial, P.; Chen, S.; Clark, E.; Clements, D. L.; Clerc, L.; Coker, J.; Communal, D.; Conversi, L.; Cox, P.; Crumb, D.; Cunningham, C.; Daly, F.; Davis, G. R.; De Antoni, P.; Delderfield, J.; Devin, N.; Di Giorgio, A.; Didschuns, I.; Dohlen, K.; Donati, M.; Dowell, A.; Dowell, C. D.; Duband, L.; Dumaye, L.; Emery, R. J.; Ferlet, M.; Ferrand, D.; Fontignie, J.; Fox, M.; Franceschini, A.; Frerking, M.; Fulton, T.; Garcia, J.; Gastaud, R.; Gear, W. K.; Glenn, J.; Goizel, A.; Griffin, D. K.; Grundy, T.; Guest, S.; Guillemet, L.; Hargrave, P. C.; Harwit, M.; Hastings, P.; Hatziminaoglou, E.; Herman, M.; Hinde, B.; Hristov, V.; Huang, M.; Imhof, P.; Isaak, K. J.; Israelsson, U.; Ivison, R. J.; Jennings, D.; Kiernan, B.; King, K. J.; Lange, A. E.; Latter, W.; Laurent, G.; Laurent, P.; Leeks, S. J.; Lellouch, E.; Levenson, L.; Li, B.; Li, J.; Lilienthal, J.; Lim, T.; Liu, J.; Lu, N.; Madden, S.; Mainetti, G.; Marliani, P.; McKay, D.; Mercier, K.; Molinari, S.; Morris, H.; Moseley, H.; Mulder, J.; Mur, M.; Naylor, D. A.; Nguyen, H.; O'Halloran, B.; Oliver, S.; Olofsson, G.; Olofsson, H. -G.; Orfei, R.; Page, M. J.; Pain, I.; Panuzzo, P.; Papageorgiou, A.; Parks, G.; Parr-Burman, P.; Pearce, A.; Pearson, C.; Pérez-Fournon, I.; Pinsard, F.; Pisano, G.; Podosek, J.; Pohlen, M.; Polehampton, E. T.; Pouliquen, D.; Rigopoulou, D.; Rizzo, D.; Roseboom, I. G.; Roussel, H.; Rowan-Robinson, M.; Rownd, B.; Saraceno, P.; Sauvage, M.; Savage, R.; Savini, G.; Sawyer, E.; Scharmberg, C.; Schmitt, D.; Schneider, N.; Schulz, B.; Schwartz, A.; Shafer, R.; Shupe, D. L.; Sibthorpe, B.; Sidher, S.; Smith, A.; Smith, A. J.; Smith, D.; Spencer, L.; Stobie, B.; Sudiwala, R.; Sukhatme, K.; Surace, C.; Stevens, J. A.; Swinyard, B. M.; Trichas, M.; Tourette, T.; Triou, H.; Tseng, S.; Tucker, C.; Turner, A.; Vaccari, M.; Valtchanov, I.; Vigroux, L.; Virique, E.; Voellmer, G.; Walker, H.; Ward, R.; Waskett, T.; Weilert, M.; Wesson, R.; White, G. J.; Whitehouse, N.; Wilson, C. D.; Winter, B.; Woodcraft, A. L.; Wright, G. S.; Xu, C. K.; Zavagno, A.; Zemcov, M.; Zhang, L.; Zonca, E.
Comments: Accepted for publication in Astronomy & Astrophyics (Herschel first results special issue)
Submitted: 2010-05-27
The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.
[148]  [pdf] - 1026751
Herschel deep far-infrared counts through the Abell 2218 cluster-lens
Comments: Accepted for publication on the A&A Herschel Special Issue, final version
Submitted: 2010-05-10, last modified: 2010-05-24
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 microns 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 microns and 2 mJy at 160 microns, 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 microns. 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 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) 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 microns, disregarding the high uncertainties of DIRBE absolute values.
[149]  [pdf] - 1032651
First detection of the Sunyaev Zel'dovich effect increment at lambda < 650 um
Comments: Accepted for publication in Astronomy and Astrophysics (Herschel special issue)
Submitted: 2010-05-20
The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy clusters. This distortion is a decrement in the CMB intensity for lambda > 1.3 mm, an increment at shorter wavelengths, and small again by lambda ~250 um. As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 and 500 um and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are Delta I_{0} = 0.097 +- 0.019 MJy sr^{-1} at 350 um and Delta I_{0} = 0.268 +- 0.031 MJy sr^{-1} at 500 um, consistent with the SZ effect spectrum derived from previous measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE data at a significance of 2.1 sigma, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE.
[150]  [pdf] - 1032650
Improving the identification of high-z Herschel sources with position priors and optical/NIR and FIR/mm photometric redshifts
Comments: Accepted for publication in Astronomy and Astrophysics (Herschel special issue)
Submitted: 2010-05-20
We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 um sources as priors. Our catalogs are able to reliably (5 sigma) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 um channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.
[151]  [pdf] - 1032647
The Herschel Lensing Survey (HLS): overview
Comments: Accepted for publication in Astronomy and Astrophysics (Herschel special issue)
Submitted: 2010-05-20
The Herschel Lensing Survey (HLS) will conduct deep PACS and SPIRE imaging of ~40 massive clusters of galaxies. The strong gravitational lensing power of these clusters will enable us to penetrate through the confusion noise, which sets the ultimate limit on our ability to probe the Universe with Herschel. Here, we present an overview of our survey and a summary of the major results from our Science Demonstration Phase (SDP) observations of the Bullet Cluster (z=0.297). The SDP data are rich, allowing us to study not only the background high-redshift galaxies (e.g., strongly lensed and distorted galaxies at z=2.8 and 3.2) but also the properties of cluster-member galaxies. Our preliminary analysis shows a great diversity of far-infrared/submillimeter spectral energy distributions (SEDs), indicating that we have much to learn with Herschel about the properties of galaxy SEDs. We have also detected the Sunyaev-Zel'dovich (SZ) effect increment with the SPIRE data. The success of this SDP program demonstrates the great potential of the Herschel Lensing Survey to produce exciting results in a variety of science areas.
[152]  [pdf] - 1032649
Deep Herschel view of obscured star formation in the Bullet cluster
Comments: Accepted for publication in Astronomy and Astrophysics (Herschel Special Issue)
Submitted: 2010-05-20
We use deep, five band (100-500um) data from the Herschel Lensing Survey (HLS) to fully constrain the obscured star formation rate, SFR_FIR, of galaxies in the Bullet cluster (z=0.296), and a smaller background system (z=0.35) in the same field. Herschel detects 23 Bullet cluster members with a total SFR_FIR = 144 +/- 14 M_sun yr^-1. On average, the background system contains brighter far-infrared (FIR) galaxies, with ~50% higher SFR_FIR (21 galaxies; 207 +/- 9 M_sun yr^-1). SFRs extrapolated from 24um flux via recent templates (SFR_24) agree well with SFR_FIR for ~60% of the cluster galaxies. In the remaining ~40%, SFR_24 underestimates SFR_FIR due to a significant excess in observed S_100/S_24 (rest frame S_75/S_18) compared to templates of the same FIR luminosity.
[153]  [pdf] - 1032648
The far-infrared/submillimeter properties of galaxies located behind the Bullet cluster
Comments: 8 pages with 4 figures. Accepted for publication in Astronomy and Astrophysics (Herschel special issue)
Submitted: 2010-05-20
The Herschel Lensing Survey (HLS) takes advantage of gravitational lensing by massive galaxy clusters to sample a population of high-redshift galaxies which are too faint to be detected above the confusion limit of current far-infrared/submillimeter telescopes. Measurements from 100-500 micron bracket the peaks of the far-infrared spectral energy distributions of these galaxies, characterizing their infrared luminosities and star formation rates. We introduce initial results from our science demonstration phase observations, directed toward the Bullet cluster (1E0657-56). By combining our observations with LABOCA 870 micron and AzTEC 1.1 mm data we fully constrain the spectral energy distributions of 19 MIPS 24 micron selected galaxies which are located behind the cluster. We find that their colors are best fit using templates based on local galaxies with systematically lower infrared luminosities.This suggests that our sources are not like local ultra-luminous infrared galaxies in which vigorous star formation is contained in a compact highly dust-obscured region. Instead, they appear to be scaled up versions of lower luminosity local galaxies with star formation occurring on larger physical scales.
[154]  [pdf] - 1032583
HerMES: Halo Occupation Number and Bias Properties of Dusty Galaxies from Angular Clustering Measurements
Comments: 5 pages; accepted for publication in A&A Special Issue on Herschel First Science Results (July 2010). HerMES information available at
Submitted: 2010-05-18
We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.
[155]  [pdf] - 214381
The Herschel-ATLAS: Evolution of the 250um luminosity function out to z=0.5
Comments: Accepted for publication in A&A
Submitted: 2010-05-13, last modified: 2010-05-17
We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.
[156]  [pdf] - 1032505
Herschel unveils a puzzling uniformity of distant dusty galaxies
Comments: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages
Submitted: 2010-05-17
The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.
[157]  [pdf] - 1032531
A search for debris disks in the Herschel ATLAS
Comments: 5 pages, 3 figures, accepted by A&A for the Herschel Special Issue
Submitted: 2010-05-17
Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.
[158]  [pdf] - 1026707
Herschel PEP: The star-formation rates of 1.5<z<2.5 massive galaxies
Comments: Accepted for publication in the A&A Herschel first results special issue. v2 Correction to the meta data only
Submitted: 2010-05-06, last modified: 2010-05-17
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<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron 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.
[159]  [pdf] - 1026840
The HerMES SPIRE submillimeter local luminosity function
Comments: Accepted for publication in the A&A Herschel Special Issue
Submitted: 2010-05-12, last modified: 2010-05-14
Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.
[160]  [pdf] - 1032480
Star formation in AGN hosts in GOODS-N
Comments: Accepted for publication in A&A Herschel issue. 4 pages, 6 figures
Submitted: 2010-05-14
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 160um observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% of X-ray AGN individually at >3sig. The host far-infrared luminosity of AGN with L2-10~10^43erg/s 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<~10^44erg/s AGN at z>~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<10^44erg/s 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.
[161]  [pdf] - 1032460
Herschel-ATLAS: Dust temperature and redshift distribution of SPIRE and PACS detected sources using submillimetre colours
Comments: 5 pages; accepted for publication in A&A Special Issue on Herschel First Science Results (July 2010). H-ATLAS information available at
Submitted: 2010-05-13
We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns 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 \pm 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 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.
[162]  [pdf] - 1032453
HerMES: The Submillimeter Spectral Energy Distributions of Herschel/SPIRE-Detected Galaxies
Comments: 5 pages, 8 figures, for special Herschel issue
Submitted: 2010-05-13
We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.
[163]  [pdf] - 955274
Herschel-ATLAS: blazars in the SDP field
Comments: 5 pages, 2 figures, accepted for publication in A&A (Herschel special issue)
Submitted: 2010-05-13
To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.
[164]  [pdf] - 1032458
The Herschel-ATLAS: Extragalactic Number Counts from 250 to 500 Microns
Comments: Accepted for publication in Astronomy and Astrophysics
Submitted: 2010-05-13
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 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. 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-200mJy 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.
[165]  [pdf] - 1032459
Herschel ATLAS: The cosmic star formation history of quasar host galaxies
Comments: Accepted by Astronomy and Astrophysics, Herschel special issue
Submitted: 2010-05-13
We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.
[166]  [pdf] - 1032422
First results from HerMES on the evolution of the submillimetre luminosity function
Comments: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages
Submitted: 2010-05-12
We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.
[167]  [pdf] - 1026838
HerMES: SPIRE galaxy number counts at 250, 350 and 500 microns
Comments: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages
Submitted: 2010-05-12
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 at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.
[168]  [pdf] - 1032429
HerMES: The SPIRE confusion limit
Comments: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages
Submitted: 2010-05-12
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.
[169]  [pdf] - 1032424
HerMES: Far-infrared properties of known AGN in the HerMES fields
Comments: Astronomy and Astrophysics, Herschel Special Issue, in press as a Letter; 5 pages
Submitted: 2010-05-12
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.
[170]  [pdf] - 1026737
PEP: first Herschel probe of dusty galaxy evolution up to z~3
Comments: Accepted for publication in the A&A Herschel first results Special Issue
Submitted: 2010-05-10
We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.
[171]  [pdf] - 1026715
Far-infrared properties of submillimeter and optically faint radio galaxies
Comments: Accepted for publication in the A&A Herschel special issue. 5 pages, 3 figures
Submitted: 2010-05-07
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 submm 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 um and 160 um multi-wavelength catalogs. About half of the galaxies in our sample are detected with PACS. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 um (only the upper limits for the OFRGs) flux densities with a single modified (beta=1.5) black body function. The median dust temperature of our SMG sample is T=36+/-8K while for our OFRG sample it is T=47+/-3K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005) found a dust temperature of T=36+/-7K for a large sample of SMGs assuming the validity of the FIR/radio correlation. The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we find <q>=2.17+/-0.19, a slightly lower value than that observed in local systems. The median IR luminosities of SMGs and OFRGs are 4.6*10^12 Lsun and 2.6*10^12 Lsun, respectively. We note that for both samples the IR luminosity estimates from the radio part of the spectral energy distribution are accurate, while estimates from the mid-IR are considerably (x3) more uncertain. Our observations confirm the remarkably high luminosities of SMGs and thus imply median star-formation rates of 960Msun yr^-1 for SMGs with S(850um)>5mJy and 460Msun yr^-1 for SMGs with S(850um)>2mJy, assuming a Chabrier IMF and no dominant AGN contribution to the far-infrared luminosity.
[172]  [pdf] - 1026710
The first Herschel view of the mass-SFR link in high-z galaxies
Comments: Accepted for publication in the A&A Herschel Special issue
Submitted: 2010-05-06
We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the 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 star formation rate (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 alpha=-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 massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.
[173]  [pdf] - 1026706
Dissecting the cosmic infra-red background with Herschel/PEP
Comments: Accepted for publication on the A&A Herschel Special Issue
Submitted: 2010-05-06
The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 um) 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 ~3.0 mJy at 100 um and ~5.7 mJy at 160 um. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 um. 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 (L(IR)~1e11 Lsun), whose SEDs resemble late-spiral galaxies, peaking at ~130 um 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 um, resolving ~45% and ~52% of the CIB, respectively. When stacking 24 um 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<=1.0, with 160 um sources found at higher redshift than 100 um ones.
[174]  [pdf] - 1026705
The far-infrared/radio correlation as probed by Herschel
Comments: A&A Herschel Special Issue, in press as a Letter. 5 pages
Submitted: 2010-05-06
We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.
[175]  [pdf] - 1026704
Herschel and SCUBA-2 imaging and spectroscopy of a bright, lensed submillimetre galaxy at z = 2.3
Comments: A&A Herschel Special Issue, in press as a Letter. 5 pages, 3 colour figures
Submitted: 2010-05-06
We present a detailed analysis of the far-IR properties of the bright, lensed, z = 2.3, SMG, SMM J2135-0102, using new observations with Herschel, SCUBA-2 and the VLA. These data allow us to constrain the galaxy's SED and show that it has an intrinsic rest-frame 8-1000um luminosity, L(bol), of (2.3 +/- 0.2) x 10^12 L(sun) and a likely SFR of ~400 M(sun)/yr. The galaxy sits on the far-IR/radio correlation for far-IR-selected galaxies. At ~>70um, the SED can be described adequately by dust components with T(d) ~ 30 and 60K. Using SPIRE's Fourier Transform Spectrometer we report a detection of the [CII] 158um cooling line. If the [CII], CO and far-IR continuum arise in photo-dissociation regions, we derive a characteristic gas density, n ~ 10^3 cm^-3, and a far-UV radiation field, G_0, 10^3x stronger than the Milky Way. L([CII])/L(bol) is significantly higher than in local ULIRGs but similar to the values found in local star-forming galaxies and starburst nuclei. This is consistent with SMM J2135-0102 being powered by starburst clumps distributed across ~2 kpc, evidence that SMGs are not simply scaled-up ULIRGs. Our results show that SPIRE's FTS has the ability to measure the redshifts of distant, obscured galaxies via the blind detection of atomic cooling lines, but it will not be competitive with ground-based CO-line searches. It will, however, allow detailed study of the integrated properties of high-redshift galaxies, as well as the chemistry of their ISM, once more suitably bright candidates have been found.
[176]  [pdf] - 1017495
Spectroscopic Identifications of SWIRE sources in ELAIS-N1
Comments: 11 pages, 8 figures, accepted for publication in MNRAS
Submitted: 2009-09-11, last modified: 2010-03-04
We present the largest spectroscopic follow-up performed in SWIRE ELAIS-N1. We were able to determine redshifts for 289 extragalactic sources. The values of spectroscopic redshifts of the latter have been compared with the estimated values from our photometric redshift code with very good agreement between the two for both galaxies and quasars. Six of the quasars are hyperluminous infrared galaxies all of which are broad line AGN. We have performed emission line diagnostics for 30 sources in order to classify them into star-forming, Seyferts, composite and LINER and compare the results to the predictions from our SED template fitting methods and mid-IR selection methods.
[177]  [pdf] - 1018048
The Herschel ATLAS
Comments: submitted to PASP
Submitted: 2009-10-22
The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.
[178]  [pdf] - 250728
AGN dust tori: the X-ray-infrared connection
Comments: 15 pages, 18 figures, accepted for publication MNRAS
Submitted: 2009-05-27, last modified: 2009-06-16
We have combined the CLASXS Chandra survey in Lockman with the Spitzer SWIRE survey data to study the X-ray-infrared connection for AGN. The sample consists of 401 X-ray-sources, of which 306 are detected by Spitzer, and a further 257 AGN candidates detected through their dust torus, but not by Chandra. For X-ray sources the X-ray hardness ratio has been modelled in terms of a power-law with absorption N(H). The optical and infrared data have been modelled in terms of our well-established optical galaxy and QSO templates, and infrared templates based on radiative transfer models. Our estimate of the N(H) distribution is consistent with other studies, but we do find a higher proportion of low absorption objects at z < 0.5 than at z > 0.5. While we find only one X-ray AGN with N(H) > 10^{24} cm^{-2}, we argue that 10 objects with torus luminosity apparently exceeding the bolometric X-ray to 3 \mu m luminosity are strong candidates for being heavily absorbed in X-rays. We also estimate that at least half of the infrared-detected AGN dust tori which are undetected in X-rays are likely to be Compton thick. Our estimate of the total number of Compton-thick objects corresponds to > 20%$ of the combined SWIRE-CLASXS sample (and with an upper limit of 39 %). The range of dust covering factors is 1-100 %, with a mean of 40 %, ie a Type 2 fraction of 40 %. Measured by the ratio of dust torus luminosity to X-ray or (for Type 1 objects) optical luminosity, the covering factor appears to decrease towards intermediate AGN luminosity, in contradiction to estimates based on ratios of narrow-line and broad-line spectra, but may increase again at low AGN luminosity.
[179]  [pdf] - 21557
Mid-infrared spectroscopy of infrared-luminous galaxies at z~0.5-3
Comments: 32 pages, 24 figures, 15 tables, accepted for publication in MNRAS
Submitted: 2009-02-19
We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.
[180]  [pdf] - 1326507
Testing the galaxy cluster mass-observable relations at z = 1 with XMM-Newton and Chandra observations of XLSSJ022403.9-041328
Comments: 11 pages, 11 figures. Accepted for publication in MNRAS. Revised to include treatment of XMM PSF. Conclusions unchanged
Submitted: 2007-09-14, last modified: 2008-04-17
We present an analysis of deep XMM-Newton and Chandra observations of the z=1.05 galaxy cluster XLSSJ022403.9-041328 (hereafter XLSSC 029), detected in the XMM-Newton large scale structure survey. Density and temperature profiles of the X-ray emitting gas were used to perform a hydrostatic mass analysis of the system. This allowed us to measure the total mass and gas fraction in the cluster and define overdensity radii R500 and R2500. The global properties of XLSSC 029 were measured within these radii and compared with those of the local population. The gas mass fraction was found to be consistent with local clusters. The mean metal abundance was 0.18 +0.17 -0.15 Zsol, with the cluster core regions excluded, consistent with the predicted and observed evolution. The properties of XLSSC 029 were then used to investigate the position of the cluster on the M-kT, YX-M, and LX-M scaling relations. In all cases the observed properties of XLSSC 029 agreed well with the simple self-similar evolution of the scaling relations. This is the first test of the evolution of these relations at z > 1 and supports the use of the scaling relations in cosmological studies with distant galaxy clusters.
[181]  [pdf] - 900202
The XMM-LSS survey: the Class 1 cluster sample over the initial 5 square degrees and its cosmological modelling
Comments: Accepted for publication by MNRAS. Full resolution images as well as additional cluster data are available through a dedicated database at
Submitted: 2007-09-13
We present a sample of 29 galaxy clusters from the XMM-LSS survey over an area of some 5deg2 out to a redshift of z=1.05. The sample clusters, which represent about half of the X-ray clusters identified in the region, follow well defined X-ray selection criteria and are all spectroscopically confirmed. For all clusters, we provide X-ray luminosities and temperatures as well as masses. The cluster distribution peaks around z=0.3 and T =1.5 keV, half of the objects being groups with a temperature below 2 keV. Our L-T(z) relation points toward self-similar evolution, but does not exclude other physically plausible models. Assuming that cluster scaling laws follow self-similar evolution, our number density estimates up to z=1 are compatible with the predictions of the concordance cosmology and with the findings of previous ROSAT surveys. Our well monitored selection function allowed us to demonstrate that the inclusion of selection effects is essential for the correct determination of the evolution of the L-T relation, which may explain the contradictory results from previous studies. Extensive simulations show that extending the survey area to 10deg2 has the potential to exclude the non-evolution hypothesis, but that constraints on more refined ICM models will probably be limited by the large intrinsic dispersion of the L-T relation. We further demonstrate that increasing the dispersion in the scaling laws increases the number of detectable clusters, hence generating further degeneracy [in addition to sigma8, Omega_m, L(M,z) and T(M,z)] in the cosmological interpretation of the cluster number counts. We provide useful empirical formulae for the cluster mass-flux and mass-count-rate relations as well as a comparison between the XMM-LSS mass sensitivity and that of forthcoming SZ surveys.
[182]  [pdf] - 1326505
The XMM-LSS catalogue: X-ray sources and associated optical data. Version I
Comments: 13 pages, 7 figures and 11 tables (fig. 1 and 6 are enclosed in reduced resolution), MNRAS Latex, accepted by MNRAS
Submitted: 2007-08-24
Following the presentation of the XMM-LSS X-ray source detection package by Pacaud et al., we provide the source lists for the first 5.5 surveyed square degrees. The catalogues pertain to the [0.5-2] and [2-10] keV bands and contain in total 3385 point-like or extended sources above a detection likelihood of 15 in either band. The agreement with deep logN-logS is excellent. The main parameters considered are position, countrate, source extent with associated likelihood values. A set of additional quantities such as astrometric corrections and fluxes are further calculated while errors on the position and countrate are deduced from simulations. We describe the construction of the band-merged catalogue allowing rapid sub-sample selection and easy cross-correlation with external multi-wavelength catalogues. A small optical CFHTLS multi-band subset of objects is associated wich each source along with an X-ray/optical overlay. We make the full X-ray images available in FITS format. The data are available at CDS and, in a more extended form, at the Milan XMM-LSS database.
[183]  [pdf] - 1327005
Structure detection in the D1 CFHTLS deep field using accurate photometric redshifts: a benchmark
Comments: 14 pages, 22 additionnal jpeg figures, accepted in A&A
Submitted: 2007-02-16
We investigate structures in the D1 CFHTLS deep field in order to test the method that will be applied to generate homogeneous samples of clusters and groups of galaxies in order to constrain cosmology and detailed physics of groups and clusters. Adaptive kernel technique is applied on galaxy catalogues. This technique needs none of the usual a-priori assumptions (luminosity function, density profile, colour of galaxies) made with other methods. Its main drawback (decrease of efficiency with increasing background) is overcame by the use of narrow slices in photometric redshift space. There are two main concerns in structure detection. One is false detection and the second, the evaluation of the selection function in particular if one wants "complete" samples. We deal here with the first concern using random distributions. For the second, comparison with detailed simulations is foreseen but we use here a pragmatic approach with comparing our results to GalICS simulations to check that our detection number is not totally at odds compared to cosmological simulations. We use XMM-LSS survey and secured VVDS redshifts up to z~1 to check individual detections. We show that our detection method is basically capable to recover (in the regions in common) 100% of the C1 XMM-LSS X-ray detections in the correct redshift range plus several other candidates. Moreover when spectroscopic data are available, we confirm our detections, even those without X-ray data.
[184]  [pdf] - 86167
A 15um selected sample of high-z starbursts and AGNs
Comments: 4 pages, 3 figures, to appear in proceedings of 'Spitzer IR Diagnostics Conference, Nov 14-16, 2005'
Submitted: 2006-10-25
We report results from a Spitzer GO-1 program of IRS spectroscopy of a large sample of Luminous Infrared Galaxies and quasars selected from the European Large Area ISO Survey (ELAIS) which have a wide multiwavelength coverage, including ISOCAM, ISOPHOT, IRAC and MIPS (from SWIRE), and optical photometry. We present the sample selection and results from the IRS spectroscopy.
[185]  [pdf] - 83628
XMM-LSS discovery of a z=1.22 galaxy cluster
Comments: 8 pages, 6 figures, MNRAS accepted
Submitted: 2006-07-18
We present details of the discovery of XLSSJ022303.0-043622, a z=1.2 cluster of galaxies. This cluster was identified from its X-ray properties and selected as a z>1 candidate from its optical/near-IR characteristics in the XMM Large-Scale Structure Survey (XMM-LSS). It is the most distant system discovered in the survey to date. We present ground-based optical and near IR observations of the system carried out as part of the XMM-LSS survey. The cluster has a bolometric X-ray luminosity of 1.1 +/- 0.7 x 10^44 erg/s, fainter than most other known z>1 X-ray selected clusters. In the optical it has a remarkably compact core, with at least a dozen galaxies inside a 125 kpc radius circle centred on the X-ray position. Most of the galaxies within the core, and those spectroscopically confirmed to be cluster members, have stellar masses similar to those of massive cluster galaxies at low redshift. They have colours comparable to those of galaxies in other z>1 clusters, consistent with showing little sign of strong ongoing star formation. The bulk of the star formation within the galaxies appears to have ceased at least 1.5 Gyr before the observed epoch. Our results are consistent with massive cluster galaxies forming at z>1 and passively evolving thereafter. We also show that the system is straightforwardly identified in Spitzer/IRAC 3.6 and 4.5 micron data obtained by the SWIRE survey emphasising the power and utility of joint XMM and Spitzer searches for the most distant clusters.
[186]  [pdf] - 83373
The XMM-LSS Survey: A well controlled X-ray cluster sample over the D1 CFHTLS area
Comments: 12 pages 5 figures, MNRAS accepted. The paper with full resolution cluster images is available at
Submitted: 2006-07-10
We present the XMM-LSS cluster catalogue corresponding to the CFHTLS D1 area. The list contains 13 spectroscopically confirmed, X-ray selected galaxy clusters over 0.8 deg2 to a redshift of unity and so constitutes the highest density sample of clusters to date. Cluster X-ray bolometric luminosities range from 0.03 to 5x10^{44} erg/s. In this study, we describe our catalogue construction procedure: from the detection of X-ray cluster candidates to the compilation of a spectroscopically confirmed cluster sample with an explicit selection function. The procedure further provides basic X-ray products such as cluster temperature, flux and luminosity. We detected slightly more clusters with a (0.5-2.0 keV) X-ray fluxes of >2x10^{-14} erg/s/cm^{-2} than we expected based on expectations from deep ROSAT surveys. We also present the Luminosity-Temperature relation for our 9 brightest objects possessing a reliable temperature determination. The slope is in good agreement with the local relation, yet compatible with a luminosity enhancement for the 0.15 < z< 0.35 objects having 1 < T < 2 keV, a population that the XMM-LSS is identifying systematically for the first time. The present study permits the compilation of cluster samples from XMM images whose selection biases are understood. This allows, in addition to studies of large-scale structure, the systematic investigation of cluster scaling law evolution, especially for low mass X-ray groups which constitute the bulk of our observed cluster population. All cluster ancillary data (images, profiles, spectra) are made available in electronic form via the XMM-LSS cluster database.
[187]  [pdf] - 83380
The XMM Large Scale Structure survey: The X-ray pipeline and survey selection function
Comments: 16 pages, 11 figures, MNRAS accepted. The paper with full resolution cluster images is available at
Submitted: 2006-07-10
We present the X-ray pipeline developed for the purpose of the cluster search in the XMM-LSS survey. It is based on a two-stage procedure via a dedicated handling of the Poisson nature of the signal: (1) source detection on multi-resolution wavelet filtered images; (2) source analysis by means of a maximum likelihood fit to the photon images. The source detection efficiency and characterisation are studied through extensive Monte-Carlo simulations. This led us to define two samples of extended sources: the C1 class that is uncontaminated, and the less restrictive C2 class that allows for 50% contamination. The resulting predicted selection function is presented and the comparison to the current XMM-LSS confirmed cluster sample shows very good agreement. We arrive at average predicted source densities of about 7 C1 and 12 C2 per deg2, which is higher than any available wide field X-ray survey. We finally notice a substantial deviation of the predicted redshift distribution for our samples from the one obtained using the usual assumption of a flux limited sample.
[188]  [pdf] - 880911
The XMM Large Scale Structure Survey: Properties and Two-Point Angular Correlations of Point-like Sources
Comments: A&A in press; High resolution version at
Submitted: 2006-07-07
We analyze X-ray sources detected over 4.2 pseudo-contiguous sq. deg. in the 0.5-2 keV and 2-10 keV bands down to fluxes of 2x10^{-15} and 8x10^{-15} erg/s/cm^2 respectively, as part of the XMM Large Scale Structure Survey. The logN-logS in both bands shows a steep slope at bright fluxes, but agrees well with other determinations below ~2x10^{-14} erg/s/cm^2. The detected sources resolve close to 30 per cent of the X-ray background in the 2-10 keV band. We study the two-point angular clustering of point sources using nearest neighbours and correlation function statistics and find a weak, positive signal for ~1130 sources in the 0.5-2 keV band, but no correlation for ~400 sources in the 2-10 keV band below scales of 100 arcsec. A sub-sample of ~200 faint sources with hard X-ray count ratios, that is likely to be dominated by obscured AGN, does show a positive signal with the data allowing for a large scaling of the angular correlation length, but only at the ~2 (3) sigma level, based on re-sampling (Poisson) statistics. We discuss possible implications and emphasize the importance of wider, complete surveys in order to fully understand the large scale structure of the X-ray sky.
[189]  [pdf] - 74892
The XMM Large-Scale Structure Survey: An initial sample of galaxy groups and clusters to a redshift z<0.6
Comments: 18 pages, MNRAS accepted
Submitted: 2005-07-29
We present X-ray and optical spectroscopic observations of twelve galaxy groups and clusters identified within the XMM Large-Scale Structure (LSS) survey. Groups and clusters are selected as extended X-ray sources from a 3.5 deg2 XMM image mosaic above a flux limit 8e-15 ergs/s/cm2 in the [0.5-2] keV energy band. Deep BVRI images and multi-object spectroscopy confirm each source as a galaxy concentration located within the redshift interval 0.29<z<0.56. We combine line-of-sight velocity dispersions with the X-ray properties of each structure computed from a two-dimensional surface brightness model and a single temperature fit to the XMM spectral data. The resulting distribution of X-ray luminosity, temperature and velocity dispersion indicate that the XMM-LSS survey is detecting low-mass clusters and galaxy groups to redshifts z < 0.6. Confirmed systems display little or no evidence for X-ray luminosity evolution at a given X-ray temperature compared to lower redshift X-ray group and cluster samples. A more complete understanding of these trends will be possible with the compilation of a statistically complete sample of galaxy groups and clusters anticipated within the continuing XMM-LSS survey.
[190]  [pdf] - 72861
Rotation of the cluster of galaxies A2107
Comments: 7 pages, 9 figures gzipped tar file. to be published in MNRAS
Submitted: 2005-05-05
We present indications of rotation in the galaxy cluster A2107 by a method that searches for the maximum gradient in the velocity field in a flat disk-like model of a cluster. Galaxies from cumulative sub-samples containing more and more distant members from the cluster centre, are projected onto an axis passing through the centre and we apply a linear regression model on the projected distances $x$ and the line-of-sight velocities $V$. The axis with the maximum linear correlation coefficient $r_{max} = \max [ r(V,x) ]$ defines the direction of the maximum velocity gradient, and consequently it presents the major axis of the apparently elliptical cluster. Because the effects of rotation are subtle, we put strong emphasis on the estimation of the uncertainties of the results by implementing different bootstrap techniques. We have found the rotational effects are more strongly expressed from distances $0.26 \div 0.54$ Mpc from the cluster centre. The total virial mass of the cluster is $(3.2\pm0.6)\times10^{14} {\cal M}_{\sun}$, while the virial mass, corrected for the rotation, is $(2.8\pm0.5) \times 10^{14}{\cal M}_{\sun}$.
[191]  [pdf] - 71618
Batch discovery of nine z~1 clusters using x-ray and K or R,z' images
Comments: MNRAS, in press. Paper with full-resolution figures can be found at
Submitted: 2005-03-10
We present results of an initial search for clusters of galaxies at z~1 and above, using data from 2.9 square degrees of XMM-Newton images. By selecting weak potentially extended X-ray sources with faint or no identifications in deep, ground-based optical imaging, we have constructed a starting sample of 19 high redshift cluster candidates. Near-IR and R,z' imaging of these fields identified nine of them as high redshift systems. Six of these were confirmed spectroscopically, three at z~1.0 and the other three in 0.8<z<0.92 range. The remaining three systems have solid photometric evidence to be at z_phot ~ 0.8, 1.0 and 1.3. The present sample significantly increases the number of such clusters. The measured density of z>~1 clusters, after discarding "low" redshift systems at z <~ 0.92 is about 1.7 /deg2 (with 68 % confidence interval equal to [1.0,2.9]) for f_X >~2.5 10^{-15} ergs/cm2/s ([0.5-2] keV) and this is a lower limit, having screened not all potential z~1 candidate clusters. Coordinates, x-ray measures and evidence for nine x-ray selected high redshift clusters is given.
[192]  [pdf] - 56664
The XMM-LSS survey. Survey design and first results
Comments: Published in Journal of Cosmology and Astroparticle Physics
Submitted: 2003-05-12, last modified: 2004-11-01
We have designed a medium deep large area X-ray survey with XMM - the XMM Large Scale Structure survey, XMM-LSS - with the scope of extending the cosmological tests attempted using ROSAT cluster samples to two redshift bins between 0<z<1 while maintaining the precision of earlier studies. Two main goals have constrained the survey design: the evolutionary study of the cluster-cluster correlation function and of the cluster number density. The results are promising and, so far, in accordance with our predictions as to the survey sensitivity and cluster number density. The feasibility of the programme is demonstrated and further X-ray coverage is awaited in order to proceed with a truly significant statistical analysis. (Abridged)
[193]  [pdf] - 880735
Galaxy evolution in clusters up to z=1.0
Comments: MNRAS, in press. 17 pages, 8 figures. Data available at
Submitted: 2004-05-31
We present a combined study of the colour-magnitude relation, colour distribution and luminosity function of a sample of 24 clusters at redshifts 0.3<z<1. The sample is largely composed of X-ray selected/detected clusters. Most of the clusters at redshifts z<0.6 display X-ray luminosity or richness typical of poor clusters or groups, rather than the more typical, massive clusters studied in literature at redshifts z>~0.3. All our clusters, including groups, display a colour-magnitude relation consistent with a passively evolving stellar population formed at a redshift z_f>~2, in accordance with observed galaxy populations in more massive clusters studied at comparable redshifts. Colours and luminosity functions show that the cluster galaxy population is consistent with the presence of at least two components: old systems formed at high redshift that have evolved passively from that epoch, together with a galaxy population displaying more recent star formation. The former population forms at ~2< z_f <~5, the latter at redshifts z<1. A model in which stars do not evolve is clearly rejected by both by the colour of reddest galaxies and by the characteristic luminosity m* measures. All clusters (with one possible exception) are detected independently by an almost three dimensional optical search employing sky position and colour - this despite the primary X-ray selection and low X-ray flux/optical richness displayed by the majority of the sample.
[194]  [pdf] - 65094
Colour-detected clusters in the XMM Large Scale Structure Survey
Comments: Contributed to XXXIX Rencontres de Moriond "Exploring the Universe. Contents and Structure of the Universe", La Thuile, Aosta, Italy, March 28 - April 4, 2004
Submitted: 2004-05-28
We present first results on the use of a colour space filter for detecting galaxy clusters at cosmological redshifts in the XMM Large Scale Structure Survey. All clusters studied, but one, are successfully colour--detected in spite of their large redshift (0.3<z<1.0), X-ray selection and intrinsic low richness (R=0 or below). We experimentally show that the cluster redshift can be derived, with good accuracy, from the colour of the red sequence, at 0.06<z<1 using a sample of about 160 clusters.
[195]  [pdf] - 64904
Galaxy evolution in clusters from z=1 to z=0
Comments: To appear in the proceedings of IAU Colloquium 195: "Outskirts of galaxy clusters: intense life in the suburbs", Torino, Italy, March 12-16, 2004
Submitted: 2004-05-20
The XMM--LSS project is detecting distant clusters of low mass, quite comparable in mass to the ones in the local universe. This allows a direct comparison of galaxy properties at different redshifts in ``similar'' clusters. We present here first results on the evolution of the reddest galaxies in ~25 clusters/groups at ~0.3< z <~1.0 and for the whole galaxy population in the same clusters. The emerging picture from the current study is that the counterparts of present day clusters tend to show two or more distinct populations : a relatively old (z_f>2-5) population evolving passively together with a younger population, ostensibly responsible for the apparent brightening of the characteristic magnitudes, m*.
[196]  [pdf] - 56665
The XMM-LSS Survey. First high redshift galaxy clusters: relaxed and collapsing systems
Comments: 11 pages, 8 figures; revised version accepted in A&A. The PDF article with full resolution figures is available from
Submitted: 2003-05-12, last modified: 2004-05-11
We present five newly found galaxy clusters at z>0.6 from the XMM Large-Scale Structure Survey (XMM-LSS). All five objects are extended X-ray sources in the XMM images. For three of them we have sufficient spectroscopically confirmed member galaxies that an estimate of the velocity dispersion is possible: XLSSC 001 at z=0.613 and \sigma_V=867^{+80}_{-60} km/s, XLSSC 002 at z=0.772 and \sigma_V=524^{+267}_{-116} km/s and XLSSC 003 at z=0.839 and \sigma_V=780^{+137}_{-75} km/s. These three clusters have X-ray bolometric luminosities L_X \sim 1-3 \times 10^{44} erg/s and temperatures 2-4 keV, and consequently are less massive than previously known clusters at similar redshifts, but nevertheless they follow the low redshift scaling relations between L_X, T and \sigma_V, within the limits of the measurement errors. One of the clusters, XLSSC 004, is detected independently as an overdensity of galaxies of a colour R-z'=1.4 that matches the redshift of the central galaxy z=0.87, although it cannot unambiguously be confirmed by the spectroscopic observations alone. The highest redshift candidate cluster pertaining to this paper, XLSSC 005, is most likely a double cluster complex at a redshift around unity, associated with an extended X-ray source with probable substructure.
[197]  [pdf] - 50046
Abell 1451 and 1RXS J131423.6-251521: a multi-wavelength study of two dynamically perturbed clusters of galaxies
Comments: 16 pages, 12 figures, some in colour; A&A accepted
Submitted: 2002-06-24
We present results from optical, X-ray and radio observations of two X-ray bright (L_X ~ 10^{45} erg/s) galaxy clusters. A1451 is at redshift z=0.1989 and has line-of-sight velocity dispersion sigma_v = 1330 km/s as measured from 57 cluster galaxies. It has regular X-ray emission without signs of substructure, a Gaussian velocity distribution, lack of a cooling flow region and significant deviations from the observed scaling laws between luminosity, temperature and velocity dispersion, indicating a possible merging shock. There is only one spectroscopically confirmed cluster radio galaxy, which is close to the X-ray peak. 1RXS J131423.6-251521 has z=0.2474 and sigma_v = 1100 km/s from 37 galaxies. There are two distinct galaxy groups with a projected separation of \~700 kpc. The velocity histogram is bi-modal with a redshift-space separation of ~1700 km/s, and the X-ray emission is double peaked. Although there are no spectroscopically confirmed cluster radio galaxies, we have identified a plausible relic source candidate.
[198]  [pdf] - 45028
Cosmology with Galaxy Clusters in the XMM Large-Scale Structure Survey
Comments: 11 LaTex pages, including 14 figures and 1 table. Accepted version to appear in A&A. Contains new section and figure on systematics. Also available at
Submitted: 2001-09-28, last modified: 2002-05-29
The upcoming XMM Large Scale Structure Survey (XMM-LSS) will ultimately provide a unique mapping of the distribution of X-ray sources in a contiguous 64 sq. deg. region. In particular, it will provide the 3-dimensional location of about 900 galaxy clusters out to a redshift of about 1. We study the prospects that this cluster catalogue offers for measuring cosmological parameters. We use the Press-Schechter formalism to predict the counts of clusters and their X-ray properties in several CDM models. We compute the detection efficiency of clusters, using realistic simulations of XMM X-ray images, and study how it differs from a conventional flux limit. We compute the expected correlation function of clusters using the extended halo model, and show that it is expected to evolve very little out to z~2, once the selection function of the survey is taken into account. The shape and the amplitude of the correlation function can be used to brake degeneracies present when cluster counts alone are considered. Ignoring systematic uncertainties, the combination of cluster counts evolution and of the correlation function yields measurements of Omega_m, sigma_8 and Gamma with a precision of about 15%, 10% and 35%, respectively, in a LCDM model. This combination will also provide a consistency check for the LCDM model, and a discrimination between this model and the OCDM model. The XMM-LSS will therefore provide important constraints on cosmological parameters, complementing that from other methods such as the Cosmic Microwave Background. We discuss how these constraints are affected by instrumental systematics and by the uncertainties in the scaling relations of clusters.
[199]  [pdf] - 47611
The XMM Large Scale Structure Survey: Scientific Motivation and First Observations
Comments: 7 pages, 6 figures (2 in colour). Proceedings of the Conference "New Visions of the X-ray Universe in the XMM-Newton and Chandra Era" 26-30 November 2001, ESTEC, NL
Submitted: 2002-02-06, last modified: 2002-04-19
Thanks to its unrivalled sensitivity and large field of view, XMM potentially occupies a leading position as a survey instrument. We present cosmological arguments in favour of a medium-sensitivity, large-scale structure survey with XMM, using galaxy clusters as tracers of the cosmic network. We show how this has motivated the definition of a concrete survey, the XMM Large-Scale Structure Survey (XMM-LSS), which will cover 64 square degrees with a sensitivity about 1000 times better than that of the ROSAT All-Sky Survey. We present our predictions for cluster counts based on the Press-Schechter formalism and detailed X-ray image simulations, and show how they agree with the cluster statistics from recent ROSAT cluster surveys. We also present the extensive multi-wavelength follow-up associated with XMM-LSS, as well as the first observations from the programme.
[200]  [pdf] - 1232750
In.XS: project for a future spaceborne hard X-ray all-sky survey
Comments: 10 pages, presented at 46th SPIE meeting (proc. vol. 4497), official website:
Submitted: 2001-09-27
The latest all-sky survey in hard X-ray band was performed by the HEAO-1 satellite (13-80 keV) with an angular resolution of 24x48 arcmin. A diffuse hard X-Ray background (HXB) was detected between 3 and 50 keV. The main scientific goal of In.XS is to resolve a large fraction of this HXB into individual sources. As no distortion by Compton up-scattering is seen in the spectrum of the microwave background, the hard X-ray background is believed to be mainly due to point sources. Type I Active Galactic Nuclei (AGN) have softer X-ray spectra than the hard X-ray background, so other sources must be considered, like faint Type II or absorbed AGN. These could be distinguished through hard X-ray spectroscopic or hardness ratio observations. Here we present In.XS - a mission concept designed to conduct the first imaging all-sky hard X-ray (2-80 keV) survey. The angular resolution of nearly 1 arcmin and good sensitivity at high energies is provided by the latest multilayer focussing mirrors, with semiconductor-based (GaAs) arrays of detectors. We also describe the mission operations, and how the all-sky survey will be complemented by follow-up pointed observations of selected fields. The good angular resolution will allow correlations and identification with objects seen at other wavelengths. In addition, since a large fraction of the Type II AGN luminosity is emitted in the hard X-ray band, this survey will provide a large unbiased sample of the AGN population. This may provide constraints on AGN evolution through the possible observation of a turnover in deep field source statistics.
[201]  [pdf] - 41008
Comparison of source detection procedures for XMM-Newton images
Comments: 20 pages, 16 figures. Accepted for publication in A&A
Submitted: 2001-02-16
Procedures based on current methods to detect sources in X-ray images are applied to simulated XMM images. All significant instrumental effects are taken into account, and two kinds of sources are considered -- unresolved sources represented by the telescope PSF and extended ones represented by a b-profile model. Different sets of test cases with controlled and realistic input configurations are constructed in order to analyze the influence of confusion on the source analysis and also to choose the best methods and strategies to resolve the difficulties. In the general case of point-like and extended objects the mixed approach of multiresolution (wavelet) filtering and subsequent detection by SExtractor gives the best results. In ideal cases of isolated sources, flux errors are within 15-20%. The maximum likelihood technique outperforms the others for point-like sources when the PSF model used in the fit is the same as in the images. However, the number of spurious detections is quite large. The classification using the half-light radius and SExtractor stellarity index is succesful in more than 98% of the cases. This suggests that average luminosity clusters of galaxies (L_[2-10] ~ 3x10^{44} erg/s) can be detected at redshifts greater than 1.5 for moderate exposure times in the energy band below 5 keV, provided that there is no confusion or blending by nearby sources. We find also that with the best current available packages, confusion and completeness problems start to appear at fluxes around 6x10^{-16} erg/s/cm^2 in [0.5-2] keV band for XMM deep surveys.
[202]  [pdf] - 38644
Probing the gravitational potential of a nearby lensing cluster Abell 2104
Comments: 13 pages, 10 postscript figures, use aa.sty, to appear in Astronomy and Astrophysics
Submitted: 2000-10-16
The cluster Abell 2104 is one of the lowest redshift clusters (z=0.153) known to have a gravitational lensing arc. We present detailed analysis of the cluster properties such as the gravitational potential using the X-ray data from ROSAT (HRI) and ASCA, as well as optical imaging and spectroscopic data from the CFHT. The cluster is highly luminous in the X-ray with a bolometric luminosity of Lx ~3x10^{45} ergs/s and a high gas temperature of ~10.4 keV. The X-ray emission extending out to at least a radius of 1.46 Mpc, displays significant substructure. The total mass deduced from the X-ray data under the assumption of hydrostatic equilibrium and isothermal gas, is found to be M_tot(r<1.46Mpc) ~(8.0 +/- 0.8)x10^{14} solar mass. The gas fraction within a radius of 1.46 Mpc is ~5-10%. The cluster galaxy velocity distribution has a dispersion of 1200 +/- 200 km/s with no obvious evidence for substructure. The total mass within 1.46 Mpc, deduced from Jean's equation using the observed galaxy number density distribution and velocity dispersion, is found to be \~6.8x10^{14} solar mass to ~2.6x10^{15} solar mass marginally consistent with the X-ray deduced total mass.
[203]  [pdf] - 103954
References for Galaxy Clusters Database
Comments: Presented at ADASSVIII conference, 1-5 Nov 1998, Urbana-Champaign, IL
Submitted: 1998-11-20
A bibliographic database will be constructed with the purpose to be a general tool for searching references for galaxy clusters. The structure of the database will be completely different from the available now databases as NED, SIMBAD, LEDA. Search based on hierarchical keyword system will be performed through web interfaces from numerous bibliographic sources -- journal articles, preprints, unpublished results and papers, theses, scientific reports. Data from the very beginning of the extragalactic research will be included as well. References for galaxy clusters (RGC) is continuation of a previous project for collecting all published information for galaxy clusters. There are neither previous attempts nor projects for compiling such database for galaxy clusters and our effort will be to include into the database all of the available bibliography through our system of keywords until the end of 1999. Now over 3000 entries are included into the preliminary version of the database.
[204]  [pdf] - 101730
Correlation Function of Superclusters of Galaxies
Comments: 19 pages, 3 figures, 4 tables. To appear in 1998 ApJ, 506, No. 2 (Oct 20)
Submitted: 1998-06-10
We present a study of the two-point correlation function of superclusters of galaxies. The largest catalogs are used. The results show negligible correlation less than 0.1-0.2 for separations up to 500-600 h^{-1} Mpc. Small correlations are obtained using various estimates and samples. Seemingly there are no structures of superclusters of galaxies.