We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) 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_{\lambda}$(Ly$\alpha$) across the range $-58 \unicode{xC5} \lesssim$ $W_{\lambda}$(Ly$\alpha$) $\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_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=-18\unicode{xC5}$) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) $\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_{\lambda}$(Ly$\alpha$) 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.

We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies.

We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$\alpha$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$\alpha$ star formation rate surface density $\Sigma_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. S\'ersic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium.

As part of our KMOS AGN Survey at High-redshift (KASHz), we present spatially-resolved VLT/KMOS and VLT/SINFONI spectroscopic data and ALMA 870$\mu$m continuum imaging of eight $z$=1.4--2.6 moderate AGN ($L_{\rm 2-10 \rm kev}$ = $10^{42} - 10^{45}$ ergs s$^{-1}$). We map [OIII], H$\alpha$ and rest-frame FIR emission to search for any spatial anti-correlation between ionised outflows (traced by the [OIII] line) and star formation (SF; traced by H$\alpha$ and FIR), that has previously been claimed for some high-z AGN and used as evidence for negative and/or positive AGN feedback. Firstly, we conclude that H$\alpha$ is unreliable to map SF inside our AGN host galaxies based on: (i) SF rates inferred from attenuation-corrected H$\alpha$ can lie below those inferred from FIR; (ii) the FIR continuum is more compact than the H$\alpha$ emission by a factor of $\sim 2$ on average; (iii) in half of our sample, we observe significant spatial offsets between the FIR and H$\alpha$ emission, with an average offset of $1.4\pm0.6$ kpc. Secondly, for the five targets with outflows we find no evidence for a spatial anti-correlation between outflows and SF using either H$\alpha$ or FIR as a tracer. This holds for our re-analysis of a famous $z$=1.6 X-ray AGN (`XID 2028') where positive and negative feedback has been previously claimed. Based on our results, any impact on SF by ionised outflows must be subtle, either occurring on scales below our resolution, or on long timescales.

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

HATLAS J084933.4+021443 was identified as a dusty starburst via its rest-frame far-IR emission. Multi-frequency imaging and spectroscopy revealed a cluster of four dusty galaxies at $z = 2.41$, covering 80 kpc. Here, we use ALMA to confirm a more distant, fifth protocluster member, and present X-ray and rest-frame optical imaging spectroscopy of the brightest, an unlensed hyperluminous IR galaxy (HyLIRG). The data reveal broad H \alpha and bright [N II] lines, and bright X-ray emission, characteristics that betray a Type-1 active galactic nucleus (AGN), strengthening evidence that AGN are ubiquitous amongst HyLIRGs. The accreting black hole is super massive, $M_{\rm bh}\approx 2\times 10^9$ M$_\odot$, with little intrinsic absorption, $N_{\rm H}\approx 5\times 10^{21}$ cm$^{-2}$. The X-ray properties suggest the accretion luminosity rivals that of the starburst, yet it is not obvious where this emerges in its panchromatic spectral energy distribution (SED). We outline three scenarios that could give rise to the observed characteristics, and how we might distinguish between them. In the first, we see the AGN through the host galaxy because of the cavity it excavates. In the others, the AGN is not co-spatial with the starburst, having been ejected via asymmetric gravitational radiation, or having evolved towards the naked quasar phase in an unseen companion.

We present an analysis of the gas outflow energetics of 529 main-sequence star-forming galaxies at z~1 using KMOS observations of the broad, underlying H-alpha and forbidden lines of [N II] and [S II]. Based on the stacked spectra for a sample with median star-formation rates and stellar masses of SFR ~ 7 Mo/yr and M* = (1.0+/-0.1)x10^10 Mo respectively, we derive a typical mass outflow rate of dM/dt = 1-4 Mo/yr and a mass loading of dM/dt/SFR = 0.2--0.4. The mass loading of the wind does not show a strong trend with star-formation rate over the range SFR ~ 2--20 Mo/yr, although we identify a trend with stellar mass such that dM/dt/SFR ~ M*^(0.26+/-0.07). Finally, we find that the line width of the broad H-alpha increases with disk circular velocity with a sub-linear scaling relation FWHM_broad ~ v^(0.21+/-0.05). As a result of this behavior, in the lowest mass galaxies (M* < 10^10 Mo), a significant fraction of the outflowing gas should have sufficient velocity to escape the gravitational potential of the halo whilst in the highest mass galaxies (M* > 10^10 Mo) most of the gas will be retained, flowing back on to the galaxy disk at later times.

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

We analyse maps of the spatially-resolved nebular emission of $\approx$1500 star-forming galaxies at $z\approx0.6$-$2.2$ from deep KMOS and MUSE observations to measure the average shape of their rotation curves. We use these to test claims for declining rotation curves at large radii in galaxies at $z\approx1$-$2$ that have been interpreted as evidence for an absence of dark matter. We show that the shape of the average rotation curves, and the extent to which they decline beyond their peak velocities, depends upon the normalisation prescription used to construct the average curve. Normalising in size by the galaxy stellar disk-scale length after accounting for seeing effects ($R_{\rm{d}}^{\prime}$), we construct stacked position-velocity diagrams that trace the average galaxy rotation curve out to $6R_{\rm{d}}^{\prime}$ ($\approx$13 kpc, on average). Combining these curves with average HI rotation curves for local systems, we investigate how the shapes of galaxy rotation curves evolve over $\approx$10 Gyr. The average rotation curve for galaxies binned in stellar mass, stellar surface mass density and/or redshift is approximately flat, or continues to rise, out to at least $6R_{\rm{d}}^{\prime}$. We find a trend between the outer slopes of galaxies' rotation curves and their stellar mass surface densities, with the higher surface density systems exhibiting flatter rotation curves. Drawing comparisons with hydrodynamical simulations, we show that the average shapes of the rotation curves for our sample of massive, star-forming galaxies at $z\approx0$-$2.2$ are consistent with those expected from $\Lambda$CDM theory and imply dark matter fractions within $6R_{\rm{d}}$ of at least $\approx60$ percent.

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.

We present an X-ray point-source catalog from the XMM-Large Scale Structure survey region (XMM-LSS), one of the XMM-Spitzer Extragalactic Representative Volume Survey (XMM-SERVS) fields. We target the XMM-LSS region with $1.3$ Ms of new XMM-Newton AO-15 observations, transforming the archival X-ray coverage in this region into a 5.3 deg$^2$ contiguous field with uniform X-ray coverage totaling $2.7$ Ms of flare-filtered exposure, with a $46$ ks median PN exposure time. We provide an X-ray catalog of 5242 sources detected in the soft (0.5-2 keV), hard (2-10 keV), and/or full (0.5-10 keV) bands with a 1% expected spurious fraction determined from simulations. A total of 2381 new X-ray sources are detected compared to previous source catalogs in the same area. Our survey has flux limits of $1.7\times10^{-15}$, $1.3\times10^{-14}$, and $6.5\times10^{-15}$ erg cm$^{-2}$ s$^{-1}$ over 90% of its area in the soft, hard, and full bands, respectively, which is comparable to those of the XMM-COSMOS survey. We identify multiwavelength counterpart candidates for 99.9% of the X-ray sources, of which 93% are considered as reliable based on their matching likelihood ratios. The reliabilities of these high-likelihood-ratio counterparts are further confirmed to be $\approx 97$% reliable based on deep Chandra coverage over $\approx 5$% of the XMM-LSS region. Results of multiwavelength identifications are also included in the source catalog, along with basic optical-to-infrared photometry and spectroscopic redshifts from publicly available surveys. We compute photometric redshifts for X-ray sources in 4.5 deg$^2$ of our field where forced-aperture multi-band photometry is available; $>70$% of the X-ray sources in this subfield have either spectroscopic or high-quality photometric redshifts.

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<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

We present the results of a new study of dust attenuation at redshifts $3 < z < 4$ based on a sample of $236$ star-forming galaxies from the VANDELS spectroscopic survey. Motivated by results from the First Billion Years (FiBY) simulation project, we argue that the intrinsic spectral energy distributions (SEDs) of star-forming galaxies at these redshifts have a self-similar shape across the mass range $8.2 \leq$ log$(M_{\star}/M_{\odot}) \leq 10.6$ probed by our sample. Using FiBY data, we construct a set of intrinsic SED templates which incorporate both detailed star formation and chemical abundance histories, and a variety of stellar population synthesis (SPS) model assumptions. With this set of intrinsic SEDs, we present a novel approach for directly recovering the shape and normalization of the dust attenuation curve. We find, across all of the intrinsic templates considered, that the average attenuation curve for star-forming galaxies at $z\simeq3.5$ is similar in shape to the commonly-adopted Calzetti starburst law, with an average total-to-selective attenuation ratio of $R_{V}=4.18\pm0.29$. We show that the optical attenuation ($A_V$) versus stellar mass ($M_{\star}$) relation predicted using our method is consistent with recent ALMA observations of galaxies at $2<z<3$ in the \emph{Hubble} \emph{Ultra} \emph{Deep} \emph{Field} (HUDF), as well as empirical $A_V - M_{\star}$ relations predicted by a Calzetti-like law. Our results, combined with other literature data, suggest that the $A_V - M_{\star}$ relation does not evolve over the redshift range $0<z<5$, at least for galaxies with log$(M_{\star}/M_{\odot}) \gtrsim 9.5$. Finally, we present tentative evidence which suggests that the attenuation curve may become steeper at log$(M_{\star}/M_{\odot}) \lesssim 9.0$.

We present results from the KMOS lensing survey-KLENS which is exploiting gravitational lensing to study the kinematics of 24 star forming galaxies at $1.4<z<3.5$ with a median mass of $\rm log(M_\star/M_\odot)=9.6$ and median star formation rate (SFR) of $\rm 7.5\,M_\odot\,yr^{-1}$. We find that 25% of these low-mass/low-SFR galaxies are rotation dominated, while the majority of our sample shows no velocity gradient. When combining our data with other surveys, we find that the fraction of rotation dominated galaxies increases with the stellar mass, and decreases for galaxies with a positive offset from the main sequence. We also investigate the evolution of the intrinsic velocity dispersion, $\sigma_0$, as a function of the redshift, $z$, and stellar mass, $\rm M_\star$, assuming galaxies in quasi-equilibrium (Toomre Q parameter equal to 1). From the $z-\sigma_0$ relation, we find that the redshift evolution of the velocity dispersion is mostly expected for massive galaxies ($\rm log(M_\star/M_\odot)>10$). We derive a $\rm M_\star-\sigma_0$ relation, using the Tully-Fisher relation, which highlights that a different evolution of the velocity dispersion is expected depending on the stellar mass, with lower velocity dispersions for lower masses, and an increase for higher masses, stronger at higher redshift. The observed velocity dispersions from this work and from comparison samples spanning $0<z<3.5$ appear to follow this relation, except at higher redshift ($z>2$), where we observe higher velocity dispersions for low masses ($\rm log(M_\star/M_\odot)\sim 9.6$) and lower velocity dispersions for high masses ($\rm log(M_\star/M_\odot)\sim 10.9$) than expected. This discrepancy could, for instance, suggest that galaxies at high-$z$ do not satisfy the stability criterion, or that the adopted parametrisation of the specific star formation rate and molecular properties fail at high redshift.

We analyse the velocity dispersion properties of 472 z~0.9 star-forming galaxies observed as part of the KMOS Redshift One Spectroscopic Survey (KROSS). The majority of this sample is rotationally dominated (83 +/- 5% with v_C/sigma_0 > 1) but also dynamically hot and highly turbulent. After correcting for beam smearing effects, the median intrinsic velocity dispersion for the final sample is sigma_0 = 43.2 +/- 0.8 km/s with a rotational velocity to dispersion ratio of v_C/sigma_0 = 2.6 +/- 0.1. To explore the relationship between velocity dispersion, stellar mass, star formation rate and redshift we combine KROSS with data from the SAMI survey (z~0.05) and an intermediate redshift MUSE sample (z~0.5). While there is, at most, a weak trend between velocity dispersion and stellar mass, at fixed mass there is a strong increase with redshift. At all redshifts, galaxies appear to follow the same weak trend of increasing velocity dispersion with star formation rate. Our results are consistent with an evolution of galaxy dynamics driven by disks that are more gas rich, and increasingly gravitationally unstable, as a function of increasing redshift. Finally, we test two analytic models that predict turbulence is driven by either gravitational instabilities or stellar feedback. Both provide an adequate description of the data, and further observations are required to rule out either model.

We use KMOS Deep Survey (KDS) galaxies, combined with results from a range of spectroscopic studies in the literature, to investigate the evolution of the stellar-mass Tully-Fisher relation since z ~ 4. We determine the slope and normalisation of the local rotation-velocity -- stellar-mass (Vc - $M_{\star}$) relationship using a reference sample of local spiral galaxies; thereafter we fix the slope, and focus on the evolution of velocity normalisation with redshift. The rotation-dominated KDS galaxies at z ~ 3.5 have rotation velocities ~ -0.1 dex lower than local reference galaxies at fixed stellar mass. By fitting 16 distant comparison samples spanning 0 < z < 3 (containing ~ 1200 galaxies), we show that the size and sign of the inferred Vc offset depends sensitively on the fraction of the parent samples used in the Tully-Fisher analysis, and how strictly the criterion of 'rotation dominated' is enforced. Confining attention to subsamples of galaxies that are especially 'disky' results in a consistent positive offset in Vc of ~ +0.1 dex, however these galaxies are not representative of the evolving-disk population at z > 1. We investigate the addition of pressure support, traced by intrinsic velocity dispersion ($\sigma_{int}$) to the KDS dynamical mass budget by adopting a 'total' effective velocity of form $V_{tot} = (Vc^{2} + 4.0\sigma_{int}^{2})^{0.5}$. The rotation-dominated and dispersion-dominated KDS galaxies fall on the same locus in the total-velocity versus stellar-mass plane, removing the need for debate over the precise selection threshold for rotation-dominated galaxies. The comparison sample offsets are in the range +0.08 to +0.15 dex in total-velocity zero-point (-0.30 to -0.55 dex in stellar-mass zero-point) from the local Tully-Fisher relation at z > 1, consistent with steady evolution of the ratio of dynamical to stellar mass with cosmic time.

We present physical properties [redshifts (z), star-formation rates (SFRs) and stellar masses (Mstar)] of bright (S850>4mJy) submm galaxies in the ~2deg2 COSMOS and UDS fields selected with SCUBA-2/JCMT. We complete the galaxy identification process for all (~2000) S/N>3.5 850-um sources, but focus our scientific analysis on a high-quality sub-sample of 651 S/N>4 sources with complete multi-wavelength coverage including 1.1-mm imaging. We check the reliability of our identifications, and the robustness of the SCUBA-2 fluxes by revisiting the recent ALMA follow-up of 29 sources in our sample. Considering >4mJy ALMA sources, our identification method has a completeness of ~86 per cent with a reliability of ~92 per cent, and only ~15-20 per cent of sources are significantly affected by multiplicity (when a secondary component contributes >1/3 of the primary source flux). The impact of source blending on the 850-um source counts as determined with SCUBA-2 is modest; scaling the single-dish fluxes by ~0.9 reproduces the ALMA source counts. For our final SCUBA-2 sample we find median z=2.40+0.10-0.04, SFR=287+-6Moyr-1, and log(Mstar/Mo)=11.12+-0.02 (the latter for 349/651 sources with optical identifications). These properties clearly locate bright submm galaxies on the high-mass end of the 'main sequence' of star-forming galaxies out to z~6, suggesting that major mergers are not a dominant driver of the high-redshift submm-selected population. Their number densities are also consistent with the evolving galaxy stellar mass function. Hence, the submm galaxy population is as expected, albeit reproducing the evolution of the main sequence of star-forming galaxies remains a challenge for theoretical models/simulations.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

We present dynamical measurements from the KMOS (K-band Multi-Object Spectrograph) Deep Survey (KDS), which is comprised of 77 typical star-forming galaxies at z $\simeq$ 3.5 in the mass range 9.0 < log(M$_{\star}$/M$_{\odot}$) < 10.5. These measurements constrain the internal dynamics, the intrinsic velocity dispersions ({\sigma}$_{int}$) and rotation velocities (V$_{C}$) of galaxies in the high redshift Universe. The mean velocity dispersion of the galaxies in our sample is {\sigma}$_{int}$ = $70.8^{+3.3}_{-3.1}$ km s$^{-1}$, revealing that the increasing average {\sigma}$_{int}$ with increasing redshift, reported for z $\lesssim2$, continues out to z $\simeq$ 3.5. Only 34 $\pm$ 8% of our galaxies are rotation-dominated (V$_{C}$/{\sigma}$_{int}$ > 1), with the sample average V$_{C}$/{\sigma}$_{int}$ value much smaller than at lower redshift. After carefully selecting comparable star-forming samples at multiple epochs, we find that the rotation-dominated fraction evolves with redshift with a z$^{-0.2}$ dependence. The rotation-dominated KDS galaxies show no clear offset from the local rotation velocity-stellar mass (i.e. V$_{C}$-M$_{\star}$) relation, although a smaller fraction of the galaxies are on the relation due to the increase in the dispersion-dominated fraction. These observations are consistent with a simple equilibrium model picture, in which random motions are boosted in high redshift galaxies by a combination of the increasing gas fractions, accretion efficiency, specific star-formation rate and stellar feedback and which may provide significant pressure support against gravity on the galactic disk scale.

We present dynamical measurements for 586 H-alpha detected star-forming galaxies from the KMOS (K-band Multi-Object Spectrograph) Redshift One Spectroscopic Survey (KROSS). The sample represents typical star-forming galaxies at this redshift (z=0.6-1.0), with a median star formation rate of ~7 Msol/yr and a stellar mass range of log[M/Msol]~9-11. We find that the rotation velocity-stellar mass relationship (the inverse of the Tully-Fisher relationship) for our rotationally-dominated sources (v/sigma>1) has a consistent slope and normalisation as that observed for z=0 disks. In contrast, the specific angular momentum (j; angular momentum divided by stellar mass), is ~0.2-0.3 dex lower on average compared to z=0 disks. The specific angular momentum scales as M^[0.6+/-0.2], consistent with that expected for dark matter (i.e., proportional to M^[2/3]). We find that z~0.9 star-forming galaxies have decreasing specific angular momentum with increasing Sersic index. Visually, the sources with the highest specific angular momentum, for a given mass, have the most disk-dominated morphologies. This implies that an angular momentum-mass-morphology relationship, similar to that observed in local massive galaxies, is already in place by z~1.

Using the exquisite depth of the Hubble Ultra Deep Field (HUDF12 programme) dataset, we explore the ongoing assembly of the outermost regions of the most massive galaxies ($\rm M_{\rm stellar}\geq$ 5$\times$10$^{10}$ M$_{\odot}$) at $z \leq$ 1. The outskirts of massive objects, particularly Early-Types Galaxies (ETGs), are expected to suffer a dramatic transformation across cosmic time due to continuous accretion of small galaxies. HUDF imaging allows us to study this process at intermediate redshifts in 6 massive galaxies, exploring the individual surface brightness profiles out to $\sim$25 effective radii. We find that 5-20\% of the total stellar mass for the galaxies in our sample is contained within 10 $< R <$ 50 kpc. These values are in close agreement with numerical simulations, and higher than those reported for local late-type galaxies ($\lesssim$5\%). The fraction of stellar mass stored in the outer envelopes/haloes of Massive Early-Type Galaxies increases with decreasing redshift, being 28.7\% at $< z > =$ 0.1, 15.1\% at $< z > =$ 0.65 and 3.5\% at $< z > =$ 2. The fraction of mass in diffuse features linked with ongoing minor merger events is $>$ 1-2\%, very similar to predictions based on observed close pair counts. Therefore, the results for our small albeit meaningful sample suggest that the size and mass growth of the most massive galaxies have been solely driven by minor and major merging from $z =$ 1 to today.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We present adaptive optics (AO) assisted SINFONI integral field unit (IFU) spectroscopy of eleven H$\alpha$ emitting galaxies selected from the High-Z Emission Line Survey (HiZELS). We obtain spatially resolved dynamics on ~kpc-scales of star-forming galaxies (stellar mass M$_\star$ = 10$^{9.5-10.5}$ M$_\odot$ and star formation rate SFR = 2-30 M$_\odot$ yr$^{-1}$) near the peak of the cosmic star-formation rate history. Combining these observations with our previous SINFONI-HiZELS campaign, we construct a sample of twenty homogeneously selected galaxies with IFU AO-aided observations -- the `SHiZELS' survey, with roughly equal number of galaxies per redshift slice, at $z$ = 0.8, 1.47, and 2.23. We measure the dynamics and identify the major kinematic axis by modelling their velocity fields to extract rotational curves and infer their inclination-corrected rotational velocities. We explore the stellar mass Tully-Fisher relationship, finding that galaxies with higher velocity dispersions tend to deviate from this relation. Using kinemetry analyses we find that galaxy interactions might be the dominant mechanism controlling the star-formation activity at $z$ = 2.23 but they become gradually less important down to $z$ = 0.8. Metallicity gradients derived from the [NII]/H$\alpha$ emission line ratio show a median negative gradient for the SHiZELS survey of $\Delta$log(O/H)/$\Delta$R = -0.026$\pm$0.008 dex kpc$^{-1}$. We find that metal-rich galaxies tend to show negative gradients, whereas metal-poor galaxies tend to exhibit positive metallicity gradients. This result suggests that the accretion of pristine gas in the periphery of galaxies plays an important role in replenishing the gas in `typical' star-forming galaxies.

We present preliminary results on on-sky test of sky subtraction methods for fiber-fed spectrograph. Using dedicated observation with FLAMES/VLT in I-band, we have tested the accuracy of the sky subtraction for 4 sky subtraction methods: mean sky, closest sky, dual stare and cross-beam switching. The cross beam-switching and dual stare method reach accuracy and precision of the sky subtraction under 1%. In contrast to the commonly held view in the literature, this result points out that fiber-fed spectrographs are adapted for the observations of faint targets.

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.

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.

Despite decades of study, we still do not fully understand why some massive galaxies abruptly switch off their star formation in the early Universe, and what causes their rapid transition to the red sequence. Post-starburst galaxies provide a rare opportunity to study this transition phase, but few have currently been spectroscopically identified at high redshift ($z>1$). In this paper we present the spectroscopic verification of a new photometric technique to identify post-starbursts in high-redshift surveys. The method classifies the broad-band optical-near--infrared spectral energy distributions (SEDs) of galaxies using three spectral shape parameters (super-colours), derived from a principal component analysis of model SEDs. When applied to the multiwavelength photometric data in the UKIDSS Ultra Deep Survey (UDS), this technique identified over 900 candidate post-starbursts at redshifts $0.5<z<2.0$. In this study we present deep optical spectroscopy for a subset of these galaxies, in order to confirm their post-starburst nature. Where a spectroscopic assessment was possible, we find the majority (19/24 galaxies; ~80 per cent) exhibit the strong Balmer absorption (H $\delta$ equivalent width $W_{\lambda}$ >5 Ang.) and Balmer break, characteristic of post-starburst galaxies. We conclude that photometric methods can be used to select large samples of recently-quenched galaxies in the distant Universe.

The KMOS Redshift One Spectroscopic Survey (KROSS) is an ESO guaranteed time survey of 795 typical star-forming galaxies in the redshift range z=0.8-1.0 with the KMOS instrument on the VLT. In this paper we present resolved kinematics and star formation rates for 584 z~1 galaxies. This constitutes the largest near-infrared Integral Field Unit survey of galaxies at z~1 to date. We demonstrate the success of our selection criteria with 90% of our targets found to be Halpha emitters, of which 81% are spatially resolved. The fraction of the resolved KROSS sample with dynamics dominated by ordered rotation is found to be 83$\pm$5%. However, when compared with local samples these are turbulent discs with high gas to baryonic mass fractions, ~35%, and the majority are consistent with being marginally unstable (Toomre Q~1). There is no strong correlation between galaxy averaged velocity dispersion and the total star formation rate, suggesting that feedback from star formation is not the origin of the elevated turbulence. We postulate that it is the ubiquity of high (likely molecular) gas fractions and the associated gravitational instabilities that drive the elevated star-formation rates in these typical z~1 galaxies, leading to the ten-fold enhanced star-formation rate density. Finally, by comparing the gas masses obtained from inverting the star-formation law with the dynamical and stellar masses, we infer an average dark matter to total mass fraction within 2.2$r_e$ (9.5kpc) of 65$\pm$12%, in agreement with the results from hydrodynamic simulations of galaxy formation.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

We present new radial velocity measurements for 82 stars, members of the Galactic globular cluster NGC 6388, obtained from ESO-VLT KMOS spectra acquired during the instrument Science Verification. The accuracy of the wavelength calibration is discussed and a number of tests of the KMOS response are presented. The cluster systemic velocity obtained (81.3 +/- 1.5 km/sec) is in very good agreement with previous determinations. While a hint of ordered rotation is found between 9'' and 20'' from the cluster centre, where the distribution of radial velocities is clearly bimodal, more data are needed before drawing any firm conclusions. The acquired sample of radial velocities has been also used to determine the cluster velocity dispersion profile between ~9'' and 70'', supplementing previous measurements at r < 2'' and r > 60'' obtained with ESO-SINFONI and ESO-FLAMES spectroscopy, respectively. The new portion of the velocity dispersion profile nicely matches the previous ones, better defining the knee of the distribution. The present work clearly shows the effectiveness of a deployable Integral Field Unit in measuring the radial velocities of individual stars for determining the velocity dispersion profile of Galactic globular clusters. It represents the pilot project for an ongoing large program with KMOS and FLAMES at the ESO-VLT, aimed at determining the next generation of velocity dispersion and rotation profiles for a representative sample of globular clusters.

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

Submillimetre (submm) galaxies are among the most rapidly star-forming and most massive high-redshift galaxies; thus, their properties provide important constraints on galaxy evolution models. However, there is still a debate about their stellar masses and their nature in the context of the general galaxy population. To test the reliability of their stellar mass determinations, we used a sample of simulated submm galaxies for which we derived stellar masses via spectral energy distribution (SED) modelling (with Grasil, Magphys, Hyperz and LePhare) adopting various star formation histories (SFHs). We found that the assumption of SFHs with two independent components leads to the most accurate stellar masses. Exponentially declining SFHs (tau) lead to lower masses (albeit still consistent with the true values), while the assumption of single-burst SFHs results in a significant mass underestimation. Thus, we conclude that studies based on the higher masses inferred from fitting the SEDs of real submm galaxies with double SFHs are most likely to be correct, implying that submm galaxies lie on the high-mass end of the main sequence of star-forming galaxies. This conclusion appears robust to assumptions of whether or not submm galaxies are driven by major mergers, since the suite of simulated galaxies modelled here contains examples of both merging and isolated galaxies. We identified discrepancies between the true and inferred stellar ages (rather than the dust attenuation) as the primary determinant of the success/failure of the mass recovery. Regardless of the choice of SFH, the SED-derived stellar masses exhibit a factor of ~2 scatter around the true value; this scatter is an inherent limitation of the SED modelling due to simplified assumptions. Finally, we found that the contribution of active galactic nuclei does not have any significant impact on the derived stellar masses.

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.

This paper presents the latest optical design for the MOONS triple-arm spectrographs. MOONS will be a Multi-Object Optical and Near-infrared Spectrograph and will be installed on one of the European Southern Observatory (ESO) Very Large Telescopes (VLT). Included in this paper is a trade-off analysis of different types of collimators, cameras, dichroics and filters.

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.

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.

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.

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.

We present a new method to classify the broad band optical-NIR spectral energy distributions (SEDs) of galaxies using three shape parameters (super-colours) based on a Principal Component Analysis of model SEDs. As well as providing a compact representation of the wide variety of SED shapes, the method allows for easy visualisation of information loss and biases caused by the incomplete sampling of the rest-frame SED as a function of redshift. We apply the method to galaxies in the UKIDSS Ultra Deep Survey with 0.9<z<1.2, and confirm our classifications by stacking rest-frame optical spectra for a fraction of objects in each class. As well as cleanly separating a tight red-sequence from star-forming galaxies, three unusual populations are identifiable by their unique colours: very dusty star-forming galaxies with high metallicity and old mean stellar age; post-starburst galaxies which have formed greater than around 10% of their mass in a recent unsustained starburst event; and metal-poor quiescent dwarf galaxies. We find that quiescent galaxies account for 45% of galaxies with log(M*/Msol)>11, declining steadily to 13% at log(M*/Msol)=10. The properties and mass-function of the post-starburst galaxies are consistent with a scenario in which gas-rich mergers contribute to the growth of the low and intermediate mass range of the red sequence.

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

KMOS is a multi-object near-infrared integral field spectrometer with 24 deployable cryogenic pick-off arms. Inevitably, data processing is a complex task that requires careful calibration and quality control. In this paper we describe all the steps involved in producing science-quality data products from the raw observations. In particular, we focus on the following issues: (i) the calibration scheme which produces maps of the spatial and spectral locations of all illuminated pixels on the detectors; (ii) our concept of minimising the number of interpolations, to the limiting case of a single reconstruction that simultaneously uses raw data from multiple exposures; (iii) a comparison of the various interpolation methods implemented, and an assessment of the performance of true 3D interpolation schemes; (iv) the way in which instrumental flexure is measured and compensated. We finish by presenting some examples of data processed using the pipeline.

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

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.

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.

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.

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.

We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope.

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

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

We use the UKIDSS Ultra-Deep Survey, the deepest degree-scale near-infrared survey to date, to investigate the clustering of star-forming and passive galaxies to z ~ 3.5. Our new measurements include the first determination of the clustering for passive galaxies at z > 2, which we achieve using a cross-correlation technique. We find that passive galaxies are the most strongly clustered, typically hosted by massive dark matter halos with M_halo > 5 x 10^12 M_sun irrespective of redshift or stellar mass. Our findings are consistent with models in which a critical halo mass determines the transition from star-forming to passive galaxies. Star-forming galaxies show no strong correlation between stellar mass and halo mass, but passive galaxies show evidence for an anti-correlation; low-mass passive galaxies appear, on average, to be located in the most massive halos. These results can be understood if the termination of star formation is most efficient for galaxies of low stellar mass in very dense environments.

Fibre-fed spectrographs now have throughputs equivalent to slit spectrographs. However, the sky subtraction accuracy that can be reached has often been pinpointed as one of the major issues associated with the use of fibres. Using technical time observations with FLAMES-GIRAFFE, two observing techniques, namely dual staring and cross beam-switching, were tested and the resulting sky subtraction accuracy reached in both cases was quantified. Results indicate that an accuracy of 0.6% on sky subtraction can be reached, provided that the cross beam-switching mode is used. This is very encouraging with regard to the detection of very faint sources with future fibre-fed spectrographs, such as VLT/MOONS or E-ELT/MOSAIC.

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.

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.

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of \alpha = -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, \alpha = -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 \HST\ program aimed at extending previous WFC3/IR observations of the UDF by quadrupling the exposure time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized the universe; our observations are designed to provide a robust determination of the star formation density at $z$$\,\gtrsim\,$8, improve measurements of the ultraviolet continuum slope at $z$$\,\sim\,7\,-\,$8, facilitate the construction of new samples of $z$$\,\sim\,9\,-\,$10 candidates, and enable the detection of sources up to $z$$\,\sim\,$12. For this project we committed to combining these and other WFC3/IR imaging observations of the UDF area into a single homogeneous dataset, to provide the deepest near-infrared observations of the sky currently achievable. In this paper we present the observational overview of the project, motivated by its scientific goals, and describe the procedures used in reducing the data as well as the final products that are produced. We have used the most up up-to-date methods for calibrating and combining the images, in particular paying attention to correcting several instrumental effects. We release the full combined mosaics, comprising a single, unified set of mosaics of the UDF, providing the deepest near-infrared blank-field view of the universe obtained to date, reaching magnitudes as deep as AB$\,\sim\,$30 in the near-infrared, and yielding a legacy dataset on this field of lasting scientific value to the community.

We present the results of the deepest search to date for star-forming galaxies beyond a redshift z~8.5 utilizing a new sequence of near-infrared Wide Field Camera 3 images of the Hubble Ultra Deep Field. This `UDF12' campaign completed in September 2012 doubles the earlier exposures with WFC3/IR in this field and quadruples the exposure in the key F105W filter used to locate such distant galaxies. Combined with additional imaging in the F140W filter, the fidelity of high redshift candidates is greatly improved. Using spectral energy distribution fitting techniques on objects selected from a deep multi-band near-infrared stack we find 7 promising z>8.5 candidates. As none of the previously claimed UDF candidates with 8.5<z<10 is confirmed by our deeper multi-band imaging, our campaign has transformed the measured abundance of galaxies in this redshift range. Although we recover the candidate UDFj-39546284 (previously proposed at z=10.3), it is undetected in the newly added F140W image, implying it lies at z=11.9 or is an intense emission line galaxy at z~2.4. Although no physically-plausible model can explain the required line intensity given the lack of Lyman alpha or broad-band UV signal, without an infrared spectrum we cannot rule out an exotic interloper. Regardless, our robust z ~ 8.5 - 10 sample demonstrates a luminosity density that continues the smooth decline observed over 6 < z < 8. Such continuity has important implications for models of cosmic reionization and future searches for z>10 galaxies with JWST.

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

This paper presents new deep and wide narrow-band surveys undertaken with UKIRT, Subaru and the VLT; a combined effort to select large, robust samples of H-alpha (Ha) emitters at z=0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyrs) in a uniform manner over ~2 deg^2 in the COSMOS and UDS fields. The deep Ha surveys reach ~3M_sun/yr out to z=2.2 for the first time, while the wide area and the coverage over two independent fields allow to greatly overcome cosmic variance. A total of 1742, 637, 515 and 807 Ha emitters are homogeneously selected at z=0.40, 0.84, 1.47 and 2.23, respectively, and used to determine the Ha luminosity function and its evolution. The faint-end slope is found to be -1.60+-0.08 over z=0-2.23, showing no evolution. The characteristic luminosity of SF galaxies, L*, evolves significantly as log[L*(z)]=0.45z+log[L*(z=0)]. This is the first time Ha has been used to trace SF activity with a single homogeneous survey at z=0.4-2.23. Overall, the evolution seen in Ha is in good agreement with the evolution seen using inhomogeneous compilations of other tracers of star formation, such as FIR and UV, jointly pointing towards the bulk of the evolution in the last 11 Gyrs being driven by a strong luminosity increase from z~0 to z~2.2. Our uniform analysis reveals an Ha star formation history of the Universe which can be simply parameterised by log(SFRD)=-2.1/(1+z) for z<2.23. Both the shape and normalisation of the Ha star formation history are consistent with the measurements of the stellar mass density growth, confirming that our Ha analysis traces the bulk of the formation of stars in the Universe up to z~2.2. Star formation activity since z~2.2 is responsible for ~95% of the total stellar mass density observed locally today (with half being assembled between z=1.2-2.2, and the rest since z<1.2).

MOONS is a new conceptual design for a Multi-Object Optical and Near-infrared Spectrograph for the Very Large Telescope (VLT), selected by ESO for a Phase A study. The baseline design consists of 1000 fibers deployable over a field of view of 500 square arcmin, the largest patrol field offered by the Nasmyth focus at the VLT. The total wavelength coverage is 0.8um-1.8um and two resolution modes: medium resolution and high resolution. In the medium resolution mode (R=4,000-6,000) the entire wavelength range 0.8um-1.8um is observed simultaneously, while the high resolution mode covers simultaneously three selected spectral regions: one around the CaII triplet (at R=8,000) to measure radial velocities, and two regions at R=20,000 one in the J-band and one in the H-band, for detailed measurements of chemical abundances. The grasp of the 8.2m Very Large Telescope (VLT) combined with the large multiplex and wavelength coverage of MOONS - extending into the near-IR - will provide the observational power necessary to study galaxy formation and evolution over the entire history of the Universe, from our Milky Way, through the redshift desert and up to the epoch of re-ionization at z>8-9. At the same time, the high spectral resolution mode will allow astronomers to study chemical abundances of stars in our Galaxy, in particular in the highly obscured regions of the Bulge, and provide the necessary follow-up of the Gaia mission. Such characteristics and versatility make MOONS the long-awaited workhorse near-IR MOS for the VLT, which will perfectly complement optical spectroscopy performed by FLAMES and VIMOS.

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.

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

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.

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.

Establishing the stellar masses (M*), and hence specific star-formation rates (sSFRs) of submillimetre galaxies (SMGs) is crucial for determining their role in the cosmic galaxy/star formation. However, there is as yet no consensus over the typical M* of SMGs. Specifically, even for the same set of SMGs, the reported average M* have ranged over an order of magnitude, from ~5x10^10 Mo to ~5x10^11 Mo. Here we study how different methods of analysis can lead to such widely varying results. We find that, contrary to recent claims in the literature, potential contamination of IRAC 3-8 um photometry from hot dust associated with an active nucleus is not the origin of the published discrepancies in derived M*. Instead, we expose in detail how inferred M* depends on assumptions made in the photometric fitting, and quantify the individual and cumulative effects of different choices of initial mass function, different brands of evolutionary synthesis models, and different forms of assumed star-formation history. We review current observational evidence for and against these alternatives as well as clues from the hydrodynamical simulations, and conclude that, for the most justifiable choices of these model inputs, the average M* of SMGs is ~2x10^11 Mo. We also confirm that this number is perfectly reasonable in the light of the latest measurements of their dynamical masses, and the evolving M* function of the overall galaxy population. M* of this order imply that the average sSFR of SMGs is comparable to that of other star-forming galaxies at z>2, at 2-3 Gyr^-1. This supports the view that, while rare outliers may be found at any M*, most SMGs simply form the top end of the main-sequence of star-forming galaxies at these redshifts. Conversely, this argues strongly against the viewpoint that SMGs are extreme pathological objects, of little relevance in the cosmic history of star-formation.

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

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.

We report the detection of an extremely bright ($\sim$37 mJy at 1100 $\mu$m and $\sim$91 mJy at 880 $\mu$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $\mu$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $\mu$m and SMA 880 $\mu$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $\sigma$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

This paper presents the results from the first wide and deep dual narrow-band survey to select H-alpha (Ha) and [OII] line emitters at z=1.47+-0.02 (using matched narrow-band filters in the H and z' bands), exploiting synergies between the UKIRT and Subaru telescopes. The Ha survey at z=1.47 reaches a flux limit of ~7x10^-17 erg/s/cm^2 and detects ~200 Ha emitters over 0.7deg^2, while the much deeper [OII] survey reaches an effective flux of ~7x10^-18 erg/s/cm^2, detecting ~1400 z=1.47 [OII] emitters in a matched co-moving volume of ~2.5x10^5 Mpc^3. The combined survey results in the identification of 190 simultaneous Ha and [OII] emitters at z=1.47. Ha and [OII] luminosity functions are derived and both are shown to evolve significantly from z~0 in a consistent way. The star formation rate density of the Universe at z=1.47 is evaluated, with the Ha analysis yielding 0.16+-0.05 M_sun/yr/Mpc^3 and the [OII] analysis 0.17+-0.04 M_sun/yr/Mpc^3. The measurements are combined with other studies, providing a self-consistent measurement of the star formation history of the Universe over the last ~11Gyrs. By using a large comparison sample at z~0.1 (from the SDSS), [OII]/Ha line ratios are calibrated as probes of dust-extinction. Ha emitters at z~1.47 show on average 1 mag of extinction at Ha, similar to the SDSS sources at z~0. Although we find that dust extinction correlates with SFR, the relation evolves by about ~0.5 mag from z~1.5 to z~0, with z~0 relations over-predicting the dust extinction corrections at high-z by that amount. Stellar mass is found to be a much more fundamental extinction predictor, with the relation between mass and extinction being valid at both z~0 and z~1.5. Dust extinction corrections as a function of optical colours are also derived, offering simpler mechanisms for estimating extinction in moderately star-forming systems over the last ~9Gyrs [Abridged].

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

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.

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.

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.

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.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

(Abridged) This paper presents an environment and stellar mass study of a large sample of star-forming (SF) galaxies at z=0.84 from the HiZELS survey, over 1.3 deg^2 in the COSMOS and UKIDSS UDS fields. By taking advantage of a truly panoramic coverage, from the field to a rich cluster, it is shown that both mass and environment play crucial roles in determining the properties of SF galaxies. The median specific SFR declines with mass in all environments, and the fraction of galaxies forming stars declines from ~40%, for M~10^10M_sun to effectively zero at M>10^11.5M_sun, confirming that mass-downsizing is generally in place by z~1. The fraction of SF galaxies also falls as a function of local environmental density from ~40% in the field to approaching zero at rich group/cluster densities. When SF does occur in high density regions, it is merger-dominated and, if only non-merging SF galaxies are considered, then the environment and mass trends are even stronger and largely independent, as in the local Universe. The median SFR of SF galaxies is found to increase with density up to intermediate (group or cluster outskirts) densities; this is clearly seen as a change in the faint-end slope of the H-alpha LF from steep (-1.9), in poor fields, to shallow (-1.1) in groups and clusters. Interestingly, the relation between median SFR and environment is only found for low to moderate-mass galaxies (below ~10^10.6M_sun), and is not seen for massive SF galaxies. Overall, these observations provide a detailed view over a sufficiently large range of mass and environment to reconcile previous observational claims: mass is the primary predictor of SF activity at z~1, but the environment, while enhancing the median SFR of (lower-mass) SF galaxies, is ultimately responsible for suppressing SF activity in all galaxies above surface densities of 10-30 Mpc^-2 (groups and clusters).

We use the UKIDSS Ultra-Deep Survey to trace the evolution of galaxy clustering to z = 3. Using photometric redshifts derived from data covering the wavelength range 0.3 - 4.5 um we examine this clustering as a function of absolute K-band luminosity, colour and star-formation rate. Comparing the deprojected clustering amplitudes, we find that red galaxies are more strongly clustered than blue galaxies out to at least z = 1.5, irrespective of rest-frame K-band luminosity. We then construct passive and star-forming samples based on stellar age, colour and star-formation histories calculated from the best fitting templates. The clustering strength of star-forming galaxies declines steadily from r_0 ~ 7 h^-1 Mpc at z ~ 2 to r_0 ~ 3 h^-1 Mpc at z ~ 0, while passive galaxies have clustering strengths up to a factor of two higher. Within the passive and star-forming subsamples, however, we find very little dependence of galaxy clustering on K-band luminosity. Galaxy `passivity' appears to be the strongest indicator of clustering strength. We compare these clustering measurements with those predicted for dark matter halos and conclude that passive galaxies typically reside in halos of mass M > 10^13 M_sun while luminous star-forming galaxies occupy halos an order of magnitude less massive over the range 0.5 < z < 1.5. The decline in the clustering strength of star-forming galaxies with decreasing redshift indicates a decline in the hosting halo mass for galaxies of a given luminosity. We find evidence for convergence of clustering in star-forming and passive galaxies around z ~ 2, which is consistent with this being the epoch at which the red sequence of galaxies becomes distinct.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

In these proceedings we report on HiZELS, the High-z Emission Line Survey, our successful panoramic narrow-band Campaign Survey using WFCAM on UKIRT to detect and study emission line galaxies at z~1-9. HiZELS employs the H2(S1) narrow-band filter together with custom-made narrow-band filters in the J and H-bands, with the primary aim of delivering large, identically-selected samples of H-alpha emitting galaxies at redshifts of 0.84, 1.47 and 2.23. Comparisons between the luminosity function, the host galaxy properties, the clustering, and the variation with environment of these H-alpha-selected samples are yielding unique constraints on the nature and evolution of star-forming galaxies, across the peak epoch of star-formation activity in the Universe. We provide a summary of the project status, and detail the main scientific results obtained so far: the measurement of the evolution of the cosmic star-formation rate density out to z > 2 using a single star-formation indicator, determination of the morphologies, environments and dust-content of the star-forming galaxies, and a detailed investigation of the evolution of their clustering properties. We also summarise the on-going work and future goals of the project.

(Abridged) The clustering properties of a well-defined sample of 734 H-alpha emitters at z=0.84 obtained as part of the Hi-z Emission Line Survey (HiZELS) are investigated. The spatial correlation function is very well-described by (r/r_0)^-1.8, with r_0=2.7+-0.3Mpc/h. The correlation length r_0 increases strongly with H-alpha luminosity, L_H-alpha, from r_0~2Mpc/h for the most quiescent galaxies (star-formation rates of ~4M_sun/yr), up to r_0>5Mpc/h for the brightest galaxies in H-alpha. The correlation length also increases with increasing rest-frame K-band luminosity (M_K), but the r_0-L_H-alpha correlation maintains its full statistical significance at fixed M_K. At z=0.84, star-forming galaxies classified as irregulars or mergers are much more clustered than discs and non-mergers, but once the samples are matched in L_H-alpha and M_K, the differences vanish, implying that the clustering is independent of morphological type at z~1. The typical H-alpha emitters found at z=0.84 reside in dark-matter haloes of ~10^12M_sun, but those with the highest SFRs reside in more massive haloes of ~10^13M_sun. Comparing the results with those of H-alpha surveys at different redshifts, it is seen that although the break of the H-alpha luminosity function, L*, evolves by a factor of ~30 from z=0.24 to z=2.23, galaxies with the same L_H-alpha/L*(z) are found in dark matter haloes of similar masses, independently of cosmic time. This not only confirms that star-formation is more efficient at higher redshift, but also suggests a fundamental connection between the strong decrease of L* since z=2.23 and the quenching of star-formation in galaxies residing within dark-matter haloes significantly more massive than 10^12M_sun at any given epoch.

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.

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.

We present the largest area survey to date (1.4 deg2) for Lyman-alpha emitters (LAEs) at z~9, as part of the Hi-z Emission Line Survey (HiZELS). The survey, which primarily targets H-alpha emitters at z < 3, uses the Wide Field CAMera on the United Kingdom Infrared Telescope and a custom narrow-band filter in the J band to reach a Lyman-alpha luminosity limit of ~10^43.8 erg/s over a co-moving volume of 1.12x10^6 Mpc^3 at z = 8.96+-0.06. Two candidates were found out of 1517 line emitters, but those were rejected as LAEs after follow-up observations. This improves the limit on the space density of bright Lyman-alpha emitters by 3 orders of magnitude and is consistent with suppression of the bright end of the Lyman-alpha luminosity function beyond z~6. Combined with upper limits from smaller but deeper surveys, this rules out some of the most extreme models for high-redshift Lyman-alpha emitters. The potential contamination of narrow-band Lyman-alpha surveys at z>7 by Galactic brown dwarf stars is also examined, leading to the conclusion that such contamination may well be significant for searches at 7.7 < z < 8.0, 9.1 < z < 9.5 and 11.7 < z < 12.2.

[Abridged] We compare H-alpha, ultraviolet (UV) and infrared (IR) indicators of star formation rate (SFR) for a sample of z = 0.84 galaxies from the High-z Emission Line Survey (HiZELS). Using multi-wavelength diagnostics, we estimate that 5 - 11 per cent of H-alpha emitters at this redshift are active galactic nuclei. We detect 35 per cent of the H-alpha emitters individually at 24 microns, and stack the star-forming emitters on deep 24-micron images in order to calculate the typical SFRs of our galaxies. These are compared to the observed H-alpha line fluxes in order to estimate the extinction at z = 0.84, and we find a significant increase in dust extinction for galaxies with higher SFRs. We demonstrate that the local relationship between SFR and extinction is suitable for our sample, and attribute the overall increase in dust extinction for z = 0.84 galaxies to an increase in the average SFR, rather than to a change in dust properties. We calculate the UV extinction, and find a similar dependence on SFR, but no evidence for a 2175 angstrom UV bump in the dust attenuation law. We calculate the conversion between the dust attenuation of nebular and stellar radiation, gamma, and show that gamma = 0.50 +/- 0.14. The extinction / SFR relationship is shown to be applicable to galaxies with a range of morphologies and bulge-to-disk ratios, to both merging and non-merging galaxies, and to galaxies within high- and low-density environments, implying that it is a fundamental property of star-forming regions. In order to allow future studies to correct for a SFR-dependent amount of dust extinction, we present an equation to predict the extinction of a galaxy, based solely upon its observed H-alpha luminosity, and use this to recalculate the H-alpha luminosity function and star formation rate density at z = 0.84.

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.

We report the redshift of a distant, highly obscured submm galaxy (SMG), based entirely on the detection of its CO line emission. We have used the newly commissioned Eight-MIxer Receiver (EMIR) at the IRAM 30m telescope, with its 8 GHz of instantaneous dual-polarization bandwidth, to search the 3-mm atmospheric window for CO emission from SMMJ14009+0252, a bright SMG detected in the SCUBA Lens Survey. A detection of the CO(3--2) line in the 3-mm window was confirmed via observations of CO(5--4) in the 2-mm window. Both lines constrain the redshift of SMMJ14009+0252 to z=2.9344, with high precision (dz=2 10^{-4}). Such observations will become routine in determining redshifts in the era of the Atacama Large Millimeter/submillimeter Array (ALMA).

New results are presented, as part of the Hi-z Emission Line Survey (HiZELS), from the largest area survey to date (1.4 sq.deg) for Lyman-alpha emitters (LAEs) at z~9. The survey, which is primarily targeting H-alpha emitters at z<3, uses the Wide Field CAMera on the United Kingdom Infrared Telescope and a custom narrow-band filter in the J band and reaches a Lyman-alpha luminosity limit of ~10^43.8 erg/s over a co-moving volume of 1.12x10^6 Mpc^3 at z=8.96+-0.06. Only 2 candidates were found out of 1517 line emitters and those were rejected as LAEs after follow-up observations. The limit on the space density of bright LAEs is improved by 3 orders of magnitude, consistent with suppression of the bright end of the Lyman-alpha luminosity function beyond z~6. Combined with upper limits from smaller but deeper surveys, this rules out some of the most extreme models for high-redshift LAEs. The potential contamination of future narrow-band Lyman-alpha surveys at z>7 by Galactic brown dwarf stars is also examined, leading to the conclusion that such contamination may well be significant for searches at 7.7<z<8.0, 9.1<z<9.5 and 11.7 < z < 12.2.

New results from a large survey of H-alpha emission-line galaxies at z=0.84 using WFCAM/UKIRT and a custom narrow-band filter in the J band are presented as part of the HiZELS survey. Reaching an effective flux limit of 1e-16 erg/s/cm^2 in a comoving volume of 1.8e5 Mpc^3, this represents the largest and deepest survey of its kind ever done at z~1. There are 1517 potential line emitters detected across 1.4 sq.deg of the COSMOS and UKIDSS UDS fields, of which 743 are selected as H-alpha emitters. These are used to calculate the H-alpha luminosity function, which is well-fitted by a Schechter function with phi*=10^(-1.92+-0.10) Mpc^-3, L*=10^(42.26+-0.05)erg/s, and alpha=-1.65+-0.15. The integrated star formation rate density (SFRD) at z=0.845 is 0.15+-0.01 M_sun/yr/Mpc^3. The results robustly confirm a strong evolution of SFRD from the present day out to z~1 and then flattening to z~2, using a single star-formation indicator. Out to z~1, both the characteristic luminosity and space density of the H-alpha emitters increase significantly; at higher redshifts, L* continues to increase, but phi* decreases. The z=0.84 H-alpha emitters are mostly disk galaxies (82+-3%), while 28+-4% of the sample show signs of merger activity and contribute ~20% to the total SFRD. Irregulars and mergers dominate the H-alpha luminosity function above L*, while disks are dominant at fainter luminosities. These results demonstrate that it is the evolution of 'normal' disk galaxies that drives the strong increase in the SFRD from the current epoch to z~1, although the continued strong evolution of L* beyond z=1 suggests an increasing importance of merger activity at higher redshifts.

We have studied the uJy radio properties of K-selected galaxies detected in the Ultra-Deep Survey portion of UKIDSS using 610- and 1,400-MHz images from the VLA and GMRT. These deep radio mosaics, combined with the largest and deepest K image currently available, allow high-S/N detections of many K-selected sub-populations, including sBzK and pBzK star-forming and passive galaxies. We find a strong correlation between the radio and K fluxes and a linear relationship between SFR and K luminosity. We find no evidence, from either radio spectral indices or a comparison with submm-derived SFRs, that the full sample is strongly contaminated by AGN. The sBzK and pBzK galaxies have similar levels of radio flux, SFR and specific SFR (SSFR) at z < 1.4, suggesting there is strong contamination of the pBzK sample by star-forming galaxies. At z > 1.4, pBzK galaxies become difficult to detect in the radio stack, though the implied SFRs are still much higher than expected for passively evolving galaxies. Their radio emission may come from low-luminosity AGN. EROs straddle the passive and star-forming regions of the BzK diagram and also straddle the two groups in terms of their radio properties. K-bright ERO samples are dominated by passive galaxies and faint ERO samples contain more star-forming galaxies. The star-formation history (SFH) from stacking all K sources in the UDS agrees well with that derived for other wavebands and other radio surveys, at least out to z ~ 2. The radio-derived SFH then appears to fall more steeply than that measured at other wavelengths. The SSFR for K-selected sources rises strongly with redshift at all stellar masses, and shows a weak dependence on stellar mass. High- and low-mass galaxies show a similar decline in SSFR since z ~ 2 (abridged).

We present the results of a study of a large sample of luminous (z'{AB}<26) Lyman break galaxies (LBGs) in the redshift interval 4.7<z<6.3, selected from a contiguous 0.63 square degree area covered by the UKIDSS Ultra Deep Survey (UDS) and the Subaru XMM-Newton Survey (SXDS). Utilising the large area coverage and the excellent available optical+nearIR data, we use a photometric redshift analysis to derive a new, robust, measurement of the bright end (L>L*) of the UV-selected luminosity function at high redshift. When combined with literature studies of the fainter LBG population, our new sample provides improved constraints on the luminosity function of redshift 5<z<6 LBGs over the luminosity range 0.1L*<L<10L*. A maximum likelihood analysis returns best-fitting Schechter function parameters of M*_1500=-20.73, phi*=0.0009 Mpc^-3 and alpha=-1.66 for the luminosity function at z=5, and M*_1500 = -20.04, phi*=0.0018 Mpc^-3 and alpha=-1.71 at z=6. In addition, an analysis of the angular clustering properties of our LBG sample demonstrates that luminous 5<z<6 LBGs are strongly clustered (r_0 = 8.1 Mpc), and are consistent with the occupation of dark matter halos with masses of ~10^{11.5-12.0} Msun. Moreover, by stacking the available multi-wavelength imaging data for the high-redshift LBGs it is possible to place useful constraints on their typical stellar mass. The results of this analysis suggest that luminous LBGs at 5<z<6 have an average stellar mass of ~10^10 Msun, consistent with the results of the clustering analysis assuming plausible values for the ratio of stellar to dark matter. Finally, by combining our luminosity function results with those of the stacking analysis we derive estimates of ~1x10^7 Msun Mpc^-3 and 4x10^6 Msun Mpc^-3 for the stellar mass density at z~5 and z~6 respectively.

We report the discovery of a giant Lya emitter (LAE) with a Spitzer/IRAC counterpart near the reionization epoch at z=6.595. The giant LAE is found from the extensive 1 deg^2 Subaru narrow-band survey for z=6.6 LAEs in the Subaru/XMM-Newton Deep Survey (SXDS) field, and subsequently identified by deep spectroscopy of Keck/DEIMOS and Magellan/IMACS. Among our 207 LAE candidates, this LAE is not only the brightest narrow-band object with L(Lya) = 3.9+/-0.2 x 10^43 erg/s in our survey volume of 10^6 Mpc^3, but also a spatially extended Lya nebula with the largest isophotal area whose major axis is at least ~3''. This object is more likely to be a large Lya nebula with a size of >~17 kpc than to be a strongly-lensed galaxy by a foreground object. Our Keck spectrum with medium-high spectral and spatial resolutions suggests that the velocity width is v(FWHM)=251+/-21 km/s, and that the line-center velocity changes by ~60 km/s in a 10-kpc range. The stellar mass and star-formation rate are estimated to be 0.9-5.0 x 10^10 Mo and >34 Mo/yr, respectively, from the combination of deep optical to infrared images of Subaru, UKIDSS-Ultra Deep Survey, and Spitzer/IRAC. Although the nature of this object is not yet clearly understood, this could be an important object for studying cooling clouds accreting onto a massive halo, or forming-massive galaxies with significant outflows contributing to cosmic reionization and metal enrichment of inter-galactic medium.

We review our current knowledge of the population of high-redshift sub-mm/mm galaxies, with particular emphasis on recent results from the SCUBA HAlf Degree Extragalactic Survey (SHADES). All available evidence indicates that these objects form the high-redshift, high-luminosity, high-mass tail of the dusty starforming galaxy population revealed at lower redshifts and luminosities by Spitzer. Current theoretical models of galaxy formation struggle to reproduce these extreme objects in the numbers indicated by current surveys.

We use the UKIDSS Ultra-deep survey (UDS), currently the deepest panoramic near infra-red survey, together with deep Subaru optical imaging to measure the clustering, number counts and luminosity function of galaxies at $z\sim 2$ selected using the BzK selection technique. We find that both star-forming (sBzK) and passive (pBzK) galaxies, to a magnitude limit of $K_{AB} < 23$, are strongly clustered. The passive galaxies are the most strongly clustered population, with scale lengths of $r_0 = 15.0^{+1.9}_{-2.2}$h$^{-1}$Mpc compared with $r_0 = 6.75^{+0.34}_{-0.37}$h$^{-1}$Mpc for star-forming galaxies. The direct implication is that passive galaxies inhabit the most massive dark-matter halos, and are thus identified as the progenitors of the most massive galaxies at the present day. In addition, the pBzKs exhibit a sharp flattening and potential turn-over in their number counts, in agreement with other recent studies. This plateau cannot be explained by the effects of incompleteness. We conclude that only very massive galaxies are undergoing passive evolution at this early epoch, consistent with the downsizing scenario for galaxy evolution. Assuming a purely passive evolution for the pBzKs from their median redshift to the present day, their luminosity function suggests that only $\sim 2.5 %$ of present day massive ellipticals had a pBzK as a main progenitor.

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.

We present a comparison between the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 and 850 microns in the Lockman Hole East with a deep Spitzer Space Telescope survey at 3.6-24 microns conducted in Guaranteed Time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer-identified submm galaxies at 1<z<1.5 appear to reside in overdensities when the density field is smoothed at 0.5-2 Mpc comoving diameters, supporting the high-redshift reversal of the local star formation - galaxy density relation. We derive the dust-shrouded cosmic star formation history of galaxies as a function of assembled stellar masses. For model stellar masses <10^11 Msun, this peaks at lower redshifts than the ostensible z~2.2 maximum for submm point sources, adding to the growing consensus for ``downsizing'' in star formation. Our surveys are also consistent with ``downsizing'' in mass assembly. Both the mean star formation rates <dM/dt> and specific star formation rates <(1/M)dM/dt> are in striking disagreement with some semi-analytic predictions from the Millenium simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field.

We have analysed the 24um properties of a radio-selected sample in the Subaru-XMM/Newton Deep Field in order to explore the behaviour of the FIR/radio relation at high redshifts. Statistically, the correlation is described by q24, the ratio between the observed flux densities at 24um and 1.4GHz, respectively. Using 24um data results in considerably more scatter in the correlation than previous work using data at 60-70um. Nevertheless, we do observe a steady correlation as a function of redshift, up to z~3.5, suggesting its validity back to primeval times. We find q24 = 0.30 +/- 0.56 for the observed and q24 = 0.71 +/- 0.47 for the k-corrected radio sample, based on sources with 300uJy < S(1.4GHz) < 3.2mJy and 24um detections. A suitable k-correction given by a M82-like mid-IR template suggests no extreme silicate absorption in the bulk of our radio sample. Using thresholds in q24 to identify radio-excess sources, we have been able to characterise the transition from radio-loud AGN to star-forming galaxies and radio-quiet AGN at faint (<1mJy) radio flux densities. Our results are in broad agreement with previous studies which show a dominant radio-loud AGN population at >1mJy. The rest-frame U-B colours of the expected radio-excess population have redder distribution than those that follow the correlation. This is therefore a promising way to select obscured Type-2 AGN, with a radio loud nature, missed by deep X-ray observations. Spectroscopic follow-up of these sources is required to fully test this method.

We present estimates of the photometric redshifts, stellar masses and star formation histories of sources in the SCUBA HAlf Degree Extragalactic Survey (SHADES). This paper describes the 60 SCUBA sources detected in the Lockman Hole covering an area of ~320 square arcmin. Using photometry spanning the B band to 8um, we find that the average SCUBA source forms a significant fraction of its stars in an early period of star formation and that most of the remainder forms in a shorter more intense burst around the redshift it is observed. This trend does not vary significantly with source redshift but the exact ratio of early to late mass is quite sensitive to the way extinction is treated in the modelling. However, the sources show a clear increase in stellar mass with redshift, consistent with downsizing. In terms of SED types, only two out of the 51 sources we have obtained photometric redshifts for are best fit by a quasar-like spectral energy distribution, with approximately 80 per cent of the sources being best fit with late-type spectra (Sc, Im and starburst). By including photometry at 850um, we conclude that the average SCUBA source is forming stars at a rate somewhere between 6 and 30 times the rate implied from the rest-frame optical in a dust obscured burst and that this burst creates 15-65 per cent of the total stellar mass. Using a simplistic calculation, we estimate from the average star formation history that between one in five and one in 15 bright (L_* +2 mag < L_optical < L_* -1 mag) galaxies in the field over the interval 0 < z < 3 will at some point in their lifetime experience a similar energetic dusty burst of star formation. Finally, we compute the evolution of the star formation rate density and find it peaks around z=2.

We constrain the mean kinetic efficiency of radio-loud active galactic nuclei by using an optically selected sample for which both the optical and the radio luminosity functions (LFs) have been determined; the former traces the bolometric luminosity L, while the latter traces the kinetic power L_k, empirically correlated to the radio emission. Thus in terms of the ratio g_k=L_k/L, we can convert the optical LF of the sample into a radio one. This computed LF is shown to match the directly observed LF for the same sample if g_k=0.10^{+0.05}_{-0.01} holds, with a scatter \sigma=0.38^{+0.04}_{-0.09} dex; with these values we also match a number of independent correlations between L_k, L and radio emission, that we derive through Monte Carlo simulations. We proceed to translate the value of g_k into a constraint on the kinetic efficiency for the production of radio jets or winds, namely, \epsilon_k=L_k/(Mdot*c^2)~0.01 in terms of the rate Mdot of mass accretion onto the central black hole. Then, on assuming that on average the radio sources share the same kinetic efficiency, we compute a solid lower limit of about 25% on the contribution of radio sources to the local black hole mass density.

This paper investigates the clustering properties of a complete sample of 1041 24um-selected sources brighter than F[24um]=400 uJy in the overlapping region between the SWIRE and UKIDSS UDS surveys. We have concentrated on the two (photometric) interval ranges z=[0.6-1.2] (low-z sample) and z>1.6 (high-z sample) as it is in these regions were we expect the mid-IR population to be dominated by intense dust-enshrouded activity such as star formation and black hole accretion. Investigations of the angular correlation function produce a correlation length are r0~15.9 Mpc for the high-z sample and r0~8.5 Mpc for the low-z one. Comparisons with physical models reveal that the high-z sources are exclusively associated with very massive (M>~10^{13} M_sun)haloes, comparable to those which locally host groups-to-clusters of galaxies, and are very common within such (rare) structures. Conversely, lower-z galaxies are found to reside in smaller halos (M_min~10^{12} M_sun) and to be very rare in such systems. While recent studies have determined a strong evolution of the 24um luminosity function between z~2 and z~0, they cannot provide information on the physical nature of such an evolution. Our clustering results instead indicate that this is due to the presence of different populations of objects inhabiting different structures, as active systems at z<~1.5 are found to be exclusively associated with low-mass galaxies, while very massive sources appear to have concluded their active phase before this epoch. Finally, we note that the small-scale clustering data seem to require steep profiles for the distribution of galaxies within their halos. This is suggestive of close encounters and/or mergers which could strongly favour both AGN and star-formation activity.

We present new results on the cosmological evolution of the near-infrared galaxy luminosity function, derived from the analysis of a new sample of \~22,000 K(AB) < 22.5 galaxies selected over an area of 0.6 square degrees from the Early Data Release of the UKIDSS Ultra Deep Survey (UDS). Our study has exploited the multi-wavelength coverage of the UDS field provided by the new UKIDSS WFCAM K and J-band imaging, the Subaru/XMM-Newton Deep Survey and the Spitzer-SWIRE Survey. The unique combination of large area and depth provided by this new survey minimises the complicating effect of cosmic variance and has allowed us, for the first time, to trace the evolution of the brightest sources out to z~2 with good statistical accuracy. In agreement with previous studies we find that the characteristic luminosity of the near-infrared luminosity function brightens by 1 magnitude between z=0 and z~2, while the total density decreases by a factor ~2. Using the rest-frame (U-B) colour to split the sample into red and blue galaxies, we confirm the classic luminosity-dependent colour bimodality at z<1. However, the strength of the colour bimodality is found to be a decreasing function of redshift, and seems to disappear by z>1.5. Due to the large size of our sample we are able to investigate the differing cosmological evolution of the red and blue galaxy populations. It is found that the space density of the brightest red galaxies (M_K < -23) stays approximately constant with redshift, and that these sources dominate the bright-end of the luminosity function at redshifts z<1. In contrast, the brightening of the characteristic luminosity and mild decrease in space density displayed by the blue galaxy population leads them to dominate the bright-end of the luminosity function at redshifts z>1.

We investigate colour selection techniques for high redshift galaxies in the UKIDSS Ultra Deep Survey Early Data Release (UDS EDR). Combined with very deep Subaru optical photometry, the depth (K_AB = 22.5) and area (0.62 deg^2) of the UDS EDR allows us to investigate optical/near-IR selection using a large sample of over 30,000 objects. By using the B-z, z-K colour-colour diagram (the BzK technique) we identify over 7500 candidate galaxies at z > 1.4, which can be further separated into passive and starforming systems (pBzK and sBzK respectively). Our unique sample allows us to identify a new feature not previously seen in BzK diagrams, consistent with the passively evolving track of early type galaxies at z < 1.4. We also compare the BzK technique with the R-K colour selection of Extremely Red Objects (EROs) and the J-K selection of Distant Red Galaxies (DRGs), and quantify the overlap between these populations. We find that the majority of DRGs, at these relatively bright magnitudes are also EROs. Since previous studies have found that DRGs at these magnitudes have redshifts of z ~ 1 we determine that these DRG/ERO galaxies have SEDs consistent with being dusty star-forming galaxies or AGN at z < 2. Finally we observe a flattening in the number counts of pBzK galaxies, similar to other studies, which may indicate that we are sampling the luminosity function of passive z > 1 galaxies over a narrow redshift range.

In a previous paper (Gavignaud et al. 2006), we presented the type-1 Active Galactic Nuclei (AGN) sample obtained from the first epoch data of the VIMOS-VLT Deep Survey (VVDS). The sample consists of 130 faint, broad-line AGN with redshift up to z=5 and 17.5< I <24.0, selected on the basis of their spectra. In this paper we present the measurement of the Optical Luminosity Function up to z=3.6 derived from this sample, we compare our results with previous results from brighter samples both at low and at high redshift. Our data, more than one magnitude fainter than previous optical surveys, allow us to constrain the faint part of the luminosity function up to high redshift. By combining our faint VVDS sample with the large sample of bright AGN extracted from the SDSS DR3 (Richards et al., 2006b) and testing a number of different evolutionary models, we find that the model which better represents the combined luminosity functions, over a wide range of redshift and luminosity, is a luminosity dependent density evolution (LDDE) model, similar to those derived from the major X-surveys. Such a parameterization allows the redshift of the AGN space density peak to change as a function of luminosity and explains the excess of faint AGN that we find at 1.0< z <1.5. On the basis of this model we find, for the first time from the analysis of optically selected samples, that the peak of the AGN space density shifts significantly towards lower redshift going to lower luminosity objects. This result, already found in a number of X-ray selected samples of AGN, is consistent with a scenario of "AGN cosmic downsizing", in which the density of more luminous AGN, possibly associated to more massive black holes, peaks earlier in the history of the Universe, than that of low luminosity ones.

We describe the number counts and spatial distribution of 239 Distant Red Galaxies (DRGs), selected from the Early Data Release of the UKIDSS Ultra Deep Survey. The DRGs are identified by their very red infrared colours with (J-K)AB>1.3, selected over 0.62 sq degree to a 90% completeness limit of KAB~20.7. This is the first time a large sample of bright DRGs has been studied within a contiguous area, and we provide the first measurements of their number counts and clustering. The population shows strong angular clustering, intermediate between those of K-selected field galaxies and optical/infrared-selected Extremely Red Galaxies. Adopting the redshift distributions determined from other recent studies, we infer a high correlation length of r0~11 h-1 Mpc. Such strong clustering could imply that our galaxies are hosted by very massive dark matter halos, consistent with the progenitors of present-day L>L* elliptical galaxies.

We present the first cluster catalogue extracted from the UKIRT Infrared Deep Sky Survey Early Data Release. The catalogue is created using UKIDSS Ultra Deep Survey infrared J and K data combined with 3.6 micro-m and 4.5 micro-m Spitzer bands and optical BVRi'z' imaging from the Subaru Telescope over 0.5 square degrees in the Subaru XMM-Newton Deep Field. We have created a new cluster-detection algorithm, based on the Friends-Of-Friends and Voronoi Tessellation methods, which utilises probability distribution functions derived from a photometric redshift analysis. We employ mock catalogues to understand the selection effects and contamination associated with the algorithm. The cluster catalogue contains 13 clusters at redshifts 0.61 <= z <= 1.39 with luminosities 10 L* <~ L_tot <~ 50 L*, corresponding to masses 5 x 10^13 M_sun <~ M_cluster <~ 3 x 10^14 M_sun for (M/M_sun) / (L/L_sun) = 75h. The measured sky surface density of ~ 10 deg^-2 for high-redshift (z=0.5-1.5), massive (>10^14 M_sun) clusters is precisely in line with theoretical predictions presented by Kneissl et al. (2001).

We construct a sample of extremely red objects (EROs) within the UKIDSS Ultra Deep Survey by combining the Early Data Release with optical data from the Subaru/XMM-Newton Deep Field. We find a total of 3715 objects over 2013 sq. arcmin with R-K>5.3 and K<=20.3, which is a higher surface density than found by previous studies. This is partly due to our ability to use a small aperture in which to measure colours, but is also the result of a genuine overdensity of objects compared to other fields. We separate our sample into passively-evolving and dusty star-forming galaxies using their RJK colours and investigate their radio properties using a deep radio map. The dusty population has a higher fraction of individually-detected radio sources and a higher mean radio flux density among the undetected objects, but the passive population has a higher fraction of bright radio sources, suggesting that AGNs are more prevalent among the passive ERO population.

We have exploited the large area coverage of the combined UKIDSS Ultra Deep Survey (UDS) and Subaru/XMM-Newton Deep Survey (SXDS) to search for bright Lyman-break galaxies (LBGs) at z >= 5. Using the available optical+near-infrared photometry to efficiently exclude low-redshift contaminants, we identify nine z >= 5 LBG candidates brighter than z'=25(AB) within the 0.6 square degree overlap region between the UDS early data release (EDR) and the optical coverage of the SXDS. Accounting for selection incompleteness, we estimate the corresponding surface density of z >= 5 LBGs with z'<=25(AB) to be 0.005+/-0.002 per square arcmin. Modelling of the optical+near-infrared photometry constrains the candidates' redshifts to lie in the range 5.1 < z < 5.9, and provides estimates for their stellar masses. Although the stellar mass estimates are individually uncertain, a stacking analysis suggests that the typical stellar mass of the LBG candidates is >~5x10^10 Msun which, if confirmed, places them amongst the most massive galaxies currently known at z >= 5. It is found that Lambda CDM structure formation can produce sufficient numbers of dark matter halos at z >= 5 to accommodate our estimated number density of massive LBGs for plausible values of sigma_8 and the ratio of stellar to dark matter. Moreover, it is found that recent galaxy formation models can also account for the existence of such massive galaxies at z >= 5. Finally, no evidence is found for the existence of LBGs with stellar masses in excess of ~3x10^11 Msun at this epoch, despite the large co-moving volume surveyed.

We present deep radio observations of the most distant complete quasar sample drawn from the Sloan Digital Sky Survey. Combining our new data with those from literature we obtain a sample which is ~100 per cent complete down to S_1.4GHz = 60 \mu Jy over the redshift range 3.8 < z < 5. The fraction of radio detections is relatively high (~43 per cent), similar to what observed locally in bright optical surveys. Even though the combined radio and optical properties of quasars remain overall unchanged from z ~ 5 to the local Universe, there is some evidence for a slight over-abundance of radio-loud sources at the highest redshifts when compared with the lower-z regime. Exploiting the deep radio VLA observations we present the first attempt to directly derive the radio luminosity function of bright quasars at z ~ 4. The unique depth -- both in radio and optical -- allows us to thoroughly explore the population of optically bright FR~II quasars up to z ~ 5 and opens a window on the behaviour of the brightest FR~I sources. A close investigation of the space density of radio loud quasars also suggests a differential evolution, with the more luminous sources showing a less pronounced cut-off at high z when compared with the less luminous ones.

Motivated by the claimed discovery of a very massive galaxy (HUDF-JD2; M~5x10^11 Msun) at extreme redshift (z = 6.5) within the Hubble Ultra Deep Field (HUDF) (Mobasher et al. 2005), we have completed a systematic search for comparably massive galaxies with z > 4 among the 2688 galaxies in our K < 23.5 (AB) catalogue within the CDFS/GOODS-South field. This search was conducted using redshift estimates based on the recently-completed, uniquely-deep 11-band imaging in this 125 square arcmin field, ~25 times larger than the NICMOS HUDF. Initial spectral fitting, based on published catalogue SExtractor photometry, led us to conclude that at least 2669 of the galaxies in our sample lie at z < 4. We carried out a detailed investigation of the 19 remaining z > 4 candidates, performing aperture photometry on all images, and including marginal detections and formal non-detections in the fitting process. This led to the rejection of a further 13 galaxies to lower redshift. Moreover, subjecting HUDF-JD2 to the same analysis, we find that it lies at z ~ 2.2, rather than the extreme redshift favoured by Mobasher et al. (2005). The 6 remaining candidates appear to be credible examples of galaxies in the redshift range z = 4 - 6, with plausible stellar ages. However, refitting with allowance for extreme values of extinction we find that, even for these objects, statistically acceptable solutions can be found at z < 3. Moreover, the recently-released Spitzer MIPS imaging in GOODS-South has revealed that 5 of our 6 final z > 4 candidates are detected at 24 microns. We conclude that there is no convincing evidence for any galaxy with M > 3 x 10^11 Msun, and z > 4 within the GOODS-South field (abridged).

At variance with most semi-analytic models, in the Anti-hierarchical Baryon Collapse scenario (Granato et al. 2001, 2004) the main driver of the galaxy formation and evolution is not the merging sequence but are baryon processes. This approach emphasizes, still in the framework of the hierarchical clustering paradigm for dark matter halos, feedback processes from supernova explosions and from active nuclei, that tie together star formation in spheroidal galaxies and the growth of black holes at their centers. We review some recent results showing the remarkably successful predictive power of this scenario, which allows us to account for the evolution with cosmic time of a broad variety of properties of galaxies and active nuclei, which proved to be very challenging for competing models.

We have selected and analysed the properties of a sample of 2905 Ks<21.5 galaxies in ~ 131 sq.arcmin of the Great Observatories Origins Deep Survey (GOODS) Chandra Deep Field South (CDFS), to obtain further constraints on the evolution of Ks-selected galaxies with respect to the results already obtained in previous studies. We made use of the public deep multiwavelength imaging from the optical B through the infrared (IR) 4.5 micron bands, in conjunction with available spectroscopic and COMBO17 data in the CDFS, to construct an optimised redshift catalogue for our galaxy sample. We computed the Ks-band LF and determined that its characteristic magnitude has a substantial brightening and a decreasing total density from z=0 to <z>=2.5. We also analysed the colours and number density evolution of galaxies with different stellar masses. Within our sample, and in contrast to what is observed for less massive systems, the vast majority (~ 85-90%) of the most massive (M>2.5x10^11 Msun) local galaxies appear to be in place before redshift z ~1. Around 65-70% of the total assemble between redshifts z=1 and z=3 and most of them display extremely red colours, suggesting that plausible star formation in these very massive systems should mainly proceed in obscured, short-timescale bursts. The remaining fraction (up to ~ 20%) could be in place at even higher redshifts z=3-4, pushing the first epoch of formation of massive galaxies beyond the limits of current near-IR surveys.

We show that the observed Velocity Dispersion Function of E/S0 galaxies matches strikingly well the distribution function of virial velocities of massive halos virializing at z > 1.5, as predicted by the standard hierarchical clustering scenario in a \LambdaCDM cosmology, for a constant ratio sigma/V_vir = 0.55 \pm 0.05, close to the value expected at virialization if it typically occurred at z > 3. This strongly suggests that dissipative processes and later merging events had little impact on the matter density profile. Adopting the above sigma/V_vir ratio, the observed relationships between photometric and dynamical properties which define the fundamental plane of elliptical galaxies, such as the luminosity-sigma (Faber-Jackson) and the luminosity-effective radius relations, as well as the M_BH-sigma relation, are nicely reproduced. Their shapes turn out to be determined by the mutual feedback of star-formation (and supernova explosions)and nuclear activity, along the lines discussed by Granato et al. (2004). To our knowledge, this is the first semi-analytic model for which simultaneous fits of the fundamental plane relations and of the epoch-dependent luminosity function of spheroidal galaxies have been presented.

The quasar sample selected by cross-correlating the FIRST and the 2dF Quasar Redshift Surveys allows us to explore, for the first time, the faint end of the radio and optical luminosity functions up to z = 2.2. We find indications (~3 \sigma) of a negative evolution for these faint sources at z > 1.8, both in radio and optical bands. This corresponds to a decrement in the space density of faint quasars of about a factor 2 at z=2.2 and confirms the presence of a differential evolution for the population of radio-active quasars. The faint end of both luminosity functions flattens and the comparison with the (optical) number density of the whole quasar population supports a dependence of the fraction of radio detected quasars on the optical luminosity. A progressive decrease in the fraction of quasars in the whole radio source population can be consistently accounted for within the `receding torus' scenario. The population of low luminosity quasars, which the FIRST-2dF detects, appears to depart from the `classical' scheme for radio-loud quasars.

We present a physical model for the coevolution of massive spheroidal galaxies and active nuclei at their centers. Supernova heating is increasingly effective in slowing down the star formation and in driving gas outflows in smaller and smaller dark matter halos. Thus the more massive protogalaxies virializing at early times are the sites of faster star formation. The correspondingly higher radiation drag causes a faster angular momentum loss by the gas and induces a larger accretion rate onto the central black hole. In turn, the kinetic energy of the outflows powered by the active nuclei can unbind the residual gas in a time shorter for larger halos. The model accounts for a broad variety of dynamical, photometric and metallicity properties of early-type galaxies, for the M_BH -- \sigma relation and for the local supermassive black-hole mass function.

We present a new approach to tackle the issue of radio loudness in quasars. We constrain a (simple) prescription for the intrinsic distribution of radio-to-optical ratios by comparing properties of Monte Carlo simulated samples with those of observed optically selected quasars. We find strong evidence for a dependence of the radio luminosity on the optical one, even though with a large scatter. The intrinsic distribution of the radio-to-optical ratios shows a peak at R*{1.4} ~ 0.3, with only ~5 per cent of objects being included in a high R*{1.4} tail which identifies the radio loud regime.

We present a new approach to tackle the issue of radio loudness in quasars. We constrain a (simple) prescription for the intrinsic distribution of radio-to-optical ratios by comparing properties of Monte Carlo simulated samples with those of observed optically selected quasars. We find strong evidence for a dependence of the radio luminosity on the optical one, even though with a large scatter. The dependence of the fraction of radio loud quasars on apparent and absolute optical magnitudes results in a selection effect related to the radio and optical limits of current surveys. The intrinsic distribution of the radio-to-optical ratios shows a peak at $R^*_{1.4} \sim 0.3$, with only $\simlt 5$ per cent of objects being included in a high $R^*_{1.4}$ tail which identifies the radio loud regime. No lack or deficit of sources -- but only a steep transition region -- is present between the radio loud and radio quiet populations at any $R^*_{1.4}$. We briefly discuss possible origins for this behaviour (e.g. absence of jets in radio quiet sources, large range of radiative radio efficiency, different life-times for the accretion and jet ejection phenomena, ...).

We present a detailed analysis of a sample of radio-detected quasars, obtained by matching together objects from the FIRST and 2dF Quasar Redshift Surveys. The dataset consists of 113 sources, spanning a redshift range 0.3 < z < 2.2, with optical magnitudes 18.25 < b_J < 20.85 and radio fluxes S_{1.4 GHz} < 1 mJy. These objects exhibit properties such as redshift and colour distribution in full agreement with those derived for the whole quasar population, suggestive of an independence of the mechanism(s) controlling the birth and life-time of quasars of their level of radio emission. The long debated question of radio-loud (RL)/radio-quiet (RQ) dichotomy is then investigated for the combined FIRST-2dF and FIRST-LBQS sample, since they present similar selection criteria. We find the fraction of radio detections to increase with magnitude from < 3% at the faintest levels up to 20% for the brightest sources. The classical RL/RQ dichotomy, in which the distribution of radio-to-optical ratios and/or radio luminosities shows a lack of sources, is ruled out by our analysis. We also find no tight relationship between optical and radio luminosities for sources in the considered sample, result that tends to exclude the mass of the quasar black hole as the physical quantity associated to the level of radio emission.