Source author record

J. S. Dunlop

J. S. Dunlop 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

121works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

121 published item(s)

preprint2023arXiv

The ultraviolet continuum slopes ($\mathbfβ$) of galaxies at $\mathbf{z\simeq8-16}$ from JWST and ground-based near-infrared imaging

We study the rest-frame ultraviolet (UV) continuum slopes ($β$) of galaxies at redshifts $8 < z < 16$, using a combination of JWST ERO and ERS NIRcam imaging and ground-based near-infrared imaging of the COSMOS field. The combination of JWST and ground-based imaging provides a wide baseline in both redshift and absolute UV magnitude ($-22.6 < M_{\rm UV} < 17.9$), sufficient to allow a meaningful comparison to previous results at lower redshift. Using a power-law fitting technique, we find that our full sample (median $M_{\rm UV}=-19.3\pm 1.3$) returns an inverse-variance weighted mean value of $\langle β\rangle = -2.10 \pm 0.05$, with a corresponding median value of $β=-2.29\pm 0.09$. These values imply that the UV colours of galaxies at $z>8$ are, on average, no bluer than the bluest galaxies in the local Universe. Moreover, we find evidence for a $β-M_{\rm UV}$ relation, such that brighter UV galaxies display redder UV slopes ($\rm{d}β/ \rm{d} M_{\rm UV} = -0.17 \pm 0.05$). Comparing to results at lower redshift, we find that the slope of our $β-M_{\rm UV}$ relation is consistent with the slope observed at $z\simeq 5$ and that, at a given $M_{\rm UV}$, our $8<z<16$ galaxies are somewhat bluer than their $z\simeq 5$ counterparts, with an inverse-variance weighted mean offset of $\langle Δβ\rangle = -0.38 \pm 0.09$. We do not find strong evidence that any objects in our sample display ultra-blue UV continuum slopes (i.e., $β\lesssim-3$) that would require their UV emission to be dominated by ultra-young, dust-free stellar populations with high Lyman-continuum escape fractions. Comparing our results to the predictions of theoretical galaxy formation models, we find that the galaxies in our sample are consistent with the young, metal-poor and moderately dust-reddened galaxies expected at $z>8$.

preprint2022arXiv

A combined VANDELS and LEGA-C study: the evolution of quiescent galaxy size, stellar mass and age from $\mathbf{\textit{z} = 0.6}$ to $\mathbf{\textit{z} = 1.3}$

We study the relationships between stellar mass, size and age within the quiescent population, using two mass-complete spectroscopic samples with $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})>10.3$, taken from VANDELS at $1.0<z<1.3$, and LEGA-C at $0.6<z<0.8$. Using robust D$_{n}$4000 values, we demonstrate that the well-known 'downsizing' signature is already in place by $z\simeq1.1$, with D$_{n}$4000 increasing by $\simeq0.1$ across a $\simeq$ 1 dex mass interval for both VANDELS and LEGA-C. We then proceed to investigate the evolution of the quiescent galaxy stellar mass-size relation from $z\simeq1.1$ to $z\simeq0.7$. We find the median size increases by a factor of $1.9\pm{0.1}$ at $\mathrm{log_{10}}(M_{\star}/\mathrm{M_{\odot}})=10.5$, and see tentative evidence for flattening of the relation, finding slopes of $α=0.72\pm0.06$ and $α=$ $0.56\pm0.04$ for VANDELS and LEGA-C respectively. We finally split our sample into galaxies above and below our fitted mass-size relations, to investigate how size and D$_{n}$4000 correlate. For LEGA-C, we see a clear difference, with larger galaxies found to have smaller D$_{n}$4000 at fixed stellar mass. Due to the faintness and smaller numbers of the VANDELS sample, we cannot confirm whether a similar relation exists at $z\simeq1.1$. We consider whether differences in stellar age or metallicity are most likely to drive this size-D$_{n}$4000 relation, finding that any metallicity differences are unlikely to fully explain the observed offset, meaning smaller galaxies must be older than their larger counterparts. We find the observed evolution in size, mass and D$_{n}$4000 across the $\simeq2$ Gyr from $z\sim1.1$ to $z\sim0.7$ can be explained by a simple toy model in which VANDELS galaxies evolve passively, whilst experiencing a series of minor mergers.

preprint2022arXiv

The discovery of rest-frame UV colour gradients and a diversity of dust morphologies in bright z ~ 7 Lyman-break galaxies

We present deep ALMA dust continuum observations for a sample of luminous ($M_{\rm UV} < -22$) star-forming galaxies at $z \simeq 7$. We detect five of the six sources in the far-infrared (FIR), providing key constraints on the obscured star-formation rate (SFR) and the infrared-excess-$β$ (IRX-$β$) relation without the need for stacking. Despite the galaxies showing blue rest-frame UV slopes ($β\simeq -2$) we find that 35-75 percent of the total SFR is obscured. We find the IRX-$β$ relation derived for these $z \simeq 7$ sources is consistent with that found for local star-burst galaxies. Using our relatively high-resolution (FWHM $\simeq 0.7\,{\rm arcsec}$) observations we identify a diversity of dust morphologies in the sample. We find both compact emission that appears offset relative to the unobscured components and extended dust emission that is co-spatial with the rest-frame UV light. In the majority of the sources we detect strong rest-frame UV colour gradients (with up to $Δβ\simeq 0.7$-$1.4$) as probed by the multi-band UltraVISTA ground-based data. The observed redder colours are spatially correlated with the location of the FIR detection. Our results show that even in bright Lyman-break galaxies at $z \simeq 7$ the peak of the star-formation is typically hosted by the fainter, redder, regions in the rest-frame UV, which have an obscured fraction of $f_{\rm obs} \ge 0.8$. As well as demonstrating the importance of dust obscured star-formation within the Epoch of Reionization, these observations provide an exciting taster of the rich spatially resolved datasets that will be obtained from JWST and high-resolution ALMA follow-up at these redshifts.

preprint2022arXiv

The stellar metallicities of massive quiescent galaxies at 1.0 < z < 1.3 from KMOS+VANDELS

We present a rest-frame UV-optical stacked spectrum representative of massive quiescent galaxies at $1.0<z<1.3$ with log$(M_*/\rm{M_\odot})>10.8$. The stack is constructed using VANDELS survey data, combined with new KMOS observations. We apply two independent full-spectral-fitting approaches, measuring a total metallicity, [Z/H]=$-0.13\pm0.08$ with Bagpipes, and [Z/H]=$0.04\pm0.14$ with Alf, a fall of $\sim0.2-0.3$ dex compared with the local Universe. We also measure an iron abundance, [Fe/H] =$-0.18\pm0.08$, a fall of $\sim0.15$ dex compared with the the local Universe. We measure the alpha enhancement via the magnesium abundance, obtaining [Mg/Fe]=$0.23\pm$0.12, consistent with similar-mass galaxies in the local Universe, indicating no evolution in the average alpha enhancement of log$(M_*/\rm{M_\odot})=11$ quiescent galaxies over the last $\sim8$ Gyr. This suggests the very high alpha enhancements recently reported for several bright $z\sim1-2$ quiescent galaxies are due to their extreme masses, log$(M_*/\rm{M_\odot})\gtrsim11.5$, rather than being typical of the $z\gtrsim1$ population. The metallicity evolution we observe with redshift (falling [Z/H], [Fe/H], constant [Mg/Fe]) is consistent with recent studies. We recover a mean stellar age of $2.5^{+0.6}_{-0.4}$ Gyr, corresponding to a formation redshift, $z_\rm{form}=2.4^{+0.6}_{-0.3}$. Recent studies have obtained varying average formation redshifts for $z\gtrsim1$ massive quiescent galaxies, and, as these studies report consistent metallicities, we identify different star-formation-history models as the most likely cause. Larger spectroscopic samples from upcoming ground-based instruments will provide precise constraints on ages and metallicities at $z\gtrsim1$. Combining these with precise JWST $z>2$ quiescent-galaxy stellar-mass functions will provide an independent test of formation redshifts derived from spectral fitting.

preprint2022arXiv

The VANDELS survey: a measurement of the average Lyman-continuum escape fraction of star-forming galaxies at z=3.5

We present a study designed to measure the average LyC escape fraction ($\langle f_{\rm esc}\rangle$) of star-forming galaxies at z=3.5. We assemble a sample of 148 galaxies from the VANDELS survey at $3.35\leq z_{\rm spec}\leq3.95$, selected to minimize line-of-sight contamination of their photometry. For this sample, we use ultra-deep, ground-based, $U-$band imaging and HST $V-$band imaging to robustly measure the distribution of $\mathcal{R_{\rm obs}}$ $=(L_{\rm LyC}/L_{\rm UV})_{\rm obs}$. We then model the distribution as a function of $\langle f_{\rm esc}\rangle$, carefully accounting for attenuation by dust, and the IGM (and CGM). A maximum likelihood fit to the $\mathcal{R_{\rm obs}}$ distribution returns a best-fitting value of $\langle f_{\rm esc}\rangle =0.07\pm0.02$, a result confirmed using an alternative Bayesian inference technique (both exclude $\langle f_{\rm esc}\rangle=0.0$ at $> 3σ$). By splitting our sample in two, we find evidence that $\langle f_{\rm esc}\rangle$ is positively correlated with Ly$α$ equivalent width, with high and low sub-samples returning best fits of $\langle f_{\rm esc}\rangle=0.12^{+0.06}_{-0.04}$ and $\langle f_{\rm esc} \rangle=0.02^{+0.02}_{-0.01}$, respectively. In contrast, we find evidence that $\langle f_{\rm esc}\rangle$ is anti-correlated with intrinsic UV luminosity and UV dust attenuation; with low UV luminosity and dust attenuation sub-samples returning best fits in the range $0.10 \leq \langle f_{\rm esc}\rangle \leq 0.22$. We do not find evidence for a clear correlation between $f_{\rm esc}$ and galaxy stellar mass, suggesting it is not a primary indicator of leakage. Although larger samples are needed to further explore these trends, they suggest that it is entirely plausible that the low dust and metallicity galaxies found at z > 6 will display the $\langle f_{\rm esc}\rangle\geq0.1$ required to drive reionization.

preprint2021arXiv

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

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

preprint2020arXiv

A lack of evolution in the very bright-end of the galaxy luminosity function from z = 8-10

We utilize deep near-infrared survey data from the UltraVISTA fourth data release (DR4) and the VIDEO survey, in combination with overlapping optical and Spitzer data, to search for bright star-forming galaxies at $z \gtrsim 7.5$. Using a full photometric redshift fitting analysis applied to the $\sim 6\,{\rm deg}^2$ of imaging searched, we find 27 Lyman-break galaxies (LBGs), including 20 new sources, with best-fitting photometric redshifts in the range $7.4 < z < 9.1$. From this sample we derive the rest-frame UV luminosity function (LF) at $z = 8$ and $z = 9$ out to extremely bright UV magnitudes ($M_{\rm UV} \simeq -23$) for the first time. We find an excess in the number density of bright galaxies in comparison to the typically assumed Schechter functional form derived from fainter samples. Combined with previous studies at lower redshift, our results show that there is little evolution in the number density of very bright ($M_{\rm UV} \sim -23$) LBGs between $z \simeq 5$ and $z\simeq 9$. The tentative detection of an LBG with best-fit photometric redshift of $z = 10.9 \pm 1.0$ in our data is consistent with the derived evolution. We show that a double power-law fit with a brightening characteristic magnitude ($ΔM^*/Δz \simeq -0.5$) and a steadily steepening bright-end slope ($Δβ/Δz \simeq -0.5$) provides a good description of the $z > 5$ data over a wide range in absolute UV magnitude ($-23 < M_{\rm UV} < -17$). We postulate that the observed evolution can be explained by a lack of mass quenching at very high redshifts in combination with increasing dust obscuration within the first $\sim 1 \,{\rm Gyr}$ of galaxy evolution.

preprint2020arXiv

ALMA 26 arcmin$^2$ survey of GOODS-S at one-millimeter (ASAGAO): millimeter properties of stellar mass selected galaxies

We make use of the ASAGAO, deep 1.2 mm continuum observations of a 26 arcmin$^2$ region in the GOODS-South field obtained with ALMA, to probe dust-enshrouded star formation in $K$-band selected (i.e., stellar mass selected) galaxies, which are drawn from the ZFOURGE catalog. Based on the ASAGAO combined map, which was created by combining ASAGAO and ALMA archival data in the GOODS-South field, we find that 24 ZFOURGE sources have 1.2 mm counterparts with a signal-to-noise ratio $>$ 4.5 (1$σ\simeq$ 30 - 70 $μ$Jy beam$^{-1}$ at 1.2 mm). Their median redshift is estimated to be $z_\mathrm{median}=$ 2.38 $\pm$ 0.14. They generally follow the tight relationship of the stellar mass versus star formation rate (i.e., the main sequence of star-forming galaxies). ALMA-detected ZFOURGE sources exhibit systematically larger infrared (IR) excess (IRX $\equiv L_\mathrm{IR}/L_\mathrm{UV}$) compared to ZFOURGE galaxies without ALMA detections even though they have similar redshifts, stellar masses, and star formation rates. This implies the consensus stellar-mass versus IRX relation, which is known to be tight among rest-frame-UV-selected galaxies, can not fully predict the ALMA detectability of stellar-mass-selected galaxies. We find that ALMA-detected ZFOURGE sources are the main contributors to the cosmic IR star formation rate density at $z$ = 2 - 3.

preprint2020arXiv

An ALMA survey of the SCUBA-2 Cosmology Legacy Survey UKIDSS/UDS field: Dust attenuation in high-redshift Lyman break Galaxies

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

preprint2020arXiv

Extended H$α$ over compact far-infrared continuum in dusty submillimeter galaxies -- Insights into dust distributions and star-formation rates at $z\sim2$

Using data from ALMA and near-infrared (NIR) integral field spectrographs including both SINFONI and KMOS on the VLT, we investigate the two-dimensional distributions of H$α$ and rest-frame far-infrared (FIR) continuum in six submillimeter galaxies at $z\sim2$. At a similar spatial resolution ($\sim$0.5" FWHM; $\sim$4.5 kpc at $z=2$), we find that the half-light radius of H$α$ is significantly larger than that of the FIR continuum in half of the sample, and on average H$α$ is a median factor of $2.0\pm0.4$ larger. Having explored various ways to correct for the attenuation, we find that the attenuation-corrected H$α$-based SFRs are systematically lower than the IR-based SFRs by at least a median factor of $3\pm1$, which cannot be explained by the difference in half-light radius alone. In addition, we find that in 40% of cases the total $V$-band attenuation ($A_V$) derived from energy balance modeling of the full ultraviolet(UV)-to-FIR spectral energy distributions (SEDs) is significantly higher than that derived from SED modeling using only the UV-to-NIR part of the SEDs, and the discrepancy appears to increase with increasing total infrared luminosity. Finally, considering all our findings along with the studies in the literature, we postulate that the dust distributions in SMGs, and possibly also in less IR luminous $z\sim2$ massive star-forming galaxies, can be decomposed into three main components; the diffuse dust heated by older stellar populations, the more obscured and extended young star-forming HII regions, and the heavily obscured central regions that have a low filling factor but dominate the infrared luminosity in which the majority of attenuation cannot be probed via UV-to-NIR emissions.

preprint2020arXiv

The rest-frame UV luminosity function at $z \simeq 4$: a significant contribution of AGN to the bright-end of the galaxy population

We measure the rest-frame UV luminosity function (LF) at $z \sim 4$ self-consistently over a wide range in absolute magnitude ($-27 \lesssim M_{\rm UV} \lesssim -20$). The LF is measured with 46,904 sources selected using a photometric redshift approach over $ \sim 6$ deg$^2$ of the combined COSMOS and XMM-LSS fields. We simultaneously fit for both AGN and galaxy LFs using a combination of Schechter or Double Power Law (DPL) functions alongside a single power law for the faint-end slope of the AGN LF. We find a lack of evolution in the shape of the bright-end of the LBG component when compared to other studies at $z \simeq 5$ and evolutionary recipes for the UV LF. Regardless of whether the LBG LF is fit with a Schechter function or DPL, AGN are found to dominate at $M_{\rm UV} < -23.5$. We measure a steep faint-end slope of the AGN LF with $α_{AGN} = -2.09^{+0.35}_{-0.38}$ ($-1.66^{+0.29}_{-0.58}$) when fit alongside a Schechter function (DPL) for the galaxies. Our results suggest that if AGN are morphologically selected it results in a bias to lower number densities. Only by considering the full galaxy population over the transition region from AGN to LBG domination can an accurate measurement of the total LF be attained.

preprint2020arXiv

The VANDELS survey: A strong correlation between Ly$α$ equivalent width and stellar metallicity at $\mathbf{3\leq z \leq 5}$

We present the results of a new study investigating the relationship between observed Ly$α$ equivalent width ($W_λ$(Ly$α$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_λ$(Ly$α$) across the range $-58 \unicode{xC5} \lesssim$ $W_λ$(Ly$α$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_λ$(Ly$α$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_λ$(Ly$α$) quartile ($\langle$$W_λ$(Ly$α$)$\rangle=-18\unicode{xC5}$) and the highest $W_λ$(Ly$α$) quartile ($\langle$$W_λ$(Ly$α$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_λ$(Ly$α$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_λ$(Ly$α$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_λ$(Ly$α$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

preprint2020arXiv

Timing the earliest quenching events with a robust sample of massive quiescent galaxies at 2 < z < 5

We present a sample of 151 massive ($M_* > 10^{10}\mathrm{M_\odot}$) quiescent galaxies at $2 < z < 5$, based on a sophisticated Bayesian spectral energy distribution fitting analysis of the CANDELS UDS and GOODS-South fields. Our sample includes a robust sub-sample of 61 objects for which we confidently exclude low-redshift and star-forming solutions. We identify 10 robust objects at $z>3$, of which 2 are at $z>4$. We report formation redshifts, demonstrating that the oldest objects formed at $z > 6$, however individual ages from our photometric data have significant uncertainties, typically $\sim0.5$ Gyr. We demonstrate that the UVJ colours of the quiescent population evolve with redshift at $z>3$, becoming bluer and more similar to post-starburst galaxies at lower redshift. Based upon this we construct a model for the time-evolution of quiescent galaxy UVJ colours, concluding that the oldest objects are consistent with forming the bulk of their stellar mass at $z\sim6-7$ and quenching at $z\sim5$. We report spectroscopic redshifts for two of our objects at $z=3.440$ and $3.396$, which exhibit extremely weak Ly$α$ emission in ultra-deep VANDELS spectra. We calculate star-formation rates based on these line fluxes, finding that these galaxies are consistent with our quiescent selection criteria, provided their Ly$α$ escape fractions are $>3$ and $>10$ per cent respectively. We finally report that our highest-redshift robust object exhibits a continuum break at $λ\sim7000$A in a spectrum from VUDS, consistent with our photometric redshift of $z_\mathrm{phot}=4.72^{+0.06}_{-0.04}$. If confirmed as quiescent this object would be the highest-redshift known quiescent galaxy. To obtain stronger constraints on the times of the earliest quenching events, high-SNR spectroscopy must be extended to $z\gtrsim3$ quiescent objects.

preprint2016arXiv

Chandra counterparts of CANDELS GOODS-S sources

Improving the capabilities of detecting faint X-ray sources is fundamental to increase the statistics on faint high-z AGN and star-forming galaxies. We performed a simultaneous Maximum Likelihood PSF fit in the [0.5-2] keV and [2-7] keV energy bands of the 4 Ms{\em Chandra} Deep Field South (CDFS) data at the position of the 34930 CANDELS H-band selected galaxies. For each detected source we provide X-ray photometry and optical counterpart validation. We validated this technique by means of a raytracing simulation. We detected a total of 698 X-ray point-sources with a likelihood $\mathcal{L}$$>$4.98 (i.e. $>$2.7$σ$). We show that the prior knowledge of a deep sample of Optical-NIR galaxies leads to a significant increase of the detection of faint (i.e. $\sim$10$^{-17}$ cgs in the [0.5-2] keV band) sources with respect to "blind" X-ray detections. By including previous X-ray catalogs, this work increases the total number of X-ray sources detected in the 4 Ms CDFS, CANDELS area to 793, which represents the largest sample of extremely faint X-ray sources assembled to date. Our results suggest that a large fraction of the optical counterparts of our X-ray sources determined by likelihood ratio actually coincides with the priors used for the source detection. Most of the new detected sources are likely star-forming galaxies or faint absorbed AGN. We identified a few sources sources with putative photometric redshift z$>$4. Despite the low number statistics and the uncertainties on the photo-z, this sample significantly increases the number of X--ray selected candidate high-z AGN.

preprint2016arXiv

Changing physical conditions in star-forming galaxies between redshifts 0 < z < 4: [OIII]/Hb evolution

We investigate the redshift evolution of the [OIII]/Hb nebular emission line ratio for a sample of galaxies spanning the redshift range 0 < z < 4. We compare the observed evolution to a set of theoretical models which account for the independent evolution of chemical abundance, ionization parameter and interstellar-medium (ISM) pressure in star-forming galaxies with redshift. Accounting for selection effects in the combined datasets, we show that the evolution to higher [OIII]/Hb ratios with redshift is a real physical effect which is best accounted for by a model in which the ionization parameter is elevated from the average values typical of local star-forming galaxies, with a possible simultaneous increase in the ISM pressure. We rule out the possibility that the observed [OIII]/Hb evolution is purely due to metallicity evolution. We discuss the implications of these results for using local empirical metallicity calibrations to measure metallicities at high redshift, and briefly discuss possible theoretical implications of our results.

preprint2016arXiv

Non-parametric analysis of the rest-frame UV sizes and morphological disturbance amongst L* galaxies at 4<z<8

We present the results of a study investigating the sizes and morphologies of redshift 4 < z < 8 galaxies in the CANDELS GOODS-S, HUDF and HUDF parallel fields. Based on non-parametric measurements and incorporating a careful treatment of measurement biases, we quantify the typical size of galaxies at each redshift as the peak of the log-normal size distribution, rather than the arithmetic mean size. Parameterizing the evolution of galaxy half-light radius as $r_{50} \propto (1+z)^n$, we find $n = -0.20 \pm 0.26$ at bright UV-luminosities ($0.3L_{*(z=3)} < L < L_*$) and $n = -0.47 \pm 0.62$ at faint luminosities ($0.12L_* < L < 0.3L_*$). Furthermore, simulations based on artificially redshifting our z~4 galaxy sample show that we cannot reject the null hypothesis of no size evolution. We show that this result is caused by a combination of the size-dependent completeness of high-redshift galaxy samples and the underestimation of the sizes of the largest galaxies at a given epoch. To explore the evolution of galaxy morphology we first compare asymmetry measurements to those from a large sample of simulated single Sérsic profiles, in order to robustly categorise galaxies as either `smooth' or `disturbed'. Comparing the disturbed fraction amongst bright ($M_{UV} \leq -20$) galaxies at each redshift to that obtained by artificially redshifting our z~4 galaxy sample, while carefully matching the size and UV-luminosity distributions, we find no clear evidence for evolution in galaxy morphology over the redshift interval 4 < z < 8. Therefore, based on our results, a bright ($M_{UV} \leq -20$) galaxy at z~6 is no more likely to be measured as `disturbed' than a comparable galaxy at z~4, given the current observational constraints.

preprint2016arXiv

SXDF-UDS-CANDELS-ALMA 1.5 arcmin$^2$ deep survey

We have conducted 1.1 mm ALMA observations of a contiguous $105'' \times 50''$ or 1.5 arcmin$^2$ window in the SXDF-UDS-CANDELS. We achieved a 5$σ$ sensitivity of 0.28 mJy, providing a flat sensus of dusty star-forming galaxies with $L_{\rm IR} \sim6\times10^{11}$ $L_\odot$ (for $T_{\rm dust}$ =40K) up to $z\sim10$ thanks to the negative K-correction at this wavelength. We detected 5 brightest sources (S/N$>$6) and 18 low-significance sources (5$>$S/N$>$4; these may contain spurious detections, though). One of the 5 brightest ALMA sources ($S_{\rm 1.1mm} = 0.84 \pm 0.09$ mJy) is extremely faint in the WFC3 and VLT/HAWK-I images, demonstrating that a contiguous ALMA imaging survey is able to uncover a faint dust-obscured population that is invisible in deep optical/near-infrared surveys. We found a possible [CII]-line emitter at $z=5.955$ or a low-$z$ CO emitting galaxy within the field, which may allow us to constrain the [CII] and/or the CO luminosity functions across the history of the universe.

preprint2016arXiv

T-PHOT version 2.0: improved algorithms for background subtraction, local convolution, kernel registration, and new options

We present the new release v2.0 of T-PHOT, a publicly available software package developed to perform PSF-matched, prior-based, multiwavelength deconfusion photometry of extragalactic fields. New features included in the code are presented and discussed: background estimation, fitting using position dependent kernels, flux prioring, diagnostical statistics on the residual image, exclusion of selected sources from the model and residual images, individual registration of fitted objects. These new options improve on the performance of the code, allowing for more accurate results and providing useful aids for diagnostics.

preprint2016arXiv

The ASTRODEEP Frontier Fields Catalogues: I - Multiwavelength photometry of Abell-2744 and MACS-J0416

We present multiwavelength photometric catalogues (HST, Spitzer and Hawk-I K band) for the first two of the Frontier Fields, Abell2744 and MACSJ0416 (plus their parallel fields). To detect faint sources even in the central regions of the clusters, we develop a robust and repeatable procedure that uses the public codes Galapagos and Galfit to model and remove most of the light contribution from both the brightest cluster members as well as the ICL. We perform the detection on the HST H160 processed image to obtain a pure H-selected sample. We also add a sample of sources which are undetected in the H160 image but appear on a stacked infrared image. Photometry in the other HST bands is obtained using SExtractor, performed again on residual images after the Galfit procedure for foreground light removal. Photometry on the Hawk-I and IRAC bands has been obtained using our PSF-matching deconfusion code T-PHOT. A similar procedure, but without the need for the foreground light removal, is adopted for the Parallel fields. The procedure allows for the detection and the photometric measurements of ~2500 sources per field. We deliver and release complete photometric H-detected catalogues, with the addition of a complementary sample of infrared-detected sources. All objects have multiwavelength coverage including B to H HST bands, plus K band from Hawk-I, and 3.6 - 4.5 μm from Spitzer. Full and detailed treatment of photometric errors is included. We perform basic sanity checks on the reliability of our results. The multiwavelength catalogues are publicly available and are ready to be used for scientific purposes. Our procedures allows for the detection of outshined objects near the bright galaxies, which, coupled with the magnification effect of the clusters, can reveal extremely faint high redshift sources. Full analysis on photometric redshifts is presented in a companion Paper II. [abridged]

preprint2016arXiv

The ASTRODEEP Frontier Fields Catalogues: II - Photometric redshifts and rest-frame properties in Abell-2744 and MACS-J0416

We present the first public release of photometric redshifts, galaxy rest-frame properties and associated magnification values in the cluster and parallel pointings of the first two Frontier Fields, Abell-2744 and MACS-J0416. We exploit a multi-wavelength catalogue ranging from HST to ground-based K and Spitzer IRAC which is specifically designed to enable detection and measurement of accurate fluxes in crowded cluster regions. The multi-band information is used to derive photometric redshifts and physical properties of sources detected either in the H-band image alone or from a stack of four WFC3 bands. To minimize systematics median photometric redshifts are assembled from six different approaches to photo-z estimates. Their reliability is assessed through a comparison with available spectroscopic samples. State of the art lensing models are used to derive magnification values on an object-by-object basis by taking into account sources positions and redshifts. We show that photometric redshifts reach a remarkable ~3-5% accuracy. After accounting for magnification the H band number counts are found in agreement at bright magnitudes with number counts from the CANDELS fields, while extending the presently available samples to galaxies intrinsically as faint as H160~32-33 thanks to strong gravitational lensing. The Frontier Fields allow to probe the galaxy stellar mass distribution at 0.5-1.5 dex lower masses, depending on magnification, with respect to extragalactic wide fields, including sources at Mstar~ 10^7-10^8 Msun at z>5. Similarly, they allow the detection of objects with intrinsic SFRs>1dex lower than in the CANDELS fields reaching 0.1-1 Msun/yr at z~6-10. [abridged]

preprint2016arXiv

The evolution of the equivalent width of the Ha emission line and specific star-formation rate in star-forming galaxies at 1<z<5

We present the results of a study which uses spectral energy distribution (SED) fitting to investigate the evolution of the equivalent width (EW) of the Halpha emission line in star-forming galaxies over the redshift interval 1<z<5. After first demonstrating the ability of our SED-fitting technique to recover EW(Ha) using a sample of galaxies at z~1.3 with EW(Ha) measurements from 3D-HST grism spectroscopy, we proceed to apply our technique to samples of spectroscopically confirmed and photometric-redshift selected star-forming galaxies at z>=1 in the CANDELS UDS and GOODS-S fields. Confining our analysis to a constant stellar mass range (9.5<log(M/Msun)<10.5), we find that the median EW(Ha) evolves only modestly with redshift, reaching a rest-frame value of EW(Ha)=301+/-30 Angs by redshift z~4.5. Furthermore, using estimates of star-formation rate (SFR) based on both UV luminosity and Ha line flux, we use our galaxy samples to compare the evolution of EW(Ha) and specific star-formation rate (sSFR). Our results indicate that over the redshift range 1<z<5, the evolution displayed by EW(Ha) and sSFR is consistent, and can be adequately parameterized as: propto (1+z)^(1.0+/-0.2). As a consequence, over this redshift range we find that the sSFR and rest-frame EW(Ha) of star-forming galaxies with stellar masses M~10^(10) Msun are related by: EW(Ha)/Ang=(63+/-7)sSFR/Gyr^(-1). Given the current uncertainties in measuring the SFRs of high-redshift galaxies, we conclude that EW(Ha) provides a useful independent tracer of sSFR for star-forming galaxies out to redshifts of z=5.

preprint2016arXiv

The Herschel-ATLAS Data Release 1 Paper II: Multi-wavelength counterparts to submillimetre sources

This paper is the second in a pair of articles presenting data release 1 (DR1) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), the largest single open-time key project carried out with the Herschel Space Observatory. The H-ATLAS is a wide-area imaging survey carried out in five photometric bands at 100, 160, 250, 350 and 500$μ$m covering a total area of 600deg$^2$. In this paper we describe the identification of optical counterparts to submillimetre sources in DR1, comprising an area of 161 deg$^2$ over three equatorial fields of roughly 12$^\circ$x4.5$^\circ$ centred at 9$^h$, 12$^h$ and 14.5$^h$ respectively. Of all the H-ATLAS fields, the equatorial regions benefit from the greatest overlap with current multi-wavelength surveys spanning ultraviolet (UV) to mid-infrared regimes, as well as extensive spectroscopic coverage. We use a likelihood-ratio technique to identify SDSS counterparts at r<22.4 for 250-$μ$m-selected sources detected at $\geq$ 4$σ$ ($\approx$28mJy). We find `reliable' counterparts (reliability R$\geq$0.8) for 44,835 sources (39 per cent), with an estimated completeness of 73.0 per cent and contamination rate of 4.7 per cent. Using redshifts and multi-wavelength photometry from GAMA and other public catalogues, we show that H-ATLAS-selected galaxies at $z<0.5$ span a wide range of optical colours, total infrared (IR) luminosities, and IR/UV ratios, with no strong disposition towards mid-IR-classified AGN in comparison with optical selection. The data described herein, together with all maps and catalogues described in the companion paper (Valiante et al. 2016), are available from the H-ATLAS website at www.h-atlas.org.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: galaxies in the deep 850-micron survey, and the star-forming `main sequence'

We investigate the properties of the galaxies selected from the deepest 850-micron survey undertaken to date with SCUBA-2 on the JCMT. This deep 850-micron imaging was taken in parallel with deep 450-micron imaging in the very best observing conditions as part of the SCUBA-2 Cosmology Legacy Survey. A total of 106 sources were uncovered at 850 microns from ~150, sq. arcmin in the centre of the COSMOS/UltraVISTA/CANDELS field, imaged to a typical rms depth of ~0.25 mJy. We utilise the wealth of available deep multi-frequency data to establish the complete redshift distribution for this sample, yielding <z> = 2.38 +- 0.09, a mean redshift comparable with that derived for all but the brightest previous sub-mm samples. We have also been able to establish the stellar masses of the majority of the galaxy identifications, enabling us to explore their location on the star-formation-rate:stellar-mass (SFR:M*) plane. Crucially, our new deep sample reaches flux densities equivalent to SFR ~ 100 Msun/yr, enabling us to confirm that sub-mm galaxies form the high-mass end of the `main sequence' (MS) of star-forming galaxies at z > 1.5 (with a mean specific SFR of sSFR = 2.25 +- 0.19 /Gyr at z ~ 2.5). Our results are consistent with no significant flattening of the MS towards high masses at these redshifts, suggesting that reports of such flattening possibly arise from under-estimates of dust-enshrouded star-formation activity in massive star-forming galaxies. However, our findings add to the growing evidence that average sSFR rises only slowly at high redshift, resulting in log(sSFR) being an apparently simple linear function of the age of the Universe.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: The EGS deep field I - Deep number counts and the redshift distribution of the recovered Cosmic Infrared Background at 450 and 850 um

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $σ_{450}=1.2$ mJy/beam and $σ_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

preprint2016arXiv

The z = 9-10 galaxy population in the Hubble Frontier Fields and CLASH surveys: The z=9 LF and further evidence for a smooth decline in UV luminosity at z >= 8

We present the results of a search for z=9-10 galaxies within the first 8 pointings of the Hubble Frontier Fields (HFF) (4 clusters plus 4 parallel fields) and 20 cluster fields from the CLASH survey. Combined with our previous analysis of the Hubble Ultra-Deep field (HUDF), we have now completed a search for z=9-10 galaxies over ~130 sq. arcmin, across 29 HST WFC3/IR pointings. As in our recent study of the first two HFF fields, we confine our primary search for high-redshift candidates in the HFF imaging to the uniformly deep (i.e. sigma_160>30 AB mag in 0.5-arcsec diameter apertures), relatively low magnification regions. In the CLASH fields our search was confined to uniformly deep regions where sigma_160>28.8 AB mag. Our SED fitting analysis unveils a sample of 33 galaxy candidates at z_phot>=8.4, five of which have primary photometric redshift solutions in the range 9.6<z_phot<11.2. By calculating a de-lensed effective volume for each candidate, the improved statistics and reduced cosmic variance provided by our new sample allows a more accurate determination of the UV-selected galaxy luminosity function (LF) at z~9. Our new results strengthen our previous conclusion that the LF appears to evolve smoothly from z=8 to z=9, an evolution which can be equally well modelled by a factor of ~2 drop in density, or a dimming of ~0.5 mag in M*. Moreover, based on our new sample, we are able to place initial constraints on the z=10 LF, finding that the number density at M_1500 ~ -19.7 is log(phi) = -4.1 (+0.2,-0.3), a factor of ~2 lower than at z=9. Finally, we use our new results to re-visit the issue of the decline in UV luminosity density at z>=8. We conclude that the data continue to support a smooth decline in rho_UV over the redshift interval 6<z<10, in agreement with simple models of early galaxy evolution driven by the growth in the underlying dark matter halo mass function.

preprint2016arXiv

VLA and ALMA Imaging of Intense, Galaxy-Wide Star Formation in z ~ 2 Galaxies

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $λ= 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

preprint2015arXiv

Cold Dust Emission from X-ray AGN in the SCUBA-2 Cosmology Legacy Survey: Dependence on Luminosity, Obscuration & AGN Activity

We study the 850um emission in X-ray selected AGN in the 2 sq-deg COSMOS field using new data from the SCUBA-2 Cosmology Legacy Survey. We find 19 850um bright X-ray AGN in a high-sensitivity region covering 0.89 sq-deg with flux densities of S850=4-10 mJy. The 19 AGN span the full range in redshift and hard X-ray luminosity covered by the sample - 0.7<z<3.5 and 43.2<log10(LX) <45. We report a highly significant stacked 850um detection of a hard X-ray flux-limited population of 699 z>1 X-ray AGN - S850=0.71+/-0.08mJy. We explore trends in the stacked 850um flux densities with redshift, finding no evolution in the average cold dust emission over the redshift range probed. For Type 1 AGN, there is no significant correlation between the stacked 850um flux and hard X-ray luminosity. However, in Type 2 AGN the stacked submm flux is a factor of 2 higher at high luminosities. When averaging over all X-ray luminosities, no significant differences are found in the stacked submm fluxes of Type 1 and Type 2 AGN as well as AGN separated on the basis of X-ray hardness ratios and optical-to-infrared colours. However, at log10(LX) >44.4, dependences in average submm flux on the optical-to-infrared colours become more pronounced. We argue that these high luminosity AGN represent a transition from a secular to a merger-driven evolutionary phase where the star formation rates and accretion luminosities are more tightly coupled. Stacked AGN 850um fluxes are compared to the stacked fluxes of a mass-matched sample of K-band selected non-AGN galaxies. We find that at 10.5<log10(M*/M0)<11.5, the non-AGN 850um fluxes are 1.5-2x higher than in Type 2 AGN of equivalent mass. We suggest these differences are due to the presence of massive dusty, red starburst galaxies in the K-band selected non-AGN sample, which are not present in optically selected catalogues covering a smaller area.

preprint2015arXiv

Early Science with the Large Millimeter Telescope: Dust constraints in a z~9.6 galaxy

Recent observations with the GISMO (Goddard-IRAM Superconducting 2 Millimeter Observer) 2 mm camera revealed a detection 8 arcsec away from the lensed galaxy MACS1149-JD1 at z=9.6. Within the 17.5 arcsec FWHM GISMO beam, this detection is consistent with the position of the high-redshift galaxy and therefore, if confirmed, this object could be claimed to be the youngest galaxy producing significant quantities of dust. We present higher resolution (8.5 arcsec) observations of this system taken with the AzTEC 1.1 mm camera mounted on the Large Millimeter Telescope Alfonso Serrano. Dust continuum emission at the position of MACS1149-JD1 is not detected with an r.m.s. of 0.17 mJy/beam. However, we find a detection ~ 11 arcsec away from MACS1149-JD1, still within the GISMO beam which is consistent with an association to the GISMO source. Combining the AzTEC and GISMO photometry, together with Herschel ancillary data, we derive a z_phot= 0.7-1.6 for the dusty galaxy. We conclude therefore that the GISMO and AzTEC detections are not associated with MACS1149-JD1. From the non-detection of MACS1149-JD1 we derive the following (3 σ) upper limits corrected for gravitational lensing magnification and for cosmic microwave background effects: dust mass < 1.6 x 10^7 M_sun, IR luminosity < 8 x 10^10 L_sun, star formation rate < 14 M_sun/yr, and UV attenuation < 2.7 mag. These limits are comparable to those derived for other high-redshift galaxies from deep Atacama Large Millimeter/submillimeter Array (ALMA) observations.

preprint2015arXiv

HerMES: Current Cosmic Infrared Background Estimates Can be Explained by Known Galaxies and their Faint Companions at z < 4

We report contributions to cosmic infrared background (CIB) intensities originating from known galaxies and their faint companions at submillimeter wavelengths. Using the publicly-available UltraVISTA catalog, and maps at 250, 350, and 500 μm from the \emph{Herschel} Multi-tiered Extragalactic Survey (HerMES), we perform a novel measurement that exploits the fact that uncatalogued sources may bias stacked flux densities --- particularly if the resolution of the image is poor --- and intentionally smooth the images before stacking and summing intensities. By smoothing the maps we are capturing the contribution of faint (undetected in K_S ~ 23.4) sources that are physically associated, or correlated, with the detected sources. We find that the cumulative CIB increases with increased smoothing, reaching 9.82 +- 0.78, 5.77 +- 0.43, and 2.32 +- 0.19$\, \rm nW m^{-2} sr^{-1}$ at 250, 350, and 500 μm at 300 arcsec FWHM. This corresponds to a fraction of the fiducial CIB of 0.94 +- 0.23, 1.07 +- 0.31, and 0.97 +- 0.26 at 250, 350, and 500 μm, where the uncertainties are dominated by those of the absolute CIB. We then propose, with a simple model combining parametric descriptions for stacked flux densities and stellar mass functions, that emission from galaxies with log(M/Msun) > 8.5 can account for the most of the measured total intensities, and argue against contributions from extended, diffuse emission. Finally, we discuss prospects for future survey instruments to improve the estimates of the absolute CIB levels, and observe any potentially remaining emission at z > 4.

preprint2015arXiv

Identification of z~>2 Herschel 500 micron sources using color-deconfusion

We present a new method to search for candidate z~>2 Herschel 500μm sources in the GOODS-North field, using a S500μm/S24μm "color deconfusion" technique. Potential high-z sources are selected against low-redshift ones from their large 500μm to 24μm flux density ratios. By effectively reducing the contribution from low-redshift populations to the observed 500μm emission, we are able to identify counterparts to high-z 500μm sources whose 24μm fluxes are relatively faint. The recovery of known z~4 starbursts confirms the efficiency of this approach in selecting high-z Herschel sources. The resulting sample consists of 34 dusty star-forming galaxies at z~>2. The inferred infrared luminosities are in the range 1.5x10^12-1.8x10^13 Lsun, corresponding to dust-obscured star formation rates (SFRs) of ~260-3100 Msun/yr for a Salpeter IMF. Comparison with previous SCUBA 850μm-selected galaxy samples shows that our method is more efficient at selecting high-z dusty galaxies with a median redshift of z=3.07+/-0.83 and 10 of the sources at z~>4. We find that at a fixed luminosity, the dust temperature is ~5K cooler than that expected from the Td-LIR relation at z<1, though different temperature selection effects should be taken into account. The radio-detected subsample (excluding three strong AGN) follows the far-infrared/radio correlation at lower redshifts, and no evolution with redshift is observed out to z~5, suggesting that the far-infrared emission is star formation dominated. The contribution of the high-z Herschel 500μm sources to the cosmic SFR density is comparable to that of SMG populations at z~2.5 and at least 40% of the extinction-corrected UV samples at z~4 (abridged).

preprint2015arXiv

S-CANDELS: The Spitzer-Cosmic Assembly Near-Infrared Deep Extragalactic Survey. Survey Design, Photometry, and Deep IRAC Source Counts

The Spitzer-Cosmic Assembly Deep Near-Infrared Extragalactic Legacy Survey (S-CANDELS; PI G. Fazio) is a Cycle 8 Exploration Program designed to detect galaxies at very high redshifts (z > 5). To mitigate the effects of cosmic variance and also to take advantage of deep coextensive coverage in multiple bands by the Hubble Space Telescope Multi-Cycle Treasury Program CANDELS, S-CANDELS was carried out within five widely separated extragalactic fields: the UKIDSS Ultra-Deep Survey, the Extended Chandra Deep Field South, COSMOS, the HST Deep Field North, and the Extended Groth Strip. S-CANDELS builds upon the existing coverage of these fields from the Spitzer Extended Deep Survey (SEDS) by increasing the integration time from 12 hours to a total of 50 hours but within a smaller area, 0.16 square degrees. The additional depth significantly increases the survey completeness at faint magnitudes. This paper describes the S-CANDELS survey design, processing, and publicly-available data products. We present IRAC dual-band 3.6+4.5 micron catalogs reaching to a depth of 26.5 AB mag. Deep IRAC counts for the roughly 135,000 galaxies detected by S-CANDELS are consistent with models based on known galaxy populations. The increase in depth beyond earlier Spitzer/IRAC surveys does not reveal a significant additional contribution from discrete sources to the diffuse Cosmic Infrared Background (CIB). Thus it remains true that only roughly half of the estimated CIB flux from COBE/DIRBE is resolved.

preprint2015arXiv

T-PHOT: A new code for PSF-matched, prior-based, multiwavelength extragalactic deconfusion photometry

We present T-PHOT, a publicly available software aimed at extracting accurate photometry from low-resolution images of deep extragalactic fields, where the blending of sources can be a serious problem for the accurate and unbiased measurement of fluxes and colours. T-PHOT has been developed within the ASTRODEEP project and it can be considered as the next generation to TFIT, providing significant improvements above it and other similar codes. T-PHOT gathers data from a high-resolution image of a region of the sky, and uses it to obtain priors for the photometric analysis of a lower resolution image of the same field. It can handle different types of datasets as input priors: i) a list of objects that will be used to obtain cutouts from the real high-resolution image; ii) a set of analytical models; iii) a list of unresolved, point-like sources, useful e.g. for far-infrared wavelength domains. We show that T-PHOT yields accurate estimations of fluxes within the intrinsic uncertainties of the method, when systematic errors are taken into account (which can be done thanks to a flagging code given in the output). T-PHOT is many times faster than similar codes like TFIT and CONVPHOT (up to hundreds, depending on the problem and the method adopted), whilst at the same time being more robust and more versatile. This makes it an optimal choice for the analysis of large datasets. In addition we show how the use of different settings and methods significantly enhances the performance. Given its versatility and robustness, T-PHOT can be considered the preferred choice for combined photometric analysis of current and forthcoming extragalactic optical to far-infrared imaging surveys. [abridged]

preprint2015arXiv

The bulge-disk decomposition of AGN host galaxies

We present the results from a study of the morphologies of moderate luminosity X-ray selected AGN host galaxies in comparison to a carefully mass-matched control sample at 0.5 < z < 3 in the CANDELS GOODS-S field. We apply a multi-wavelength morphological decomposition analysis to these two samples and report on the differences between the morphologies as fitted from single Sersic and multiple Sersic models, and models which include an additional nuclear point-source component. Thus, we are able to compare the widely adopted single Sersic fits from previous studies to the results from a full morphological decomposition, and address the issue of how biased the inferred properties of AGN hosts are by a potential nuclear contribution from the AGN itself. We find that the AGN hosts are mixed systems which have higher bulge fractions than the control sample in our highest redshift bins at the >99.7% confidence level, according to all model fits even those which adopt a point-source component. This serves to alleviate concerns that previous, purely single Sersic, analyses of AGN hosts could have been spuriously biased towards higher bulge fractions. This dataset allows us to further probe the physical nature of these point-source components; we find no strong correlation between the point-source component and AGN activity, and that these point-source components are best modelled physically by nuclear starbursts. Our analysis of the bulge and disk fractions of these AGN hosts in comparison to a mass-matched control sample reveals a similar morphological evolutionary track for both the active and non-active populations, providing further evidence in favour of a model where AGN activity is triggered by secular processes.

preprint2015arXiv

The galaxy luminosity function at z ~ 6 and evidence for rapid evolution in the bright end from z ~ 7 to 5

We present the results of a search for bright (-22.7 < M_UV < -20.5) Lyman-break galaxies at z ~ 6 within a total of 1.65 square degrees of imaging in the UltraVISTA/COSMOS and UKIDSS UDS/SXDS fields. The deep near-infrared imaging available in the two independent fields, in addition to deep optical (including z'-band) data, enables the sample of z ~ 6 star-forming galaxies to be securely detected long-ward of the break (in contrast to several previous studies). We show that the expected contamination rate of our initial sample by cool galactic brown dwarfs is < 3 per cent and demonstrate that they can be effectively removed by fitting brown dwarf spectral templates to the photometry. At z ~ 6 the galaxy surface density in the UltraVISTA field exceeds that in the UDS by a factor of ~ 1.8, indicating strong cosmic variance even between degree-scale fields at z > 5. We calculate the bright end of the rest-frame Ultra-Violet (UV) luminosity function (LF) at z ~ 6. The galaxy number counts are a factor of ~1.7 lower than predicted by the recent LF determination by Bouwens et al.. In comparison to other smaller area studies, we find an evolution in the characteristic magnitude between z ~ 5 and z ~ 7 of dM* ~ 0.4 mag, and show that a double power-law or a Schechter function can equally well describe the LF at z = 6. Furthermore, the bright-end of the LF appears to steepen from z ~ 7 to z ~ 5, which could indicate the onset of mass quenching or the rise of dust obscuration, a conclusion supported by comparing the observed LFs to a range of theoretical model predictions.

preprint2014arXiv

A mass threshold in the number density of passive galaxies at z$\sim$2

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

preprint2014arXiv

A reassessment of the redshift distribution and physical properties of luminous (sub-)millimetre galaxies

Motivated by the current controversy over the redshift distribution and physical properties of luminous (sub-)mm sources, we have undertaken a new study of the brightest sample of unlensed (sub-)mm sources with pre-ALMA interferometric follow-up in the COSMOS field. Exploiting the very latest multi-frequency supporting data, we find that this sample displays a redshift distribution indistinguishable from that of the lensed sources uncovered with the South Pole Telescope (SPT), with z_median ~ 3.5. We also find that, over the redshift range z ~ 2-6 the median stellar mass of the most luminous (sub-)mm sources is M_* ~ 3 * 10^11 M_odot, yielding a typical specific star-formation rate sSFR ~ 3 Gyr^-1. Consistent with recent ALMA and SMA studies, we confirm that source blending is not a serious issue in the study of luminous (sub-)mm sources uncovered by ground-based, single-dish surveys; only ~ 10-15% of bright (S_850 ~ 5 - 10 mJy) (sub-)mm sources arise from significant (i.e. >20%) blends, and so our conclusions are largely unaffected by whether we adopt the original single-dish mm/sub-mm flux-densities/positions, or the interferometric data. Our results suggest that apparent disagreements over the redshift distribution of (sub-)mm sources are a result of "down-sizing" in dust-enshrouded star-formation, consistent with existing knowledge of the star-formation histories of massive galaxies. They also indicate that extreme star-forming galaxies at high redshift are, on average, subject to the same star-formation rate-limiting processes as less luminous objects, and lie on the "main sequence" of star-forming galaxies at z>3.

preprint2014arXiv

A submillimeter galaxy illuminating its circumgalactic medium: Ly-alpha scattering in a cold, clumpy outflow

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

preprint2014arXiv

Constraints on the star-formation rate of z~3 LBGs with measured metallicity in the CANDELS GOODS-South field

We analyse 14 LBGs at z~2.8-3.8 constituting the only sample where both a spectroscopic measurement of their metallicity and deep IR observations (CANDELS+HUGS survey) are available. Fixing the metallicity of population synthesis models to the observed values, we determine best-fit physical parameters under different assumptions about the star-formation history and also consider the effect of nebular emission. For comparison we determine the UV slope of the objects, and use it to estimate their SFR_UV99 by correcting the UV luminosity following Meurer et al. (1999). A comparison between SFR obtained through SED-fitting (SFR_fit) and the SFR_UV99 shows that the latter are underestimated by a factor 2-10, regardless of the assumed SFH. Other SFR indicators (radio, far-IR, X-ray, recombination lines) coherently indicate SFRs a factor of 2-4 larger than SFR_UV99 and in closer agreement with SFR_fit. This discrepancy is due to the solar metallicity implied by the usual beta-A1600 conversion factor. We propose a refined relation, appropriate for sub-solar metallicity LBGs: A1600 = 5.32+1.99beta. This relation reconciles the dust-corrected UV with the SED-fitting and the other SFR indicators. We show that the fact that z~3 galaxies have sub-solar metallicity implies an upward revision by a factor of ~1.5-2 of the global SFRD, depending on the assumptions about the age of the stellar populations. We find very young best-fit ages (10-500 Myrs) for all our objects. From a careful examination of the uncertainties in the fit and the amplitude of the Balmer break we conclude that there is little evidence of the presence of old stellar population in at least half of the LBGs in our sample, suggesting that these objects are probably caught during a huge star-formation burst, rather than being the result of a smooth evolution.

preprint2014arXiv

Imaging the environment of a z = 6.3 submillimeter galaxy with SCUBA-2

We describe a search for submillimeter emission in the vicinity of one of the most distant, luminous galaxies known, HerMES FLS3 at z=6.34, exploiting it as a signpost to a potentially biased region of the early Universe, as might be expected in hierarchical structure formation models. Imaging to the confusion limit with the innovative, wide-field submillimeter bolometer camera, SCUBA-2, we are sensitive to colder and/or less luminous galaxies in the surroundings of HFLS3. We use the Millennium Simulation to illustrate that HFLS3 may be expected to have companions if it is as massive as claimed, but find no significant evidence from the surface density of SCUBA-2 galaxies in its vicinity, or their colors, that HFLS3 marks an over-density of dusty, star-forming galaxies. We cannot rule out the presence of dusty neighbours with confidence, but deeper 450-um imaging has the potential to more tightly constrain the redshifts of nearby galaxies, at least one of which likely lies at z>~5. If associations with HFLS3 can be ruled out, this could be taken as evidence that HFLS3 is less biased than a simple extrapolation of the Millennium Simulation may imply. This could suggest either that it represents a rare short-lived, but highly luminous, phase in the evolution of an otherwise typical galaxy, or that this system has suffered amplification due to a foreground gravitational lens and so is not as intrinsically luminous as claimed.

preprint2014arXiv

The bright end of the galaxy luminosity function at z ~ 7: before the onset of mass quenching?

We present the results of a new search for bright star-forming galaxies at z ~ 7 within the UltraVISTA DR2 and UKIDSS UDS DR10 data, which together provide 1.65 sq deg of near-infrared imaging with overlapping optical and Spitzer data. Using a full photo-z analysis to identify high-z galaxies and reject contaminants, we have selected a sample of 34 luminous (-22.7 < M_UV < -21.2) galaxies with 6.5 < z < 7.5. Crucially, the deeper imaging provided by UltraVISTA DR2 confirms all of the robust objects previously uncovered by Bowler et al. (2012), validating our selection technique. Our sample includes the most massive galaxies known at z ~ 7, with M_* ~ 10^{10} M_sun, and the majority are resolved, consistent with larger sizes (r_{1/2} ~ 1 - 1.5 kpc) than displayed by less massive galaxies. From our final sample, we determine the form of the bright end of the rest-frame UV galaxy luminosity function (LF) at z ~ 7, providing strong evidence that the bright end of the z = 7 LF does not decline as steeply as predicted by the Schechter function fitted to fainter data. We consider carefully, and exclude the possibility that this is due to either gravitational lensing, or significant contamination of our galaxy sample by AGN. Rather, our results favour a double power-law form for the galaxy LF at high z, or, more interestingly, a LF which simply follows the form of the dark-matter halo mass function at bright magnitudes. This suggests that the physical mechanism which inhibits star-formation activity in massive galaxies (i.e. AGN feedback or some other form of `mass quenching') has yet to impact on the observable galaxy LF at z ~ 7, a conclusion supported by the estimated masses of our brightest galaxies which have only just reached a mass comparable to the critical `quenching mass' of M_* = 10 ^{10.2} M_sun derived from studies of the mass function of star-forming galaxies at lower z.

preprint2014arXiv

The Bulge-Disk Decomposed Evolution of Massive Galaxies at 1<z<3 in CANDELS

We present the results of a new and improved study of the morphological and spectral evolution of massive galaxies over the redshift range 1<z<3. Our analysis is based on a bulge-disk decomposition of 396 galaxies with Mstar>10^11 Msolar from the CANDELS WFC3/IR imaging within the COSMOS and UKIDSS UDS survey fields. We find that, by modelling the H(160) image of each galaxy with a combination of a de Vaucouleurs bulge (Sersic index n=4) and an exponential disk (n=1), we can then lock all derived morphological parameters for the bulge and disk components, and successfully reproduce the shorter-wavelength J(125), i(814), v(606) HST images simply by floating the magnitudes of the two components. This then yields sub-divided 4-band HST photometry for the bulge and disk components which, with no additional priors, is well described by spectrophotometric models of galaxy evolution. Armed with this information we are able to properly determine the masses and star-formation rates for the bulge and disk components, and find that: i) from z=3 to z=1 the galaxies move from disk-dominated to increasingly bulge-dominated, but very few galaxies are pure bulges/ellipticals by z=1; ii) while most passive galaxies are bulge-dominated, and most star-forming galaxies disk-dominated, 18+/-5% of passive galaxies are disk-dominated, and 11+/-3% of star-forming galaxies are bulge-dominated, a result which needs to be explained by any model purporting to connect star-formation quenching with morphological transformations; iii) there exists a small but significant population of pure passive disks, which are generally flatter than their star-forming counterparts (whose axial ratio distribution peaks at b/a~0.7); iv) flatter/larger disks re-emerge at the highest star-formation rates, consistent with recent studies of sub-mm galaxies, and with the concept of a maximum surface-density for star-formation activity.

preprint2014arXiv

The colour distribution of galaxies at redshift five

We present the results of a study investigating the rest-frame ultra-violet (UV) spectral slopes of redshift z~5 Lyman-break galaxies (LBGs). By combining deep Hubble Space Telescope imaging of the CANDELS and HUDF fields with ground-based imaging from the UKIDSS Ultra Deep Survey (UDS), we have produced a large sample of z~5 LBGs spanning an unprecedented factor of >100 in UV luminosity. Based on this sample we find a clear colour-magnitude relation (CMR) at z~5, such that the rest-frame UV slopes (beta) of brighter galaxies are notably redder than their fainter counterparts. We determine that the z~5 CMR is well described by a linear relationship of the form: d beta = (-0.12 +/- 0.02) d Muv, with no clear evidence for a change in CMR slope at faint magnitudes (i.e. Muv > -18.9). Using the results of detailed simulations we are able, for the first time, to infer the intrinsic (i.e. free from noise) variation of galaxy colours around the CMR at z~5. We find significant (12 sigma) evidence for intrinsic colour variation in the sample as a whole. Our results also demonstrate that the width of the intrinsic UV slope distribution of z~5 galaxies increases from Delta(beta)=0.1 at Muv=-18 to Delta(beta)=0.4 at Muv=-21. We suggest that the increasing width of the intrinsic galaxy colour distribution and the CMR itself are both plausibly explained by a luminosity independent lower limit of beta=-2.1, combined with an increase in the fraction of red galaxies in brighter UV-luminosity bins.

preprint2014arXiv

The Decomposed Bulge and Disk Size-Mass Relations of Massive Galaxies at 1<z<3 in CANDELS

We have constructed a mass-selected sample of Mstar>10^11Msolar galaxies at 1<z<3 in the CANDELS UDS and COSMOS fields and have decomposed these systems into their separate bulge and disk components according to their H(160)-band morphologies. By extending this analysis to multiple bands we have been able to conduct individual bulge and disk component SED fitting which has provided us with stellar-mass and star-formation rate estimates for the separate bulge and disk components. These have been combined with size measurements to explore the evolution of these massive high-redshift galaxies. By utilising the new decomposed stellar-mass estimates, we confirm that the bulge components display a stronger size evolution than the disks. This can be seen from both the fraction of bulge components which lie below the local relation and the median sizes of the bulge components, where the bulges are a median factor of 2.93+/-0.32 times smaller than similarly massive local galaxies at 1<z<2 and 3.41+/-0.58 smaller at 2<z<3; for the disks the corresponding factors are 1.65+/-0.14 and 1.99+/-0.25. Moreover, by splitting our sample into the passive and star-forming bulge and disk sub-populations and examining their sizes as a fraction of their present-day counter-parts, we find that the star-forming and passive bulges are equally compact, star-forming disks are larger, while the passive disks have intermediate sizes. This trend is not evident when classifying galaxy morphology on the basis of single-Sersic fits and adopting the overall star-formation rates. Finally, by evolving the star-formation histories of the passive disks back to the redshifts when the passive disks were last active, we show that the passive and star-forming disks have consistent sizes at the relevant epoch. These trends need to be reproduced by any mechanisms which attempt to explain the morphological evolution of galaxies.

preprint2014arXiv

The galaxy stellar mass function at 3.5<z<7.5 in the CANDELS/UDS, GOODS-South, and HUDF fields

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

preprint2014arXiv

The Hawk-I UDS and GOODS Survey (HUGS): Survey design and deep K-band number counts

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

preprint2014arXiv

The Mass-Metallicity-SFR Relation at z >~ 2 with 3D-HST

We present new accurate measurements of the physical properties of a statistically significant sample of 103 galaxies at z~2 using near-infrared spectroscopy taken as part of the 3D-HST survey. We derive redshifts, metallicities and star formation rates (SFRs) from the [OII], [OIII] and Hbeta nebular emission lines and exploit the multi-wavelength photometry available in CANDELS to measure stellar masses. We find the mass-metallicity relation (MZR) derived from our data to have the same trend as previous determinations in the range 0<z<3, with lower mass galaxies having lower metallicities. However we find an offset in the relation compared to the previous determination of the z~2 MZR by Erb et al. 2006b, who measure metallicities using the [NII]/Halpha ratio, with metallicities lower at a given mass. Incorporating our SFR information we find that our galaxies are offset from the Fundamental Metallicity Relation (FMR) by ~0.3 dex. We investigate the photoionization conditions and find that our galaxies are consistent with the elevated ionization parameter previously reported in high-redshift galaxies. Using the BPT diagram we argue that, if this is the case, metallicity indicators based on [NII] and Halpha may not be consistent with the ones obtained via oxygen lines and Hbeta. Using a recent determination of the theoretical evolution of the star forming sequence in the BPT diagram we convert our measured [OIII]/Hbeta line ratios to [NII]/Halpha ratios. From the [NII]/Halpha ratio we infer systematically higher metallicities in better agreement with the FMR. Our results thus suggest the evolution of the FMR previously reported at z~2-3 may be an artifact of the differential evolution in metallicity indicators, and caution against using locally calibrated metallicity relations at high redshift which do not account for evolution in the physical conditions of star-forming regions.

preprint2014arXiv

The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the rest-frame far-infrared emission of sub-millimeter galaxies

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

preprint2014arXiv

The SCUBA-2 Cosmology Legacy Survey: the submillimetre properties of Lyman break galaxies at z=3-5

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

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

A near-infrared morphological comparison of high-redshift submm and radio galaxies: massive star-forming discs vs relaxed spheroids

We present deep, high-quality K-band images of complete subsamples of powerful radio and sub-mm galaxies at z=2. The data were obtained in the best available seeing at UKIRT and Gemini North, with integration times scaled to ensure that comparable rest-frame surface brightness levels are reached for all galaxies. We fit two-dimensional axi-symmetric galaxy models to determine galaxy morphologies at rest-frame optical wavelengths > 4000A, varying luminosity, axial ratio, half-light radius, and Sersic index. We find that, while some images show evidence of galaxy interactions, >95% of the rest-frame optical light in all galaxies is well-described by these simple models. We also find a clear difference in morphology between these two classes of galaxy; fits to the individual images and image stacks reveal that the radio galaxies are moderately large (<r{1/2}>=8.4+-1.1kpc; median r{1/2}=7.8), de Vaucouleurs spheroids (<n> = 4.07+-0.27; median n=3.87), while the sub-mm galaxies appear to be moderately compact (<r{1/2}>=3.4+-0.3kpc; median r{1/2}=3.1kpc) exponential discs (<n>=1.44+-0.16; median n=1.08). We show that the z=2 radio galaxies display a well-defined Kormendy relation but that, while larger than other recently-studied high-z massive galaxy populations, they are still ~1.5 times smaller than their local counterparts. The scalelengths of the starlight in the sub-mm galaxies are comparable to those reported for the molecular gas. Their sizes are also similar to those of comparably massive quiescent galaxies at z>1.5. In terms of stellar mass surface density, the majority of the radio galaxies lie within the locus defined by local ellipticals. In contrast, while best modelled as discs, most of the sub-mm galaxies have higher stellar mass densities than local galaxies, and appear destined to evolve into present-day massive ellipticals.

preprint2013arXiv

A new multi-field determination of the galaxy luminosity function at z=7-9 incorporating the 2012 Hubble Ultra Deep Field imaging

We present a new determination of the UV galaxy luminosity function (LF) at redshift z ~ 7 and z ~ 8, and a first estimate at z ~ 9. An accurate determination of the form and evolution of the LF at high z is crucial for improving our knowledge of early galaxy evolution and cosmic reionization. Our analysis exploits fully the new, deepest WFC3/IR imaging from our HST UDF12 campaign, and includes a new, consistent analysis of all appropriate, shallower/wider-area HST data. Our new measurement of the evolving LF at z ~ 7-8 is based on a final catalogue of ~600 galaxies, and involves a step-wise maximum likelihood determination based on the redshift probability distribution for each object; this makes full use of the 11-band imaging now available in the HUDF, including the new UDF12 F140W data, and the deep Spitzer IRAC imaging. The final result is a determination of the z ~ 7 LF extending down to M_UV = -16.75, and the z ~ 8 LF down to M_UV = -17.00. Fitting a Schechter function, we find M* = -19.90 (+0.23/-0.28), log phi* = -2.96 (+0.18/-0.23), and a faint-end slope alpha=-1.90 (+0.14/-0.15) at z~7, and M* = -20.12 (+0.37/-0.48), log phi* = -3.35 (+0.28/-0.47), alpha=-2.02 (+0.22/-0.23) at z~8. These results strengthen suggestions that the evolution at z > 7 is more akin to `density evolution' than the apparent `luminosity evolution' seen at z ~ 5-7. We also provide the first meaningful information on the LF at z ~ 9, explore alternative extrapolations to higher z, and consider the implications for the evolution of UV luminosity density. Finally, we provide catalogues (including z_phot, M_UV and all photometry) for the 100 most robust z~6.5-11.9 galaxies in the HUDF used in this analysis. We discuss our results in the context of earlier work and the results of an independent analysis of the UDF12 data based on colour-colour selection (Schenker et al. 2013).

preprint2013arXiv

Chronos: A NIR spectroscopic galaxy survey. From the formation of galaxies to the peak of activity

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

preprint2013arXiv

High velocity outflows from young star-forming galaxies in the UKIDSS Ultra-Deep Survey

We investigate galactic-scale outflows in the redshift range 0.71 < z < 1.63, using 413 K-band selected galaxies observed in the spectroscopic follow-up of the UKIDSS Ultra-Deep Survey (UDSz). The galaxies have an average stellar mass of ~10^9.5 solar masses and span a wide range in rest-frame colours, representing typical star-forming galaxies at this epoch. We stack the spectra by various galaxy properties, including stellar mass, [OII] equivalent width, star-formation rate, specific star-formation rate and rest-frame spectral indices. We find that outflows are present in virtually all spectral stacks, with velocities ranging from 100-1000 km s^-1, indicating that large-scale outflowing winds are a common property at these redshifts. The highest velocity outflows (>500 km s^-1) are found in galaxies with the highest stellar masses and the youngest stellar populations. Our findings suggest that high velocity galactic outflows are mostly driven by star-forming processes rather than AGN, with implied mass outflow rates comparable to the rates of star formation. Such behaviour is consistent with models required to reproduce the high-redshift mass-metallicity relation.

preprint2013arXiv

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

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

preprint2013arXiv

On-sky characterisation of the VISTA NB118 narrow-band filters at 1.19 micron

Observations of the high redshift Universe through narrow-band filters have proven very successful in the last decade. The 4-meter VISTA telescope, equipped with the wide-field camera VIRCAM, offers a major step forward in wide-field near-infrared imaging, and in order to utilise VISTA's large field-of-view and sensitivity, the Dark Cosmology Centre provided a set of 16 narrow-band filters for VIRCAM. These NB118 filters are centered at a wavelength near 1.19 micron in a region with few airglow emission lines. The filters allow the detection of Halpha emitters at z = 0.8, Hbeta and [OIII] emitters at z ~ 1.4, [OII] emitters at z = 2.2, and Ly-alpha emitters at z = 8.8. Based on guaranteed time observations of the COSMOS field we here present a detailed description and characterization of the filters and their performance. In particular we provide sky-brightness levels and depths for each of the 16 detector/filter sets and find that some of the filters show signs of some red-leak. We identify a sample of 2 x 10^3 candidate emission-line objects in the data. Cross-correlating this sample with a large set of galaxies with known spectroscopic redshifts we determine the "in situ" passbands of the filters and find that they are shifted by about 3.5-4 nm (corresponding to 30% of the filter width) to the red compared to the expectation based on the laboratory measurements. Finally, we present an algorithm to mask out persistence in VIRCAM data. Scientific results extracted from the data will be presented separately.

preprint2013arXiv

The Morphologies of Massive Galaxies at 1 < z < 3 in the CANDELS-UDS Field

We have used high-resolution, HST WFC3/IR, near-infrared imaging to conduct a detailed bulge-disk decomposition of the morphologies of ~200 of the most massive (M_star > 10^11 M_solar) galaxies at 1<z<3 in the CANDELS-UDS field. We find that, while such massive galaxies at low redshift are generally bulge-dominated, at redshifts 1<z<2 they are predominantly mixed bulge+disk systems, and by z>2 they are mostly disk-dominated. Interestingly, we find that while most of the quiescent galaxies are bulge-dominated, a significant fraction (25-40%) of the most quiescent galaxies, have disk-dominated morphologies. Thus, our results suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies.

preprint2013arXiv

The properties of (sub)millimetre-selected galaxies as revealed by CANDELS HST WFC3/IR imaging in GOODS-South

We have exploited the HST CANDELS WFC3/IR imaging to study the properties of (sub-)mm galaxies in GOODS-South. After using the deep radio and Spitzer imaging to identify galaxy counterparts for the (sub-)mm sources, we have used the new CANDELS data in two ways. First, we have derived improved photometric redshifts and stellar masses, confirming that the (sub-)mm galaxies are massive (<M*>=2.2x10^11 M_solar) galaxies at z=1-3. Second, we have exploited the depth and resolution of the WFC3/IR imaging to determine the sizes and morphologies of the galaxies at rest-frame optical wavelengths, fitting two-dimensional axi-symmetric Sersic models. Crucially, the WFC3/IR H-band imaging enables modelling of the mass-dominant galaxy, rather than the blue high-surface brightness features which often dominate optical (rest-frame UV) images of (sub-)mm galaxies, and can confuse visual morphological classification. As a result of this analysis we find that >95% of the rest-frame optical light in almost all of the (sub-)mm galaxies is well-described by either a single exponential disk, or a multiple-component system in which the dominant constituent is disk-like. We demonstrate that this conclusion is consistent with the results of high-quality ground-based K-band imaging, and explain why. The massive disk galaxies which host luminous (sub-)mm emission are reasonably extended (r_e=4 kpc), consistent with the sizes of other massive star-forming disks at z~2. In many cases we find evidence of blue clumps within the sources, with the mass-dominant disk becoming more significant at longer wavelengths. Finally, only a minority of the sources show evidence for a major galaxy-galaxy interaction. Taken together, these results support the view that most (sub-)mm galaxies at z~2 are simply the most extreme examples of normal star-forming galaxies at that era.

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 SCUBA-2 Cosmology Legacy Survey: Ultraluminous star-forming galaxies in a z=1.6 cluster

We analyse new SCUBA-2 submillimeter and archival SPIRE far-infrared imaging of a z=1.62 cluster, Cl0218.3-0510, which lies in the UKIDSS/UDS field of the SCUBA-2 Cosmology Legacy Survey. Combining these tracers of obscured star formation activity with the extensive photometric and spectroscopic information available for this field, we identify 31 far-infrared/submillimeter-detected probable cluster members with bolometric luminosities >1e12 Lo and show that by virtue of their dust content and activity, these represent some of the reddest and brightest galaxies in this structure. We exploit Cycle-1 ALMA submillimeter continuum imaging which covers one of these sources to confirm the identification of a SCUBA-2-detected ultraluminous star-forming galaxy in this structure. Integrating the total star-formation activity in the central region of the structure, we estimate that it is an order of magnitude higher (in a mass-normalised sense) than clusters at z~0.5-1. However, we also find that the most active cluster members do not reside in the densest regions of the structure, which instead host a population of passive and massive, red galaxies. We suggest that while the passive and active populations have comparable near-infrared luminosities at z=1.6, M(H)~-23, the subsequent stronger fading of the more active galaxies means that they will evolve into passive systems at the present-day which are less luminous than the descendants of those galaxies which were already passive at z~1.6 (M(H)~-20.5 and M(H)~-21.5 respectively at z~0). We conclude that the massive galaxy population in the dense cores of present-day clusters were already in place at z=1.6 and that in Cl0218.3-0510 we are seeing continuing infall of less extreme, but still ultraluminous, star-forming galaxies onto a pre-existing structure.

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.

preprint2013arXiv

X-ray selected AGN Hosts are Similar to Inactive Galaxies out to z=3: Results from CANDELS/CDF-S

We use multi-band spatially resolved photometry from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) in the 4 MSec Chandra Deep Field-South (CDF-S) to explore the nuclear and extended colors, color gradients and stellar populations of X-ray selected AGN host galaxies out to z=3. Based on a study of their central light, we develop X-ray based criteria to exclude objects with strong AGN contamination. We use stellar masses from the FIREWORKS database to understand and account for stellar mass selection effects, and carefully study, for the first time, the resolved host galaxy properties of AGNs at z~2 in their rest-frame optical light without substantial nuclear contamination. AGN hosts span a sizable range of stellar masses, colors and color gradients at these redshifts. Their colors, color gradients and stellar population properties are very similar to inactive galaxies of the same stellar mass. At z~1, we find a slightly narrower range in host colors compared to inactive galaxies, as well as hints of more recent star-formation. These differences are weaker or non-existent among AGN hosts at z~2. We discuss the importance of AGN driven feedback in the quenching of galaxies at z>~1 and speculate on possible evolution in the relationship between black hole accretion and the host galaxy towards high redshifts.

preprint2012arXiv

A critical analysis of the UV-continuum slopes of high-redshift galaxies; no evidence (yet) for extreme stellar populations at z > 6

It has recently been reported that the faintest galaxies at z~6-7 display extremely blue UV continuum slopes, with a UV power-law index beta ~ -3. Such slopes are bluer than previously reported for any other galaxy population, and imply extinction-free, young, very low-metallicity stellar populations with a high ionizing photon escape fraction. Here we undertake a critical study of the evidence for such extreme values of beta, combining three new WFC3/IR-selected samples of galaxies spanning ~2 decades in UV luminosity over the redshift range z~4.5-8. We explore the impact of inclusion/exclusion of less robust high-z candidates, and use the varying depths of the samples to explore the effects of noise and selection bias. Simple data-consistency arguments suggest that artificially blue average values of beta can result when the analysis is extended into the deepest ~ 0.5-mag bin of these WFC3/IR-selected samples, regardless of the actual luminosity or z range probed. By confining attention to robust, well-detected high-z galaxy candidates, we find that the average value of beta is consistent with -2.05 +/- 0.10 for z=5-7, and -22 < M_UV < -18. We create and analyse a set of simulations which demonstrate that a bias towards artifically low/blue average values of beta is indeed expected when the UV slope analysis is extended towards the source detection threshold, and conclude that there is as yet no clear evidence for UV slopes significantly bluer than beta ~ -2, the typical value displayed by the bluest star-forming galaxies at more modest z. A robust measurement of beta for the faintest galaxies at z~7-8 remains a key observational goal, as it provides a fundamental test for high escape fractions from a potentially abundant source of reionizing photons. This goal is achievable with HST, but requires still deeper WFC3/IR imaging in the HUDF.

preprint2012arXiv

A Herschel view of the far-infrared properties of submillimetre galaxies

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

preprint2012arXiv

A remarkably high fraction of strong Ly_alpha emitters amongst luminous redshift 6.0<z<6.5 Lyman break galaxies in the UKIDSS Ultra-Deep Survey

We present spectroscopic confirmation of ten highly luminous (L >= 2L*) Lyman alpha emitters in the redshift range 6.01<z<6.49 (nine galaxies and one AGN), initially drawn from a sample of fourteen z_phot >= 6 Lyman break galaxies (LBGs) selected from an area of 0.25 square degrees within the UKIDSS Ultra-deep Survey (UDS). Overall, our high rate of spectroscopic confirmation (>= 71%) and low rate of contamination provides a strong vindication of the photometric redshift analysis used to define the original sample. By considering star-formation rate estimates based on the Ly_alpha and UV continuum luminosity we conclude that our sample is consistent with a Ly_alpha escape fraction of ~25%. Moreover, after careful consideration of the potential uncertainties and biases, we find that 40%-50% of our sample of L >= 2L* galaxies at 6.0<z<6.5 display strong Ly_alpha emission (rest-frame equivalent width >= 25 Angs), a fraction which is a factor of ~2 higher than previously reported for L <= L* galaxies at z~6. Our results suggest that, as the epoch of reionization is approached, it is plausible that the Ly_alpha emitter fraction amongst luminous (L >=2 L*) LBGs shows a similarly sharp increase to that observed in their lower-luminosity (L <= L*) counterparts.

preprint2012arXiv

AzTEC half square degree survey of the SHADES fields - II. Identifications, redshifts, and evidence for large-scale structure

The AzTEC 1.1 mm survey of the SCUBA HAlf Degree Extragalactic Survey (SHADES) fields is the largest (0.7 deg2) blank-field millimetre-wavelength survey undertaken to date at a resolution of ~18" and a depth of ~1 mJy. We have used the deep optical-to-radio multi-wavelength data in the SHADES Lockman Hole East and SXDF/UDS fields to obtain galaxy identifications for ~64% (~80% with tentative identifications) of the 148 AzTEC-SHADES 1.1 mm sources, exploiting deep radio and 24 um data complemented by methods based on 8 um flux-density and red optical-infrared (i-K) colour. This unusually high identification rate can be attributed to the relatively bright millimetre-wavelength flux-density threshold, combined with the relatively deep supporting multi-frequency data. We have further exploited the optical-mid-infrared-radio data to derive a ~60% (~75% with tentative identifications) complete redshift distribution for the AzTEC-SHADES sources, yielding a median redshift of z~2.2, with a high-redshift tail extending to at least z~4. Despite the larger area probed by the AzTEC survey relative to the original SCUBA SHADES imaging, the redshift distribution of the AzTEC sources is consistent with that displayed by the SCUBA sources, and reinforces tentative evidence that the redshift distribution of mm/sub-mm sources in the Lockman Hole field is significantly different from that found in the SXDF/UDS field. Comparison with simulated surveys of similar scale extracted from semi-analytic models based on the Millennium simulation indicates that this is as expected if the mm/sub-mm sources are massive (M>10^11 Mo) star-forming galaxies tracing large-scale structures over scales of 10-20 Mpc. This confirms the importance of surveys covering several square degrees (as now underway with SCUBA2) to obtain representative samples of bright (sub)mm-selected galaxies.

preprint2012arXiv

CANDELS Observations of the Structural Properties and Evolution of Galaxies in a Cluster at z=1.62

We discuss the structural and morphological properties of galaxies in a z=1.62 proto-cluster using near-IR imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: galaxies with colors of quiescent stellar populations generally have morphologies consistent with spheroids, and galaxies with colors consistent with ongoing star formation have disk-like and irregular morphologies. The size distribution of the quiescent cluster galaxies shows a deficit of compact (< 1kpc), massive galaxies compared to CANDELS field galaxies at z=1.6. As a result the cluster quiescent galaxies have larger average effective sizes compared to field galaxies at fixed mass at greater than 90% significance. Combined with data from the literature, the size evolution of quiescent cluster galaxies is relatively slow from z~1.6 to the present, growing as (1+z)^(-0.6+/-0.1). If this result is generalizable, then it implies that physical processes associated with the denser cluster region seems to have caused accelerated size growth in quiescent galaxies prior to z=1.6 and slower subsequent growth at z<1.6 compared to galaxies in the lower density field. The quiescent cluster galaxies at z=1.6 have higher ellipticities compared to lower redshift samples at fixed mass, and their surface-brightness profiles suggest that they contain extended stellar disks. We argue the cluster galaxies require dissipationless (i.e., gas-poor or "dry") mergers to reorganize the disk material and to match the relations for ellipticity, stellar mass, size, and color of early-type galaxies in z<1 clusters.

preprint2012arXiv

Discovery of bright z ~ 7 galaxies in the UltraVISTA survey

We have exploited the new, deep, near-infrared UltraVISTA imaging of the COSMOS field, in tandem with deep optical and mid-infrared imaging, to conduct a new search for luminous galaxies at redshifts z ~ 7. The unique multi-wavelength dataset provided by VISTA, CFHT, Subaru, HST and Spitzer over a common area of 1 deg^2 has allowed us to select galaxy candidates at z > 6.5 by searching first for Y+J-detected (< 25 AB mag) objects which are undetected in the CFHT+HST optical data. This sample was then refined using a photometric redshift fitting code, enabling the rejection of lower-redshift galaxy contaminants and cool galactic M,L,T dwarf stars.The final result of this process is a small sample of (at most) ten credible galaxy candidates at z > 6.5 which we present in this paper. The first four of these appear to be robust galaxies at z > 6.5, and fitting to their stacked SED yields z = 6.98+-0.05 with a stellar mass M* = 5x10^9 Msun, and rest-frame UV spectral slope beta = -2.0+-0.2. The next three are also good candidates for z > 6.5 galaxies, but the possibility that they are low-redshift galaxies or dwarf stars cannot be excluded. Our final subset of three additional candidates is afflicted not only by potential dwarf-star contamination, but also contains objects likely to lie at redshifts just below z = 6.5. We show that the three even-brighter z > 7 galaxy candidates reported in the COSMOS field by Capak et al. (2011) in fact all lie at z ~ 1.5-3.5. Consequently the new z ~ 7 galaxies reported here are the first credible z ~ 7 Lyman-break galaxies discovered in the COSMOS field and, as the most UV-luminous discovered to date at these redshifts, are prime targets for deep follow-up spectroscopy. We explore their physical properties, and briefly consider the implications of their inferred number density for the form of the galaxy luminosity function at z = 7.

preprint2012arXiv

FMOS near-IR spectroscopy of Herschel selected galaxies: star formation rates, metallicity and dust attenuation at z~1

We investigate the properties (e.g. star formation rate, dust attentuation, stellar mass and metallicity) of a sample of infrared luminous galaxies at z \sim 1 via near-IR spectroscopy with Subaru-FMOS. Our sample consists of Herschel SPIRE and Spitzer MIPS selected sources in the COSMOS field with photometric redshifts in the range 0.7 < z-phot < 1.8, which have been targeted in 2 pointings (0.5 sq. deg.) with FMOS. We find a modest success rate for emission line detections, with candidate Hα emission lines detected for 57 of 168 SPIRE sources (34 per cent). By stacking the near-IR spectra we directly measure the mean Balmer decrement for the Hα and Hβ lines, finding a value of <E(B-V)> = 0.51\pm0.27 for <LIR> = 10^12 Lsol sources at <z> = 1.36. By comparing star formation rates estimated from the IR and from the dust uncorrected Hα line we find a strong relationship between dust attenuation and star formation rate. This relation is broadly consistent with that previously seen in star-forming galaxies at z ~ 0.1. Finally, we investigate the metallicity via the N2 ratio, finding that z ~ 1 IR-selected sources are indistinguishable from the local mass-metallicity relation. We also find a strong correlation between dust attentuation and metallicity, with the most metal-rich IR-sources experiencing the largest levels of dust attenuation.

preprint2012arXiv

Herschel-ATLAS: Multi-wavelength SEDs and physical properties of 250 micron-selected galaxies at z < 0.5

We present a pan-chromatic analysis of an unprecedented sample of 1402 250 micron-selected galaxies at z < 0.5 (mean z = 0.24) from the Herschel-ATLAS survey. We complement our Herschel 100-500 micron data with UV-K-band photometry from the Galaxy And Mass Assembly (GAMA) survey and apply the MAGPHYS energy-balance technique to produce pan-chromatic SEDs for a representative sample of 250 micron selected galaxies spanning the most recent 5 Gyr of cosmic history. We derive estimates of physical parameters, including star formation rates, stellar masses, dust masses and infrared luminosities. The typical H-ATLAS galaxy at z < 0.5 has a far-infrared luminosity in the range 10^10 - 10^12 Lsolar (SFR: 1-50 Msolar/yr) thus is broadly representative of normal star forming galaxies over this redshift range. We show that 250 micron-selected galaxies contain a larger mass of dust at a given infra-red luminosity or star formation rate than previous samples selected at 60 micron from IRAS. We derive typical SEDs for H-ATLAS galaxies, and show that the emergent SED shape is most sensitive to specific star formation rate. The optical-UV SEDs also become more reddened due to dust at higher redshifts. Our template SEDs are significantly cooler than existing infra-red templates. They may therefore be most appropriate for inferring total IR luminosities from moderate redshift submillimetre selected samples and for inclusion in models of the lower redshift submillimetre galaxy populations.

preprint2012arXiv

Herschel-ATLAS/GAMA: spatial clustering of low-redshift sub-mm galaxies

We have measured the clustering properties of low-redshift (z < 0.3) sub-mm galaxies detected at 250 micron in the Herschel-ATLAS Science Demonstration Phase (SDP) field. We selected a sample for which we have high-quality spectroscopic redshifts, obtained from reliably matching the 250-micron sources to a complete (for r < 19.4) sample of galaxies from the GAMA database. Both the angular and spatial clustering strength are measured for all z < 0.3 sources as well as for five redshift slices with thickness delta z=0.05 in the range 0.05 < z < 0.3. Our measured spatial clustering length r_0 is comparable to that of optically-selected, moderately star-forming (blue) galaxies: we find values around 5 Mpc. One of the redshift bins contains an interesting structure, at z = 0.164.

preprint2012arXiv

The ages, masses and star-formation rates of spectroscopically confirmed z~6 galaxies in CANDELS

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5<z<6.5, all with WFC3/IR and IRAC imaging from the HST/CANDELS and Spitzer/SEDS surveys. Based on fitting the observed photometry with SED templates covering a wide range of different star-formation histories, and a self consistent treatment of Lyman-alpha emission, we find that the derived stellar masses lie within the range 10^9 Msun < M*< 10^10 Msun and are robust to within a factor of two. In contrast, we confirm previous reports that the ages of the stellar populations are poorly constrained. Although the best-fitting models for three objects have ages >= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged).

preprint2012arXiv

The Morphologies of Massive Galaxies at 1<z<3 in the CANDELS-UDS Field: Compact Bulges, and the Rise and Fall of Massive Disks

We have used deep, HST, near-IR imaging to study the morphological properties of the most massive galaxies at high z, modelling the WFC3/IR H-band images of the ~200 galaxies in the CANDELS-UDS field with 1 < z_phot < 3, and stellar masses M_star > 10^11 M_sun. We have used both single-Sersic and bulge+disk models, have investigated the errors/biases introduced by uncertainties in the background and the PSF, and have obtained formally-acceptable model fits to >90% of the galaxies. Our results indicate that these massive galaxies at 1 < z < 3 lie both on and below the local size-mass relation, with a median R_e~2.6 kpc, a factor of ~2.3 smaller than comparably-massive local galaxies. Moreover, we find that bulge-dominated objects in particular show evidence for a growing bimodality in the size-mass relation with increasing z, and by z > 2 the compact bulges display effective radii a factor ~4 smaller than local ellipticals of comparable mass. These trends appear to extend to the bulge components of disk-dominated galaxies, and vice versa. We also find that, while such massive galaxies at low z are bulge-dominated, at 1 < z < 2 they are predominantly mixed bulge+disk systems, and by z > 2 they are mostly disk-dominated. The majority of the disk-dominated galaxies are actively forming stars, but this is also true for many of the bulge-dominated systems. Interestingly, however, while most of the quiescent galaxies are bulge-dominated, we find that a significant fraction (25-40%) of the most quiescent galaxies have disk-dominated morphologies. Thus, while our results show that the massive galaxy population is undergoing dramatic changes at this crucial epoch, they also suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies into present-day giant ellipticals.

preprint2012arXiv

The optically unbiased GRB host (TOUGH) survey. VI. Radio observations at z<1 and consistency with typical star-forming galaxies

The objective of this paper is to determine the level of obscured star formation activity and dust attenuation in a sample of gamma-ray burst (GRB) hosts; and to test the hypothesis that GRB hosts have properties consistent with those of the general star-forming galaxy populations. We present a radio continuum survey of all z<1 GRB hosts in The Optically Unbiased GRB Host (TOUGH) sample supplemented with radio data for all (mostly pre-Swift) GRB-SN hosts discovered before October 2006. We present new radio data for 22 objects and have obtained a detection for three of them (GRB 980425, 021211, 031203; none in the TOUGH sample), increasing the number of radio-detected GRB hosts from two to five. The star formation rate (SFR) for the GRB 021211 host of ~825 Mo yr^-1, the highest ever reported for a GRB host, places it in the category of ultraluminous infrared galaxies. We found that at least 63% of GRB hosts have SFR < 100 Mo yr^-1 and at most 8% can have SFR > 500 Mo yr^-1. For the undetected hosts the mean radio flux (<35 uJy 3sigma) corresponds to an average SFR < 15 Mo yr^-1. Moreover, ~88% of the z<1 GRB hosts have ultraviolet dust attenuation A_UV < 6.7 mag (visual attenuation A_V < 3 mag). Hence we did not find evidence for large dust obscuration in a majority of GRB hosts. Finally, we found that the distributions of SFRs and A_UV of GRB hosts are consistent with those of Lyman break galaxies, Halpha emitters at similar redshifts and of galaxies from cosmological simulations. The similarity of the GRB population with other star-forming galaxies is consistent with the hypothesis that GRBs, a least at z<1, trace a large fraction of all star formation, and are therefore less biased indicators than once thought.

preprint2012arXiv

The SCUBA-2 Cosmology Legacy Survey: blank-field number counts of 450um-selected galaxies and their contribution to the cosmic infrared background

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

preprint2012arXiv

The size-luminosity relation at z=7 in CANDELS and its implication on reionization

The exploration of the relation between galaxy sizes and other physical parameters has provided important clues for understanding galaxy formation. We use the CANDELS Deep+Wide surveys in the GOODS-South, UDS and EGS fields, complemented by data from the HUDF09 program, to address the relation between size and luminosity at z\sim7. We select 153 z-band drop-out galaxies in six different fields characterized by a wide combination of depth and areal coverage, ideally suited to sample without biases the observed size-magnitude plane. Detailed simulations allow us to derive the completeness as a function of size and magnitude and to quantify measurements errors/biases. We find that the half light radius distribution function of z\sim7 galaxies fainter than J=26.6 is peaked at <0.1 arcsec (or equivalently 0.5 kpc proper), while at brighter magnitudes high-z galaxies are typically larger than ~0.15 arcsec. We also find a well defined size-luminosity relation, Rh\simL^1/2. We compute the Luminosity Function in the HUDF and P12HUDF fields, finding large spatial variation on the number density of faint galaxies. Adopting the size distribution and the size-luminosity relation found for faint galaxies at z=7, we derive a mean slope of -1.7\pm0.1 for the luminosity function of LBGs at this redshift. Using this LF, we find that the amount of ionizing photons cannot keep the Universe re-ionized if the IGM is clumpy (C_HII>3) and the Lyman continuum escape fraction of high-z LBGs is relatively low (f_esc<0.3). With future CANDELS data, we can put severe limits to the role of galaxies in the reionization of the Universe.

preprint2012arXiv

The sizes, masses and specific star-formation rates of massive galaxies at 1.3<z<1.5: strong evidence in favour of evolution via minor mergers

We report the results of a comprehensive study of the relationship between galaxy size, stellar mass and specific star-formation rate (sSFR) at redshifts 1.3<z<1.5. Based on a mass complete (M_star >= 6x10^10 Msun), spectroscopic sample from the UKIDSS Ultra-deep Survey (UDS), with accurate stellar-mass measurements derived from spectro photometric fitting, we find that at z~1.4 the location of massive galaxies on the size-mass plane is determined primarily by their sSFR. At this epoch we find that massive galaxies which are passive (sSFR <= 0.1 Gyr^-1) follow a tight size-mass relation, with half-light radii a factor f=2.4+/-0.2 smaller than their local counterparts. Moreover, amongst the passive sub-sample we find no evidence that the off-set from the local size-mass relation is a function of stellar population age. Based on a sub-sample with dynamical mass estimates we also derive an independent estimate of f=2.3+/-0.3 for the typical growth in half-light radius between z~1.4 and the present day. Focusing on the passive sub-sample, we conclude that to produce the necessary evolution predominantly via major mergers would require an unfeasible number of merger events and over populate the high-mass end of the local stellar mass function. In contrast, we find that a scenario in which mass accretion is dominated by minor mergers can produce the necessary evolution, whereby an increase in stellar mass by a factor of ~2, accompanied by an increase in size by a factor of ~3.5, is sufficient to reconcile the size-mass relation at z~1.4 with that observed locally. Finally, we note that a significant fraction (44+/-12%) of the passive galaxies in our sample have a disk-like morphology, providing additional evidence that separate physical processes are responsible for the quenching of star-formation and the morphological transformation of massive galaxies (abridged).

preprint2012arXiv

The source counts of submillimetre galaxies detected at 1.1 mm

The source counts of galaxies discovered at sub-millimetre and millimetre wavelengths provide important information on the evolution of infrared-bright galaxies. We combine the data from six blank-field surveys carried out at 1.1 mm with AzTEC, totalling 1.6 square degrees in area with root-mean-square depths ranging from 0.4 to 1.7 mJy, and derive the strongest constraints to date on the 1.1 mm source counts at flux densities S(1100) = 1-12 mJy. Using additional data from the AzTEC Cluster Environment Survey to extend the counts to S(1100) ~ 20 mJy, we see tentative evidence for an enhancement relative to the exponential drop in the counts at S(1100) ~ 13 mJy and a smooth connection to the bright source counts at >20 mJy measured by the South Pole Telescope; this excess may be due to strong lensing effects. We compare these counts to predictions from several semi-analytical and phenomenological models and find that for most the agreement is quite good at flux densities > 4 mJy; however, we find significant discrepancies (>3sigma) between the models and the observed 1.1 mm counts at lower flux densities, and none of them are consistent with the observed turnover in the Euclidean-normalised counts at S(1100) < 2 mJy. Our new results therefore may require modifications to existing evolutionary models for low luminosity galaxies. Alternatively, the discrepancy between the measured counts at the faint end and predictions from phenomenological models could arise from limited knowledge of the spectral energy distributions of faint galaxies in the local Universe.

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 unbiased measurement of UV spectral slopes in low luminosity galaxies at z=7

The Ultraviolet (UV) continuum slope beta, typically observed at z=7 in Hubble Space Telescope (HST) WFC3/IR bands via the J-H colour, is a useful indicator of the age, metallicity, and dust content of high-redshift stellar populations. Recent studies have shown that the redward evolution of beta with cosmic time from redshift 7 to 4 can be largely explained by a build up of dust. However, initial claims that faint z=7 galaxies in the Hubble Ultra Deep Field WFC3/IR imaging (HUDF09) were blue enough to require stellar populations of zero reddening, low metallicity and young ages, hitherto unseen in star-forming galaxies, have since been refuted and revised. Here we revisit the question of how best to measure the UV slope of z=7 galaxies through source recovery simulations, within the context of present and future ultra-deep imaging from HST. We consider how source detection, selection and colour measurement have each biased the measurement of beta in previous studies. After finding a robust method for measuring beta in the simulations (via a power law fit to all the available photometry), we remeasure the UV slopes of a sample of previously published low luminosity z=7 galaxy candidates. The mean UV slope of faint galaxies in this sample appears consistent with an intrinsic distribution of normal star-forming galaxies with beta=-2, although properly decoding the underlying distribution will require further imaging from the ongoing HUDF12 programme. We therefore go on to consider strategies for obtaining better constraints on the underlying distribution of UV slopes at z=7 from these new data, which will benefit particularly from the addition of imaging in a second J-band filter: F140W. We find that a precise and unbiased measurement of beta should then be possible.

preprint2012arXiv

The VISTA Deep Extragalactic Observations (VIDEO) Survey

In this paper we describe the first data release of the the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is a ~12degree^2 survey in the near-infrared Z,Y,J,H and K_s bands, specifically designed to enable the evolution of galaxies and large structures to be traced as a function of both epoch and environment from the present day out to z=4, and active galactic nuclei (AGN) and the most massive galaxies up to and into the epoch of reionization. With its depth and area, VIDEO will be able to fully explore the period in the Universe where AGN and starburst activity were at their peak and the first galaxy clusters were beginning to virialize. VIDEO therefore offers a unique data set with which to investigate the interplay between AGN, starbursts and environment, and the role of feedback at a time when it was potentially most crucial. We provide data over the VIDEO-XMM3 tile, which also covers the Canada-France-Hawaii-Telescope Legacy Survey Deep-1 field (CFHTLS-D1). The released VIDEO data reach a 5-sigma AB-magnitude depth of Z=25.7, Y=24.5, J=24.4, H=24.1 and K_s=23.8 in 2 arcsec diameter apertures (the full depth of Y=24.6 will be reached within the full integration time in future releases). The data are compared to previous surveys over this field and we find good astrometric agreement with the Two-Micron All Sky Survey, and source counts in agreement with the recently released UltraVISTA survey data. The addition of the VIDEO data to the CFHTLS-D1 optical data increases the accuracy of photometric redshifts and significantly reduces the fraction of catastrophic outliers over the redshift range 0<z<1 from 5.8 to 3.1 per cent in the absence of an i-band luminosity prior. (Truncated Abstract)

preprint2011arXiv

A robust sample of galaxies at redshifts 6.0<z<8.7: stellar populations, star-formation rates and stellar masses

We present the results of a photometric redshift analysis designed to identify z>6 galaxies from the near-IR HST imaging in three deep fields (HUDF, HUDF09-2 & ERS). By adopting a rigorous set of criteria for rejecting low-z interlopers, and by employing a deconfusion technique to allow the available IRAC imaging to be included in the candidate selection process, we have derived a robust sample of 70 Lyman-break galaxies (LBGs) spanning the redshift range 6.0<z<8.7. Based on our final sample we investigate the distribution of UV spectral slopes (beta), finding a variance-weighted mean value of <beta>=-2.05 +/- 0.09 which, contrary to some previous results, is not significantly bluer than displayed by lower-redshift starburst galaxies. We confirm the correlation between UV luminosity and stellar mass reported elsewhere, but based on fitting galaxy templates featuring a range of star-formation histories, metallicities and reddening we find that, at z>=6, the range in mass-to-light ratio (M*/L_UV) at a given UV luminosity could span a factor of ~50. Focusing on a sub-sample of twenty-one candidates with IRAC detections at 3.6-microns we find that L* LBGs at z~6.5 have a median stellar mass of M* = (2.1 +/- 1.1) x 10^9 Msun and a median specific star-formation rate of 1.9 +/- 0.8 Gyr^-1. Using the same sub-sample we have investigated the influence of nebular continuum and line emission, finding that for the majority of candidates (16 out of 21) the best-fitting stellar-mass estimates are reduced by less than a factor of 2.5. Finally, a detailed comparison of our final sample with the results of previous studies suggests that, at faint magnitudes, several high-redshift galaxy samples in the literature are significantly contaminated by low-redshift interlopers (abridged).

preprint2011arXiv

AGN environments at z<1.5 in the UKIDSS Ultra-Deep Survey

We investigate the environments of both X-ray and radio-loud AGN within the UKIDSS Ultra-deep Survey (UDS) using deep infrared selection to sample the galaxy density field in the redshift range 1.0 <= z <= 1.5. Using angular cross-correlation techniques we find that both X-ray and radio-loud AGN preferentially reside in overdense environments. We also find that both types of AGN cluster more strongly with those galaxies classified as `passive' rather than those that are actively star-forming. We infer clustering scale lengths comparable to those of passive red galaxies, suggesting that typical AGN at these epochs reside in dark-matter halos of mass M >~ 10^13 M_sun. A closer look at the small-scale environments of the AGN reveals that the neighbouring galaxies of radio-loud AGN have U-B colours more skewed towards the `green-valley' and the red sequence, whereas the neighbours of X-ray AGN show no difference to the general galaxy population. This suggests that although both AGN types live in overdense environments, the radio-loud AGN may be preferentially located in more evolved cluster cores, in a similar environment to low-powered radio AGN in the local Universe.

preprint2011arXiv

Extreme Emission Line Galaxies in CANDELS: Broad-Band Selected, Star-Bursting Dwarf Galaxies at z>1

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z~1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. Supported by spectroscopic confirmation of strong [OIII] emission lines -- with rest-frame equivalent widths ~1000Å-- in the four candidates that have HST/WFC3 grism observations, we conclude that these objects are galaxies with 10^8 Msol in stellar mass, undergoing an enormous starburst phase with M_*/(dM_*/dt) of only ~15 Myr. These bursts may cause outflows that are strong enough to produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the co-moving number density (3.7x10^-4 Mpc^-3) can produce in ~4 Gyr much of the stellar mass density that is presently contained in 10^8-10^9 Msol dwarf galaxies. Therefore, our observations provide a strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, short-lived bursts, mostly at z>1.

preprint2011arXiv

Galaxy Environments in the UKIDSS Ultra Deep Survey (UDS)

We present a study of galaxy environments to z~2, based on a sample of over 33,000 K-band selected galaxies detected in the UKIDSS Ultra Deep Survey (UDS). The combination of infrared depth and area in the UDS allows us to extend previous studies of galaxy environment to z>1 without the strong biases associated with optical galaxy selection. We study the environments of galaxies divided by rest frame (U-B) colours, in addition to `passive' and `star-forming' subsets based on template fitting. We find that galaxy colour is strongly correlated with galaxy overdensity on small scales (<1Mpc diameter), with red/passive galaxies residing in significantly denser environments than blue/star-forming galaxies to z~1.5. On smaller scales (<0.5Mpc diameter) we also find a relationship between galaxy luminosity and environment, with the most luminous blue galaxies at z~1 inhabiting environments comparable to red, passive systems at the same redshift. Monte Carlo simulations demonstrate that these conclusions are robust to the uncertainties introduced by photometric redshift errors.

preprint2011arXiv

Herschel-ATLAS: counterparts from the UV--NIR in the science demonstration phase catalogue

We present a technique to identify optical counterparts of 250 um-selected sources from the Herschel-ATLAS survey. Of the 6621 250 um > 32 mJy sources in our science demonstration catalogue we find that ~60 percent have counterparts brighter than r=22.4 mag in the Sloan Digital Sky Survey. Applying a likelihood ratio technique we are able to identify 2423 of the counterparts with a reliability R > 0.8. This is approximately 37 percent of the full 250 micron catalogue. We have estimated photometric redshifts for each of these 2423 reliable counterparts, while 1099 also have spectroscopic redshifts collated from several different sources, including the GAMA survey. We estimate the completeness of identifying counterparts as a function of redshift, and present evidence that 250 um-selected Herschel-ATLAS galaxies have a bimodal redshift distribution. Those with reliable optical identifications have a redshift distribution peaking at z ~ 0.25 +/- 0.05, while sub-mm colours suggest that a significant fraction with no counterpart above the r-band limit have z > 1. We also suggest a method for selecting populations of strongly-lensed high redshift galaxies. Our identifications are matched to UV--NIR photometry from the GAMA survey, and these data are available as part of the Herschel-ATLAS public data release.

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

Keck Spectroscopy of Faint 3<z<8 Lyman Break Galaxies:- Evidence for a Declining Fraction of Emission Line Sources In the Redshift Range 6<z<8

Using deep Keck spectroscopy of Lyman break galaxies selected from infrared imaging data taken with WFC3/IR onboard the Hubble Space Telescope, we present new evidence for a reversal in the redshift-dependent fraction of star forming galaxies with detectable Lyman alpha emission in the redshift range 6.3 < z < 8.8. Our earlier surveys with the DEIMOS spectrograph demonstrated a significant increase with redshift in the fraction of line emitting galaxies over the interval 4 < z < 6, particularly for intrinsically faint systems which dominate the luminosity density. Using the longer wavelength sensitivities of LRIS and NIRSPEC, we have targeted 19 Lyman break galaxies selected using recent WFC3/IR data whose photometric redshifts are in the range 6.3 < z < 8.8 and which span a wide range of intrinsic luminosities. Our spectroscopic exposures typically reach a 5-sigma sensitivity of < 50 A for the rest-frame equivalent width (EW) of Lyman alpha emission. Despite the high fraction of emitters seen only a few hundred million years later, we find only 2 convincing and 1 possible line emitter in our more distant sample. Combining with published data on a further 7 sources obtained using FORS2 on the ESO VLT, and assuming continuity in the trends found at lower redshift, we discuss the significance of this apparent reversal in the redshift-dependent Lyman alpha fraction in the context of our range in continuum luminosity. Assuming all the targeted sources are at their photometric redshift and our assumptions about the Lyman alpha EW distribution are correct, we would expect to find so few emitters in less than 1% of the realizations drawn from our lower redshift samples. Our new results provide further support for the suggestion that, at the redshifts now being probed spectroscopically, we are entering the era where the intergalactic medium is partially neutral.

preprint2011arXiv

Luminous Starbursts in the Redshift Desert at z~1-2: Star Formation Rates, Masses & Evidence for Outflows

We present a spectroscopic catalogue of 40 luminous starburst galaxies at z=0.7--1.7 (median z=1.3). 19 of these are submillimetre galaxies (SMGs) and 21 are submillimetre-faint radio galaxies (SFRGs). This sample helps to fill in the redshift desert at z=1.2--1.7 in previous studies as well as probing a lower luminosity population of galaxies. Radio fluxes are used to determine star-formation rates for our sample which range from around 50 to 500 M$_\odot$ yr$^{-1}$ and are generally lower than those in z$\sim$2 SMGs. We identify nebular [OII] 3727 emission in the rest-UV spectra and use the linewidths to show that SMGs and SFRGs in our sample have larger linewidths and therefore dynamical masses than optically selected star-forming galaxies at similar redshifts. The linewidths are indistinguishable from those measured in the z$\sim$2 SMG populations suggesting little evolution in the dynamical masses of the galaxies between redshift 1--2. [NeV] and [NeIII] emission lines are identified in a subset of the spectra indicating the presence of an active galactic nucleus (AGN). In addition, a host of interstellar absorption lines corresponding to transitions of MgII and FeII ions are also detected. These features show up prominently in composite spectra and we use these composites to demonstrate that the absorption lines are present at an average blueshift of $-240\pm$50 kms$^{-1}$ relative to the systemic velocities of the galaxies derived from [OII]. This indicates the presence of large-scale outflowing interstellar gas in these systems (Abridged)

preprint2011arXiv

The BLAST 250-micron selected galaxy population in GOODS-South

We identify and investigate the nature of the 20 brightest 250-micron sources detected by the BLAST within the central 150 sq. arcmin of the GOODS-South field. Aided by the available deep VLA radio imaging, reaching S_1.4 = 30 micro-Jy, we have identified radio counterparts for 17/20 of the 250-micron sources. The resulting enhanced positional accuracy of ~1 arcsec has then allowed us to exploit the deep multi-frequency imaging of GOODS-South to establish secure galaxy counterparts for the 17 radio-identified sources, and plausible galaxy candidates for the 3 radio-unidentified sources. Confusion is a serious issue for this deep BLAST 250-micron survey, due to the large size of the beam. Nevertheless, we argue that our chosen counterparts are significant, and often dominant contributors to the measured BLAST flux densities. For all of these 20 galaxies we have been able to determine spectroscopic (8) or photometric (12) redshifts. The result is the first `complete' redshift distribution for a deep 250-micron selected galaxy sample. This reveals that 250-micron surveys reaching detection limits of ~30 mJy contain not only low-redshift spirals/LIRGs, but also the extreme z~2 dust-enshrouded starburst galaxies previously discovered at sub-millimetre wavelengths. Based on their IRAC colours, we find that virtually all of the BLAST galaxy identifications appear better described as analogues of the M82 starburst galaxy, or Sc star-forming discs rather than highly obscured ULIRGs. Inspection of the LABOCA 870-micron imaging of the GOODS-South field yields detections of 7/11 of the z>1 BLAST sources, and reveals 250/870 flux-density ratios consistent with a standard 40K modified black-body fit with a dust emissivity index beta=1.5.

preprint2011arXiv

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

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 μm by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth τ = (ν/ν0)^β where β = +1.8 and ν0 = c/100 μm. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 μm selected samples.

preprint2011arXiv

The LABOCA Survey of the Extended Chandra Deep Field South: Clustering of submillimetre galaxies

We present a measurement of the spatial clustering of submillimetre galaxies (SMGs) at z = 1-3. Using data from the 870 micron LESS survey, we employ a novel technique to measure the cross-correlation between SMGs and galaxies, accounting for the full probability distributions for photometric redshifts of the galaxies. From the observed projected two-point cross-correlation function we derive the linear bias and characteristic dark matter (DM) halo masses for the SMGs. We detect clustering in the cross-correlation between SMGs and galaxies at the > 4 sigma level. For the SMG autocorrelation we obtain r_0 = 7.7 (+1.8,-2.3) h^-1 Mpc, and derive a corresponding DM halo mass of log(M_halo [h^-1 M_sun]) = 12.8 (+0.3,-0.5). Based on the evolution of DM haloes derived from simulations, we show that that the z = 0 descendants of SMGs are typically massive (~2-3 L*) elliptical galaxies residing in moderate- to high-mass groups (log(M_halo [h^-1 M_sun]) = 13.3 (+0.3,-0.5). From the observed clustering we estimate an SMG lifetime of ~100 Myr, consistent with lifetimes derived from gas consumption times and star-formation timescales, although with considerable uncertainties. The clustering of SMGs at z ~ 2 is consistent with measurements for optically-selected quasi-stellar objects (QSOs), supporting evolutionary scenarios linking starbursts and QSOs. Given that SMGs reside in haloes of characteristic mass ~ 6 x 10^12 h^-1 M_sun, we demonstrate that the redshift distribution of SMGs can be described remarkably well by the combination of two effects: the cosmological growth of structure and the evolution of the molecular gas fraction in galaxies. We conclude that the powerful starbursts in SMGs likely represent a short-lived but universal phase in massive galaxy evolution, associated with the transition between cold gas-rich, star-forming galaxies and passively evolving systems. [Abridged]

preprint2011arXiv

The luminosity-dependent high-redshift turnover in the steep spectrum radio luminosity function: clear evidence for downsizing in the radio-AGN population

This paper presents a new grid-based method for investigating the evolution of the steep-spectrum radio luminosity function, with the aim of quantifying the high-redshift cut-off suggested by previous work. To achieve this, the Combined EIS-NVSS Survey of Radio Sources (CENSORS) has been developed; this is a 1.4 GHz radio survey, containing 135 sources complete to a flux density of 7.2 mJy, selected from the NRAO VLA Sky Survey (NVSS) over 6 deg^2 of the ESO Imaging Survey (EIS) Patch D. The sample is currently 73% spectroscopically complete, with the remaining redshifts estimated via the K-z or I-z magnitude-redshift relation. CENSORS is combined with additional radio data from the Parkes All-Sky, Parkes Selected Regions, Hercules and VLA COSMOS samples to provide comprehensive coverage of the radio power vs. redshift plane. The redshift distributions of these samples, together with radio source count determinations, and measurements of the local luminosity function, provide the input to the fitting process. The modelling reveals clear declines, at > 3sigma significance, in comoving density at z > 0.7 for lower luminosity sources (log P = 25-26); these turnovers are still present at log P > 27, but move to z > 3, suggesting a luminosity-dependent evolution of the redshift turnover, similar to the `cosmic downsizing' seen for other AGN populations. These results are shown to be robust to the estimated redshift errors and to increases in the spectral index for the highest redshift sources. Analytic fits to the best-fitting steep spectrum grid are provided so that the results presented here can be easily accessed by the reader, as well as allowing plausible extrapolations outside of the regions covered by the input datasets

preprint2011arXiv

The zCOSMOS redshift survey : Influence of luminosity, mass and environment on the galaxy merger rate

The contribution of major mergers to galaxy mass assembly along cosmic time is an important ingredient to the galaxy evolution scenario. We aim to measure the evolution of the merger rate for both luminosity/mass selected galaxy samples and investigate its dependence with the local environment. We use a sample of 10644 spectroscopically observed galaxies from the zCOSMOS redshift survey to identify pairs of galaxies destined to merge, using only pairs for which the velocity difference and projected separation of both components with a confirmed spectroscopic redshift indicate a high probability of merging. We have identified 263 spectroscopically confirmed pairs with r_p^{max} = 100 h^{-1} kpc. We find that the density of mergers depends on luminosity/mass, being higher for fainter/less massive galaxies, while the number of mergers a galaxy will experience does not depends significantly on its intrinsic luminosity but rather on its stellar mass. We find that the pair fraction and merger rate increase with local galaxy density, a property observed up to redshift z=1. We find that the dependence of the merger rate on the luminosity or mass of galaxies is already present up to redshifts z=1, and that the evolution of the volumetric merger rate of bright (massive) galaxies is relatively flat with redshift with a mean value of 3*10^{-4} (8*10^{-5} respectively) mergers h^3 Mpc^{-3} Gyr^{-1}. The dependence of the merger rate with environment indicates that dense environments favors major merger events as can be expected from the hierarchical scenario. The environment therefore has a direct impact in shapping-up the mass function and its evolution therefore plays an important role on the mass growth of galaxies along cosmic time.

preprint2010arXiv

A Spitzer-Selected Galaxy Cluster at z=1.62

We report the discovery of a galaxy cluster at z=1.62 located in the Spitzer Wide-Area Infrared Extragalactic survey XMM-LSS field. This structure was selected solely as an overdensity of galaxies with red Spitzer/IRAC colors, satisfying [3.6]-[4.5] > -0.1 AB mag. Photometric redshifts derived from Subaru XMM Deep Survey (BViz-bands), UKIRT Infrared Deep Survey-Ultra-Deep Survey (UKIDSS-UDS, JK-bands), and from the Spitzer Public UDS survey (3.6-8.0 micron) show that this cluster corresponds to a surface density of galaxies at z ~ 1.6 that is more than 20 sigma above the mean at this redshift. We obtained optical spectroscopic observations of galaxies in the cluster region using IMACS on the Magellan telescope. We measured redshifts for seven galaxies in the range z=1.62-1.63 within 2.8 arcmin (<1.4 Mpc) of the astrometric center of the cluster. A posteriori analysis of the XMM data in this field reveal a weak (4 sigma) detection in the [0.5-2 keV] band compatible with the expected thermal emission from such a cluster. The color-magnitude diagram of the galaxies in this cluster shows a prominent red-sequence, dominated by a population of red galaxies with (z-J) > 1.7 mag. The photometric redshift probability distributions for the red galaxies are strongly peaked at z=1.62, coincident with the spectroscopically confirmed galaxies. The rest-frame (U-B) color and scatter of galaxies on the red-sequence are consistent with a mean luminosity-weighted age of 1.2 +/- 0.1 Gyr, yielding a formation redshift z_f = 2.35 +/- 0.10, and corresponding to the last significant star-formation period in these galaxies.

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

GBT Zpectrometer CO(1-0) Observations of the Strongly-Lensed Submillimeter Galaxies from the Herschel ATLAS

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) has uncovered a population of strongly-lensed submillimeter galaxies (SMGs). The Zpectrometer instrument on the Green Bank Telescope (GBT) was used to measure the redshifts and constrain the masses of the cold molecular gas reservoirs for two candidate high-redshift lensed sources. We derive CO(1-0) redshifts of z=3.042+/-0.001 and z=2.625+/-0.001, and measure molecular gas masses of (1--3)x10^{10}Msun, corrected for lens amplification and assuming a conversion factor of alhpa=0.8 Msun(K km/s pc^2)^{-1}. We find typical L(IR)/L'(CO) ratios of 120+/-40 and 140+/-50 Lsun (K km/s pc^2)^{-1}, which are consistent with those found for local ULIRGs and other high-redshift SMGs. From analysis of published data, we find no evidence for enhanced L(IR)/L'(CO(1-0)) ratios for the SMG population in comparison to local ULIRGs. The GBT results highlight the power of using the CO lines to derive blind redshifts, which is challenging for the SMG population at optical wavelengths given their high obscuration.

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: The far-infrared--radio correlation at z<0.5

We use data from the Herschel-ATLAS to investigate the evolution of the far-infrared--radio correlation over the redshift range 0<z<0.5. Using the total far-infrared luminosity of all >5sigma sources in the Herschel-ATLAS Science Demonstration Field and cross-matching these data with radio data from the Faint Images of the Radio Sky at Twenty-Centimetres (FIRST) survey and the NRAO VLA Northern Sky Survey (NVSS), we obtain 104 radio counterparts to the Herschel sources. With these data we find no evidence for evolution in the far-infrared--radio correlation over the redshift range 0<z<0.5, where the median value for the ratio between far-infrared and radio luminosity, $q_{\rm IR}$, over this range is $q_{\rm IR} = 2.40\pm 0.12$ (and a mean of $q_{\rm IR}=2.52 \pm 0.03$ accounting for the lower limits), consistent with both the local value determined from {\em IRAS} and values derived from surveys targeting the high-redshift Universe. By comparing the radio fluxes of our sample measured from both FIRST and NVSS we show that previous results suggesting an increase in the value of $q_{\rm IR}$ from high to low redshift may be the result of resolving out extended emission of the low-redshift sources with relatively high-resolution interferometric data, although AGN contamination could still play a significant role. We also find tentative evidence that the longer wavelength, cooler dust is heated by an evolved stellar population which does not trace the star-formation rate as closely as the shorter wavelength $\ltsim 250~μ$m emission or the radio emission, supporting suggestions based on detailed models of individual galaxies.

preprint2010arXiv

Mid-infrared spectroscopy of candidate AGN-dominated submillimeter galaxies

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

preprint2010arXiv

The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

preprint2010arXiv

The Herschel-ATLAS: Evolution of the 250um luminosity function out to z=0.5

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

preprint2010arXiv

The stellar mass function of the most massive galaxies at 3<=z<5 in the UKIDSS Ultra Deep Survey

We have analysed a sample of 1292 4.5 micron-selected galaxies at z>=3, over 0.6 square degrees of the UKIRT Infrared Deep Survey (UKIDSS) Ultra Deep Survey (UDS). Using photometry from the U band through 4.5 microns, we have obtained photometric redshifts and derived stellar masses for our sources. Only two of our galaxies potentially lie at z>5. We have studied the galaxy stellar mass function at 3<=z<5, based on the 1213 galaxies in our catalogue with [4.5]<= 24.0. We find that: i) the number density of M > 10^11 Msun galaxies increased by a factor > 10 between z=5 and 3, indicating that the assembly rate of these galaxies proceeded > 20 times faster at these redshifts than at 0<z<2; ii) the Schechter function slope alpha is significantly steeper than that displayed by the local stellar mass function, which is both a consequence of the steeper faint end and the absence of a pure exponential decline at the high-mass end; iii) the evolution of the comoving stellar mass density from z=0 to 5 can be modelled as log10 (rho_M) =-(0.05 +/- 0.09) z^2 - (0.22 -/+ 0.32) z + 8.69. At 3<=z<4, more than 30% of the M > 10^11 Msun galaxies would be missed by optical surveys with R<27 or z<26. Thus, our study demonstrates the importance of deep mid-IR surveys over large areas to perform a complete census of massive galaxies at high z and trace the early stages of massive galaxy assembly.

preprint2010arXiv

Unveiling the Nature of Submillimeter Galaxy SXDF850.6

We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emission that is found ~8" north of the coordinates determined from SCUBA. Despite the lack of detectable counterparts in deep UV/optical images, we find a source at the SMA position in near-infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry (u, B, V, R, i', z', J, H, K, 3.6 micron, 4.5 micron, 5.8 micron, and 8.0 micron) and to submm-radio photometry (850 micron, 880 micron, 1100 micron, and 21 cm) independently, and we find both are well described by starburst templates at a redshift of z ~= 2.2 (+/- 0.3). The best-fit parameters from the UV-optical-IR SED fit are a redshift of z = 1.87 (+0.15/-0.07), a stellar mass of M_star = 2.5 +2.2/-0.3 x 10^11 M_sun, an extinction of A_V = 3.0 (+0.3/-1.0) mag, and an age of 720 (+1880/-210) Myr. The submm-radio SED fit provides a consistent redshift of z ~ 1.8-2.5, an IR luminosity of L_IR = (7-26) x 10^12 L_sun, and a star formation rate of 1300-4500 M_sun/yr. These results suggest that SXDF850.6 is a mature system already having a massive amount of old stellar population constructed before its submm bright phase and is experiencing a dusty starburst, possibly induced by major mergers.

preprint2009arXiv

AzTEC Half Square Degree Survey of the SHADES Fields -- I. Maps, Catalogues, and Source Counts

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

preprint2009arXiv

Deep multi-frequency radio imaging in the Lockman Hole using the GMRT and VLA: II. The spectral index of submillimetre galaxies

We have employed the GMRT and the VLA to map the Lockman Hole. At 610 and 1,400 MHz, we reach noise levels of 15 and 6 uJy/beam, respectively, with well-matched resolutions (~5"). At this depth we obtained reliable detections for about half of the known submm galaxies (SMGs) in the field. For radio-identified SMGs, which are typically at z ~ 2, we measure a mean radio spectral index of alpha = -0.75 +/- 0.06 and standard deviation of 0.29, between rest-frame ~1.8 and ~4.2 GHz. The slope of their continuum emission is indistinguishable from that of local star-forming galaxies and suggests that extended optically-thin synchrotron emission dominates the radio output of SMGs. Cooling effects by synchrotron emission and Inverse Compton scattering off the CMB do not seem to affect their radio SEDs. For those SMGs judged by Spitzer mid-IR colours and spectroscopy to host obscured AGN, we find a clear deviation from the rest of the sample - they typically have steeper radio spectral indices, alpha ~< -1.0. These findings suggest these mid-IR-/AGN-selected SMGs may have an intrinsically different injection mechanism for relativistic particles, or they might reside in denser environments. This work provides a reliable spectral template for the estimation of far-IR/radio photometric redshifts, and will enable accurate statistical K-corrections for the large samples of SMGs expected with SCUBA-2 and Herschel.

preprint2009arXiv

Galaxies at z = 6 - 9 from the WFC3/IR imaging of the HUDF

We present the results of a systematic search for galaxies in the redshift range z = 6 - 9, within the new, deep, near-infrared imaging of the Hubble Ultra Deep Field provided by the Wide Field Camera 3 (WFC3) on HST. We have performed full SED fitting to the optical+infrared photometry of all high-redshift galaxy candidates detected at greater than 5-sigma in at least one of the WFC3/IR broad-band filters. After rejection of contaminants, the result is a sample of 49 galaxies with primary redshift solutions z > 5.9. Our sample, selected without recourse to specific colour cuts, re-selects all but the faintest one of the 16 z-drops selected by Oesch et al. (2009), recovers all 5 of the Y-drops reported by Bouwens et al. (2009), and adds a further 29 galaxy candidates, of which 12 lie beyond z = 6.3, and 4 lie beyond z = 7. We also present confidence intervals on our photometric redshift estimates, and caution that acceptable low-redshift (z < 2) solutions exist for 28 out of the 37 galaxies at z > 6.3, and for all 8 galaxy candidates at z > 7.5. Nevertheless, the very highest redshift candidates appear to be strongly clustered in the field. We derive new estimates of the ultraviolet galaxy luminosity function at z = 7 and z = 8. Where our results are most robust, at a characteristic luminosity M(1500) ~ -19.5 (AB), we find that the comoving number density of galaxies declines by a factor of ~ 2.5 between z = 6 and z = 7, and by a further factor of ~ 2 by z = 8. These results suggest that it is difficult for the observed population of high-redshift star-forming galaxies to achieve reionisation by z ~ 6 without a significant contribution from galaxies well below the detection limits, plus alterations in the escape fraction of ionising photons and/or continued vigorous star formation at z > 15.

preprint2009arXiv

The LABOCA Survey of the Extended Chandra Deep Field South

We present a sensitive 870 micron survey of the Extended Chandra Deep Field South (ECDFS) using LABOCA on the APEX telescope. The LABOCA ECDFS Submillimetre Survey (LESS) covers the full 30' x 30' field size of the ECDFS and has a uniform noise level of 1.2 mJy/beam. LESS is thus the largest contiguous deep submillimetre survey undertaken to date. The noise properties of our map show clear evidence that we are beginning to be affected by confusion noise. We present a catalog of 126 SMGs detected with a significance level above 3.7 sigma. The ECDFS exhibits a deficit of bright SMGs relative to previously studied blank fields but not of normal star-forming galaxies that dominate the extragalactic background light (EBL). This is in line with the underdensities observed for optically defined high redshift source populations in the ECDFS (BzKs, DRGs,optically bright AGN and massive K-band selected galaxies). The differential source counts in the full field are well described by a power law with a slope of alpha=-3.2, comparable to the results from other fields. We show that the shape of the source counts is not uniform across the field. The integrated 870 micron flux densities of our source-count models account for >65% of the estimated EBL from COBE measurements. We have investigated the clustering of SMGs in the ECDFS by means of a two-point correlation function and find evidence for strong clustering on angular scales <1'. Assuming a power law dependence for the correlation function and a typical redshift distribution for the SMGs we derive a spatial correlation length of r_0=13+/-6 h^-1 Mpc.

preprint2009arXiv

X-ray groups and clusters of galaxies in the Subaru-XMM Deep Field

We present the results of a search for galaxy clusters in Subaru-XMM Deep Field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 2x10^{-15} ergs cm^{-2} s^{-1} over one of the widest XMM-Newton contiguous raster surveys, covering an area of 1.3 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. The red sequence technique allows us to identify 57 cluster candidates. We report on the progress with the cluster spectroscopic follow-up and derive their properties based on the X-ray luminosity and cluster scaling relations. In addition, 3 sources are identified as X-ray counterparts of radio lobes, and in 3 further sources, X-ray counterpart of radio lobes provides a significant fraction of the total flux of the source. In the area covered by NIR data, our identification success rate achieves 86%. We detect a number of radio galaxies within our groups and for a luminosity-limited sample of radio galaxies we compute halo occupation statistics using a marked cluster mass function. We compare the cluster detection statistics in the SXDF with the predictions of concordance cosmology and current knowledge of the X-ray cluster properties, concluding that a reduction of concordance sigma_8 value by 5% is required in order to match the prediction of the model and the data. This conclusion still needs verification through the completion of cluster follow-up.

preprint2008arXiv

A new measurement of the evolving near-infrared galaxy luminosity function out to z~4: a continuing challenge to theoretical models of galaxy formation

We present the most accurate measurement to date of cosmological evolution of the near-infrared galaxy luminosity function, from the local Universe out to z~4. The analysis is based on a large and highly complete sample of galaxies selected from the first data release of the UKIDSS Ultra Deep Survey. Exploiting a master catalogue of K- and z-band selected galaxies over an area of 0.7 square degrees, we analyse a sample of ~50,000 galaxies, all with reliable photometry in 16-bands from the far-ultraviolet to the mid-infrared. The unique combination of large area and depth provided by the Ultra Deep Survey allows us to trace the evolution of the K-band luminosity function with unprecedented accuracy. In particular, via a maximum likelihood analysis we obtain a simple parameterization for the luminosity function and its cosmological evolution, including both luminosity and density evolution, which provides an excellent description of the data from z =0 up to z~4. We find differential evolution for galaxies dependent on galaxy luminosity, revealing once again the ``down-sizing behaviour'' of galaxy formation. Finally, we compare our results with the predictions of the latest theoretical models of galaxy formation, based both on semi-analytical prescriptions, and on full hydrodynamical simulations.

preprint2007arXiv

The UKIRT Infrared Deep Sky Survey First Data Release

The First Data Release (DR1) of the UKIRT Infrared Deep Sky Survey (UKIDSS) took place on 2006 July 21. UKIDSS is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. DR1 is the first large release of survey-quality data from UKIDSS and includes 320 sq degs of multicolour data to (Vega) K=18, complete (depending on the survey) in three to five bands from the set ZYJHK, together with 4 sq degs of deep JK data to an average depth K=21. In addition the release includes a similar quantity of data with incomplete filter coverage. In JHK, in regions of low extinction, the photometric uniformity of the calibration is better than 0.02 mag. in each band. The accuracy of the calibration in ZY remains to be quantified, and the same is true of JHK in regions of high extinction. The median image FWHM across the dataset is 0.82 arcsec. We describe changes since the Early Data Release in the implementation, pipeline and calibration, quality control, and archive procedures. We provide maps of the areas surveyed, and summarise the contents of each of the five surveys in terms of filters, areas, and depths. DR1 marks completion of 7 per cent of the UKIDSS 7-year goals.

preprint2006arXiv

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

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

preprint2005arXiv

A robust sample of submillimetre galaxies: constraints on the prevalence of dusty, high-redshift starbursts

The modest significance of most sources detected in current submm surveys can potentially compromise some analyses due to the inclusion of spurious sources in catalogues typically selected at >3.0-3.5sigma. Here, we develop and apply a dual-survey extraction technique to SCUBA and MAMBO images of the Lockman Hole. Cut above 5sigma, our catalogue of SMGs is more robust than previous samples, with a reduced likelihood of real, but faint SMGs (beneath and around the confusion limit) entering via superposition with noise. The effective flux limit of the survey is well matched to our deep 1.4-GHz image. The former is sensitive to luminous, dusty galaxies at extreme redshifts whilst the latter probes the z<3 regime. A high fraction of our robust SMGs (~80%) have radio counterparts which, given the ~10% contamination by spurious sources, suggests that very distant SMGs (z>>3) are unlikely to make up more than ~10% of the bright SMG population. This implies that almost all of the bright SMG population is amenable to study via the deepest current radio imaging. We use these radio counterparts to provide an empirical calibration of the positional uncertainty in SMG catalogues. We then go on to outline the acquisition of redshifts for radio-identified SMGs, from sample selection in the submm, to counterpart selection in the radio and optical/IR, to slit placement on spectrograph masks. We determine a median of z=2.05 +/- 0.41 from a sample of six secure redshifts for unambigious radio-identified submm sources and z=2.14 +/- 0.27 when we include submm sources with multiple radio counterparts and/or less reliable redshifts. These figures are consistent with previous estimates, suggesting that our knowledge of the median redshift of bright SMGs population has not been biased by the low significance of the source catalogues employed.

preprint2005arXiv

The evolution of Ks-selected galaxies in the GOODS/CDFS deep ISAAC field

We present estimated redshifts and derived properties for a sample of 1663 galaxies with Ks <= 22 (Vega), selected from 50.4 sq.arcmin of the GOODS/CDFS field with deep ISAAC imaging, and make an extensive comparison of their properties with those of the extremely red galaxies (ERGs) selected in the same field. We study in detail the evolution of Ks-selected galaxies up to redshifts z ~ 4 and clarify the role of ERGs within the total Ks-band galaxy population. We compute the total Ks-band luminosity function (LF) and compare its evolution with the ERG LF. Up to <z_phot>=2.5, the bright end of the Ks-band LF shows no sign of decline, and is progressively well reproduced by the ERGs with increasing redshift. We also explore the evolution of massive systems present in our sample: up to 20%-25% of the population of local galaxies with assembled stellar mass M>1x10^11 Msun has been formed before redshift z ~ 4, and contains ~ 45% to 70% of the stellar mass density of the Universe at that redshift. Within our sample, the comoving number density of these massive systems is then essentially constant down to redshift z ~ 1.5, by which point most of them have apparently evolved into (I-Ks)-selected ERGs. The remaining massive systems observed in the local Universe are assembled later, at redshifts z <= 1.5. Our results therefore suggest a two-fold assembly history for massive galaxies, in which galaxy/star formation proceeds very efficiently in high mass haloes at very high redshift.

preprint2002arXiv

Discovery of the host galaxy of HDF850.1, the brightest sub-mm source in the Hubble Deep Field

Despite extensive observational efforts, the brightest sub--mm source in the Hubble Deep Field, HDF850.1, has failed to yield a convincing optical/infrared identification almost 4 years after its discovery. This failure is all the more notable given the availability of supporting multi-frequency data of unparalleled depth, and sub-arcsec positional accuracy for the sub-mm/mm source. Consequently, HDF850.1 has become a test case of the possibility that the most violently star-forming objects in the universe are too red and/or distant to be seen in the deepest optical images. Here we report the discovery of the host galaxy of HDF850.1. This object has been revealed by careful analysis of a new, deep K-prime image of the HDF obtained with the Subaru 8.2-m telescope. Its reality is confirmed by a similar analysis of the HST NICMOS F160W image of the same region. This object is extremely faint (K=23.5), clumpy (on sub-arcsec scales) and very red (I-K > 5.2; H-K = 1.4 +- 0.35). The likelihood that it is the correct identification is strongly reinforced by a reanalysis of the combined MERLIN+VLA 1.4-GHz map of the field which yields a new radio detection of HDF850.1 only 0.1 arcsec from the new near-ir counterpart, and with sufficient positional accuracy to exclude all previously considered alternative optical candidates. We have calculated new confidence limits on the estimated redshift of HDF850.1 and find z = 4.1 +- 0.5. We also calculate that the flux density of HDF850.1 has been boosted by a factor of ~3 through lensing by the intervening elliptical 3-586.0, consistent with predictions that a small but significant fraction of blank-field sub-mm sources are lensed by foreground galaxies. We discuss the wider implications of these results for the sub-mm population and cosmic star-formation history.

preprint2001arXiv

A NICMOS imaging study of high-z quasar host galaxies

We present the first results from a major Hubble Space Telescope program designed to investigate the cosmological evolution of quasar host galaxies from z~2 to the present day. Here we describe J and H-band NICMOS imaging of two quasar samples at redshifts of 0.9 and 1.9 respectively. Each sample contains equal numbers of radio-loud and radio-quiet quasars, selected to lie within the same narrow range of optical absolute magnitude (-24 > M_V > -25). Filter and target selection were designed to ensure that at each redshift the images sample the same part of the object's rest-frame spectrum, avoiding potential contamination by [OIII]lambda5007 and H-alpha emission lines. At z=1 the hosts of both radio-loud and radio-quiet quasars lie on the same Kormendy relation described by 3CR radio galaxies at comparable redshift. There is some evidence for a gap of ~1 mag between the host luminosities of RLQs and RQQs, a difference that cannot be due to emission-line contamination given the design of our study. However, within current uncertainties, simple passive stellar evolution is sufficient to link these galaxies with the elliptical hosts of low-redshift quasars of comparable nuclear output, implying that the hosts are virtually fully assembled by z=1. At z=2 the luminosity gap appears to have widened further to ~1.5 mag. Thus while the hosts of radio-loud quasars remain consistent with a formation epoch of z>3, allowing for passive evolution implies that the hosts of radio-quiet quasars are ~2-4 times less massive at z=2 than at low z.

preprint2001arXiv

On the black hole - bulge mass relation in active and inactive galaxies

New black-hole mass estimates are presented for a sample of 72 AGN covering three decades in optical luminosity. Using a sub-sample of Seyfert galaxies, which have black-hole mass estimates from both reverberation mapping and stellar velocity dispersions, we investigate the geometry of the AGN broad-line region (BLR). It is demonstrated that a model in which the orbits of the line-emitting material have a flattened geometry is favoured over randomly orientated orbits. Using this model we investigate the M_{bh}-L_{bulge} relation for a combined 90-object sample, consisting of the AGN plus a sample of 18 nearby inactive elliptical galaxies with dynamical black-hole mass measurements. It is found that, for all reasonable mass-to-light ratios, the M_{bh}-L_{bulge} relation is equivalent to a linear scaling between bulge and black-hole mass. The best-fitting normalization of the M_{bh}-M_{bulge} relation is found to be M_{bh}=0.0012M_{bulge}, in agreement with recent black-hole mass studies based on stellar velocity dispersions. Furthermore, the scatter around the M_{bh}-L_{bulge} relation for the full sample is found to be significantly smaller than has been previously reported (Delta \log M_{bh}=0.39 dex). Finally, using the nearby inactive elliptical galaxy sample alone, it is shown that the scatter in the M_{bh}-L_{bulge} relation is only 0.33 dex, comparable to that of the M_{bh}-sigma relation. These results indicate that reliable black-hole mass estimates can be obtained for high redshift galaxies.

preprint2000arXiv

Starburst galaxies and structure in the submillimetre background towards the Hubble Deep Field

We use an 850-micron SCUBA map of the HDF to study the dust properties of optically-selected starburst galaxies at high redshift. Although the 850-micron map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S_{850}=0.2 mJy for an apparent UV star-formation rate of 1 h^{-2} msun yr^{-1}. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-micron background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-micron background originates over an extremely wide range of redshifts -- perhaps 1 < z < 6.

preprint1998arXiv

The radio properties of radio-quiet quasars

Although radio-quiet quasars (RQQs) constitute >90% of optically-identified quasar samples their radio properties are only poorly understood. In this paper we present the results of a multi-frequency VLA study of 27 low-redshift RQQs. We detect radio emission from 20 objects, half of which are unresolved (< 0.24''). In cases where significant structure can be resolved, double, triple and linear radio sources on scales of a few kpc are found. The radio emission (typically) has a steep spectrum (alpha ~ 0.7), and high brightness temperatures (T_B > 10^5 K) are measured in some of the radio components. The RQQs form a natural extension to the radio luminosity - absolute magnitude distribution of nearby Seyfert 1s. We conclude that a significant fraction of the radio emission in RQQs originates in a compact nuclear source directly associated with the quasar. There are no significant differences between the radio properties of RQQs with elliptical hosts and those in disc galaxies within the current sample.

preprint1997arXiv

High-redshift radio galaxies and quasars at sub-millimetre wavelengths: assessing their evolutionary status

We present new results of a study of the sub-millimetre continuum emission from a sample of 9 radio galaxies and 4 radio-quiet quasars at redshifts $z = 0.75 - 4.26$. The observations were made at 800$μm$, using the single-element bolometer UKT14 on the James Clerk Maxwell Telescope (JCMT), reaching a typical r.m.s. sensitivity of $σ_{rms} \sim 4$ mJy and represent some of the deepest submillimetre extragalactic measurements made to date. Three detections were achieved, of which two are secure (4C41.17, Dunlop et al. (1994) and H1413$+$117) and one (53W002) is tentative, whilst comparable upper-limits were obtained for 7 of the 10 remaining sources. We use these data as the motivation for a detailed discussion of the conversion from submillimetre and millimetre continuum fluxes to dust/gas masses and star formation rates at high-redshift, and determine these quantities from our own and other data on high-redshift radio-galaxies and quasars. In particular we have investigated the impact of the four main sources of uncertainty in deriving physical quantities from such data, namely i) potential contamination by galactic cirrus, ii) uncertainty in the value of the dust rest-frequency mass-absorption coefficient, iii) difficulty in constraining the dust temperature, and iv) estimation of the appropriate gas:dust ratio in high-redshift objects. Our discussion emphasises how important it will be to quantify and, where possible, minimize such uncertainties (via, for example, appropriate observational strategies) in order to fully capitalize on the ten-fold improvement in sensitivity offered by the imminent arrival of the next generation of bolometer arrays, such as SCUBA on the JCMT.

preprint1997arXiv

Using New Submillimetre Surveys to Identify the Evolutionary Status of High-z Galaxies

This paper describes a key submillimetre survey which we are currently conducting to address some of the outstanding questions in cosmology - how, at what epoch and over what period of time did massive galaxies form at high redshift? A summary of the technical feasibility of future submillimetre observations with new ground-based, airborne and satellite telescopes is also presented.