A search for debris disks in the Herschel-ATLAS

Original article can be found at: http://www.aanda.org/ Copyright The European Southern Observatory (ESO)Aims. We aim to demonstrate that the Herschel-ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods. We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results. We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.Peer reviewe

ALMA observations of lensed Herschel sources: testing the dark matter halo paradigm

With the advent of wide-area submillimetre surveys, a large number of high-redshift gravitationally lensed dusty star-forming galaxies have been revealed. Because of the simplicity of the selection criteria for candidate lensed sources in such surveys, identified as those with S500 μm > 100 mJy, uncertainties associated with the modelling of the selection function are expunged. The combination of these attributes makes submillimetre surveys ideal for the study of strong lens statistics. We carried out a pilot study of the lensing statistics of submillimetre-selected sources by making observations with the Atacama Large Millimeter Array (ALMA) of a sample of strongly lensed sources selected from surveys carried out with the Herschel Space Observatory. We attempted to reproduce the distribution of image separations for the lensed sources using a halo mass function taken from a numerical simulation that contains both dark matter and baryons. We used three different density distributions, one based on analytical fits to the haloes formed in the EAGLE simulation and two density distributions [Singular Isothermal Sphere (SIS) and SISSA] that have been used before in lensing studies. We found that we could reproduce the observed distribution with all three density distributions, as long as we imposed an upper mass transition of ∼1013 M⊙ for the SIS and SISSA models, above which we assumed that the density distribution could be represented by a Navarro–Frenk–White profile. We show that we would need a sample of ∼500 lensed sources to distinguish between the density distributions, which is practical given the predicted number of lensed sources in the Herschel surveys

Extremely red galaxies: dust attenuation and classification

We re-address the classification criterion for extremely red galaxies (ERGs) of Pozzetti and Mannucci (2000 -- PM00), which aims to separate, in the Ic-K (or Rc-K) vs. J-K colour--colour diagram, passively evolving, old (> 1 Gyr) stellar populations in a dust-free environment, associated with ellipticals (Es), from dusty starburst galaxies (DSGs), both at 1 < z < 2. We explore a category of objects not considered previously, i.e., galaxies forming in this redshift range on short (0.1 Gyr) timescales and observed also in their early, dusty post-starburst phase. We also investigate the impact of structure of the dusty medium and dust amount on the observed optical/near-IR colours of high-z DSGs/DPSGs, through multiple-scattering radiative transfer calculations for a dust/stars configuration and an extinction function calibrated with nearby dusty starbursts. As a main result, we find that dusty post-starburst galaxies (DPSGs), with ages between 0.2 and 1 Gyr, at 1.3 < z < 2 mix with Es at 1 < z < 2 for a large range in dust amount. This ``intrusion'' is a source of concern for the present two-colour classification of ERGs. On the other hand, we confirm, in agreement with PM00, that DSGs are well separated from Es, both at 1 < z < 2, in the Ic-K vs. J-K colour--colour diagram, whatever the structure (two-phase clumpy or homogeneous) of their dusty medium and their dust amount are. This result holds under the new hypothesis of high-z Es being as dusty as nearby ones. Thus the interpretation of the optical/near-IR colours of high-z Es may suffer from a multiple degeneracy among age, metallicity, dust and redshift. We also find that DPSGs at z around 1 mix with DSGs at 1 < z < 2, as a function of dust amount and structure of the dusty medium. All these results help explaining the complexity of the ERG classification... (Abridged)Comment: 17 pages, 19 figures, accepted for publication in MNRA

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.Comment: 13 pages, 12 figures, accepted by MNRA

Herschel-ATLAS: Far-infrared properties of radio-loud and radio-quiet quasars

This is pre-copyedited, author-produced pdf of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record, E. Kalfountzou, et al., ‘Herschel-ATLAS: Far-infrared properties of radio-loud and radio-quiet quasars’, MNRAS, Vol 42(2): 1181-1196, first published online June 11, 2014, is available online via doi: 10.1093/mnras/stu782 Published by Oxford University Press on behalf of the Royal Astronomical Society.We have constructed a sample of radio-loud and radio-quiet quasars from the Faint Images Radio Sky at Twenty-one centimetres and the Sloan Digital Sky Survey Data Release 7, over the Herschel-ATLAS Phase 1 area (9h, 12h and 14 h . 5 ). Using a stacking analysis, we find a significant correlation between the far-infrared (FIR) luminosity and 1.4-GHz luminosity for radio-loud quasars. Partial correlation analysis confirms the intrinsic correlation after removing the redshift contribution, while for radio-quiet quasars, no partial correlation is found. Using a single-temperature grey-body model, we find a general trend of lower dust temperatures in case of radio-loud quasars compared to radio-quiet quasars. Also, radio-loud quasars are found to have almost constant mean values of dust mass along redshift and optical luminosity bins. In addition, we find that radio-loud quasars at lower optical luminosities tend to have on average higher FIR and 250-μm luminosity with respect to radio-quiet quasars with the same optical luminosites. Even if we use a two-temperature grey-body model to describe the FIR data, the FIR luminosity excess remains at lower optical luminosities. These results suggest that powerful radio jets are associated with star formation especially at lower accretion ratesPeer reviewe

The Spitzer High Redshift Radio Galaxy Survey

We present results from a comprehensive imaging survey of 70 radio galaxies at redshifts 1<z<5.2 using all three cameras onboard the Spitzer Space Telescope. The resulting spectral energy distributions unambiguously show a stellar population in 46 sources and hot dust emission associated with the active nucleus in 59. Using a new restframe S_3um/S_1.6um versus S_um/S_3um criterion, we identify 42 sources where the restframe 1.6um emission from the stellar population can be measured. For these radio galaxies, the median stellar mass is high, 2x10^11 M_sun, and remarkably constant within the range 13, there is tentative evidence for a factor of two decrease in stellar mass. This suggests that radio galaxies have assembled the bulk of their stellar mass by z~3, but confirmation by more detailed decomposition of stellar and AGN emission is needed. The restframe 500 MHz radio luminosities are only marginally correlated with stellar mass but are strongly correlated with the restframe 5um hot dust luminosity. This suggests that the radio galaxies have a large range of Eddington ratios. We also present new Very Large Array 4.86 and 8.46 GHz imaging of 14 radio galaxies and find that radio core dominance --- an indicator of jet orientation --- is strongly correlated with hot dust luminosity. While all of our targets were selected as narrow-lined, type 2 AGNs, this result can be understood in the context of orientation-dependent models if there is a continuous distribution of orientations from obscured type 2 to unobscured type 1 AGNs rather than a clear dichotomy. Finally, four radio galaxies have nearby (<6") companions whose mid-IR colors are suggestive of their being AGNs. This may indicate an association between radio galaxy activity and major mergers.Comment: 31 pages, 125 figures. Accepted for publication in the Astrophysical Journa

A search for debris disks in the Herschel -ATLAS

Aims. We aim to demonstrate that the Herschel-ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods. We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results. We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The UKIDSS project is defined in Lawrence et al. (2007). UKIDSS uses the UKIRT Wide Field Camera (WFCAM; Casali et al. 2007). The photometric system is described in Hambly et al. (2008), and the calibration is described in Hodgkin et al. (2009). The pipeline processing and science archive are described in Hambly et al. (2008). M.A.T. would like to thank two of our undergraduate project students, Sam Richards and Max Podger, who carried out initial database searches and also David Pinfield and Ralf Napiwotski for discussions on low mass stars

A new search for distant radio galaxies in the Southern hemisphere -- III. Optical spectroscopy and analysis of the MRCR--SUMSS sample

We have compiled a sample of 234 ultra-steep-spectrum(USS)-selected radio sources in order to find high-redshift radio galaxies (HzRGs). The sample is in the southern sky at -40 deg < DEC < -30 deg which is the overlap region of the 408-MHz Revised Molonglo Reference Catalogue, 843-MHz Sydney University Molonglo Sky Survey (the MRCR--SUMSS sample) and the 1400-MHz NRAO VLA Sky Survey. This is the third in a series of papers on the MRCR--SUMSS sample. Here we present optical spectra from the ANU 2.3-m telescope, ESO New Technology Telescope and ESO Very Large Telescope for 52 of the identifications from Bryant et al. (2009, Paper II), yielding redshifts for 36 galaxies, 13 of which have z>2. We analyse the K-z distribution and compare 4-arcsec-aperture magnitudes with 64-kpc aperture magnitudes in several surveys from the literature; the MRCR--SUMSS sample is found to be consistent with models for 10^{11}-10^{12} solar mass galaxies. Dispersions about the fits in the K-z plot support passive evolution of radio galaxy hosts since z>3. By comparing USS-selected samples in the literature, we find that the resultant median redshift of the samples shown is not dependent on the flux density distribution or selection frequency of each sample. In addition, our finding that the majority of the radio spectral energy distributions remain straight over a wide frequency range suggests that a k-correction is not responsible for the success of USS-selection in identifying high redshift radio galaxies and therefore the steep radio spectra may be intrinsic to the source or a product of the environment. Two galaxies have been found to have both compact radio structures and strong self-absorption in the Ly-alpha line, suggesting they are surrounded by a dense medium...abridged.Comment: Accepted for MNRAS. 25 page

Evolution of the Far-Infrared-Radio Correlation and Infrared SEDs of Massive Galaxies over z = 0 - 2

We investigate the far-infrared-radio correlation (FRC) of stellar-mass-selected galaxies in the Extended Chandra Deep Field South using far-infrared imaging from Spitzer and radio imaging from the Very Large Array and Giant Metre-Wave Radio Telescope. We stack in redshift bins to probe galaxies below the noise and confusion limits. Radio fluxes are K-corrected using observed flux ratios, leading to tentative evidence for an evolution in spectral index. We compare spectral energy distribution (SED) templates of local galaxies for K-correcting FIR fluxes, and show that the data are best fit by a quiescent spiral template (M51) rather than a warm starburst (M82) or ULIRG (Arp220), implying a predominance of cold dust in massive galaxies at high redshift. In contrast we measure total infrared luminosities that are consistent with high star-formation rates. We observe that the FRC index (q) does not evolve significantly over z=0-2 when computed from K-corrected 24 or 160-mum photometry, but that using 70-mum fluxes leads to an apparent decline in q beyond z~1. This suggests some change in the SED at high redshift, either a steepening of the spectrum at rest-frame ~25-35mum or a deficiency at ~70mum leading to a drop in the total infrared/radio ratios. We compare our results to other work in the literature and find synergies with recent findings on the high-redshift FRC, high specific star-formation rates of massive galaxies and the cold dust temperatures in these galaxies.Comment: 21 pages, 13 figures, accepted for publication in MNRA

Dissecting the roles of mass and environment quenching in galaxy evolution with EAGLE

We exploit the pioneering cosmological hydrodynamical simulation, EAGLE, to study how the connection between halo mass (M_halo), stellar mass (M*) and star-formation rate (SFR) evolves across redshift. Using Principal Component Analysis we identify the key axes of correlation between these physical quantities, for the full galaxy sample and split by satellite/central and low/high halo mass. The first principal component of the z=0 EAGLE galaxy population is a positive correlation between M_halo, M* and SFR. This component is particularly dominant for central galaxies in low mass haloes. The second principal component, most significant in high mass haloes, is a negative correlation between M_halo and SFR, indicative of environmental quenching. For galaxies above M*~10^10M_solar, however, the SFR is seen to decouple from the M_halo-M* correlation; this result is found to be independent of environment, suggesting that mass quenching effects are also in operation. We find extremely good agreement between the EAGLE principal components and those of SDSS galaxies; this lends confidence to our conclusions. Extending our study to EAGLE galaxies in the range z=0-4, we find that, although the relative numbers of galaxies in the different subsamples change, their principal components do not change significantly with redshift. This indicates that the physical processes that govern the evolution of galaxies within their dark matter haloes act similarly throughout cosmic time. Finally, we present halo occupation distribution model fits to EAGLE galaxies and show that one flexible 6-parameter functional form is capable of fitting a wide range of different mass- and SFR-selected subsamples.Comment: 17 pages, 5 figures; accepted for publication in MNRA