Source author record

M. Vaccari

M. Vaccari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

120works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

120 published item(s)

preprint2025arXiv

Atomic hydrogen reservoirs in quiescent galaxies at z = 0.4

Context. Based on Local Universe observations, quiescent galaxies (QGs) host lower to no HI compared to star-forming galaxies (SFGs), but no constraints have been derived so far at higher redshift (z>0.1). Understanding whether QGs can retain significant HI reservoirs at higher z is crucial to refine quenching and gas accretion models and to constrain overall star formation efficiency at different epochs. Aims. We aim to probe HI in candidate QGs at intermediate redshifts (z=0.36) and to understand whether there exists a class of QGs retaining consistent HI reservoirs and which parameters (dust content, stellar mass, Dn4000, morphology, environment) effectively capture HI-rich QGs. Methods. We perform 21-cm spectral line stacking on MIGHTEE-HI data at z=0.36, targeting two different samples of QGs, defined by means of a color-selection criterion and a spectroscopic criterion based on Dn4000, respectively. We also perform stacking on subsamples of the spectroscopically-selected quiescent sample to investigate the correlation between the HI content and other galaxy properties. Results. We find that QGs with an IR counterpart (i.e., dusty galaxies) are found to host a substantial HI content, on average just 40% lower than SFGs. In contrast, color-selected QGs still hold HI, but lower than SFGs by a factor 3. Among dusty objects, we find morphology to have a mild impact on the atomic gas content, with spirals hosting approximately 15-30% more HI than spheroids. Environmental effects are also present, with low-density regions hosting galaxies that are HI-richer than in high-density ones, by approximately 30% for spirals and 60% for spheroids. We suggest that, in general, HI content is driven by several factors, including quenching mechanisms and ISM enrichment processes. Also, quiescent galaxies - and especially dusty systems - seem to yield HI more consistently than in the Local Universe.

preprint2022arXiv

A structure function analysis of VST-COSMOS AGN

We present our sixth work in a series dedicated to variability studies of active galactic nuclei (AGN) based on the survey of the COSMOS field by the VLT Survey Telescope (VST). Its 54 r-band visits over 3.3 yr and single-visit depth of 24.6 r-band mag make this dataset a valuable scaled-down version that can help forecast the performance of the Rubin Observatory Legacy Survey of Space and Time (LSST). This work is centered on the analysis of the structure function (SF) of VST-COSMOS AGN, investigating possible differences in its shape and slope related to how the AGN were selected, and explores possible connections between the ensemble variability of AGN and black-hole mass, accretion rate, bolometric luminosity, redshift, and obscuration of the source. Given its features, our dataset opens up the exploration of samples ~2 mag fainter than most of the literature to date. We identify several samples of AGN - 677 in total - obtained by a variety of selection techniques which partly overlap. Our analysis compares results for the various samples. We split each sample in two based on the median of the physical property of interest, and analyze differences in the shape and slope of the SF, and possible causes. While the shape of the SF does not change with depth, it is highly affected by the type of AGN (unobscured/obscured) included in the sample. Where a linear region can be identified, we find that the variability amplitude anticorrelates with accretion rate and bolometric luminosity, consistent with previous literature on the topic, while no dependence on black-hole mass emerges from this study. With its longer baseline and denser and more regular sampling, the LSST will allow an improved characterization of the SF and its dependencies on the mentioned physical properties over much larger AGN samples.

preprint2022arXiv

Deep ASKAP EMU Survey of the GAMA23 field: Properties of radio sources

We present the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the Galaxy and Mass Assembly (GAMA)-23h field. The survey was carried out at 887.5 MHz and covers a 83 square degree field. We imaged the calibrated visibility data, taken as part of the Evolutionary Mapping of Universe (EMU) Early Science Programme, using the latest version of the ASKAPSoft pipeline. The final mosaic has an angular resolution of 10 arcsec and a central rms noise of around 38 $μ$Jy beam$^{-1}$. The derived radio source catalogue has 39812 entries above a peak flux density threshold of 5$σ$. We searched for the radio source host galaxy counterparts using the GAMA spectroscopic (with an i-band magnitude limit of 19.2 mag) and multi-wavelength catalogues that are available as part of the collaboration. We identified hosts with GAMA spectroscopic redshifts for 5934 radio sources. We describe the data reduction, imaging, and source identification process, and present the source counts. Thanks to the wide area covered by our survey, we obtain very robust counts down to 0.2 mJy. ASKAP's exceptional survey speed, providing efficient, sensitive and high resolution mapping of large regions of the sky in conjunction with the multi-wavelength data available for the GAMA23 field, allowed us to discover 63 giant radio galaxies. The data presented here demonstrate the excellent capabilities of ASKAP in the pre-SKA era.

preprint2022arXiv

MIGHTEE: the nature of the radio-loud AGN population

We study the nature of the faint radio source population detected in the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Early Science data in the COSMOS field, focusing on the properties of the radio-loud active galactic nuclei (AGN). Using the extensive multi-wavelength data available in the field, we are able to classify 88 per cent of the 5223 radio sources in the field with host galaxy identifications as AGN (35 per cent) or star-forming galaxies (54 per cent). We select a sample of radio-loud AGN with redshifts out to $z \sim 6$ and radio luminosities $10^{20} < \textrm{L}_{1.4~\textrm{GHz}} / \textrm{W Hz}^{-1} < 10^{27}$ and classify them as high-excitation and low-excitation radio galaxies (HERGs and LERGs). The classification catalogue is released with this work. We find no significant difference in the host galaxy properties of the HERGs and LERGs in our sample. In contrast to previous work, we find that the HERGs and LERGs have very similar Eddington-scaled accretion rates; in particular we identify a population of very slowly accreting AGN that are formally classified as HERGs at these low radio luminosities, where separating into HERGs and LERGs possibly becomes redundant. We investigate how black hole mass affects jet power, and find that a black hole mass $\gtrsim 10^{7.8}~\textrm{M}_\odot$ is required to power a jet with mechanical power greater than the radiative luminosity of the AGN ($L_\textrm{mech}/L_\textrm{bol} > 1$). We discuss that both a high black hole mass and black hole spin may be necessary to launch and sustain a dominant radio jet.

preprint2022arXiv

The LOFAR Two-metre Sky Survey -- V. Second data release

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data.

preprint2022arXiv

The star-formation rates of QSOs

We examine the far-IR properties of a sample of 5391 optically selected QSOs in the 0.5<z<2.65 redshift range down to log[nuLnu,2500 (erg/s)]>44.7, using SPIRE data from Herschel-ATLAS. We split the sample in a grid of 74 luminosity-redshift bins and compute the average optical-infrared spectral energy distribution (SED) in each bin. By normalising an intrinsic AGN template to the AGN optical power (at 5100A) we decompose the total infrared emission (L_IR; 8-1000um) into an AGN (L_IR,AGN) and star-forming component (L_IR,SF). We find that the AGN contribution to L_IR increases as a function of AGN power which manifests as a reduction of the `far-IR bump' in the average QSO SEDs. We note that L_IR,SF does not correlate with AGN power; the mean star formation rates (SFRs) of AGN host galaxies are a function of redshift only and they range from ~6 Msun/yr at z~0 to a plateau of <200 Msun/yr at z~2.6. Our results indicate that the accuracy of far-IR emission as a proxy for SFR decreases with increasing AGN luminosity. We show that, at any given redshift, observed trends between infrared luminosity (whether monochromatic or total) and AGN power (in the optical or X-rays) can be explained by a simple model which is the sum of two components: (A) the infrared emission from star-formation, uncorrelated with AGN power and (B) the infrared emission from AGN, directly proportional to AGN power in the optical or X-rays.

preprint2021arXiv

COSMOS2020: A panchromatic view of the Universe to $z\sim10$ from two complementary catalogs

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25<i<27$ reach a precision of $5\,\%$. Finally, these results are discussed in the context of previous, current, and future surveys in the COSMOS field. Compared to COSMOS2015, reaches the same photometric redshift precision at almost one magnitude deeper. Both photometric catalogs and their photometric redshift solutions and physical parameters will be made available through the usual astronomical archive systems (ESO Phase 3, IPAC IRSA, and CDS).

preprint2020arXiv

A random forest-based selection of optically variable AGN in the VST-COSMOS field

The survey of the COSMOS field by the VLT Survey Telescope is an appealing testing ground for variability studies of active galactic nuclei (AGN). With 54 r-band visits over 3.3 yr and a single-visit depth of 24.6 r-band mag, the dataset is also particularly interesting in the context of performance forecasting for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This work is the fifth in a series dedicated to the development of an automated, robust, and efficient methodology to identify optically variable AGN, aimed at deploying it on future LSST data. We test the performance of a random forest (RF) algorithm in selecting optically variable AGN candidates, investigating how the use of different AGN labeled sets (LSs) and features sets affects this performance. We define a heterogeneous AGN LS and choose a set of variability features and optical and near-infrared colors based on what can be extracted from LSST data. We find that an AGN LS that includes only Type I sources allows for the selection of a highly pure (91%) sample of AGN candidates, obtaining a completeness with respect to spectroscopically confirmed AGN of 69% (vs. 59% in our previous work). The addition of colors to variability features mildly improves the performance of the RF classifier, while colors alone prove less effective than variability in selecting AGN as they return contaminated samples of candidates and fail to identify most host-dominated AGN. We observe that a bright (r < 21 mag) AGN LS is able to retrieve candidate samples not affected by the magnitude cut, which is of great importance as faint AGN LSs for LSST-related studies will be hard to find and likely imbalanced. We estimate a sky density of 6.2 million AGN for the LSST main survey down to our current magnitude limit.

preprint2020arXiv

A Wide-area GMRT 610~MHz survey of ELAIS N1 field

In this paper we present a wide-area 610 MHz survey of the ELAIS\,N1 field with the GMRT, covering an area of 12.8 deg$^2$ at a resolution of 6 arcsec and with an rms noise of $\sim 40$ $μ$Jy beam$^{-1}$. This is equivalent to $\sim 20$ $μ$Jy beam$^{-1}$ rms noise at 1.4 GHz for a spectral index of $-0.75$. The primary goal of the survey was to study the polarised sky at sub-mJy flux densities at $<$ GHz frequencies. In addition, a range of other science goals, such as investigations in to the nature of the low-frequency $μ$Jy source populations and alignments of radio jets. A total of 6,400 sources were found in this region, the vast majority of them compact. The sample jointly detected by GMRT at 610 MHz and by VLA FIRST at 1.4\,GHz has a median spectral index of $-0.85 \pm 0.05$ and a median 610 MHz flux density of 4.5 mJy. This region has a wealth of ancillary data which is useful to characterize the detected sources. The multi-wavelength cross matching resulted optical/IR counterparts to $\sim 90$ per~cent of the radio sources, with a significant fraction having at least photometric redshift. Due to the improved sensitivity of this survey over preceding ones, we have discovered six giant radio sources (GRS), with three of them at $z \sim 1$ or higher. This implies that the population of GRS may be more abundant and common than known to date and if true this has implications for the luminosity function and the evolution of radio sources. We have also identified several candidate extended relic sources.

preprint2020arXiv

Extending the variability selection of active galactic nuclei in the W-CDF-S and SERVS/SWIRE region

Variability has proven to be a powerful tool to detect active galactic nuclei (AGN) in multi-epoch surveys. The new-generation facilities expected to become operational in the next few years will mark a new era in time-domain astronomy and their wide-field multi-epoch campaigns will favor extensive variability studies. We present our analysis of AGN variability in the second half of the VST survey of the Wide Chandra Deep Field South (W-CDF-S), performed in the r band and covering a 2 sq. deg. area. The analysis complements a previous work, in which the first half of the area was investigated. We provide a reliable catalog of variable AGN candidates, which will be critical targets in future variability studies. We selected a sample of optically variable sources and made use of infrared data from the Spitzer mission to validate their nature by means of color-based diagnostics. We obtain a sample of 782 AGN candidates among which 12 are classified as supernovae, 54 as stars, and 232 as AGN. We estimate a contamination $\lesssim 20\%$ and a completeness $\sim 38\%$ with respect to mid-infrared selected samples.

preprint2020arXiv

SCUBA-2 overdensities associated with candidate protoclusters selected from Planck data

We measure the 850-$μ$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$μ$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$μ$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$μ$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters.

preprint2019arXiv

Cosmic evolution of star-forming galaxies to $z \simeq 1.8$ in the faint low-frequency radio source population

We study the properties of star-forming galaxies selected at 610 MHz with the GMRT in a survey covering $\sim$1.86 deg$^2$ down to a noise of $\sim$7.1\,$μ$Jy / beam. These were identified by combining multiple classification diagnostics: optical, X-ray, infrared and radio data. Of the 1685 SFGs from the GMRT sample, 496 have spectroscopic redshifts whereas 1189 have photometric redshifts. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4 GHz radio luminosity ratio $\rm{q_{IR}}$, decreases with increasing redshift: $\rm{q_{IR}\,=\,2.86\pm0.04(1\,+\,z)^{-0.20\pm0.02}}$ out to $z \sim 1.8$. We use the $\rm{V/V_{max}}$ statistic to quantify the evolution of the co-moving space density of the SFG sample. Averaged over luminosity our results indicate $\rm{\langle V/V_{max} \rangle}$ to be $\rm{0.51\,\pm\, 0.06}$, which is consistent with no evolution in overall space density. However we find $\rm V/V_{max}$ to be a function of radio luminosity, indicating strong luminosity evolution with redshift. We explore the evolution of the SFGs radio luminosity function by separating the source into five redshift bins and comparing to theoretical model predictions. We find a strong redshift trend that can be fitted with a pure luminosity evolution of the form $\rm{L_{610\,MHz}\,\propto\,(\,1+\,z)^{(2.95\pm0.19)-(0.50\pm0.15)z}}$. We calculate the cosmic SFR density since $\rm{z \sim 1.5}$ by integrating the parametric fits of the evolved 610\,MHz luminosity function. Our sample reproduces the expected steep decline in the star formation rate density since $\rm{z\,\sim\,1}$.

preprint2019arXiv

The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$σ$ goal of 70 $μ$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

preprint2016arXiv

Candidate Gravitationally Lensed Dusty Star-forming Galaxies in the Herschel Wide Area Surveys

We present a list of candidate gravitationally lensed dusty star-forming galaxies (DSFGs) from the HerMES Large Mode Survey (HeLMS) and the Herschel Stripe 82 Survey (HerS). Together, these partially overlapping surveys cover 372 deg$^{2}$ on the sky. After removing local spiral galaxies and known radio-loud blazars, our candidate list of lensed DSFGs is composed of 77 sources with 500 $μ$m flux densities ($S_{500}$) greater than 100 mJy. Such sources are dusty starburst galaxies similar to the first bright Sub Millimeter Galaxies (SMGs) discovered with SCUBA. We expect a large fraction of this list to be strongly lensed, with a small fraction made up of bright SMG-SMG mergers that appear as Hyper-Luminous Infrared Galaxies (HyLIRGs; $\rm L_{IR}>10^{13} L_{\odot}$). Thirteen of the 77 candidates have spectroscopic redshifts from CO spectroscopy with ground-based interferometers, putting them at $z>1$ and well above the redshift of the foreground lensing galaxies. The surface density of our sample is 0.21 $\pm$ 0.03 deg$^{-2}$. We present follow-up imaging of a few of the candidates confirming their lensing nature. The sample presented here is an ideal tool for higher resolution imaging and spectroscopic observations to understand detailed properties of starburst phenomena in distant galaxies.

preprint2016arXiv

HELP: star formation as function of galaxy environmentwith Herschel

The Herschel Extragalactic Legacy Project (HELP) brings together a vast range of data from many astronomical observatories. Its main focus is on the Herschel data, which maps dust obscured star formation over 1300 deg$^2$. With this unprecedented combination of data sets, it is possible to investigate how the star formation vs stellar mass relation (main-sequence) of star-forming galaxies depends on environment. In this pilot study we explore this question between 0.1 < z < 3.2 using data in the COSMOS field. We estimate the local environment from a smoothed galaxy density field using the full photometric redshift probability distribution. We estimate star formation rates by stacking the SPIRE data from the Herschel Multi-tiered Extragalactic Survey (HerMES). Our analysis rules out the hypothesis that the main-sequence for star-forming systems is independent of environment at 1.5 < z < 2, while a simple model in which the mean specific star formation rate declines with increasing environmental density gives a better description. However, we cannot exclude a simple hypothesis in which the main-sequence for star-forming systems is independent of environment at z < 1.5 and z > 2. We also estimate the evolution of the star formation rate density in the COSMOS field and our results are consistent with previous measurements at z < 1.5 and z > 2 but we find a $1.4^{+0.3}_{-0.2}$ times higher peak value of the star formation rate density at $z \sim 1.9$.

preprint2016arXiv

HELP: XID+, The Probabilistic De-blender for Herschel SPIRE maps

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{μ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 https://github.com/H-E-L-P/XID_plus. 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.

preprint2016arXiv

The galaxy-halo connection in the VIDEO Survey at 0.5<z<1.7

We present a series of results from a clustering analysis of the first data release of the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is the only survey currently capable of probing the bulk of stellar mass in galaxies at redshifts corresponding to the peak of star formation on degree scales. Galaxy clustering is measured with the two-point correlation function, which is calculated using a non parametric kernel based density estimator. We use our measurements to investigate the connection between the galaxies and the host dark matter halo using a halo occupation distribution methodology, deriving bias, satellite fractions, and typical host halo masses for stellar masses between $10^{9.35}M_{\odot}$ and $10^{10.85}M_{\odot}$, at redshifts $0.5<z<1.7$. Our results show typical halo mass increasing with stellar mass (with moderate scatter) and bias increasing with stellar mass and redshift consistent with previous studies. We find the satellite fraction increased towards low redshifts, increasing from $\sim 5\%$ at $z\sim 1.5$, to $\sim 20\%$ at $z\sim 0.6$, also increasing for lower mass galaxies. We combine our results to derive the stellar mass to halo mass ratio for both satellites and centrals over a range of halo masses and find the peak corresponding to the halo mass with maximum star formation efficiency to be $ \sim 2 \times10^{12} M_{\odot}$ over cosmic time, finding no evidence for evolution.

preprint2016arXiv

The Herschel-ATLAS: a sample of 500μm-selected lensed galaxies over 600 square degrees

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500μm extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500μm and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

preprint2016arXiv

The Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole field

We present the Spitzer-IRAC/MIPS Extragalactic survey (SIMES) in the South Ecliptic Pole (SEP) field. The large area covered (7.7 deg$^2$), together with one of the lowest Galactic cirrus emissions in the entire sky and a very extensive coverage by Spitzer, Herschel, Akari, and GALEX, make the SIMES field ideal for extragalactic studies. The elongated geometry of the SIMES area ($\approx$4:1), allowing for a significant cosmic variance reduction, further improves the quality of statistical studies in this field. Here we present the reduction and photometric measurements of the Spitzer/IRAC data. The survey reaches a depth of 1.93 and 1.75 $μ$Jy (1$σ$) at 3.6 and 4.5 $μ$m, respectively. We discuss the multiwavelength IRAC--based catalog, completed with optical, mid-- and far--IR observations. We detect 341,000 sources with F$_{3.6μm} \geq 3σ$. Of these, 10% have an associated 24 $μ$m counterpart, while 2.7% have an associated SPIRE source. We release the catalog through the NASA/IPAC Infrared Science Archive (IRSA). Two scientific applications of these IRAC data are presented in this paper: first we compute integral number counts at 3.6 $μ$m. Second, we use the [3.6]--[4.5] color index to identify galaxy clusters at z$>$1.3. We select 27 clusters in the full area, a result consistent with previous studies at similar depth.

preprint2015arXiv

A new search for variability-selected active galaxies within the VST SUDARE-VOICE survey: the Chandra Deep Field South and the SERVS-SWIRE area

This work makes use of the VST observations to select variable sources. We use also the IR photometry, SED fitting and X-ray information where available to confirm the nature of the AGN candidates. The IR data, available over the full survey area, allow to confirm the consistency of the variability selection with the IR color selection method, while the detection of variability may prove useful to detect the presence of an AGN in IR selected starburst galaxies.

preprint2015arXiv

HerMES: Disentangling active galactic nuclei and star formation in the radio source population

We separate the extragalactic radio source population above ~50 uJy into active galactic nuclei (AGN) and star-forming sources. The primary method of our approach is to fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IRAC and MIPS and Herschel/SPIRE photometry, of 380 radio sources in the Extended Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN and star-forming contributions to their infrared emission. With the inclusion of other AGN diagnostics such as X-ray luminosity, Spitzer/IRAC colours, radio spectral index and the ratio of star-forming total infrared flux to k-corrected 1.4 GHz flux density, qIR, we determine whether the radio emission in these sources is powered by star formation or by an AGN. The majority of these radio sources (60 per cent) show the signature of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent are hybrid systems for which the radio emission is being powered by star formation. This implies that radio sources which have likely been selected on their star formation have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding a strong contribution to the source counts from pure star-forming sources, we find that hybrid sources constitute 20-65 per cent of the sources. This result suggests that hybrid sources have a significant contribution, along with sources that do not host a detectable AGN, to the observed flattening of the source counts at ~1mJy for the extragalactic radio source population.

preprint2015arXiv

Spitzer Imaging of Strongly-Lensed Herschel-Selected Dusty Star Forming Galaxies

We present the rest-frame optical spectral energy distribution and stellar masses of six Herschel- selected gravitationally lensed dusty, star-forming galaxies (DSFGs) at 1 < z < 3. These galaxies were first identified with Herschel/SPIRE imaging data from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). The targets were observed with Spitzer/IRAC at 3.6 and 4.5um. Due to the spatial resolution of the IRAC observations at the level of 2 arcseconds, the lensing features of a background DSFG in the near-infrared are blended with the flux from the foreground lensing galaxy in the IRAC imaging data. We make use of higher resolution Hubble/WFC3 or Keck/NIRC2 Adaptive Optics imaging data to fit light profiles of the foreground lensing galaxy (or galaxies) as a way to model the foreground components, in order to successfully disentangle the foreground lens and background source flux densities in the IRAC images. The flux density measurements at 3.6 and 4.5um, once combined with Hubble/WFC3 and Keck/NIRC2 data, provide important constraints on the rest-frame optical spectral energy distribution of the Herschel-selected lensed DSFGs. We model the combined UV- to millimeter-wavelength SEDs to establish the stellar mass, dust mass, star-formation rate, visual extinction, and other parameters for each of these Herschel-selected DSFGs. These systems have inferred stellar masses in the range 8 x 10^10 to 4 x 10^11 Msun and star-formation rates of around 100 Msun yr-1. This puts these lensed sub-millimeter systems well above the SFR-M* relation observed for normal star-forming galaxies at similar redshifts. The high values of SFR inferred for these systems are consistent with a major merger-driven scenario for star formation.

preprint2015arXiv

Star Formation in Herschel's Monsters versus Semi-Analytic Models

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.

preprint2015arXiv

SUDARE-VOICE variability-selection of Active Galaxies in the Chandra Deep Field South and the SERVS/SWIRE region

One of the most peculiar characteristics of Active Galactic Nuclei (AGN) is their variability over all wavelengths. This property has been used in the past to select AGN samples and is foreseen to be one of the detection techniques applied in future multi-epoch surveys, complementing photometric and spectroscopic methods. In this paper, we aim to construct and characterise an AGN sample using a multi-epoch dataset in the r band from the SUDARE-VOICE survey. Our work makes use of the VST monitoring program of an area surrounding the Chandra Deep Field South to select variable sources. We use data spanning a six month period over an area of 2 square degrees, to identify AGN based on their photometric variability. The selected sample includes 175 AGN candidates with magnitude r < 23 mag. We distinguish different classes of variable sources through their lightcurves, as well as X-ray, spectroscopic, SED, optical and IR information overlapping with our survey. We find that 12% of the sample (21/175) is represented by SN. Of the remaining sources, 4% (6/154) are stars, while 66% (102/154) are likely AGNs based on the available diagnostics. We estimate an upper limit to the contamination of the variability selected AGN sample of about 34%, but we point out that restricting the analysis to the sources with available multi-wavelength ancillary information, the purity of our sample is close to 80% (102 AGN out of 128 non-SN sources with multi-wavelength diagnostics). Our work thus confirms the efficiency of the variability selection method in agreement with our previous work on the COSMOS field; in addition we show that the variability approach is roughly consistent with the infrared selection.

preprint2015arXiv

The faint radio source population at 15.7 GHz - II. Multi-wavelength properties

A complete, flux density limited sample of 96 faint ($> 0.5$ mJy) radio sources is selected from the 10C survey at 15.7 GHz in the Lockman Hole. We have matched this sample to a range of multi-wavelength catalogues, including SERVS, SWIRE, UKIDSS and optical data; multi-wavelength counterparts are found for 80 of the 96 sources and spectroscopic redshifts are available for 24 sources. Photometric reshifts are estimated for the sources with multi-wavelength data available; the median redshift of the sample is 0.91 with an interquartile range of 0.84. Radio-to-optical ratios show that at least 94 per cent of the sample are radio loud, indicating that the 10C sample is dominated by radio galaxies. This is in contrast to samples selected at lower frequencies, where radio-quiet AGN and starforming galaxies are present in significant numbers at these flux density levels. All six radio-quiet sources have rising radio spectra, suggesting that they are dominated by AGN emission. These results confirm the conclusions of Paper I that the faint, flat-spectrum sources which are found to dominate the 10C sample below $\sim 1$ mJy are the cores of radio galaxies. The properties of the 10C sample are compared to the SKADS Simulated Skies; a population of low-redshift starforming galaxies predicted by the simulation is not found in the observed sample.

preprint2015arXiv

The HerMES sub-millimetre local and low-redshift luminosity functions

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$\,μ$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$\,μ$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 $Φ_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $Φ_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,μ$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$.

preprint2014arXiv

A Far-Infrared Spectroscopic Survey of Intermediate Redshift (Ultra) Luminous Infrared Galaxies

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.

preprint2014arXiv

Dust Spectral Energy Distributions of Nearby Galaxies: an Insight from the Herschel Reference Survey

We gather infrared (IR) photometric data from 8 to 500 microns (Spitzer, WISE, IRAS and Herschel) for all of the HRS galaxies. Draine & Li (2007) models are fit to the data from which the stellar contribution has been carefully removed. We find that our photometric coverage is sufficient to constrain all of the models parameters and that a strong constraint on the 20-60 microns range is mandatory to estimate the relative contribution of the photo-dissociation regions to the IR SED. The SED models tend to systematically under-estimate the observed 500 microns flux densities, especially for low mass systems. We provide the output parameters for all of the galaxies: the minimum intensity of the interstellar radiation field (ISRF), the fraction of PAH, the relative contribution of PDR and evolved stellar population to the dust heating, the $M_{dust}$ and the $L_{IR}$. For a subsample of gas-rich galaxies, we analyze the relations between these parameters and the integrated properties of galaxies, such as $M_*$, SFR, metallicity, H$α$ and H-band surface brightness, and the FUV attenuation. A good correlation between the fraction of PAH and the metallicity is found implying a weakening of the PAH emission in galaxies with low metallicities. The intensity of the IRSF and the H-band and H$α$ surface brightnesses are correlated, suggesting that the diffuse dust component is heated by both the young stars in star forming regions and the diffuse evolved population. We use these results to provide a new set of IR templates calibrated with Herschel observations on nearby galaxies and a mean SED template to provide the z=0 reference for cosmological studies. For the same purpose, we put our sample on the SFR-$M_*$ diagram. The templates are compared to the most popular IR SED libraries, enlightening a large discrepancy between all of them in the 20-100 microns range.

preprint2014arXiv

Herschel Observations of Far-Infrared Cooling Lines in intermediate Redshift (Ultra)-luminous Infrared Galaxies

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.

preprint2014arXiv

Linking the X-ray and infrared properties of star-forming galaxies at z<1.5

We present the most complete study to date of the X-ray emission from star-formation in high redshift (median z=0.7; z<1.5), IR-luminous (L_IR=10^10-10^13 L_sun) galaxies detected by Herschel's PACS and SPIRE instruments. For our purpose we take advantage of the deepest X-ray data to date, the Chandra deep fields (North and South). Sources which host AGN are removed from our analysis by means of multiple AGN indicators. We find an AGN fraction of 18+/-2 per cent amongst our sample and note that AGN entirely dominate at values of log[L_X/L_IR]>-3 in both hard and soft X-ray bands. From the sources which are star-formation dominated, only a small fraction are individually X-ray detected and for the bulk of the sample we calculate average X-ray luminosities through stacking. We find an average soft X-ray to infrared ratio of log[L_SX/L_IR]=-4.3 and an average hard X-ray to infrared ratio of log[L_HX/L_IR]=-3.8. We report that the X-ray/IR correlation is approximately linear through the entire range of L_IR and z probed and, although broadly consistent with the local (z<0.1) one, it does display some discrepancies. We suggest that these discrepancies are unlikely to be physical, i.e. due to an intrinsic change in the X-ray properties of star-forming galaxies with cosmic time, as there is no significant evidence for evolution of the L_X/L_IR ratio with redshift. Instead they are possibly due to selection effects and remaining AGN contamination. We also examine whether dust obscuration in the galaxy plays a role in attenuating X-rays from star-formation, by investigating changes in the L_X/L_IR ratio as a function of the average dust temperature. We conclude that X-rays do not suffer any measurable attenuation in the host galaxy.

preprint2014arXiv

PACS photometry of the Herschel Reference Survey - Far-infrared/sub-millimeter colours as tracers of dust properties in nearby galaxies

We present Herschel/PACS 100 and 160 micron integrated photometry for the 323 galaxies in the Herschel Reference Survey (HRS), a K-band-, volume-limited sample of galaxies in the local Universe. Once combined with the Herschel/SPIRE observations already available, these data make the HRS the largest representative sample of nearby galaxies with homogeneous coverage across the 100-500 micron wavelength range. In this paper, we take advantage of this unique dataset to investigate the properties and shape of the far-infrared/sub-millimeter spectral energy distribution in nearby galaxies. We show that, in the stellar mass range covered by the HRS (8<log(M*/Msun)<12), the far-infrared/sub-millimeter colours are inconsistent with a single modified black-body having the same dust emissivity index beta for all galaxies. In particular, either beta decreases, or multiple temperature components are needed, when moving from metal-rich/gas-poor to metal-poor/gas-rich galaxies. We thus investigate how the dust temperature and mass obtained from a single modified black-body depend on the assumptions made on beta. We show that, while the correlations between dust temperature, galaxy structure and star formation rate are strongly model dependent, the dust mass scaling relations are much more reliable, and variations of beta only change the strength of the observed trends.

preprint2014arXiv

The Far-Infrared-Radio Correlation in MS0451-03

We present a multi-wavelength analysis of star-forming galaxies in the massive cluster MS0451.6-0305 at z $\sim$ 0.54 to shed new light on the evolution of the far-infrared-radio relationship in distant rich clusters. We have derived total infrared luminosities for a spectroscopically confirmed sample of cluster and field galaxies through an empirical relation based on $Spitzer$ MIPS 24 $μ$m photometry. The radio flux densities were measured from deep Very Large Array 1.4 GHz radio continuum observations. We find the ratio of far-infrared to radio luminosity for galaxies in an intermediate redshift cluster to be $q_{\rm FIR}$ = 1.80$\pm$0.15 with a dispersion of 0.53. Due to the large intrinsic dispersion, we do not find any observable change in this value with either redshift or environment. However, a higher percentage of galaxies in this cluster show an excess in their radio fluxes when compared to low redshift clusters ($27^{+23}_{-13}\%$ to $11\%$), suggestive of a cluster enhancement of radio-excess sources at this earlier epoch. In addition, the far-infrared-radio relationship for blue galaxies, where $q_{\rm FIR}$ = 2.01$\pm$0.14 with a dispersion of 0.35, is consistent with the predicted value from the field relationship, although these results are based on a sample from a single cluster.

preprint2014arXiv

The temperature dependence of the far-infrared-radio correlation in the Herschel-ATLAS

We use 10,387 galaxies from the Herschel Astrophysical TeraHertz Large Area Survey (H-ATLAS) to probe the far-infrared radio correlation (FIRC) of star forming galaxies as a function of redshift, wavelength, and effective dust temperature. All of the sources in our 250 μm-selected sample have spectroscopic redshifts, as well as 1.4 GHz flux density estimates measured from the Faint Images of the Radio Sky at Twenty centimetres (FIRST) survey. This enables us to study not only individual sources, but also the average properties of the 250 μm selected population using median stacking techniques. We find that individual sources detected at $\geq 5σ$ in both the H-ATLAS and FIRST data have logarithmic flux ratios (i.e. FIRC $q_λ$ parameters) consistent with previous studies of the FIRC. In contrast, the stacked values show larger $q_λ$, suggesting excess far-IR flux density/luminosity in 250μm selected sources above what has been seen in previous analyses. In addition, we find evidence that 250 μm sources with warm dust SEDs have a larger 1.4 GHz luminosity than the cooler sources in our sample. Though we find no evidence for redshift evolution of the monochromatic FIRC, our analysis reveals significant temperature dependence. Whilst the FIRC is reasonably constant with temperature at 100 μm, we find increasing inverse correlation with temperature as we probe longer PACS and SPIRE wavelengths. These results may have important implications for the use of monochromatic dust luminosity as a star formation rate indicator in star-forming galaxies, and in the future, for using radio data to determine galaxy star formation rates.

preprint2014arXiv

Tracing the cosmic growth of super massive black holes to z ~ 3 with Herschel

We study a sample of Herschel-PACS selected galaxies within the GOODS-South and the COSMOS fields in the framework of the PACS Evolutionary Probe (PEP) project. Starting from the rich multi-wavelength photometric data-sets available in both fields, we perform a broad-band Spectral Energy Distribution (SED) decomposition to disentangle the possible active galactic nucleus (AGN) contribution from that related to the host galaxy. We find that 37 per cent of the Herschel-selected sample shows signatures of nuclear activity at the 99 per cent confidence level. The probability to reveal AGN activity increases for bright ($L_{\rm 1-1000} > 10^{11} \rm L_{\odot}$) star-forming galaxies at $z>0.3$, becoming about 80 per cent for the brightest ($L_{\rm 1-1000} > 10^{12} \rm L_{\odot}$) infrared (IR) galaxies at $z \geq 1$. Finally, we reconstruct the AGN bolometric luminosity function and the super-massive black hole growth rate across cosmic time up to $z \sim 3$ from a Far-Infrared (FIR) perspective. This work shows general agreement with most of the panchromatic estimates from the literature, with the global black hole growth peaking at $z \sim 2$ and reproducing the observed local black hole mass density with consistent values of the radiative efficiency $ε_{\rm rad}$ ($\sim$0.07).

preprint2014arXiv

Variability-selected Active Galactic Nuclei in the VST-SUDARE/VOICE Survey of the COSMOS Field

Optical variability has proven to be an effective way of detecting AGNs in imaging surveys, lasting from weeks to years. In the present work we test its use as a tool to identify AGNs in the VST multi-epoch survey of the COSMOS field, originally tailored to detect supernova events. We make use of the multi-wavelength data provided by other COSMOS surveys to discuss the reliability of the method and the nature of our AGN candidates. Our selection returns a sample of 83 AGN candidates; based on a number of diagnostics, we conclude that 67 of them are confirmed AGNs (81% purity), 12 are classified as supernovae, while the nature of the remaining 4 is unknown. For the subsample of AGNs with some spectroscopic classification, we find that Type 1 are prevalent (89%) compared to Type 2 AGNs (11%). Overall, our approach is able to retrieve on average 15% of all AGNs in the field identified by means of spectroscopic or X-ray classification, with a strong dependence on the source apparent magnitude. In particular, the completeness for Type 1 AGNs is 25%, while it drops to 6% for Type 2 AGNs. The rest of the X-ray selected AGN population presents on average a larger r.m.s. variability than the bulk of non variable sources, indicating that variability detection for at least some of these objects is prevented only by the photometric accuracy of the data. We show how a longer observing baseline would return a larger sample of AGN candidates. Our results allow us to assess the usefulness of this AGN selection technique in view of future wide-field surveys.

preprint2013arXiv

An Overview of the Dwarf Galaxy Survey

The Dwarf Galaxy Survey (DGS) program is studying low-metallicity galaxies using 230h of far-infrared (FIR) and submillimetre (submm) photometric and spectroscopic observations of the Herschel Space Observatory and draws to this a rich database of a wide range of wavelengths tracing the dust, gas and stars. This sample of 50 galaxies includes the largest metallicity range achievable in the local Universe including the lowest metallicity (Z) galaxies, 1/50 Zsun, and spans 4 orders of magnitude in star formation rates. The survey is designed to get a handle on the physics of the interstellar medium (ISM) of low metallicity dwarf galaxies, especially on their dust and gas properties and the ISM heating and cooling processes. The DGS produces PACS and SPIRE maps of low-metallicity galaxies observed at 70, 100, 160, 250, 350, and 500 mic with the highest sensitivity achievable to date in the FIR and submm. The FIR fine-structure lines, [CII] 158 mic, [OI] 63 mic, [OI] 145 mic, [OIII] 88 mic, [NIII] 57 mic and [NII] 122 and 205 mic have also been observed with the aim of studying the gas cooling in the neutral and ionized phases. The SPIRE FTS observations include many CO lines (J=4-3 to J=13-12), [NII] 205 mic and [CI] lines at 370 and 609 mic. This paper describes the sample selection and global properties of the galaxies, the observing strategy as well as the vast ancillary database available to complement the Herschel observations. The scientific potential of the full DGS survey is described with some example results included.

preprint2013arXiv

Far-IR Emission From Dust-Obscured Galaxies

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3σ). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

preprint2013arXiv

HerMES: Candidate High-Redshift Galaxies Discovered with Herschel/SPIRE

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $μ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 $μ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\,μ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.

preprint2013arXiv

HerMES: Clusters of Dusty Galaxies Uncovered by Herschel and Planck

The potential for Planck to detect clusters of dusty, star-forming galaxies at z greater than 1 is tested by examining the Herschel-SPIRE images of Planck Early Release Compact Source Catalog (ERCSC) sources lying in fields observed by the HerMES survey. Of the 16 Planck sources that lie in the roughly 90 sq. deg. examined, we find that twelve are associated with single bright Herschel sources. The remaining four are associated with overdensities of Herschel sources, making them candidate clusters of dusty, starforming galaxies. We use complementary optical and NIR data for these clumps to test this idea, and find evidence for the presence of galaxy clusters in all four cases. We use photometric redshifts and red sequence galaxies to estimate the redshifts of these clusters, finding that they range from 0.8 to 2.3. These redshifts imply that the Herschel sources in these clusters, which contribute to the detected Planck flux, are forming stars very rapidly, with typical total cluster star formation rates greater than 1000Msun per yr. The high redshift clusters discovered in these observations are used to constrain the epoch of cluster galaxy formation, finding that the galaxies in our clusters are 1 to 1.5 Gy old at z about 1 to 2. Prospects for the discovery of further clusters of dusty galaxies are discussed, using not only all sky Planck surveys, but also deeper, smaller area, Herschel surveys.

preprint2013arXiv

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

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 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.

preprint2013arXiv

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

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.

preprint2013arXiv

HerMES: Point source catalogues from Herschel-SPIRE observations II

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

preprint2013arXiv

HerMES: The Contribution to the Cosmic Infrared Background from Galaxies Selected by Mass and Redshift

We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin^2) K-selected sources (K_AB < 24.0), split according to their rest-frame U - V vs. V - J colors into 72,216 star-forming and 9,034 quiescent galaxies, on maps from Spitzer/MIPS (24um), Herschel/SPIRE (250, 350, 500um), Herschel/PACS (100, 160um), and AzTEC (1100um). The fraction of the CIB resolved by our catalog is (69 $\pm$ 15)% at 24um, (78 $\pm$ 17)% at 70um, (58 $\pm$ 13)% at 100um, (78 $\pm$ 18)% at 160um, (80 $\pm$ 17)% at 250um, (69 $\pm$ 14)% at 350um, (65 $\pm$ 12)% at 500um, and (45 $\pm$ 8)% at 1100um. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at $λ$ < 200um is sourced predominantly from galaxies at z < 1, while at $λ$ > 200um the bulk originates from 1 < z < 2. Galaxies with stellar masses log(M/ M_sun)=9.5-11 are responsible for the majority of the CIB, with those in the log(M/ M_sun)=9.5-10 contributing mostly at $λ$ < 250um, and those in the log(M/ M_sun)=10.5-11 bin dominating at $λ$ > 350um. The contribution from galaxies in the log(M/ M_sun)=9.0-9.5 and log(M/ M_sun)=11.0-12.0 stellar mass bins contribute the least, both of order 5%, although the highest stellar-mass bin is a significant contributor to the luminosity density at z > 2. The luminosities of the galaxies responsible for the CIB shifts from a combination of "normal" and luminous infrared galaxies (LIRGs) at $λ$ < 160um, to LIRGs at 160um < $λ$ < 500um, to finally LIRGs and ultra-luminous infrared galaxies (ULIRGs) at $λ$ > 500um. Stacking analyses were performed with SIMSTACK (available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html) which accounts for possible biases due to clustering.

preprint2013arXiv

On the redshift evolution (0 < z < 4) of dust attenuation and of the total (UV+IR) star formation rate density

Using new homogeneous LFs in the FUV and in the FIR Herschel/PEP and Herschel/HerMES, we study the evolution of the dust attenuation with redshift. With this information in hand, we are able to estimate the redshift evolution of the total (FUV + FIR) star formation rate density SFRD_TOT. By integrating SFRD_TOT, we follow the mass building and analyze the redshift evolution of the stellar mass density (SMD). This letter aims at providing a complete view of star formation from the local universe to z = 4 and, using assumptions on earlier star formation history, compares this evolution to what was known before in an attempt to draw a homogeneous picture of the global evolution of star formation in galaxies. The main conclusions of this letter are: 1) the dust attenuation A_FUV is found to increase from z = 0 to z \sim 1.2 and then starts to decrease up to our last data point at z = 3.6; 2) the estimated SFRD confirms published results up to z = 2. At z > 2, we observe either a plateau or a small increase up to z = 3 and then a likely decrease up to z = 3.6; 3) the peak of A_FUV is delayed with respect to the plateau of SFRD_TOT and a likely origin might be found in the evolution of the bright ends of the FUV and FIR LFs; 4) using assumptions (namely exponential rise and linear rise with time) for the evolution of the star formation density from z = 3.6 to z_form = 10, we integrate SFRD_TOT and find a good agreement with the published SMDs.

preprint2013arXiv

Revealing the cold dust in low-metallicity environments: I - Photometry analysis of the Dwarf Galaxy Survey with Herschel

We present new photometric data from our Herschel Key Programme, the Dwarf Galaxy Survey (DGS), dedicated to the observation of the gas and dust in 48 low-metallicity environments. They were observed with PACS and SPIRE onboard Herschel at 70,100,160,250,350, and 500 microns. We focus on a systematic comparison of the derived FIR properties (FIR luminosity, dust mass, dust temperature and emissivity index) with more metal-rich galaxies and investigate the detection of a potential submm excess. The data reduction method is adapted for each galaxy to derive the most reliable photometry from the final maps. PACS flux densities are compared with the MIPS 70 and 160 microns bands. We use colour-colour diagrams and modified blackbody fitting procedures to determine the dust properties of the DGS galaxies. We also include galaxies from the Herschel KINGFISH sample, containing more metal-rich environments, totalling 109 galaxies. The location of the DGS galaxies on Herschel colour-colour diagrams highlights the differences in global environments of low-metallicity galaxies. The dust in DGS galaxies is generally warmer than in KINGFISH galaxies (T_DGS~32 K, T_KINGFISH~23 K). The emissivity index, beta, is ~1.7 in the DGS, but metallicity does not make a strong effect on beta. The dust-to-stellar mass ratio is lower in low-metallicity galaxies: M_dust/M_star~0.02% for the DGS vs 0.1% for KINGFISH. Per unit dust mass, dwarf galaxies emit ~6 times more in the FIR than higher metallicity galaxies. Out of the 22 DGS galaxies detected at 500 micron, 41% present an excess in the submm not explained by our dust SED model. The excess mainly appears in lower metallicity galaxies (12+log(O/H) < 8.3), and the strongest excesses are detected in the most metal-poor galaxies. We stress the need for observations longwards of the Herschel wavelengths to detect any submm excess appearing beyond 500 micron.

preprint2013arXiv

The Herschel census of infrared SEDs through cosmic time

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.

preprint2013arXiv

The Herschel PEP/HerMES Luminosity Function. I: Probing the Evolution of PACS selected Galaxies to z~4

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.

preprint2013arXiv

The rapid assembly of an elliptical galaxy of 400 billion solar masses at a redshift of 2.3

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.

preprint2013arXiv

The suppression of star formation by powerful active galactic nuclei

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.

preprint2012arXiv

Clustering of star-forming galaxies detected in mid-infrared with the Spitzer wide-area survey

We discuss the clustering properties of galaxies with signs of ongoing star formation detected by the Spitzer Space Telescope at 24mum band in the SWIRE Lockman Hole field. The sample of mid-IR-selected galaxies includes ~20,000 objects detected above a flux threshold of S24mum=310muJy. We adopt optical/near-IR color selection criteria to split the sample into the lower-redshift and higher-redshift galaxy populations. We measure the angular correlation function on scales of theta=0.01-3.5 deg, from which, using the Limber inversion along with the redshift distribution established for similarly selected source populations in the GOODS fields (Rodighiero et al. 2010), we obtain comoving correlation lengths of r0=4.98+-0.28 h^-1 Mpc and r0 =8.04+-0.69 h^-1 Mpc for the low-z (<z>=0.7) and high-z (<z>=1.7) subsamples, respectively. Comparing these measurements with the correlation functions of dark matter halos identified in the Bolshoi cosmological simulation (Klypin et al. 2011}, we find that the high-redshift objects reside in progressively more massive halos reaching Mtot>3e12 h^-1 Msun, compared to Mtot>7e11 h^-1 Msun for the low-redshift population. Approximate estimates of the IR luminosities based on the catalogs of 24mum sources in the GOODS fields show that our high-z subsample represents a population of "distant ULIRGs" with LIR>10^12Lsun, while the low-z subsample mainly consists of "LIRGs", LIR~10^11Lsun. The comparison of number density of the 24mum selected galaxies and of dark matter halos with derived minimum mass Mtot shows that only 20% of such halos may host star-forming galaxies.

preprint2012arXiv

Cosmic Evolution of Star-Formation Enhancement in Close Major-Merger Galaxy Pairs Since z = 1

The infrared (IR) emission of M_* galaxies (10^{10.4} < M_{star} < 10^{11.0} M_\sun) in galaxy pairs, derived using data obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement induced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies (SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation enhancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z=0 and z=1, this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by a factor of ~10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic medium IGM in the dark matter halos (DMH) hosting these pairs.

preprint2012arXiv

Evolution of the far-infrared luminosity functions in the Spitzer Wide-area Infrared Extragalactic Legacy Survey

We present new observational determination of the evolution of the rest-frame 70 and 160 micron and total infrared (TIR) galaxy luminosity functions (LFs) using 70 micron data from the Spitzer Wide-area Infrared Extragalactic Legacy Survey (SWIRE). The LFs were constructed for sources with spectroscopic redshifts only in the XMM-LSS and Lockman Hole fields from the SWIRE photometric redshift catalogue. The 70 micron and TIR LFs were constructed in the redshift range 0<z<1.2 and the 160 micron LF was constructed in the redshift range 0<z<0.5 using a parametric Bayesian and the vmax methods. We assume in our models, that the faint-end power-law index of the LF does not evolve with redshifts. We find the the double power-law model is a better representation of the IR LF than the more commonly used power-law and Gaussian model. We model the evolution of the FIR LFs as a function of redshift where where the characteristic luminosity, $L^\ast$ evolve as $\propto(1+z)^{α_\textsc{l}}$. The rest-frame 70 micron LF shows a strong luminosity evolution out to z=1.2 with alpha_l=3.41^{+0.18}_{-0.25}. The rest-frame 160 micron LF also showed rapid luminosity evolution with alpha_l=5.53^{+0.28}_{-0.23} out to z=0.5. The rate of evolution in luminosity is consistent with values estimated from previous studies using data from IRAS, ISO and Spitzer. The TIR LF evolves in luminosity with alpha_l=3.82^{+0.28}_{-0.16} which is in agreement with previous results from Spitzer 24 micron which find strong luminosity evolution. By integrating the LF we calculated the co-moving IR luminosity density out to z=1.2, which confirm the rapid evolution in number density of LIRGs and ULIRGs which contribute ~68^{+10}_{-07} % to the co-moving star formation rate density at z=1.2. Our results based on 70 micron data confirms that the bulk of the star formation at z=1 takes place in dust obscured objects.

preprint2012arXiv

HerMES: A Statistical Measurement of the Redshift Distribution of Herschel-SPIRE Sources Using the Cross-correlation Technique

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $μ$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Boötes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $μ$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $μ$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

preprint2012arXiv

HerMES: Candidate Gravitationally Lensed Galaxies and Lensing Statistics at Submillimeter Wavelengths

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

preprint2012arXiv

Submillimetre Photometry of 323 Nearby Galaxies from the Herschel Reference Survey

The Herschel Reference Survey (HRS) is a guaranteed time Herschel key project aimed at studying the physical properties of the interstellar medium in galaxies of the nearby universe. This volume limited, K-band selected sample is composed of galaxies spanning the whole range of morphological types (from ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster). We present flux density measurements of the whole sample of 323 galaxies of the HRS in the three bands of the Spectral and Photometric Imaging Receiver (SPIRE), at 250, 350 and 500 microns. Aperture photometry is performed on extended galaxies and point spread function (PSF) fitting on timeline data for unresolved objects; we carefully estimate errors and upper limits. The flux densities are found to be in good agreement with those of the HeViCS and KINGFISH key projects in all SPIRE bands, and of the Planck consortium at 350 and 550 microns, for the galaxies in common. This submillimetre catalogue of nearby galaxies is a benchmark for the study of the dust properties in the local universe, giving the zero redshift reference for any cosmological survey.

preprint2012arXiv

The AGN content in luminous IR galaxies at z\sim2 from a global SED analysis including Herschel data

We use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

preprint2012arXiv

The complex physics of dusty star-forming galaxies at high redshifts as revealed by Herschel and Spitzer

We combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $μ$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $Δ$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $Δ$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

preprint2012arXiv

The Herschel Exploitation of Local Galaxy Andromeda (HELGA). I: Global far-infrared and sub-mm morphology

We have obtained Herschel images at five wavelengths from 100 to 500 micron of a ~5.5x2.5 degree area centred on the local galaxy M31 (Andromeda), our nearest neighbour spiral galaxy, as part of the Herschel guaranteed time project "HELGA". The main goals of HELGA are to study the characteristics of the extended dust emission, focusing on larger scales than studied in previous observations of Andromeda at an increased spatial resolution, and the obscured star formation. In this paper we present data reduction and Herschel maps, and provide a description of the far-infrared morphology, comparing it with features seen at other wavelengths. We use high--resolution maps of the atomic hydrogen, fully covering our fields, to identify dust emission features that can be associated to M31 with confidence, distinguishing them from emission coming from the foreground Galactic cirrus. Thanks to the very large extension of our maps we detect, for the first time at far-infrared wavelengths, three arc-like structures extending out to ~21, ~26 and ~31 kpc respectively, in the south-western part of M31. The presence of these features, hosting ~2.2e6 Msol of dust, is safely confirmed by their detection in HI maps. Overall, we estimate a total dust mass of ~5.8e7 Msol, about 78% of which is contained in the two main ring-like structures at 10 and 15 kpc, at an average temperature of 16.5 K. We find that the gas-to-dust ratio declines exponentially as a function of the galacto-centric distance, in agreement with the known metallicity gradient, with values ranging from 66 in the nucleus to ~275 in the outermost region. [Abridged]

preprint2012arXiv

The Herschel Multi-tiered Extragalactic Survey: HerMES

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 μm), and Herschel-PACS (at 100 and 160 μ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σ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.

preprint2012arXiv

The JCMT Nearby Galaxies Legacy Survey VIII. CO data and the L(CO3-2)-L(FIR) correlation in the SINGS sample

The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.

preprint2012arXiv

The Spitzer Extragalactic Representative Volume Survey (SERVS): survey definition and goals

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

preprint2011arXiv

Discovery of a Multiply-Lensed Submillimeter Galaxy in Early HerMES Herschel/SPIRE Data

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 $μ\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.

preprint2011arXiv

Dynamical Structure of the Molecular Interstellar Medium in an Extremely Bright, Multiply Lensed z~3 Submillimeter Galaxy Discovered with Herschel

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

preprint2011arXiv

HerMES: detection of cosmic magnification of sub-mm galaxies using angular cross-correlation

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

preprint2011arXiv

HerMES: Lyman Break Galaxies individually detected at 0.7 < z < 2.0 in GOODS-N with Herschel/SPIRE

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

preprint2011arXiv

HerMES: point source catalogues from deep Herschel-SPIRE observations

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

preprint2011arXiv

Herschel/HerMES: The X-ray - Infrared correlation for star-forming galaxies at z~1

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

preprint2011arXiv

Modeling of the HERMES J105751.1+573027 submillimeter source lensed by a dark matter dominated foreground group of galaxies

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 $μ=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.

preprint2011arXiv

Observation of H2O in a strongly lensed Herschel-ATLAS source at z=2.3

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

preprint2011arXiv

Photometry and Photometric Redshift catalogs for the Lockman Hole Deep Field

We present broad band photometry and photometric redshifts for 187611 sources located in ~0.5deg^2 in the Lockman Hole area. The catalog includes 389 X-ray detected sources identified with the very deep XMM-Newton observations available for an area of 0.2 deg^2. The source detection was performed on the Rc, z' and B band images and the available photometry is spanning from the far ultraviolet to the mid infrared, reaching in the best case scenario 21 bands. Astrometry corrections and photometric cross-calibrations over the entire dataset allowed the computation of accurate photometric redshifts. Special treatment is undertaken for the X-ray sources, the majority of which is active galactic nuclei. Comparing the photometric redshifts to the available spectroscopic redshifts we achieve for normal galaxies an accuracy of σ_{Δz/(1+z)}=0.036, with 12.7% outliers, while for the X-ray detected sources the accuracy is σ_{Δz/(1+z)}=0.069, with 18.3% outliers, where the outliers are defined as sources with |z_{phot}-z_{spec}|>0.15 (1+z_{spec})}. These results are a significant improvement over the previously available photometric redshifts for normal galaxies in the Lockman Hole, while it is the first time that photometric redshifts are computed and made public for AGN for this field.

preprint2011arXiv

Physical conditions of the interstellar medium of high-redshift, strongly lensed submillimetre galaxies from the Herschel-ATLAS

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.

preprint2011arXiv

Redshift Determination and CO Line Excitation Modeling for the Multiply-Lensed Galaxy HLSW-01

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

preprint2011arXiv

Specific star-formation and the relation to stellar mass from 0<z<2 as seen in the far-infrared at 70 and 160mu

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0<z<2 and find that the sLFIR evolves as (1+z)^alpha with alpha=4.4\pm0.3 for galaxies with 10.5 < log M*/Msun < 12. This is considerably steeper than the {Delucia2007} semi-analytic model (alpha \sim 2.5). When separating galaxies into early and late types on the basis of the optical/IR spectral energy distributions we find that the decrease in sLFIR with stellar mass is stronger in early type galaxies (beta ~ -0.46), while late type galaxies exhibit a flatter trend (beta \sim -0.15). The evolution is strong for both classes but stronger for the early type galaxies. The early types show a trend of decreasing strength of evolution as we move from lower to higher masses while the evolution of the late type galaxies has little dependence on stellar mass. We suggest that in late-type galaxies we are seeing a consistently declining sSFR..

preprint2011arXiv

Spectroscopic follow-up of 70μm sources in Spitzer Wide-area Infrared Extragalactic Legacy Survey

We present spectroscopic follow-up observations of 70μm selected galaxies from the SWIRE XMMLSS and Lockman Hole fields. We have measured spectroscopic redshifts for 293 new sources down to a 70μm flux limit of 9mJy and R < 22. The redshift distribution peaks at z ~ 0.3 and has a high redshift tail out to z = 3.5. We perform emission line diagnostics for 91 sources where [OIII], Hβ, [NII], Hα and [SII] emission lines are available to determine their power source. We find in our sample 13 QSOs, 1 Seyfert II galaxy, 33 star forming galaxies, 30 composite galaxies, 5 LINERs and 21 ambiguous galaxies. We fit single temperature dust spectral energy distributions (SEDs) to 81 70μm sources with 160μm photometry to estimate dust temperatures and masses. Assuming the dust emissivity factor (β) as 1.5, we determine temperatures in the range ~ 20-60K and dust masses with a range of 10{^6}-10{^9} M{_\odot}. Plotting these objects in the luminosity-temperature diagram suggests that these objects have lower dust temperatures than local IR luminous galaxies. The Herschel Space Observatory will be crucial in understanding the nature of these sources and to accurately determining the shape of the Rayleigh-Jeans tail of the dust SED. We then model SEDs from optical to far-IR for each source using a set of galaxy and quasar templates in the optical and near-IR (NIR) and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst and AGN dust torus) in the mid-IR (MIR) to far-IR (FIR). The number of objects fit with each dust template are: 57 Arp 220, 127 M82, 9 cirrus, 1 AGN dust torus, 70 M82 and cirrus, 26 M82 and AGN dust torus and 3 Arp 220 and AGN dust torus. We determine the total IR luminosity (LIR) in range 10{^8}-10{^{15}} L{_\odot} by integrating the SED models from 8 to 1000μm.

preprint2011arXiv

Sub-millimetre galaxies reside in dark matter halos with masses greater than 3x10^11 solar masses

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.

preprint2011arXiv

The Herschel Multi-Tiered Extragalactic Survey: SPIRE-mm Photometric Redshifts

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 μ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 τ = (ν/ν0)^β where β = +1.8 and ν0 = c/100 μ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 μm selected samples.

preprint2011arXiv

The Spitzer Extragalactic Representative Volume Survey (SERVS): The Environments of High-z SDSS Quasi-Stellar-Objects

This paper presents a study of the environments of SDSS Quasi-Stellar-Objects (QSOs) in the Spitzer Extragalactic Representative Volume Survey (SERVS). We concentrate on the high-redshift QSOs as these have not been studied in large numbers with data of this depth before. We use the IRAC 3.6-4.5μm colour of objects and ancillary r-band data to filter out as much foreground contamination as possible. This technique allows us to find a significant (> 4-σ) over-density of galaxies around QSOs in a redshift bin centred on z ~ 2.0 and a (> 2-σ) over-density of galaxies around QSOs in a redshift bin centred on z ~ 3.3. We compare our findings to the predictions of a semi-analytic galaxy formation model, based on the ΛCDM millennium simulation, and find for both redshift bins that the model predictions match well the source-density we have measured from the SERVS data.

preprint2011arXiv

Ultra Steep Spectrum radio sources in the Lockman Hole: SERVS identifications and redshift distribution at the faintest radio fluxes

Ultra Steep Spectrum (USS) radio sources have been successfully used to select powerful radio sources at high redshifts (z>~2). Typically restricted to large-sky surveys and relatively bright radio flux densities, it has gradually become possible to extend the USS search to sub-mJy levels, thanks to the recent appearance of sensitive low-frequency radio facilities. Here a first detailed analysis of the nature of the faintest USS sources is presented. By using Giant Metrewave Radio Telescope and Very Large Array radio observations of the Lockman Hole at 610 MHz and 1.4 GHz, a sample of 58 USS sources, with 610 MHz integrated fluxes above 100 microJy, is assembled. Deep infrared data at 3.6 and 4.5 micron from the Spitzer Extragalactic Representative Volume Survey (SERVS) is used to reliably identify counterparts for 48 (83%) of these sources, showing an average total magnitude of [3.6](AB)=19.8 mag. Spectroscopic redshifts for 14 USS sources, together with photometric redshift estimates, improved by the use of the deep SERVS data, for a further 19 objects, show redshifts ranging from z=0.1 to z=2.8, peaking at z~0.6 and tailing off at high redshifts. The remaining 25 USS sources, with no redshift estimate, include the faintest [3.6] magnitudes, with 10 sources undetected at 3.6 and 4.5 micron (typically [3.6]>22-23 mag, from local measurements), which suggests the likely existence of higher redshifts among the sub-mJy USS population. The comparison with the Square Kilometre Array Design Studies Simulated Skies models indicate that Fanaroff-Riley type I radio sources and radio-quiet Active Galactic Nuclei may constitute the bulk of the faintest USS population, and raises the possibility that the high efficiency of the USS technique for the selection of high redshift sources remains even at the sub-mJy level.

preprint2010arXiv

Evolution of Dust Temperature of Galaxies through Cosmic Time as seen by Herschel

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.

preprint2010arXiv

FIR colours and SEDs of nearby galaxies observed with Herschel

We present infrared colours (in the 25-500 mic spectral range) and UV to radio continuum spectral energy distributions of a sample of 51 nearby galaxies observed with SPIRE on Herschel. The observed sample includes all morphological classes, from quiescent ellipticals to active starbursts. Active galaxies have warmer colour temperatures than normal spirals. In ellipticals hosting a radio galaxy, the far-infrared (FIR) emission is dominated bynthe synchrotron nuclear emission. The colour temperature of the cold dust is higher in quiescent E-S0a than in star-forming systems probably because of the different nature of their dust heating sources (evolved stellar populations, X-ray, fast electrons) and dust grain properties. In contrast to the colour temperature of the warm dust, the f350/f500 index sensitive to the cold dust decreases with star formation and increases with metallicity, suggesting an overabundance of cold dust or an emissivity parameter beta<2 in low metallicity, active systems.

preprint2010arXiv

First results from HerMES on the evolution of the submillimetre luminosity function

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.

preprint2010arXiv

H-ATLAS: PACS imaging for the Science Demonstration Phase

We describe the reduction of data taken with the PACS instrument on board the Herschel Space Observatory in the Science Demonstration Phase of the Herschel-ATLAS (H-ATLAS) survey, specifically data obtained for a 4x4-deg^2 region using Herschel's fast-scan (60 arcsec/s) parallel mode. We describe in detail a pipeline for data reduction using customised procedures within HIPE from data retrieval to the production of science-quality images. We found that the standard procedure for removing Cosmic-Ray glitches also removed parts of bright sources and so implemented an effective two-stage process to minimise these problems. The pronounced 1/f noise is removed from the timelines using 3.4- and 2.5-arcmin boxcar high-pass filters at 100 and 160-um. Empirical measurements of the point-spread function (PSF) are used to determine the encircled energy fraction as a function of aperture size. For the 100- and 160-um bands, the effective PSFs are ~9 and ~13 arcsec (FWHM), and the 90-per-cent encircled energy radii are 13 and 18 arcsec. Astrometric accuracy is good to ~<2 arcsec. The noise in the final maps is correlated between neighbouring pixels and rather higher than advertised prior to launch. For a pair of cross-scans, the 5-sigma point-source sensitivities are 125-165 mJy for 9-13-arcsec radius apertures at 100-um and 150-240 mJy for 13-18-arcsec radius apertures at 160-um.

preprint2010arXiv

HerMES : SPIRE detection of high redshift massive compact galaxies in GOODS-N field

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μm data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24μm and 15/45 (all MIPS-24μm detections) are detected at SPIRE-250μ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).

preprint2010arXiv

HerMES: Deep Galaxy Number Counts from a P(D) Fluctuation Analysis of SPIRE Science Demonstration Phase Observations

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-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.

preprint2010arXiv

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

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

preprint2010arXiv

HerMES: Halo Occupation Number and Bias Properties of Dusty Galaxies from Angular Clustering Measurements

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.

preprint2010arXiv

HerMES: SPIRE galaxy number counts at 250, 350 and 500 microns

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.

preprint2010arXiv

HerMES: The SPIRE confusion limit

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

preprint2010arXiv

HerMES: The Submillimeter Spectral Energy Distributions of Herschel/SPIRE-Detected Galaxies

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.

preprint2010arXiv

Herschel photometric observations of the low metallicity dwarf galaxy NGC 1705

We present Herschel SPIRE and PACS photometeric observations of the low metallicity (Z ~ 0.35 solar) nearby dwarf galaxy, NGC 1705, in six wavelength bands as part of the Dwarf Galaxy Survey guaranteed time Herschel Key Program. We confirm the presence of two dominant circumnuclear IR-bright regions surrounding the central super star cluster that had been previously noted at mid-IR wavelengths and in the sub-mm by LABOCA. On constructing a global spectral energy distribution using the SPIRE and PACS photometry, in conjunction with archival IR measurements, we note the presence of an excess at sub-mm wavelengths. This excess suggests the presence of a significant cold dust component within NGC 1705 and was modeled as an additional cold component in the SED. Although alternative explanations for the sub-mm excess beyond 350 microns, such as changes to the dust emissivity cannot be ruled out, the most likely explanation for the observed submillimetre excess is that of an additional cold dust component.

preprint2010arXiv

Herschel photometric observations of the nearby low metallicity irregular galaxy NGC 6822

We present the first Herschel PACS and SPIRE images of the low-metallicity galaxy NGC6822 observed from 70 to 500 mu and clearly resolve the HII regions with PACS and SPIRE. We find that the ratio 250/500 is dependent on the 24 mu surface brightness in NGC6822, which would locally link the heating processes of the coldest phases of dust in the ISM to the star formation activity. We model the SEDs of some regions HII regions and less active regions across the galaxy and find that the SEDs of HII regions show warmer ranges of dust temperatures. We derive very high dust masses when graphite is used in our model to describe carbon dust. Using amorphous carbon, instead, requires less dust mass to account for submm emission due to its lower emissivity properties. This indicates that SED models including Herschel constraints may require different dust properties than commonly used.

preprint2010arXiv

Herschel unveils a puzzling uniformity of distant dusty galaxies

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.

preprint2010arXiv

Herschel-ATLAS: Statistical Properties of Galactic Cirrus in the GAMA-9 Hour Science Demonstration Phase Field

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $β$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = â H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

preprint2010arXiv

Herschel-SPIRE observations of the disturbed galaxy NGC4438

We present Herschel-SPIRE observations of the perturbed galaxy NGC4438 in the Virgo cluster. These images reveal the presence of extra-planar dust up to ~4-5 kpc away from the galaxy's disk. The dust closely follows the distribution of the stripped atomic and molecular hydrogen, supporting the idea that gas and dust are perturbed in a similar fashion by the cluster environment. Interestingly, the extra-planar dust lacks a warm temperature component when compared to the material still present in the disk, explaining why it was missed by previous far-infrared investigations. Our study provides evidence for dust stripping in clusters of galaxies and illustrates the potential of Herschel data for our understanding of environmental effects on galaxy evolution.

preprint2010arXiv

Herschel-SPIRE, Far-Infrared Properties of Millimetre-Bright and -Faint Radio Galaxies

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

preprint2010arXiv

Mapping the interstellar medium in galaxies with Herschel/SPIRE

The standard method of mapping the interstellar medium in a galaxy, by observing the molecular gas in the CO 1-0 line and the atomic gas in the 21-cm line, is largely limited with current telescopes to galaxies in the nearby universe. In this letter, we use SPIRE observations of the galaxies M99 and M100 to explore the alternative approach of mapping the interstellar medium using the continuum emission from the dust. We have compared the methods by measuring the relationship between the star-formation rate and the surface density of gas in the galaxies. We find the two methods give relationships with a similar dispersion, confirming that observing the continuum emission from the dust is a promising method of mapping the interstellar medium in galaxies.

preprint2010arXiv

Mid- and Far-infrared Luminosity Functions and Galaxy Evolution from Multiwavelength Spitzer Observations up to z~2.5

[Abridged]We exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent 24, 15, 12 and 8micron monochromatic and bolometric IR luminosity functions (LF) over the 0<z<2.5 redshift range. Our analysis is based on the combination of data from deep Spitzer surveys in the VVDS-SWIRE and GOODS areas. To our limiting flux of S(24)=400microJy our derived sample in VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S(24)=80microJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degs. We obtain reliable optical identifications and redshifts, providing us a rich and robust dataset for our luminosity function determination. Based on the multi-wavelength information available, we constrain the LFs at 8, 12, 15 and 24micron. We also extrapolate total IR luminosities from our best-fit to the observed SEDs of each source, and use this to derive the bolometric LF and comoving volume emissivity up to z~2.5. In the 0<z<1 interval, the bolometric IR luminosity density evolves as (1+z)^3.8+/-0.4. Although more uncertain at higher-z, our results show a flattening of the IR luminosity density at z>1. The mean redshift of the peak in the source number density shifts with luminosity: the brighest IR galaxies appear to be forming stars earlier in cosmic time (z>1.5), while the less luminous ones keep doing it at more recent epochs (z~1 for L(IR)<10^11L_sun). Our results suggest a rapid increase of the galaxy IR comoving volume emissivity back to z~1 and a constant average emissivity at z>1. We also seem to find a difference in the evolution rate of the source number densities as a function of luminosity, a downsizing evolutionary pattern similar to that reported from other samples of cosmic sources.

preprint2010arXiv

Mid-infrared spectroscopy of candidate AGN-dominated submillimeter galaxies

Spitzer spectroscopy has revealed that ~80% of submm galaxies (SMGs) are starburst (SB) dominated in the mid-infrared. Here we focus on the remaining ~20% that show signs of harboring powerful active galactic nuclei (AGN). We have obtained Spitzer-IRS spectroscopy of a sample of eight SMGs which are candidates for harboring powerful AGN on the basis of IRAC color-selection (S8/S4.5>2; i.e. likely power-law mid-infrared SEDs). SMGs with an AGN dominating (>50%) their mid-infrared emission could represent `missing link' sources in an evolutionary sequence involving a major merger. First of all, we detect PAH features in all of the SMGs, indicating redshifts from 2.5-3.4, demonstrating the power of the mid-infrared to determine redshifts for these optically faint dusty galaxies. Secondly, we see signs of both star-formation (from the PAH features) and AGN activity (from continuum emission) in our sample: 62% of the sample are AGN-dominated in the mid-infrared with a median AGN content of 56%, compared with <30% on average for typical SMGs, revealing that our IRAC color selection has successfully singled out sources with proportionately more AGN emission than typical SB-dominated SMGs. However, we find that only about 10% of these AGN dominate the bolometric emission of the SMG when the results are extrapolated to longer infrared wavelengths, implying that AGN are not a significant power source to the SMG population overall, even when there is evidence in the mid-infrared for substantial AGN activity. When existing samples of mid-infrared AGN-dominated SMGs are considered, we find that S8/S4.5>1.65 works well at selecting mid-infrared energetically dominant AGN in SMGs, implying a duty cycle of ~15% if all SMGs go through a subsequent mid-infrared AGN-dominated phase in the proposed evolutionary sequence.

preprint2010arXiv

On the origin of M81 group extended dust emission

Galactic cirrus emission at far-infrared wavelengths affects many extragalactic observations. Separating this emission from that associated with extragalactic objects is both important and difficult. In this paper we discuss a particular case, the M81 group, and the identification of diffuse structures prominent in the infrared, but also detected at optical wavelengths. The origin of these structures has previously been controversial, ranging from them being the result of a past interaction between M81 and M82 or due to more local Galactic emission. We show that over of order a few arcminute scales the far-infrared (Herschel 250 &μ&m) emission correlates spatially very well with a particular narrow velocity (2-3 km/s) component of the Galactic HI. We find no evidence that any of the far-infrared emission associated with these features actually originates in the M81 group. Thus we infer that the associated diffuse optical emission must be due to galactic light back scattered off dust in our galaxy. Ultra-violet observations pick out young stellar associations around M81, but no detectable far-infrared emission. We consider in detail one of the Galactic cirrus features, finding that the far-infrared HI relation breaks down below arc minute scales and that at smaller scales there can be quite large dust temperature variations.

preprint2010arXiv

Probing the molecular interstellar medium of M82 with Herschel-SPIRE spectroscopy

We present the observations of the starburst galaxy M82 taken with the Herschel SPIRE Fourier Transform Spectrometer. The spectrum (194-671 μm) shows a prominent CO rotational ladder from J = 4-3 to 13-12 emitted by the central region of M82. The fundamental properties of the gas are well constrained by the high J lines observed for the first time. Radiative transfer modeling of these high-S/N 12CO and 13CO lines strongly indicates a very warm molecular gas component at ~500 K and pressure of ~3x10^6 K cm^-3, in good agreement with the H_2 rotational lines measurements from Spitzer and ISO. We suggest that this warm gas is heated by dissipation of turbulence in the interstellar medium (ISM) rather than X-rays or UV flux from the straburst. This paper illustrates the promise of the SPIRE FTS for the study of the ISM of nearby galaxies.

preprint2010arXiv

Radial distribution of gas and dust in the two spiral galaxies M99 and M100

By combining Herschel-SPIRE data with archival Spitzer, HI, and CO maps, we investigate the spatial distribution of gas and dust in the two famous grand-design spirals M99 and M100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SPIRE photometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices f350/f500 and f250/f350 decrease radially consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of both galaxies.

preprint2010arXiv

SAGACE: the Spectroscopic Active Galaxies And Clusters Explorer

The SAGACE experiment consists of a mm/sub-mm telescope with a 3-m diameter primary mirror, coupled to a cryogenic multi-beam differential spectrometer. SAGACE explores the sky in the 100-760 GHz frequency range, using four diffraction-limited bolometer arrays. The instrument is designed to perform spectroscopic surveys of the Sunyaev-Zeldovich effects of thousands of galaxy clusters, of the spectral energy distribution of active galactic nuclei, and of the [CII] line of a thousand galaxies in the redshift desert. In 2008 a full phase-A study for a national small mission was completed and delivered to the Italian Space Agency (ASI). We have shown that taking advantage of the differential operation of the Fourier Transform Spectrometer, this ambitious instrument can operate from a Molniya orbit, and can be built and operated within the tight budget of a small mission.

preprint2010arXiv

SPIRE imaging of M82: cool dust in the wind and tidal streams

M82 is a unique representative of a whole class of galaxies, starbursts with superwinds, in the Very Nearby Galaxy Survey with Herschel. In addition, its interaction with the M81 group has stripped a significant portion of its interstellar medium from its disk. SPIRE maps now afford better characterization of the far-infrared emission from cool dust outside the disk, and sketch a far more complete picture of its mass distribution and energetics than previously possible. They show emission coincident in projection with the starburst wind and in a large halo, much more extended than the PAH band emission seen with Spitzer. Some complex substructures coincide with the brightest PAH filaments, and others with tidal streams seen in atomic hydrogen. We subtract the far-infrared emission of the starburst and underlying disk from the maps, and derive spatially-resolved far-infrared colors for the wind and halo. We interpret the results in terms of dust mass, dust temperature, and global physical conditions. In particular, we examine variations in the dust physical properties as a function of distance from the center and the wind polar axis, and conclude that more than two thirds of the extraplanar dust has been removed by tidal interaction, and not entrained by the starburst wind.

preprint2010arXiv

The AGN fraction of submm-selected galaxies and contributions to the submm/mm-wave extragalactic background light

We present a comparison of the SCUBA Half Degree Extragalactic Survey (SHADES) at 450, 850 and 1100 microns with deep guaranteed time 15 microns AKARI FU-HYU survey data and Spitzer guaranteed time data at 3.6-24 microns in the Lockman Hole East. The AKARI data was analysed using bespoke software based in part on the drizzling and minimum-variance matched filtering developed for SHADES, and was cross calibrated against Infrared Space Observatory (ISO) fluxes. Our stacking analyses find AKARI 15um galaxies with >~200 microJy contribute >10% of the 450 micron background, but only <4% of the 1100 micron background, suggesting that different populations contribute at mm-wavelengths. We confirm our earlier result that the ultra-deep 450 micron SCUBA-2 Cosmology Survey will be dominated by populations already detected by AKARI and Spitzer mid-infrared surveys. The superb mid-infrared wavelength coverage afforded by combining Spitzer and AKARI photometry is an excellent diagnostic of AGN contributions, and we find that (23-52)% of submm-selected galaxies have AGN bolometric fractions f_AGN>0.3.

preprint2010arXiv

The central region of spiral galaxies as seen by Herschel. M81, M99 and M100

With appropriate spatial resolution, images of spiral galaxies in thermal infrared (~10 micron and beyond) often reveal a bright central component, distinct from the stellar bulge, superimposed on a disk with prominent spiral arms. ISO and Spitzer studies have shown that much of the scatter in the mid-infrared colors of spiral galaxies is related to changes in the relative importance of these two components, rather than to other modifications, such as the morphological type or star formation rate, that affect the properties of the galaxy as a whole. With the Herschel imaging capability from 70 to 500 micron, we revisit this two-component approach at longer wavelengths, to see if it still provides a working description of the brightness distribution of galaxies, and to determine its implications on the interpretation of global far-infrared properties of galaxies.

preprint2010arXiv

The Deep SPIRE HerMES Survey: Spectral Energy Distributions and their Astrophysical Indications at High Redshift

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.

preprint2010arXiv

The dust morphology of the elliptical Galaxy M86 with SPIRE

We present Herschel-SPIRE observations at 250-500um of the giant elliptical galaxy M86 and examine the distribution of the resolved cold dust emission and its relation with other galactic tracers. The SPIRE images reveal three dust components: emission from the central region; a dust lane extending north-south; and a bright emission feature 10kpc to the south-east. We estimate that approximately 10^6 solar masses of dust is spatially coincident with atomic and ionized hydrogen, originating from stripped material from the nearby spiral NGC4438 due to recent tidal interactions with M86. The gas-to-dust ratio of the cold gas component ranges from ~20-80. We discuss the different heating mechanisms for the dust features.

preprint2010arXiv

The far-infrared/radio correlation as probed by Herschel

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.

preprint2010arXiv

The HerMES SPIRE submillimeter local luminosity function

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

preprint2010arXiv

The Herschel Multi-Tiered Extragalactic Survey: Source Extraction and Cross-Identifications in Confusion-Dominated SPIRE Images

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to 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.

preprint2010arXiv

The Herschel Reference Survey

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 micron. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K-band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, (ii) the connection between the dust content and composition and the other phases of the interstellar medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multi-frequency data to carry out an analysis of the statistical properties of the sample.

preprint2010arXiv

The Herschel Space Observatory view of dust in M81

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh-Jeans slope of dust emission observed at 70-500 microns in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 microns are primarily dependent on radius but that the ratio of 70 to 160 micron emission shows no clear dependence on surface brightness or radius. These results along with analyses of the spectral energy distributions imply that the 160-500 micron emission traces 15-30 K dust heated by evolved stars in the bulge and disc whereas the 70 micron emission includes dust heated by the active galactic nucleus and young stars in star forming regions.

preprint2010arXiv

The Herschel-SPIRE instrument and its in-flight performance

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.

preprint2010arXiv

The JCMT Nearby Galaxies Legacy Survey IV. Velocity Dispersions in the Molecular Interstellar Medium in Spiral Galaxies

An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 12 nearby spiral galaxies reveals low velocity dispersions in the molecular component of the interstellar medium. The three lowest luminosity galaxies show a relatively flat velocity dispersion as a function of radius while the remaining nine galaxies show a central peak with a radial fall-off within 0.2-0.4 r(25). Correcting for the average contribution due to the internal velocitydispersions of a population of giant molecular clouds, the average cloud-cloud velocity dispersion across the galactic disks is 6.1 +/- 1.0 km/s (standard deviation 2.9 km/s), in reasonable agreement with previous measurements for the Galaxy andM33. The cloud-cloud velocity dispersion derived from the CO data is on average two times smaller than the HI velocity dispersion measured in the same galaxies. The low cloud-cloudvelocity dispersion implies that the molecular gas is the critical component determining the stability of the galactic disk against gravitational collapse, especially in those regions of the disk which are H2 dominated. The cloud-cloud velocity dispersion shows a significant positivecorrelation with both the far-infrared luminosity, which traces the star formation activity, and the K-band absolute magnitude, which traces the total stellar mass. For three galaxies in the Virgo cluster, smoothing the data to a resolution of 4.5 kpc (to match the typical resolution of high redshift CO observations) increases the measured velocity dispersion by roughly a factor of two, comparable to the dispersion measured recently in a normal galaxy at z=1. This comparison suggests that the mass and star formation rate surface densities may be similar in galaxies from z=0-1 and that the high star formation rates seen at z=1 may be partly due to the presence of physically larger molecular gas disks.

preprint2009arXiv

Spectroscopic Cosmological Surveys in the Far-IR

We show the feasibility of spectroscopic cosmological surveys with the SAFARI instrument onboard of SPICA. The work is done through simulations that make use of both empirical methods, i.e. the use of observed luminosity functions and theoretical models for galaxy formation and evolution. The relations assumed between the line emission to trace AGN and star formation activity have been derived from the observations of local samples of galaxies. The results converge to indicate the use of blind spectroscopy with the SAFARI FTS at various resolutions to study galaxy evolution from the local to the distant (z~3) Universe. Specifically, two different and independent galaxy evolution models predict about 7-10 sources to be spectroscopically detected in more than one line in a 2'x 2'SAFARI field of view, down to the expected flux limits of SAFARI, with about 20% of sources to be detected at z>2. SPICA-SAFARI will be therefore excellent at detecting high-z sources and at assessing in a direct way their nature (e.g whether mainly AGN or Star Formation powered) thanks to blind spectroscopy.

preprint2009arXiv

Testing the starburst/AGN connection with SWIRE X-ray/70 micron sources

We explore the nature of X-ray sources with 70 micron counterparts selected in the SWIRE fields ELAIS-N1, Lockman Hole and Chandra Deep Field South, for which Chandra X-ray data are available. A total of 28 X-ray/70 micron sources in the redshift interval 0.5<z<1.3 are selected. The X-ray luminosities and the shape of the X-ray spectra show that these sources are AGN. Modelling of the optical to far-infrared Spectral Energy Distribution indicates that most of them (27/28) have a strong starburst component (>50 solar masses per year) that dominates in the infrared. It is found that the X-ray and infrared luminosities of the sample sources are broadly correlated, consistent with a link between AGN activity and star-formation. Contrary to the predictions of some models for the co-evolution of AGN and galaxies, the X-ray/70 micron sources in the sample are not more obscured at X-ray wavelengths compared to the overall X-ray population. It is also found that the X-ray/70 micron sources have lower specific star-formation rates compared to the general 70 micron population, consistent with AGN feedback moderating the star-formation in the host galaxies.

preprint2009arXiv

The Herschel ATLAS

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.

preprint2009arXiv

The JCMT Nearby Galaxies Legacy Survey III: Comparisons of cold dust, polycyclic aromatic hydrocarbons, molecular gas, and atomic gas in NGC 2403

We used 3.6, 8.0, 70, 160 micron Spitzer Space Telescope data, James Clerk Maxwell Telescope HARP-B CO J=(3-2) data, National Radio Astronomy Observatory 12 meter telescope CO J=(1-0) data, and Very Large Array HI data to investigate the relations among PAHs, cold (~20 K) dust, molecular gas, and atomic gas within NGC 2403, an SABcd galaxy at a distance of 3.13 Mpc. The dust surface density is mainly a function of the total (atomic and molecular) gas surface density and galactocentric radius. The gas-to-dust ratio monotonically increases with radius, varying from ~100 in the nucleus to ~400 at 5.5 kpc. The slope of the gas-to-dust ratio is close to that of the oxygen abundance, suggesting that metallicity strongly affects the gas-to-dust ratio within this galaxy. The exponential scale length of the radial profile for the CO J=(3-2) emission is statistically identical to the scale length for the stellar continuum-subtracted 8 micron (PAH 8 micron) emission. However, CO J=(3-2) and PAH 8 micron surface brightnesses appear uncorrelated when examining sub-kpc sized regions.

preprint2006arXiv

The SCUBA HAlf Degree Extragalactic Survey (SHADES) -- II. Submillimetre maps, catalogue and number counts

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

preprint2004arXiv

The European Large Area ISO Survey (ELAIS): The Final Band-merged Catalogue

We present the final band-merged ELAIS catalogue at 6.7, 15, 90, and 175 $μ$m, and the associated data at u,g,r,i,z,J,H,K, and 20cm. The origin of the survey, infrared and radio observations, data-reduction and optical identifications are briefly reviewed, and a summary of the area covered, and completeness limit for each infrared band is given. A detailed discussion of the band-merging and optical association strategy is given. The total catalogues consists of 2860 sources. For extragalactic sources observed in 3 or more infrared bands, colour-colour diagrams are presented and discussed in terms of the contributing infrared populations. Spectral energy distributions are shown for selected sources and compared with cirrus, M82 and Arp220 starburst, and AGN dust torus models. Spectroscopic redshifts are tabulated, where available. For the N1 and N2 areas, the INT ugriz Wide Field Survey, permits photometric redshifts to be estimated for galaxies and quasars. These agree well with the spectroscopic redshifts, within the uncertainty of the photometric method. The redshift distribution is given for selected ELAIS bands and colour-redshift diagrams are discussed. There is a high proportion of ultraluminous infrared galaxies in the ELAIS Catalogue ($> 10 %$ of 15 $μ$m sources), many with Arp220-like colours. 10 hyperluminous infrared galaxies and 10 EROs are found in the survey. The large numbers of ultraluminous galaxies imply very strong evolution in the star-formation rate between z = 0 and 1.