Source author record

E. Hatziminaoglou

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

54works
4topics
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

54 published item(s)

preprint2026arXiv

Calibrating Mid-Infrared Emission Features As Diagnostics of Star Formation in Infrared-Luminous Galaxies via Radiative Transfer Modeling

Luminous infrared galaxies are key sites of obscured stellar mass assembly at z > 0.5. Their star formation rates (SFRs) are often estimated using the luminosities of the 6.2 micron and 11.2 micron polycyclic aromatic hydrocarbon (PAH) features, or those of the [Ne II] and [Ne III] fine-structure lines, as they are minimally affected by obscuration. It is uncertain whether the calibration of these features as SFR tracers depends on the starburst bolometric luminosity or the level of Active Galactic Nucleus (AGN) activity. We here investigate the relationship between the luminosities of PAH and Neon lines with star formation rate for highly luminous objects using radiative transfer modeling and archival observations of 42 local Ultraluminous (>= 10^12 L_sun) Infrared Galaxies (ULIRGs). We find that PAH and [Ne II] features arise mainly in star-forming regions, with small contributions from the AGN or host, but that the [Ne III] line has a mixed contribution from both star formation and AGN activity. We present relations between L_PAH and L_NeII, and both starburst luminosity and SFR. We find relations for lower luminosity (L_IR ~= 10^10-10^12 L_sun) systems underestimate the SFRs in local ULIRGs by up to ~1 dex. The 6.2 micron and 11.2 micron PAH features, and the [Ne II] line, are thus good tracers of SFR in ULIRGs. We do not find that a more luminous AGN affects the relationship between SFR and PAH or Neon luminosity, but that it can make PAH emission harder to discern. Our results and derived relations are relevant to studies of star-forming and composite galaxies at z < 3 with the James Webb Space Telescope.

preprint2022arXiv

The molecular gas properties in local Seyfert 2 galaxies

We present a multi-wavelength study of the molecular gas properties of a sample of local Seyfert 2 galaxies to assess if, and to what extent, the presence of an active galactic nucleus (AGN) can affect the Interstellar Medium (ISM) properties in a sample of 33 local Seyfert 2 galaxies. We compare the molecular gas content (MH2), derived from new and archival low-J CO line measurements of a sample of AGN and a control sample of star-forming galaxies (SFGs). Both the AGN and the control sample are characterised in terms of host-galaxy properties (e.g., stellar and dust masses, Mstar and Mdust, respectively; and star formation rate, SFR). We also investigate the effect of AGN activity onto the emission of polycyclic aromatic hydrocarbon (PAH) molecules in the mid-infrared (MIR), a waveband where the dust-reprocessed emission from the obscured AGN contributes the most. The AGN hosted in less massive galaxies (i.e., log(Mstar/Msun)<10.5; log(Mdust/Msun)<7.5) show larger molecular gas content with respect to SFGs matched in stellar and dust mass. When comparing their depletion time(tdep~MH2/SFR), AGN show tdep~0.3-1.0 Gyr, similar to those observed in the control sample of SFGs. Seyfert 2 galaxies show fainter PAH luminosity at increasingly larger dominance of the nuclear activity in the MIR. We find no clear evidence for a systematic reduction of the molecular gas reservoir at galactic scale in Seyfert galaxies with respect to SFGs. This is in agreement with recent studies showing that molecular gas content only is reduced in regions of sub-kpc size, where the emission from the accreting supermassive black hole dominates. Nonetheless, we show that the impact of AGN activity on the ISM is clearly visible as suppression of the PAH luminosity.

preprint2019arXiv

A CO molecular gas wind 340 pc away from the Seyfert 2 nucleus in ESO 420-G13 probes an elusive radio jet

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 μm}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

preprint2016arXiv

Parsec-scale HI absorption structure in a low-redshift galaxy seen against a Compact Symmetric Object

We present global VLBI observations of the 21-cm transition of atomic hydrogen seen in absorption against the radio source J0855+5751. The foreground absorber (SDSS~J085519.05+575140.7) is a dwarf galaxy at $z$ = 0.026. As the background source is heavily resolved by VLBI, the data allow us to map the properties of the foreground HI gas with a spatial resolution of 2pc. The absorbing gas corresponds to a single coherent structure with an extent $>$35pc, but we also detect significant and coherent variations, including a change in the HI optical depth by a factor of five across a distance of $\leq$6pc. The large size of the structure provides support for the Heiles & Troland model of the ISM, as well as its applicability to external galaxies. The large variations in HI optical depth also suggest that caution should be applied when interpreting $T_S$ measurements from radio-detected DLAs. In addition, the distorted appearance of the background radio source is indicative of a strong jet-cloud interaction in its host galaxy. We have measured its redshift ($z$ = 0.54186) using optical spectroscopy on the William Herschel Telescope and this confirms that J0855+5751 is a FRII radio source with a physical extent of $<$1kpc and supports the previous identification of this source as a Compact Symmetric Object. These sources often show absorption associated with the host galaxy and we suggest that both HI and OH should be searched for in J0855+5751.

preprint2015arXiv

ALMA Long Baseline Observations of the Strongly Lensed Submillimeter Galaxy HATLAS J090311.6+003906 at z=3.042

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

preprint2015arXiv

An Overview of the 2014 ALMA Long Baseline Campaign

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

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

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

preprint2013arXiv

A Dust-Obscured Massive Maximum-Starburst Galaxy at a Redshift of 6.34

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

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

Mid- to far infrared properties of star-forming galaxies and active galactic nuclei

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

preprint2013arXiv

Panchromatic spectral energy distributions of Herschel sources

(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is then fitted with a modified version of the MAGPHYS code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an AGN. The defined Gaussian grouping is also used to identify rare sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z~1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 micron detected sources in PEP GOODS fields. AGN appear to be distributed in the stellar mass (M*) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M*-SFR-z space.

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 SCUBA-2 Cosmology Legacy Survey: demographics of the 450 μm-population

We investigate the multi-wavelength properties of a sample of 450-μm selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4σ in deep SCUBA-2 450-μm observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of σ450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450μm-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <β_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between β_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 μm sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

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

A Population of z>2 Far-Infrared Herschel-SPIRE selected Starbursts

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

preprint2012arXiv

A Redshift Survey of Herschel Far-Infrared Selected Starbursts and Implications for Obscured Star Formation

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

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

Probing the Nuclear and Circumnuclear Activity of NGC1365 in the Infrared

We present new far-infrared (70-500micron) Herschel PACS and SPIRE imaging observations as well as new mid-IR Gemini/T-ReCS imaging (8.7 and 18.3micron) and spectroscopy of the inner Lindblad resonance (ILR) region (R<2.5kpc) of the spiral galaxy NGC1365. We complemented these observations with archival Spitzer imaging and spectral mapping observations. The ILR region of NGC1365 contains a Seyfert 1.5 nucleus and a ring of star formation with an approximate diameter of 2kpc. The strong star formation activity in the ring is resolved by the Herschel/PACS imaging data, as well as by the Spitzer 24micron continuum emission, [NeII]12.81micron line emission, and 6.2 and 11.3micron PAH emission. The AGN is the brightest source in the central regions up to lambda~24micron, but it becomes increasingly fainter in the far-infrared when compared to the emission originating in the infrared clusters (or groups of them) located in the ring. We modeled the AGN unresolved infrared emission with the CLUMPY torus models and estimated that the AGN contributes only to a small fraction (~5%) of the infrared emission produced in the inner ~5kpc. We fitted the non-AGN 24-500micron spectral energy distribution of the ILR region and found that the dust temperatures and mass are similar to those of other nuclear and circumnuclear starburst regions. Finally we showed that within the ILR region of NGC1365 most of the on-going star formation activity is taking place in dusty regions as probed by the 24micron emission.

preprint2012arXiv

SAFIR: testing the coexistence of AGN and star formation activity and the nature of the dusty torus in the local universe

We present the Seyfert and star formation Activity in the Far-InfraRed (SAFIR) project, a small (15.1h) Herschel guaranteed time proposal performing PACS and SPIRE imaging of a small sample of nearby Seyfert galaxies. This project is aimed at studying the physical nature of the nuclear IR emission by means of multi-component spectral energy distribution (SED) fitting and the star formation properties of AGN hosts, as traced by cold dust. We summarize the results achieved so far and outline the on-going work.

preprint2012arXiv

Smooth and Clumpy Dust Distribution in AGN: a Direct Comparison of two Commonly Explored Infrared Emission Models

The geometry of the dust distribution within the inner regions of Active Galactic Nuclei (AGN) is still a debated issue and relates directly with the AGN unified scheme. Traditionally, models discussed in the literature assume one of two distinct dust distributions in what is believed to be a toroidal region around the Supermassive Black Holes: a continuous distribution, customarily referred to as smooth, and a concentration of dust in clumps or clouds, referred to as clumpy. In this paper we perform a thorough comparison between two of the most popular models in the literature, namely the smooth models by Fritz. et al. 2006 and the clumpy models by Nenkova et al. 2008a, in their common parameters space. Particular attention is paid to the silicate features at ~9.7 and ~18 micron, the width of the infrared bump, the near-infrared index and the luminosity at 12.3 micron, all previously reported as possible diagnostic tools to distinguish between the two dust distributions. We find that, due to the different dust chemical compositions used in the two models, the behaviour of the silicate features at 9.7 and 18 micron is quite distinct between the two models. The width of the infrared bump and the peak of the infrared emission can take comparable values, their distributions do, however, vary. The near-infrared index is also quite different, due partly to the primary sources adopted by the two models. Models with matched parameters do not produce similar SEDs and virtually no random parameter combinations can result in seemingly identical SEDs.

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

preprint2012arXiv

The starburst-AGN connection in the merger galaxy Mrk 938: an infrared and X-ray view

Mrk938 is a luminous infrared galaxy in the local Universe believed to be the remnant of a galaxy merger. It shows a Seyfert 2 nucleus and intense star formation according to optical spectroscopic observations. We have studied this galaxy using new Herschel far-IR imaging data in addition to archival X-ray, UV, optical, near-IR and mid-IR data. Mid- and far-IR data are crucial to characterise the starburst contribution, allowing us to shed new light on its nature and to study the coexistence of AGN and starburst activity in the local Universe. The decomposition of the mid-IR Spitzer spectrum shows that the AGN bolometric contribution to the mid-IR and total infrared luminosity is small (Lbol(AGN)/LIR~0.02), which agrees with previous estimations. We have characterised the physical nature of its strong infrared emission and constrained it to a relatively compact emitting region of <2kpc. It is in this obscured region where most of the current star formation activity is taking place as expected for LIRGs. We have used Herschel imaging data for the first time to constrain the cold dust emission with unprecedented accuracy. We have fitted the integrated far-IR spectral energy distribution and derived the properties of the dust, obtaining a dust mass of 3x10^7Msun. The far-IR is dominated by emission at 35K, consistent with dust heated by the on-going star formation activity.

preprint2011arXiv

Herschel-ATLAS: The link between accretion luminosity and star formation in quasar host galaxies

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^θ (1 + z)^ζ, we find θ = 0.22 +/- 0.08 and ζ = 1.6 +/- 0.4, although there is substantial additional uncertainty in ζ of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

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

Resolving the nuclear dust distribution of the Seyfert 2 galaxy NGC 3081

We report far-infrared (FIR) imaging of the Seyfert 2 galaxy NGC 3081 in the range 70-500 micron, obtained with an unprecedented angular resolution, using the Herschel Space Observatory instruments PACS and SPIRE. The 11 kpc (~70 arcsec) diameter star-forming ring of the galaxy appears resolved up to 250 micron. We extracted infrared (1.6-500 micron) nuclear fluxes, that is active nucleus-dominated fluxes, and fitted them with clumpy torus models, which successfully reproduce the FIR emission with small torus sizes. Adding the FIR data to the near- and mid-infrared spectral energy distribution (SED) results in a torus radial extent of Ro=4(+2/-1) pc, as well as in a flat radial distribution of the clouds (i.e. the q parameter). At wavelengths beyond 200 micron, cold dust emission at T=28+/-1 K from the circumnuclear star-forming ring of 2.3 kpc (~15 arcsec) in diameter starts making a contribution to the nuclear emission. The dust in the outer parts of the galaxy is heated by the interstellar radiation field (19+/-3 K).

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.

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

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

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: Herschel-SPIRE observations of Lyman Break Galaxies

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

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 ATLAS: The cosmic star formation history of quasar host galaxies

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

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: blazars in the SDP field

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$μ$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $γ$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

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 environments of z~1 Active Galactic Nuclei at 3.6um

We present an analysis of a large sample of AGN environments at z~1 using stacked Spitzer data at 3.6um. The sample contains type-1 and type-2 AGN in the form of quasars and radio galaxies, and spans a large range in both optical and radio luminosity. We find, on average, that 2 to 3 massive galaxies containing a substantial evolved stellar population lie within a 200-300 kpc radius of the AGN, constituting a >8-sigma excess relative to the field. Secondly, we find evidence for the environmental source density to increase with the radio luminosity of AGN, but not with black-hole mass. This is shown first by dividing the AGN into their classical AGN types, where we see more significant over-densities in the fields of the radio-loud AGN. If instead we dispense with the classical AGN definitions, we find that the source over-density as a function of radio luminosity for all our AGN exhibits a positive correlation. One interpretation of this result is that the Mpc-scale environment is in some way influencing the radio emission that we observe from AGN. This could be explained by the confinement of radio jets in dense environments leading to enhanced radio emission or, alternatively, may be linked to more rapid black-hole spin brought on by galaxy mergers.

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

preprint2009arXiv

A Deep XMM-Newton Serendipitous Survey of a middle-latitude area. II. New deeper X-ray and optical observations

The radio-quiet neutron star 1E1207.4-5209 has been the target of several XMM-Newton observations, with a total exposure of ~350 ks. The source is located at intermediate galactic latitude (b~10 degrees), i.e. in a sky region with an extremely interesting mix of both galactic and extra-galactic X-ray sources. The aim of our work is to investigate the properties of both the intermediate-latitude galactic and extra-galactic X-ray source populations in the 1E1207.4-5209 field. We performed a coherent analysis of the whole XMM-Newton observation data set to build a catalogue of serendipitous X-ray sources detected with high confidence and to derive information on the source flux, spectra, and time variability. In addition, we performed a complete multi-band (UBVRI) optical coverage of the field with the Wide Field Imager (WFI) of the ESO/MPG 2.2m telescope (La Silla) to search for candidate optical counterparts to the X-ray sources, down to a V-band limiting magnitude of ~24.5. We detected a total of 144 serendipitous X-ray sources. Thanks to the refined X-ray positions and to the WFI observations, we found candidate optical counterparts for most of the X-ray sources in our compilation. For most of the brightest ones we proposed a likely classification based on both the X-ray spectra and the optical colours. Our results indicate that at intermediate galactic latitude the X-ray source population is dominated by the extra-galactic component, but with a significant contribution from the galactic component in the soft energy band, below 2 keV.

preprint2009arXiv

Resolved Dust Emission in a Quasar at z=3.65

We present submillimetre observations of the z=3.653 quasar SDSS160705+533558 together with data in the optical and infrared. The object is unusually bright in the far-IR and submm with an IR luminosity of ~10^14 L_sun. We ascribe this luminosity to a combination of AGN and starburst emission, with the starburst forming stars at a rate of a few thousand solar masses per year. Submillimetre Array (SMA) imaging observations with a resolution ~1" show that the submm (850 micron) emission is extended on scales of 10--35kpc and is offset from the optical position by ~10 kpc. This morphology is dissimilar to that found in submm galaxies, which are generally un- or marginally resolved on arcsecond scales, or submm-luminous AGN where the AGN lies at the peak of the submm or molecular emission. The simplest explanation is that the object is in the early stages of a merger between a gas rich galaxy, which hosts the starburst, and a gas-poor AGN-host galaxy, which is responsible for the quasar emission. It is also possible that jet induced star formation might contribute to the unusual morphology.

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.

preprint2001arXiv

ESO Imaging Survey. Deep Public Survey: Multi-Color Optical Data for the Chandra Deep Field South

This paper presents multi-passband optical data obtained from observations of the Chandra Deep Field South (CDF-S), located at alpha ~ 3h 32m, delta ~ -27d 48m. The observations were conducted at the ESO/MPG 2.2m telescope at La Silla using the 8kx8k Wide-Field Imager (WFI). This data set, taken over a period of one year, represents the first field to be completed by the ongoing Deep Public Survey (DPS) being carried out by the ESO Imaging Survey (EIS) project. This paper describes the optical observations, the techniques employed for un-supervised pipeline processing and the general characteristics of the final data set. The paper includes data taken in six different filters U'UBVRI. The data cover an area of about 0.25 square degrees reaching 5 sigma limiting magnitudes of U'_AB=26.0, U_AB=25.7, B_AB=26.4$, V_AB=25.4, R_AB=25.5 and I_AB= 24.7 mag, as measured within a 2xFWHM aperture. The optical data covers the area of ~ 0.1