Search

In this work, we use 8 years of deep near-infrared imaging to select and study a new set of 601 active galaxies identified through long-term near-infrared (NIR) variability in the UKIDSS Ultra Deep Survey (UDS). These objects are compared to 710 X-ray bright AGN detected by the Chandra X-ray observatory. We show that infrared variability and X-ray emission select distinct sets of active galaxies, finding only a 37 per cent overlap of galaxies detected by both techniques and confirming NIR-variable AGN to be typically X-ray quiet. Examining the mass functions of the active galaxies shows that NIR variability detects AGN activity in galaxies over a significantly wider range of host stellar mass compared to X-ray detection. For example, at z $\sim$ 1, variable AGN are identified among approximately 1 per cent of galaxies in a roughly flat distribution above the stellar mass completeness limit (> 10$^{9}$M$_{\odot}$), while X-ray detection primarily identifies AGN in galaxies of higher mass (> 10$^{10}$M$_{\odot}$). We conclude that long-term near-infrared variability provides an important new tool for obtaining more complete samples of AGN in deep survey fields., Comment: Accepted for publication in MNRAS. 12 pages, 11 figures

We report ultra-deep, medium-resolution spectroscopic observations for 4 quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot})>11$ at $3 < z < 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first two galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by $860$ pkpc on the sky, within a larger structure for which we confirm several other members. Both formed at $z\simeq8-10$. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other two ultra-massive quiescent galaxies are previously known at $z=3.99$ and $3.19$, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $\Lambda$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $\alpha-$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these 4 galaxies are consistent with the $\Lambda$-CDM halo-mass function using an extreme value statistics approach. We find that the $z=4.62$ objects and the $z=3.19$ object are unlikely within our area under the assumption of standard stellar fractions ($f_*\simeq0.1-0.2$). However, these objects roughly align with the most massive galaxies expected under the assumption of 100 per cent conversion of baryons to stars ($f_*$=1). Our results suggest extreme galaxy formation physics during the first billion years, but no conflict with $\Lambda$-CDM cosmology., Comment: 25 pages, 10 figures, accepted for publication in MNRAS

We conduct the first study of how the relative quenching probability of galaxies depends on environment over the redshift range $0.5 < z < 3$, using data from the UKIDSS Ultra-Deep Survey. By constructing the stellar mass functions for quiescent and post-starburst (PSB) galaxies in high, medium and low density environments to $z = 3$, we find an excess of quenched galaxies in dense environments out to at least $z \sim 2$. Using the growth rate in the number of quenched galaxies, combined with the star-forming galaxy mass function, we calculate the probability that a given star-forming galaxy is quenched per unit time. We find a significantly higher quenching rate in dense environments (at a given stellar mass) at all redshifts. Massive galaxies (M$_* > 10^{10.7}$ M$_{\odot}$) are on average 1.7 $\pm$ 0.2 times more likely to quench per Gyr in the densest third of environments compared to the sparsest third. Finally, we compare the quiescent galaxy growth rate to the rate at which galaxies pass through a PSB phase. Assuming a visibility timescale of 500 Myr, we find that the PSB route can explain $\sim$ 50\% of the growth in the quiescent population at high stellar mass (M$_* > 10^{10.7}$ M$_{\odot}$) in the redshift range $0.5 < z < 3$, and potentially all of the growth at lower stellar masses., Comment: 12 pages, 8 figures. Accepted for publication in MNRAS

We present new empirical constraints on the evolution of $\rho_{\rm H_2}$, the cosmological mass density of molecular hydrogen, back to $z\approx2.5$. We employ a statistical approach measuring the average observed $850\mu{\rm m}$ flux density of near-infrared selected galaxies as a function of redshift. The redshift range considered corresponds to a span where the $850\mu{\rm m}$ band probes the Rayleigh-Jeans tail of thermal dust emission in the rest-frame, and can therefore be used as an estimate of the mass of the interstellar medium (ISM). Our sample comprises of ${\approx}150,000$ galaxies in the UKIDSS-UDS field with near-infrared magnitudes $K_{\rm AB}\leq25$ mag and photometric redshifts with corresponding probability distribution functions derived from deep 12-band photometry. With a sample approximately 2 orders of magnitude larger than in previous works we significantly reduce statistical uncertainties on $\rho_{\rm H_2}$ to $z\approx2.5$. Our measurements are in broad agreement with recent direct estimates from blank field molecular gas surveys, finding that the epoch of molecular gas coincides with the peak epoch of star formation with $\rho_{\rm H_2}\approx2\times10^7\,{\rm M_\odot}\,{\rm Mpc^{-3}}$ at $z\approx2$. We demonstrate that $\rho_{\rm H_2}$ can be broadly modelled by inverting the star-formation rate density with a fixed or weakly evolving star-formation efficiency. This 'constant efficiency' model shows a similar evolution to our statistically derived $\rho_{\rm H_2}$, indicating that the dominant factor driving the peak star formation history at $z\approx2$ is a larger supply of molecular gas in galaxies rather than a significant evolution of the star-formation rate efficiency within individual galaxies., Comment: Accepted for publication in ApJ

Information on the structure around active galactic nuclei (AGN) has long been derived from measuring lags in their varying light output at different wavelengths. In principle, infrared data would reach to larger radii, potentially even probing reprocessed radiation in any surrounding dusty torus. In practice, this has proved challenging because high quality data are required to detect such variability, and the observations must stretch over a long period to probe the likely month-scale lags in variability. In addition, large numbers of sources would need to be observed to start searching for any patterns in such lags. Here, we show that the UKIDSS Ultra Deep Survey, built up from repeated observations over almost a decade, provides an ideal data set for such a study. For 94 sources identified as strongly-varying AGN within its square-degree field, we find that the K-band light curves systematically lag the J-band light curves by an average of around a month. The lags become smaller at higher redshift, consistent with the band shift to optical rest-frame emission. The less luminous AGN also display shorter lags, as would be expected if their physical size scales with luminosity., Comment: 4 pages, 4 figures, accepted for publication as a letter in MNRAS

We present an analysis of the spatial clustering of a large sample of high-resolution, interferometically identified, submillimetre galaxies (SMGs). We measure the projected cross-correlation function of ~350 SMGs in the UKIDSS Ultra Deep-Survey Field across a redshift range of $z=1.5-3$ utilising a method that incorporates the uncertainties in the redshift measurements for both the SMGs and cross-correlated galaxies through sampling their full probability distribution functions. By measuring the absolute linear bias of the SMGs we derive halo masses of $\log_{10}(M_{\rm halo}[{h^{-1}\,\rm M_{\odot}}])\sim12.8$ with no evidence of evolution in the halo masses with redshift, contrary to some previous work. From considering models of halo mass growth rates we predict that the SMGs will reside in haloes of mass $\log_{10}(M_{\rm halo}[{h^{-1}\,\rm M_{\odot}}])\sim13.2$ at $z=0$, consistent with the expectation that the majority of $z=1.5-3$ SMGs will evolve into present-day spheroidal galaxies. Finally, comparing to models of stellar-to-halo mass ratios, we show that SMGs may correspond to systems that are maximally efficient at converting their gas reservoirs into stars. We compare them to a simple model for gas cooling in halos that suggests that the unique properties of the SMG population, including their high levels of star-formation and their redshift distribution, are a result of the SMGs being the most massive galaxies that are still able to accrete cool gas from their surrounding intragalactic medium., Comment: 13 pages, 7 figures, submitted to MNRAS

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high signal to noise, medium resolution spectra of galaxies at redshift between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We give a detailed description of sample selection, observations and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high Signal to Noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5{\AA} allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-IR photometry and physical parameters derived through SED fitting. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range 8.3< Log(M*/Msolar)<11.7. All catalogues and spectra are accessible through the survey database (http://vandels.inaf.it) where all information can be queried interactively, and via the ESO Archive (https://www.eso.org/qi/)., Comment: 16 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

We analyse a robust sample of 30 near-infrared-faint (K>25.3, 5 sigma) submillimetre galaxies selected across a 0.96 deg^2 field, to investigate their properties and the cause of their lack of detectable optical/near-infrared emission. Our analysis exploits precise identifications based on ALMA 870um continuum imaging, combined with the very deep near-infrared imaging from the UKIDSS-UDS survey. We estimate that K>25.3 submillimetre galaxies represent 15+/-2 per cent of the total population brighter than S870=3.6mJy, with an expected surface density of ~450/deg^2 above S870>1mJy. As such they pose a source of contamination in surveys for both high-redshift "quiescent" galaxies and very-high-redshift Lyman-break galaxies. We show that these K-faint submillimetre galaxies are simply the tail of the broader submillimetre population, with comparable dust and stellar masses to K<25.3 mag submillimetre galaxies, but lying at significantly higher redshifts (z=3.44+/-0.06 versus z=2.36+/-0.11) and having higher dust attenuation (Av=5.2+/-0.3 versus Av=2.9+/-0.1). We investigate the origin of the strong dust attenuation and find indications that these K-faint galaxies have smaller dust continuum sizes than the K<25.3 galaxies, as measured by ALMA, which suggests their high attenuation is related to their compact sizes. We find a correlation of dust attenuation with star-formation rate surface density (Sigma_SFR), with the K-faint submillimetre galaxies representing the higher-Sigma_SFR and highest-Av galaxies. The concentrated, intense star-formation activity in these systems is likely to be associated with the formation of spheroids in compact galaxies at high redshifts, but as a result of their high obscuration these are completely missed in UV, optical and even near-infrared surveys., Comment: 9 pages, 4 figures. Submitted to MNRAS, comments welcome

We study the radio properties of 706 sub-millimeter galaxies (SMGs) selected at 870$\mu$m with the Atacama Large Millimeter Array from the SCUBA-2 Cosmology Legacy Survey map of the Ultra Deep Survey field. We detect 273 SMGs at $>4\sigma$ in deep Karl G. Jansky Very Large Array 1.4 GHz observations, of which a subset of 45 SMGs are additionally detected in 610 MHz Giant Metre-Wave Radio Telescope imaging. We quantify the far-infrared/radio correlation through parameter $q_\text{IR}$, defined as the logarithmic ratio of the far-infrared and radio luminosity, and include the radio-undetected SMGs through a stacking analysis. We determine a median $q_\text{IR} = 2.20\pm0.03$ for the full sample, independent of redshift, which places these $z\sim2.5$ dusty star-forming galaxies $0.44\pm0.04$ dex below the local correlation for both normal star-forming galaxies and local ultra-luminous infrared galaxies (ULIRGs). Both the lack of redshift-evolution and the offset from the local correlation are likely the result of the different physical conditions in high-redshift starburst galaxies, compared to local star-forming sources. We explain the offset through a combination of strong magnetic fields ($B\gtrsim0.2$mG), high interstellar medium (ISM) densities and additional radio emission generated by secondary cosmic rays. While local ULIRGs are likely to have similar magnetic field strengths, we find that their compactness, in combination with a higher ISM density compared to SMGs, naturally explains why local and high-redshift dusty star-forming galaxies follow a different far-infrared/radio correlation. Overall, our findings paint SMGs as a homogeneous population of galaxies, as illustrated by their tight and non-evolving far-infrared/radio correlation., Comment: Accepted for publication in ApJ; 18 pages + appendices; 7 figures in main text

We present the first attempt to select AGN using long-term NIR variability. By analysing the K-band light curves of all the galaxies in the UKIDSS Ultra Deep Survey, the deepest NIR survey over ~1 sq degree, we have isolated 393 variable AGN candidates. A comparison to other selection techniques shows that only half of the variable sources are also selected using either deep Chandra X-ray imaging or IRAC colour selection, suggesting that using NIR variability can locate AGN that are missed by more standard selection techniques. In particular, we find that long-term NIR variability identifies AGN at low luminosities and in host galaxies with low stellar masses, many of which appear relatively X-ray quiet., Comment: 12 pages, 18 figures, accepted for publication in MNRAS, an error in Equation 1 has been fixed in this version

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., Comment: 17 pages, 12 figures. Resubmitted to MNRAS after referee report

We analyse the physical properties of a large, homogeneously selected sample of ALMA-located sub-mm galaxies (SMGs). This survey, AS2UDS, identified 707 SMGs across the ~1 sq.deg. field, including ~17 per cent, which are undetected at $K$>~25.7 mag. We interpret their ultraviolet-to-radio data using MAGPHYS and determine a median redshift of z=2.61+-0.08 (1$\sigma$ range of z=1.8-3.4) with just ~6 per cent at z>4. Our survey provides a sample of massive dusty galaxies at z>~1, with median dust and stellar masses of $M_d$=(6.8+-0.3)x10$^{8}$M$_\odot$ (thus, gas masses of ~10$^{11}$M$_\odot$) and $M_\ast=$(1.26+-0.05)x10$^{11}$M$_\odot$. We find no evolution in dust temperature at a constant far-infrared luminosity across z~1.5-4. The gas mass function of our sample increases to z~2-3 and then declines at z>3. The space density and masses of SMGs suggest that almost all galaxies with $M_\ast$>~3x10$^{11}$M$_\odot$ have passed through an SMG-like phase. The redshift distribution is well fit by a model combining evolution of the gas fraction in halos with the growth of halo mass past a threshold of $M_h$~6x10$^{12}$M$_\odot$, thus SMGs may represent the highly efficient collapse of gas-rich massive halos. We show that SMGs are broadly consistent with simple homologous systems in the far-infrared, consistent with a centrally illuminated starburst. Our study provides strong support for an evolutionary link between the active, gas-rich SMG population at z>1 and the formation of massive, bulge-dominated galaxies across the history of the Universe., Comment: 33 pages, 16 figures, MNRAS accepted version

We present an analysis of the morphology and profiles of the dust continuum emission in 153 bright sub-millimetre galaxies (SMGs) detected with ALMA at S/N ratios of $>8$ in high-resolution $0.18''$ ($\sim1$kpc) 870$\mu$m maps. We measure sizes, shapes and light profiles for the rest-frame far-infrared emission from these luminous star-forming systems and derive a median effective radius ($R_e$) of $0.10''\pm0.04''$ for our sample with a median flux of $S_{870}=5.6\pm0.2$mJy. We find that the apparent axial ratio ($b/a$) distribution of the SMGs peaks at $b/a\sim0.63\pm0.24$ and is best described by triaxial morphologies, while their emission profiles are best fit by a Sersic model with $n\simeq1.0\pm0.1$, similar to exponential discs. This combination of triaxiality and $n\sim1$ Sersic index are characteristic of bars and we suggest that the bulk of the 870$\mu$m dust continuum emission in the central $\sim2$kpc of these galaxies arises from bar-like structures. By stacking our 870$\mu$m maps we recover faint extended dust continuum emission on $\sim4$kpc scales which contributes $13\pm1$% of the total 870$\mu$m emission. The scale of this extended emission is similar to that seen for the molecular gas and rest-frame optical light in these systems, suggesting that it represents an extended dust and gas disc at radii larger than the more active bar component. Including this component in our estimated size of the sources we derive a typical effective radius of $\simeq0.15''\pm0.05''$ or $1.2\pm0.4$kpc. Our results suggest that kpc-scale bars are ubiquitous features of high star-formation rate systems at $z\gg1$, while these systems also contain fainter and more extended gas and stellar envelopes. We suggest that these features, seen some $10-12$Gyrs ago, represent the formation phase of the earliest galactic-scale components: stellar bulges., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

We combine high-resolution ALMA and HST/CANDELS observations of 20 submillimeter galaxies (SMGs) predominantly from the AS2UDS survey at z~2 with bright rest-frame optical counterparts (Ks < 22.9) to investigate the resolved structural properties of their dust and stellar components. We derive two-dimensional stellar-mass distributions that are inferred from spatial mass-to-light ratio (M/L) corrections based on rest-frame optical colors. Due to the high central column densities of dust in our SMGs, our mass distributions likely represent a lower limit to the true central mass density. The centroid positions between the inferred stellar-mass and the dust distributions agree within 1.1 kpc, indicating an overall good spatial agreement between the two components. The majority of our sources exhibit compact dust configurations relative to the stellar component (with a median ratio of effective radii Re,dust/Re,Mstar = 0.6). This ratio does not change with specific star-formation rate (sSFR) over the factor of 30 spanned by our targets, sampling the locus of "normal" main sequence galaxies up to the starburst regime, log(sSFR/sSFRMS) > 0.5. Our results imply that massive SMGs are experiencing centrally enhanced star formation unlike typical spiral galaxies in the local Universe. The sizes and stellar densities of our SMGs are in agreement with those of the passive population at z=1.5, consistent with these systems being the descendants of z~2 SMGs., Comment: 26 pages, 11 figures, accepted for publication in ApJ

We analyse the multi-frequency radio spectral properties of $41$ 6GHz-detected ALMA-identified, submillimetre galaxies (SMGs), observed at 610MHz, 1.4GHz, 6GHz with GMRT and the VLA. Combining high-resolution ($\sim0.5''$) 6GHz radio and ALMA $870\,\mu$m imaging (tracing rest-frame $\sim20$GHz, and $\sim250\,\mu$m dust continuum), we study the far-infrared/radio correlation via the logarithmic flux ratio $q_{\rm IR}$, measuring $\langle q_{\rm IR}\rangle=2.20\pm 0.06$ for our sample. We show that the high-frequency radio sizes of SMGs are $\sim1.9\pm 0.4\times$ ($\sim2$-$3$kpc) larger than those of the cool dust emission, and find evidence for a subset of our sources being extended on $\sim 10$kpc scales at 1.4GHz. By combining radio flux densities measured at three frequencies, we can move beyond simple linear fits to the radio spectra of high-redshift star-forming galaxies, and search for spectral curvature, which has been observed in local starburst galaxies. At least a quarter (10/41) of our sample show evidence of a spectral break, with a median $\langle\alpha^{1.4\,{\rm GHz}}_{610\,{\rm GHz}}\rangle=-0.60\pm 0.06$, but $\langle\alpha^{6\,{\rm GHz}}_{1.4\,{\rm GHz}}\rangle=-1.06\pm 0.04$ -- a high-frequency flux deficit relative to simple extrapolations from the low-frequency data. We explore this result within this subset of sources in the context of age-related synchrotron losses, showing that a combination of weak magnetic fields ($B\sim35\,\mu$G) and young ages ($t_{\rm SB}\sim40$--$80\,$Myr) for the central starburst can reproduce the observed spectral break. Assuming these represent evolved (but ongoing) starbursts and we are observing these systems roughly half-way through their current episode of star formation, this implies starburst durations of $\lesssim100$Myr, in reasonable agreement with estimates derived via gas depletion timescales., Comment: 25 pages, 10 colour figures, 2 tables; accepted for publication in ApJ

We present the catalogue and properties of sources in AS2UDS, an 870-$\mu$m continuum survey with the Atacama Large Millimetre/sub-millimetre Array (ALMA) of 716 single-dish sub-millimetre sources detected in the UKIDSS/UDS field by the SCUBA-2 Cosmology Legacy Survey. In our sensitive ALMA follow-up observations we detect 708 sub-millimetre galaxies (SMGs) at $>$\,4.3$\sigma$ significance across the $\sim$\,1-degree diameter field. We combine our precise ALMA positions with the extensive multi-wavelength coverage in the UDS field to fit the spectral energy distributions of our SMGs to derive a median redshift of $z_{\rm phot}=$\,2.61$\pm$0.09. This large sample reveals a statistically significant trend of increasing sub-millimetre flux with redshift suggestive of galaxy downsizing. 101 ALMA maps do not show a $>$\,4.3$\sigma$ SMG, but we demonstrate from stacking {\it Herschel} SPIRE observations at these positions, that the vast majority of these blank maps correspond to real single-dish sub-millimetre sources. We further show that these blank maps contain an excess of galaxies at $z_{\rm phot}=$\,1.5--4 compared to random fields, similar to the redshift range of the ALMA-detected SMGs. In addition, we combine X-ray and mid-infrared active galaxy nuclei activity (AGN) indicators to yield a likely range for the AGN fraction of 8--28\,\% in our sample. Finally, we compare the redshifts of this population of high-redshift, strongly star-forming galaxies with the inferred formation redshifts of massive, passive galaxies being found out to $z\sim$\,2, finding reasonable agreement -- in support of an evolutionary connection between these two classes of massive galaxy., Comment: Submitted to MNRAS; Comments welcome. Full catalogue will be made publicly available on acceptance of paper

We describe the application of the supervised machine-learning algorithms to identify the likely multi-wavelength counterparts to submillimeter sources detected in panoramic, single-dish submillimeter surveys. As a training set, we employ a sample of 695 ($S_{\rm 870\mu m}$ >1 mJy) submillimeter galaxies (SMGs) with precise identifications from the ALMA follow-up of the SCUBA-2 Cosmology Legacy Survey's UKIDSS-UDS field (AS2UDS). We show that radio emission, near-/mid-infrared colors, photometric redshift, and absolute $H$-band magnitude are effective predictors that can distinguish SMGs from submillimeter-faint field galaxies. Our combined radio+machine-learning method is able to successfully recover $\sim$85 percent of ALMA-identified SMGs which are detected in at least three bands from the ultraviolet to radio. We confirm the robustness of our method by dividing our training set into independent subsets and using these for training and testing respectively, as well as applying our method to an independent sample of $\sim$100 ALMA-identified SMGs from the ALMA/LABOCA ECDF-South Survey (ALESS). To further test our methodology, we stack the 870$\mu$m ALMA maps at the positions of those $K$-band galaxies that are classified as SMG counterparts by the machine-learning but do not have a $>$4.3$ \sigma$ ALMA detection. The median peak flux density of these galaxies is $S_{\rm 870\mu m}=(0.61\pm0.03)$ mJy, demonstrating that our method can recover faint and/or diffuse SMGs even when they are below the detection threshold of our ALMA observations. In future, we will apply this method to samples drawn from panoramic single-dish submillimeter surveys which currently lack interferometric follow-up observations, to address science questions which can only be tackled with large, statistical samples of SMGs., Comment: 25 pages, 10 figures, three tables, accepted for publication in ApJ

We report the first results of AS2UDS: an 870 $\mu$m continuum survey with the Atacama Large Millimeter/Submillimeter Array (ALMA) of a total area of $\sim$ 50 arcmin$^2$ comprising a complete sample of 716 submillimeter sources drawn from the SCUBA-2 Cosmology Legacy Survey (S2CLS) map of the UKIDSS/UDS field. The S2CLS parent sample covers a 0.96 degree$^2$ field at $\sigma_{850}=0.90\pm0.05$ mJy beam$^{-1}$. Our deep, high-resolution ALMA observations with $\sigma_{\rm 870}\sim$ 0.25 mJy and a 0.15"--0.30" FWHM synthesized beam, provide precise locations for 695 submillimetre galaxies (SMGs) responsible for the submillimeter emission corresponding to 606 sources in the low resolution, single-dish map. We measure the number counts of SMGs brighter than $S_{\rm 870}\geq$ 4 mJy, free from the effects of blending and show that the normalisation of the counts falls by 28 $\pm$ 2% in comparison to the SCUBA-2 parent sample, but that the shape remains unchanged. We determine that 44$^{+16}_{-14}$% of the brighter single-dish sources with $S_{850}\geq$ 9 mJy consist of a blend of two or more ALMA-detectable SMGs brighter than $S_{870}\sim$ 1 mJy (corresponding to a galaxy with a total-infrared luminosity of $L_{IR}\geq$ 10$^{12}$ L$_\odot$), in comparison to 28 $\pm$ 2% for the single-dish sources at $S_{\rm 850}\geq$ 5 mJy. Using the 46 single-dish submillimeter sources that contain two or more ALMA-detected SMGs with photometric redshifts, we show that there is a significant statistical excess of pairs of SMGs with similar redshifts ($<$1% probability of occurring by chance), suggesting that at least 30% of these blends arise from physically associated pairs of SMGs., Comment: 11 pages, 4 figures, fixed typo (sorry Rob)

We report the results of a search for serendipitous [CII] 157.74$\mu$m emitters at $z\sim4.4$-$4.7$ using the Atacama Large Millimeter/submillimeter Array (ALMA). The search exploits the AS2UDS continuum survey, which covers ~50 arcmin$^2$ of the sky towards 695 luminous ($S_{870}\gtrsim1$mJy) submillimeter galaxies (SMGs), selected from the SCUBA-2 Cosmology Legacy Survey (S2CLS) 0.96deg$^2$ Ultra Deep Survey (UDS) field. We detect ten candidate line emitters, with an expected false detection rate of ten percent. All of these line emitters correspond to 870$\mu$m continuum-detected sources in AS2UDS. The emission lines in two emitters appear to be high-J CO, but the remainder have multi-wavelength properties consistent with [CII] from $z\simeq4.5$ galaxies. Using our sample, we place a lower limit of $>5\times10^{-6}$Mpc$^{-3}$ on the space density of luminous ($L_{\rm IR} \simeq 10^{13}$Lsun) SMGs at $z=4.40$-$4.66$, suggesting $\ge7$percent of SMGs with $S_{870\mu{\rm m}}\gtrsim1$mJy lie at $4

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection., Comment: 19 pages, 9 figures, accepted for publication in MNRAS

We use a sample of 4178 Lyman break galaxies (LBGs) at z = 3, 4 and 5 in the UKIRT Infrared Deep Sky Survey (UKIDSS) Ultra Deep Survey (UDS) field to investigate the relationship between the observed slope of the stellar continuum emission in the ultraviolet, {\beta}, and the thermal dust emission, as quantified via the so-called 'infrared excess' (IRX = LIR/LUV). Through a stacking analysis we directly measure the 850-{\mu}m flux density of LBGs in our deep (0.9mJy) James Clerk Maxwell Telescope (JCMT) SCUBA-2 850-{\mu}m map, as well as deep public Herschel/SPIRE 250-, 350- and 500-{\mu}m imaging. We establish functional forms for the IRX-{\beta} relation to z ~ 5, confirming that there is no significant redshift evolution of the relation and that the resulting average IRX-{\beta} curve is consistent with a Calzetti-like attenuation law. We compare our results with recent work in the literature, finding that discrepancies in the slope of the IRX-{\beta} relation are driven by biases in the methodology used to determine the ultraviolet slopes. Consistent results are found when IRX-{\beta} is evaluated by stacking in bins of stellar mass, M, and we argue that the near-linear IRX-M relationship is a better proxy for correcting observed UV luminosities to total star formation rates, provided an accurate handle on M can be had, and also gives clues as to the physical driver of the role of dust-obscured star formation in high-redshift galaxies., Comment: 12 pages, 10 figures, submitted to MNRAS

This paper defines the UK Infra-Red Telescope (UKIRT) Hemisphere Survey (UHS) and release of the remaining ~12,700 sq.deg of J-band survey data products. The UHS will provide continuous J and K-band coverage in the northern hemisphere from a declination of 0 deg to 60 deg by combining the existing Large Area Survey, Galactic Plane Survey and Galactic Clusters Survey conducted under the UKIRT Infra-red Deep Sky Survey (UKIDSS) programme with this new additional area not covered by UKIDSS. The released data includes J-band imaging and source catalogues over the new area, which, together with UKIDSS, completes the J-band UHS coverage over the full ~17,900 sq.deg area. 98 per cent of the data in this release have passed quality control criteria, the remaining 2 per cent being scheduled for re-observation. The median 5-sigma point source sensitivity of the released data is 19.6 mag (Vega). The median full width at half-maximum of the point spread function across the dataset is 0.75 arcsec. In this paper, we outline the survey management, data acquisition, processing and calibration, quality control and archiving as well as summarising the characteristics of the released data products. The data are initially available to a limited consortium with a world-wide release scheduled for August 2018., Comment: 14 pages, 9 figures, MNRAS in press. v2 includes some changes suggested by referee

We present a multi-wavelength analysis of 52 sub-millimeter galaxies (SMGs), identified using ALMA 870$\mu$m continuum imaging in a pilot program to precisely locate bright SCUBA2-selected sub-mm sources in the UKIDSS Ultra Deep Survey (UDS) field. Using the available deep (especially near-infrared), panoramic imaging of the UDS field at optical-to-radio wavelengths we characterize key properties of the SMG population. The median photometric redshift of the bright ALMA/SCUBA-2 UDS (AS2UDS) SMGs that are detected in a sufficient number of wavebands to derive a robust photometric redshift is $z$=2.65$\pm$0.13. However, similar to previous studies, 27% of the SMGs are too faint at optical-to-near-infrared wavelengths to derive a reliable photometric redshift. Assuming that these SMGs lie at z$\gtrsim$3 raises the median redshift of the full sample to $z$=2.9$\pm$0.2. A subset of 23, unlensed, bright AS2UDS SMGs have sizes measured from resolved imaging of their rest-frame far-infrared emission. We show that the extent and luminosity of the far-infrared emission are consistent with the dust emission arising from regions that are optically thick, on average, at a wavelength of $\lambda_0$$\ge$75$\mu$m (1-$\sigma$ dispersion of 55-90$\mu$m). Using the dust masses derived from our optically-thick spectral energy distribution models we determine that these galaxies have a median hydrogen column density of $N_{H}$=9.8$_{-0.7}^{+1.4}$$\times$10$^{23}$cm$^{-2}$, or a corresponding median $V$-band obscuration of $A_\mathrm{v}$=540$^{+80}_{-40}$mag, averaged along the line of sight to the source of their restframe $\sim$200$\mu$m emission. We discuss the implications of this extreme attenuation by dust for the multiwavelength study of dusty starbursts and reddening-sensitive tracers of star formation., Comment: Accepted for publication in ApJ

We investigate the effects of dense environments on galaxy evolution by examining how the properties of galaxies in the z = 1.6 protocluster Cl 0218.3-0510 depend on their location. We determine galaxy properties using spectral energy distribution fitting to 14-band photometry, including data at three wavelengths that tightly bracket the Balmer and 4000A breaks of the protocluster galaxies. We find that two-thirds of the protocluster galaxies, which lie between several compact groups, are indistinguishable from field galaxies. The other third, which reside within the groups, differ significantly from the intergroup galaxies in both colour and specific star formation rate. We find that the fraction of red galaxies within the massive protocluster groups is twice that of the intergroup region. These excess red galaxies are due to enhanced fractions of both passive galaxies (1.7 times that of the intergroup region) and dusty star-forming galaxies (3 times that of the intergroup region). We infer that some protocluster galaxies are processed in the groups before the cluster collapses. These processes act to suppress star formation and change the mode of star formation from unobscured to obscured., Comment: Accepted for publication in MNRAS. http://mnras.oxfordjournals.org/content/early/2016/09/16/mnras.stw2359

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., Comment: Submitted to MNRAS, comments welcome. Catalogue and maps at http://dx.doi.org/10.5281/zenodo.57792

We present a comprehensive picture of the Cl0218.3-0510 protocluster at $z=1.623$ across 10 co-moving Mpc. Using filters that tightly bracket the Balmer and 4000 Angstrom breaks of the protocluster galaxies we obtain precise photometric redshifts resulting in a protocluster galaxy sample that is 89+/-5% complete and has a contamination of only 12+/-5%. Both star forming and quiescent protocluster galaxies are located allowing us to map the structure of the forming cluster for the first time. The protocluster contains 6 galaxy groups, the largest of which is the nascent cluster. Only a small minority of the protocluster galaxies are in the nascent cluster (11%) or in the other galaxy groups (22%), as most protocluster galaxies reside between the groups. Unobscured star forming galaxies predominantly reside between the protocluster's groups, whereas red galaxies make up a large fraction of the groups' galactic content, so observing the protocluster through only one of these types of galaxies results in a biased view of the protocluster's structure. The structure of the protocluster reveals how much mass is available for the future growth of the cluster and we use the Millennium Simulation, scaled to a Planck cosmology, to predict that Cl0218.3-0510 will evolve into a 2.7x 10^14 Msun cluster by the present day., Comment: Accepted for publication in MNRAS

We use a sample of 37 of the densest clusters and protoclusters across $1.3 \le z \le 3.2$ from the Clusters Around Radio-Loud AGN (CARLA) survey to study the formation of massive cluster galaxies. We use optical $i'$-band and infrared 3.6$\mu$m and 4.5$\mu$m images to statistically select sources within these protoclusters and measure their median observed colours; $\langle i'-[3.6] \rangle$. We find the abundance of massive galaxies within the protoclusters increases with decreasing redshift, suggesting these objects may form an evolutionary sequence, with the lower redshift clusters in the sample having similar properties to the descendants of the high redshift protoclusters. We find that the protocluster galaxies have an approximately unevolving observed-frame $i'-[3.6]$ colour across the examined redshift range. We compare the evolution of the $\langle i'-[3.6] \rangle$ colour of massive cluster galaxies with simplistic galaxy formation models. Taking the full cluster population into account, we show that the formation of stars within the majority of massive cluster galaxies occurs over at least 2Gyr, and peaks at $z \sim 2$-3. From the median $i'-[3.6]$ colours we cannot determine the star formation histories of individual galaxies, but their star formation must have been rapidly terminated to produce the observed red colours. Finally, we show that massive galaxies at $z>2$ must have assembled within 0.5Gyr of them forming a significant fraction of their stars. This means that few massive galaxies in $z>2$ protoclusters could have formed via dry mergers., Comment: Accepted for publication in MNRAS. 20 pages including 13 figures & 2 appendices

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored., Comment: Accepted for publication by ApJ. Replaced to match the final version (two tables and one figure added in the appendix)

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., Comment: 16 pages, 7 figures, accepted for publication in ApJ

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)., Comment: Accepted for publication on Astronomy and Astrophysics

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

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., Comment: ApJ in press, 16 pages, 8 figures, 1 table

Galaxy formation models invoke the presence of strong feedback mechanisms that regulate the growth of massive galaxies at high redshifts. In this paper we aim to: (1) confirm spectroscopically the redshifts of a sample of massive galaxies selected with photometric redshifts z > 2.5; (2) investigate the properties of their stellar and interstellar media; (3) detect the presence of outflows, and measure their velocities. To achieve this, we analysed deep, high-resolution (R~2000) FORS2 rest-frame UV spectra for 11 targets. We confirmed that 9 out of 11 have spectroscopic redshifts z > 2.5. We also serendipitously found two mask fillers at redshift z > 2.5, which originally were assigned photometric redshifts 2.0 < z < 2.5. In the four highest-quality spectra we derived outflow velocities by fitting the absorption line profiles with models including multiple dynamical components. We found strongly asymmetric, high-ionisation lines, from which we derived outflow velocities ranging from 480 to 1518 km/s. The two galaxies with highest velocity show signs of AGN. We revised the spectral energy distribution fitting U-band through 8 micron photometry, including the analysis of a power-law component subtraction to identify the possible presence of active galactic nuclei (AGN). The revised stellar masses of all but one of our targets are >1e10 Msun, with four having stellar masses > 5e10 Msun. Three galaxies have a significant power-law component in their spectral energy distributions, which indicates that they host AGN. We conclude that massive galaxies are characterised by significantly higher velocity outflows than the typical Lyman break galaxies at z ~ 3. The incidence of high-velocity outflows (~40% within our sample) is also much higher than among massive galaxies at z < 1, which is consistent with the powerful star formation and nuclear activity that most massive galaxies display at z > 2., Comment: 17 pages, 14 figures, Accepted for publication in A&A

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 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\mu m-selected sources is presented, showing a broad peak in the redshift range 1=10^12 +/- 0.8 L_sol, =42 +/- 11 K and <\beta_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 \beta_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 \mu 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., Comment: 20 pages, 12 figures, accepted for publication in MNRAS

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., Comment: 16 pages, 15 figures, accepted by MNRAS

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 =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., Comment: 9 pages, 5 figures, submitted to MNRAS

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= 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)., Comment: 23 pages, 12 figures, for submission to MNRAS

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= 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)., Comment: 21 pages, 11 figures, accepted for publication in MNRAS. Replaced to match accepted version

We present new improved constraints on the Hubble parameter H(z) in the redshift range 0.15 < z < 1.1, obtained from the differential spectroscopic evolution of early-type galaxies as a function of redshift. We extract a large sample of early-type galaxies (\sim11000) from several spectroscopic surveys, spanning almost 8 billion years of cosmic lookback time (0.15 < z < 1.42). We select the most massive, red elliptical galaxies, passively evolving and without signature of ongoing star formation. Those galaxies can be used as standard cosmic chronometers, as firstly proposed by Jimenez & Loeb (2002), whose differential age evolution as a function of cosmic time directly probes H(z). We analyze the 4000 {\AA} break (D4000) as a function of redshift, use stellar population synthesis models to theoretically calibrate the dependence of the differential age evolution on the differential D4000, and estimate the Hubble parameter taking into account both statistical and systematical errors. We provide 8 new measurements of H(z) (see Tab. 4), and determine its change in H(z) to a precision of 5-12% mapping homogeneously the redshift range up to z \sim 1.1; for the first time, we place a constraint on H(z) at z \neq 0 with a precision comparable with the one achieved for the Hubble constant (about 5-6% at z \sim 0.2), and covered a redshift range (0.5 < z < 0.8) which is crucial to distinguish many different quintessence cosmologies. These measurements have been tested to best match a \Lambda CDM model, clearly providing a statistically robust indication that the Universe is undergoing an accelerated expansion. This method shows the potentiality to open a new avenue in constrain a variety of alternative cosmologies, especially when future surveys (e.g. Euclid) will open the possibility to extend it up to z \sim 2., Comment: 34 pages, 15 figures, 6 tables, published in JCAP. It is a companion to Moresco et al. (2012b, http://arxiv.org/abs/1201.6658) and Jimenez et al. (2012, http://arxiv.org/abs/1201.3608). The H(z) data can be downloaded at http://www.physics-astronomy.unibo.it/en/research/areas/astrophysics/cosmology-with-cosmic-chronometers

We present spectroscopic confirmation of ten highly luminous (L >= 2L*) Lyman alpha emitters in the redshift range 6.01= 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= 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., Comment: accepted by MNRAS, 13 pages, 7 figures

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., Comment: 12 pages, 7 figures, Accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. The paper contains 9 pages, 11 figures and 1 table

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 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., Comment: 16 pages, 15 figures, 2 tables. Accepted for publication in MNRAS

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 (=8.4+-1.1kpc; median r{1/2}=7.8), de Vaucouleurs spheroids ( = 4.07+-0.27; median n=3.87), while the sub-mm galaxies appear to be moderately compact (=3.4+-0.3kpc; median r{1/2}=3.1kpc) exponential discs (=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., Comment: 24 pages, 9 figures

We present a study on the clustering of a stellar mass selected sample of 18,482 galaxies with stellar masses M*>10^10M(sun) at redshifts 0.410^13M(sun) we find that this ratio is <0.02, much lower than the universal baryonic mass fraction. We show that the remaining baryonic mass is included partially in stars within satellite galaxies in these haloes, and as diffuse hot and warm gas. We also find that, at a fixed stellar mass, the stellar-to-total-mass ratio increases at lower redshifts. This suggests that galaxies at a fixed stellar mass form later in lower mass dark matter haloes, and earlier in massive haloes. We interpret this as a "halo downsizing" effect, however some of this evolution could be attributed to halo assembly bias., Comment: Accepted for publication in MNRAS. 19 pages, 8 figures and 3 tables.

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., Comment: 15 pages, MNRAS subm