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

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.

Strong He II emission is produced by low-metallicity stellar populations. Here, we aim to identify and study a sample of He II $\lambda 1640$-emitting galaxies at redshifts of $z \sim 2.5-5$ in the deep VANDELS spectroscopic survey.. We identified a total of 33 Bright He II emitters (S/N > 2.5) and 17 Faint emitters (S/N < 2.5) in the VANDELS survey and used the available deep multi-wavelength data to study their physical properties. After identifying seven potential AGNs in our sample and discarding them from further analysis, we divided the sample of \emph{Bright} emitters into 20 \emph{Narrow} (FWHM < 1000 km s$^{-1}$) and 6 \emph{Broad} (FWHM > 1000 km s$^{-1}$) He II emitters. We created stacks of Faint, Narrow, and Broad emitters and measured other rest-frame UV lines such as O III] and C III] in both individual galaxies and stacks. We then compared the UV line ratios with the output of stellar population-synthesis models to study the ionising properties of He II emitters. We do not see a significant difference between the stellar masses, star-formation rates, and rest-frame UV magnitudes of galaxies with He II and no He II emission. The stellar population models reproduce the observed UV line ratios from metals in a consistent manner, however they under-predict the total number of \heii ionising photons, confirming earlier studies and suggesting that additional mechanisms capable of producing He II are needed, such as X-ray binaries or stripped stars. The models favour subsolar metallicities ($\sim0.1Z_\odot$) and young stellar ages ($10^6 - 10^7$ years) for the He II emitters. However, the metallicity measured for He II emitters is comparable to that of non-He II emitters at similar redshifts. We argue that galaxies with He II emission may have undergone a recent star-formation event, or may be powered by additional sources of He II ionisation.

We present the star formation histories of 39 galaxies with high quality rest-frame optical spectra at 0.5<z<1.3 selected to have strong Balmer absorption lines and/or Balmer break, and compare to a sample of spectroscopically selected quiescent galaxies at the same redshift. Photometric selection identifies a majority of objects that have clear evidence for a recent short-lived burst of star formation within the last 1.5 Gyr, i.e. "post-starburst" galaxies, however we show that good quality continuum spectra are required to obtain physical parameters such as burst mass fraction and burst age. Dust attenuation appears to be the primary cause for misidentification of post-starburst galaxies, leading to contamination in spectroscopic samples where only the [OII] emission line is available, as well as a small fraction of objects lost from photometric samples. The 31 confirmed post-starburst galaxies have formed 40-90% of their stellar mass in the last 1-1.5 Gyr. We use the derived star formation histories to find that the post-starburst galaxies are visible photometrically for 0.5-1 Gyr. This allows us to update a previous analysis to suggest that 25-50% of the growth of the red sequence at z~1 could be caused by a starburst followed by rapid quenching. We use the inferred maximum historical star formation rates of several 100-1000 Msun/yr and updated visibility times to confirm that sub-mm galaxies are likely progenitors of post-starburst galaxies. The short quenching timescales of 100-200 Myr are consistent with cosmological hydrodynamic models in which rapid quenching is caused by the mechanical expulsion of gas due to an AGN.

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

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

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

We present a new measurement of the evolving galaxy far-IR luminosity function (LF) extending out to redshifts z~5, with resulting implications for the level of dust-obscured star-formation density in the young Universe. To achieve this we have exploited recent advances in sub-mm/mm imaging with SCUBA-2 on the James Clerk Maxwell Telescope (JCMT) and the Atacama Large Millimeter/Submillimeter Array (ALMA), which together provide unconfused imaging with sufficient dynamic range to provide meaningful coverage of the luminosity-redshift plane out to z>4. Our results support previous indications that the faint-end slope of the far-IR LF is sufficiently flat that comoving luminosity-density is dominated by bright objects (~L*). However, we find that the number-density/luminosity of such sources at high redshifts has been severely over-estimated by studies that have attempted to push the highly-confused Herschel SPIRE surveys beyond z~2. Consequently we confirm recent reports that cosmic star-formation density is dominated by UV-visible star formation at z>4. Using both direct (1/Vmax) and maximum likelihood determinations of the LF, we find that its high-redshift evolution is well characterized by continued positive luminosity evolution coupled with negative density evolution (with increasing redshift). This explains why bright sub-mm sources continue to be found at z>5, even though their integrated contribution to cosmic star-formation density at such early times is very small. The evolution of the far-IR galaxy LF thus appears similar in form to that already established for active galactic nuclei, possibly reflecting a similar dependence on the growth of galaxy mass.

We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We firstly constrain the stellar mass vs stellar age relationship, finding a steep trend towards earlier average formation with increasing stellar mass of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* > 10^{11} M_\odot$. We show that this is consistent with other spectroscopic studies from $0 < z < 2$. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by Simba and IllustrisTNG agree well with observations at $z=0.1$, they are too shallow at $z=1$, predicting an evolution of $<0.5$ Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} M_\odot$) are consistent with formation in recent ($z<2$), intense starburst events, with timescales of $\lesssim500$ Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by $z=1$, they formed at $z\sim5$ and quenched by $z=3$. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity.

We investigate the prevalence of galactic-scale outflows in post-starburst (PSB) galaxies at high redshift ($1 < z < 1.4$), using the deep optical spectra available in the UKIDSS Ultra Deep Survey (UDS). We use a sample of $\sim40$ spectroscopically confirmed PSBs, recently identified in the UDS field, and perform a stacking analysis in order to analyse the structure of strong interstellar absorption features such as Mg ii ($\lambda2800$ Ang.). We find that for massive ($M_* > 10^{10}\rm\,M_{\odot}$) PSBs at $z > 1$, there is clear evidence for a strong blue-shifted component to the Mg ii absorption feature, indicative of high-velocity outflows ($v_{\rm out}\sim1150\pm160\rm\,km\,s^{-1}$) in the interstellar medium. We conclude that such outflows are typical in massive PSBs at this epoch, and potentially represent the residual signature of a feedback process that quenched these galaxies. Using full spectral fitting, we also obtain a typical stellar velocity dispersion $\sigma_*$ for these PSBs of $\sim200\rm\,km\,s^{-1}$, which confirms they are intrinsically massive in nature (dynamical mass $M_{\rm d}\sim10^{11}\rm\,M_{\odot}$). Given that these high-$z$ PSBs are also exceptionally compact ($r_{\rm e}\sim1$--$2\rm\,kpc$) and spheroidal (Sersic index $n\sim3$), we propose that the outflowing winds may have been launched during a recent compaction event (e.g. major merger or disc collapse) that triggered either a centralised starburst or active galactic nuclei (AGN) activity. Finally, we find no evidence for AGN signatures in the optical spectra of these PSBs, suggesting they were either quenched by stellar feedback from the starburst itself, or that if AGN feedback is responsible, the AGN episode that triggered quenching does not linger into the post-starburst phase.

We constrain the distribution of spatially offset Lyman-alpha emission (Ly$\alpha$) relative to rest-frame ultraviolet emission in $\sim300$ high redshift ($3<z<5.5$) Lyman-break galaxies (LBGs) exhibiting Ly$\alpha$ emission from VANDELS, a VLT/VIMOS slit-spectroscopic survey of the CANDELS Ultra Deep Survey and Chandra Deep Field South fields (${\simeq0.2}~\mathrm{deg}^2$ total). Because slit spectroscopy compresses two-dimensional spatial information into one spatial dimension, we use Bayesian inference to recover the underlying Ly$\alpha$ spatial offset distribution. We model the distribution using a 2D circular Gaussian, defined by a single parameter $\sigma_{r,\mathrm{Ly}\alpha}$, the standard deviation expressed in polar coordinates. Over the entire redshift range of our sample ($3<z<5.5$), we find $\sigma_{r,\mathrm{Ly}\alpha}=1.70^{+0.09}_{-0.08}$ kpc ($68\%$ conf.), corresponding to $\sim0.25$ arcsec at $\langle z\rangle=4.5$. We also find that $\sigma_{r,\mathrm{Ly}\alpha}$ decreases significantly with redshift. Because Ly$\alpha$ spatial offsets can cause slit-losses, the decrease in $\sigma_{r,\mathrm{Ly}\alpha}$ with redshift can partially explain the increase in the fraction of Ly$\alpha$ emitters observed in the literature over this same interval, although uncertainties are still too large to reach a strong conclusion. If $\sigma_{r,\mathrm{Ly}\alpha}$ continues to decrease into the reionization epoch, then the decrease in Ly$\alpha$ transmission from galaxies observed during this epoch might require an even higher neutral hydrogen fraction than what is currently inferred. Conversely, if spatial offsets increase with the increasing opacity of the IGM, slit losses may explain some of the drop in Ly$\alpha$ transmission observed at $z>6$. Spatially resolved observations of Ly$\alpha$ and UV continuum at $6<z<8$ are needed to settle the issue.

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 investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

We present a new analysis of the potential power of deep, near-infrared, imaging surveys with the James Webb Space Telescope (JWST) to improve our knowledge of galaxy evolution. In this work we properly simulate what can be achieved with realistic survey strategies, and utilise rigorous signal:noise calculations to calculate the resulting posterior constraints on the physical properties of galaxies. We explore a broad range of assumed input galaxy types (>20,000 models, including extremely dusty objects) across a wide redshift range (out to z~12), while at the same time considering a realistic mix of galaxy properties based on our current knowledge of the evolving population (as quantified through the Empirical Galaxy Generator: EGG). While our main focus is on imaging surveys with NIRCam, spanning lambda(obs) = 0.6-5.0 microns, an important goal of this work is to quantify the impact/added-value of: i) parallel imaging observations with MIRI at longer wavelengths, and ii) deeper supporting optical/UV imaging with HST (potentially prior to JWST launch) in maximising the power and robustness of a major extragalactic NIRCam survey. We show that MIRI parallel 7.7-micron imaging is of most value for better constraining the redshifts and stellar masses of the dustiest (A_V > 3) galaxies, while deep B-band imaging (reaching~28.5 AB mag) with ACS on HST is vital for determining the redshifts of the large numbers of faint/low-mass, z < 5 galaxies that will be detected in a deep JWST NIRCam survey.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Parametric models for galaxy star-formation histories (SFHs) are widely used, though they are known to impose strong priors on physical parameters. This has consequences for measurements of the galaxy stellar-mass function (GSMF), star-formation-rate density (SFRD) and star-forming main sequence (SFMS). We investigate the effects of the exponentially declining, delayed exponentially declining, lognormal and double power law SFH models using BAGPIPES. We demonstrate that each of these models imposes strong priors on specific star-formation rates (sSFRs), potentially biasing the SFMS, and also imposes a strong prior preference for young stellar populations. We show that stellar mass, SFR and mass-weighted age inferences from high-quality mock photometry vary with the choice of SFH model by at least 0.1, 0.3 and 0.2 dex respectively. However the biases with respect to the true values depend more on the true SFH shape than the choice of model. We also demonstrate that photometric data cannot discriminate between SFH models, meaning it is important to perform independent tests to find well-motivated priors. We finally fit a low-redshift, volume-complete sample of galaxies from the Galaxy and Mass Assembly (GAMA) Survey with each model. We demonstrate that our stellar masses and SFRs at redshift, $z\sim0.05$ are consistent with other analyses. However, our inferred cosmic SFRDs peak at $z\sim0.4$, approximately 6 Gyr later than direct observations suggest, meaning our mass-weighted ages are significantly underestimated. This makes the use of parametric SFH models for understanding mass assembly in galaxies challenging. In a companion paper we consider non-parametric SFH models.

We provide access to the current VANDELS data release (DR2) through a dedicated database, accessible at the link http://vandels.inaf.it/db. The data release includes the spectra for all galaxies for which the scheduled integration time was completed during the first two observing seasons and the spectra for those galaxies for which the scheduled integration time was 50\% complete at the end of season two (i.e. 20/40 hours and 40/80 hours). The total number of spectra released is 1362 (586 in CDFS and 776 in UDS). Spectra are stored as multi-extension FITS files containing the 1D extracted spectrum, the 2D linearly re-sampled spectrum, the 1D sky spectrum, the 1D noise estimate and the image thumbnail of the object.

We present Atacama Large Millimeter/Submillimeter Array observations of the rest-frame far-infrared (FIR) dust continuum emission of six bright Lyman-break galaxies (LBGs) at $z \simeq 7$. One LBG is detected ($5.2\sigma$ at peak emission), while the others remain individually undetected at the $3\sigma$ level. The average FIR luminosity of the sample is found to be $L_{\rm FIR} \simeq 2 \times 10^{11}\,{\rm L}_{\odot}$, corresponding to an obscured star-formation rate (SFR) that is comparable to that inferred from the unobscured UV emission. In comparison to the infrared excess (IRX$\,=L_{\rm FIR}/L_{\rm UV}$)-$\beta$ relation, our results are consistent with a Calzetti-like attenuation law (assuming a dust temperature of T = 40-50 K). We find a physical offset of 3 kpc between the dust continuum emission and the rest-frame UV light probed by Hubble Space Telescope imaging for galaxy ID65666 at $z = 7.17^{+0.09}_{-0.06}$. The offset is suggestive of an inhomogeneous dust distribution, where 75% of the total star formation activity (SFR$ \,\simeq 70\,{\rm M}_{\odot}/{\rm yr}$) of the galaxy is completely obscured. Our results provide direct evidence that dust obscuration plays a key role in shaping the bright-end of the observed rest-frame UV luminosity function at $z \simeq 7$, in agreement with cosmological galaxy formation simulations. The existence of a heavily-obscured component of galaxy ID65666 indicates that dusty star-forming regions, or even entire galaxies, that are "UV-dark" are significant even in the $z \simeq 7$ galaxy population.

We present Bayesian Analysis of Galaxies for Physical Inference and Parameter EStimation, or BAGPIPES, a new Python tool which can be used to rapidly generate complex model galaxy spectra and to fit these to arbitrary combinations of spectroscopic and photometric data using the MultiNest nested sampling algorithm. We extensively test our ability to recover realistic star-formation histories (SFHs) by fitting mock observations of quiescent galaxies from the MUFASA simulation. We then perform a detailed analysis of the SFHs of a sample of 9289 quiescent galaxies from UltraVISTA with stellar masses, $M_* > 10^{10}\ \mathrm{M_\odot}$ and redshifts $0.25 < z < 3.75$. The majority of our sample exhibit SFHs which rise gradually then quench relatively rapidly, over $1{-}2$ Gyr. This behaviour is consistent with recent cosmological hydrodynamic simulations, where AGN-driven feedback in the low-accretion (jet) mode is the dominant quenching mechanism. At $z > 1$ we also find a class of objects with SFHs which rise and fall very rapidly, with quenching timescales of $< 1$ Gyr, consistent with quasar-mode AGN feedback. Finally, at $z < 1$ we find a population with SFHs which quench more slowly than they rise, over $>3$ Gyr, which we speculate to be the result of diminishing overall cosmic gas supply. We confirm the mass-accelerated evolution (downsizing) trend, and a trend towards more rapid quenching at higher stellar masses. However, our results suggest that the latter is a natural consequence of mass-accelerated evolution, rather than a change in quenching physics with stellar mass. We find $61\pm8$ per cent of $z > 1.5$ massive quenched galaxies undergo significant further evolution by $z = 0.5$. BAGPIPES is available at https://bagpipes.readthedocs.io

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 the results of a new study of the relationship between infrared excess (IRX), UV spectral slope (beta) and stellar mass at redshifts 2<z<3, based on a deep Atacama Large Millimeter Array (ALMA) 1.3-mm continuum mosaic of the Hubble Ultra Deep Field (HUDF). Excluding the most heavily-obscured sources, we use a stacking analysis to show that z~2.5 star-forming galaxies in the mass range 9.25 <= log(M/Msun) <= 10.75 are fully consistent with the IRX-beta relation expected for a relatively grey attenuation curve, similar to the commonly adopted Calzetti law. Based on a large, mass complete, sample of 2 <= z <= 3 star-forming galaxies drawn from multiple surveys, we proceed to derive a new empirical relationship between beta and stellar mass, making it possible to predict UV attenuation (A_1600) and IRX as a function of stellar mass, for any assumed attenuation law. Once again, we find that z~2.5 star-forming galaxies follow A_1600-mass and IRX-mass relations consistent with a relatively grey attenuation law, and find no compelling evidence that star-forming galaxies at this epoch follow a reddening law as steep as the Small Magellanic Cloud (SMC) extinction curve. In fact, we use a simple simulation to demonstrate that previous determinations of the IRX-beta relation may have been biased toward low values of IRX at red values of beta, mimicking the signature expected for an SMC-like dust law. We show that this provides a plausible mechanism for reconciling apparently contradictory results in the literature and that, based on typical measurement uncertainties, stellar mass provides a cleaner prediction of UV attenuation than beta. Although the situation at lower stellar masses remains uncertain, we conclude that for 2<z<3 star-forming galaxies with log(M/Msun) >= 9.75, both the IRX-beta and IRX-mass relations are well described by a Calzetti-like attenuation law.

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 investigate the relationship between the quenching of star formation and the structural transformation of massive galaxies, using a large sample of photometrically-selected post-starburst galaxies in the UKIDSS UDS field. We find that post-starburst galaxies at high-redshift ($z>1$) show high S\'ersic indices, significantly higher than those of active star-forming galaxies, but with a distribution that is indistinguishable from the old quiescent population. We conclude that the morphological transformation occurs before (or during) the quenching of star formation. Recently quenched galaxies are also the most compact; we find evidence that massive post-starburst galaxies (M$_{\ast}> 10^{10.5} ~$M$_{\odot}$) at high redshift ($z>1$) are on average smaller than comparable quiescent galaxies at the same epoch. Our findings are consistent with a scenario in which massive passive galaxies are formed from three distinct phases: (1) gas-rich dissipative collapse to very high densities, forming the proto-spheroid; (2) rapid quenching of star formation, to create the "red nugget" with post-starburst features; (3) a gradual growth in size as the population ages, perhaps as a result of minor mergers.

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

We reinvestigate a claimed sample of 22 X-ray detected active galactic nuclei (AGN) at redshifts z > 4, which has reignited the debate as to whether young galaxies or AGN reionized the Universe. These sources lie within the GOODS-S/CANDELS field, and we examine both the robustness of the claimed X-ray detections (within the Chandra 4Ms imaging) and perform an independent analysis of the photometric redshifts of the optical/infrared counterparts. We confirm the reality of only 15 of the 22 reported X-ray detections, and moreover find that only 12 of the 22 optical/infrared counterpart galaxies actually lie robustly at z > 4. Combining these results we find convincing evidence for only 7 X-ray AGN at z > 4 in the GOODS-S field, of which only one lies at z > 5. We recalculate the evolving far-UV (1500 Angstrom) luminosity density produced by AGN at high redshift, and find that it declines rapidly from z = 4 to z = 6, in agreement with several other recent studies of the evolving AGN luminosity function. The associated rapid decline in inferred hydrogen-ionizing emissivity contributed by AGN falls an order-of-magnitude short of the level required to maintain hydrogen ionization at z ~ 6. We conclude that all available evidence continues to favour a scenario in which young galaxies reionized the Universe, with AGN making, at most, a very minor contribution to cosmic hydrogen reionization.

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 search for passive galaxies at $z>3$ in the GOODS-South field, using different techniques based on photometric data, and paying attention to develop methods that are sensitive to objects that have become passive shortly before the epoch of observation. We use CANDELS HST catalogues, ultra-deep $Ks$ data and new IRAC photometry, performing spectral energy distribution fitting using models with abruptly quenched star formation histories. We then single out galaxies which are best fitted by a passively evolving model, and having only low probability $( < 5 \%)$ star-forming solutions. We verify the effects of including nebular lines emission, and we consider possible solutions at different redshifts. The number of selected sources dramatically depends on the models used in the SED fitting. Without including emission lines and with photometric redshifts fixed at the CANDELS estimate, we single out 30 candidates; the inclusion of nebular lines emission reduces the sample to 10 objects; allowing for solutions at different redshifts, only 2 galaxies survive as robust candidates. Most of the candidates are not far-infrared emitters, corroborating their association with passive galaxies. Our results translate into an upper limit in the number density of $\sim 0.173$ arcmin$^2$ above the detection limit. However, we conclude that the selection of passive galaxies at $z>3$ is still subject to significant uncertainties, being sensitive to assumptions in the SED modeling adopted and to the relatively low S/N of the objects. By means of dedicated simulations, we show that JWST will greatly enhance the accuracy, allowing for a much more robust classification.

We investigate the relation between the star formation rate (SFR) and the stellar mass (M), i.e. the Main Sequence (MS) relation of star-forming galaxies, at 1.3<=z<6 in the first four HST Frontier Fields, based on rest-frame UV observations. Gravitational lensing combined with deep HST observations allows us to extend the analysis of the MS down to stellar masses as low as logM/Msun~7.5 at z<~4 and logM/Msun~8 at higher redshifts, a factor of ~10 below most previous results. We perform an accurate simulation to take into account the effect of observational uncertainties on the MS and correct for the Eddington bias. This step allows us to reliably measure the MS and in particular its slope. While the normalization increases with redshift, we fit an unevolving and approximately linear slope. We nicely extend to lower masses the results of brighter surveys. Thanks to the large dynamic range in mass and by making use of the simulation, we analyzed any possible mass dependence of the dispersion around the MS. We find tentative evidence that the scatter decreases with increasing stellar masses, suggesting a larger variety of star formation histories in low mass galaxies. This trend agrees with theoretical predictions, and is explained as either a consequence of the smaller number of progenitors of low mass galaxies in a hierarchical scenario and/or of the efficient but intermittent stellar feedback processes in low mass halos. Finally, we observe an increase in the SFR per unit stellar mass with redshift milder than predicted by theoretical models, implying a still incomplete understanding of the processes responsible for galaxy growth.

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.

We present the multiwavelength photometry of two Frontier Fields massive galaxy clusters MACS-J0717 and MACS-J1149 and their parallel fields, ranging from HST to ground based K and Spitzer IRAC bands, and the public release of photometric redshifts and rest frame properties of galaxies found in cluster and parallel pointings. This work was done within ASTRODEEP project and aims to provide a reference for future investigations of the extragalactic populations. To fully exploit the depth of the images and detect faint sources we used an accurate procedure which carefully removes the foreground light of bright cluster sources and the intra-cluster light thus enabling detection and measurement of accurate fluxes in crowded cluster regions. This same procedure has been successfully used to derive the photometric catalogue of MACS-J0416 and Abell-2744. The obtained multi-band photometry was used to derive photometric redshifts, magnification and physical properties of sources. In line with the first two FF catalogues released by ASTRODEEP, the photometric redshifts reach $\sim$4$\%$ accuracy. Moreover we extend the presently available samples to galaxies intrinsically as faint as H160$\sim$32-34 mag thanks the magnification factors induced to strong gravitational lensing. Our analysis allows us to probe galaxy masses larger then 10$^{7}$ M$\odot$ and/or SFR=0.1-1M$\odot$/yr out to redshift z$>6$.

We present the results of a study investigating the dust attenuation law at $z\simeq 5$, based on synthetic spectral energy distributions (SEDs) calculated for a sample of N=498 galaxies drawn from the First Billion Years (FiBY) simulation project. The simulated galaxies at $z\simeq 5$, which have M$_{1500} \leq -18.0$ and $7.5 \leq \rm{log(M/M}_{\odot}\rm{)} \leq 10.2$, display a mass-dependent $\alpha$-enhancement, with a median value of $[\alpha/\rm{Fe}]_{z=5}~\simeq~4~\times~[\alpha/\rm{Fe}]_{Z_{\odot}}$. The median Fe/H ratio of the simulated galaxies is $0.14\pm0.05$ which, even including the effects of nebular continuum, produces steep intrinsic UV continuum slopes; $\langle \beta_{i} \rangle = -2.4 \pm 0.05$. Using a set of simple dust attenuation models, in which the wavelength-dependent attenuation is assumed to be of the form $A(\lambda) \propto \lambda^{n}$, we explore the parameter values which best reproduce the observed $z=5$ luminosity function (LF) and colour-magnitude relation (CMR). We find that a simple model in which the absolute UV attenuation is a linearly increasing function of log stellar mass, and the dust attenuation slope ($n$) is within the range $-0.7 \leq n \leq-0.3$, can successfully reproduce the LF and CMR over a wide range of stellar population synthesis model (SPS) assumptions. This range of attenuation curves is consistent with a power-law fit to the Calzetti attenuation law in the UV ($n=-0.55$), and other similarly `grey' star-forming galaxy attenuation curves recently derived at $z\simeq2$. In contrast, attenuation curves as steep as the Small Magellanic Cloud (SMC) extinction curve ($n=-1.24$) are formally ruled out. Finally, we show that our models are consistent with recent 1.3mm ALMA observations of the Hubble Ultra Deep Field (HUDF), and predict the form of the $z\simeq5$ IRX$-\beta$ relation.

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.

The z = 6.6 Lyman-$\alpha$ emitter 'CR7' has been claimed to have a Population III-like stellar population, or alternatively, be a candidate Direct Collapse Black Hole (DCBH). In this paper we investigate the evidence for these exotic scenarios using recently available, deeper, optical, near-infrared and mid-infrared imaging. We find strong Spitzer/IRAC detections for the main component of CR7 at 3.6 and 4.5 microns, and show that it has a blue colour ([3.6] - [4.5] $= -1.2\pm 0.3$). This colour cannot be reproduced by current Pop. III or pristine DCBH models. Instead, the results suggest that the [3.6] band is contaminated by the [OIII]4959,5007 emission line with an implied rest-frame equivalent width of EW_0 (H$\beta$ + [OIII]) $\gtrsim 2000$\AA. Furthermore, we find that new near-infrared data from the UltraVISTA survey supports a weaker He II 1640 emission line than previously measured, with EW_0 $= 40 \pm 30$\AA. For the fainter components of CR7 visible in Hubble Space Telescope imaging, we find no evidence that they are particularly red as previously claimed, and show that the derived masses and ages are considerably uncertain. In light of the likely detection of strong [OIII] emission in CR7 we discuss other more standard interpretations of the system that are consistent with the data. We find that a low-mass, narrow-line AGN can reproduce the observed features of CR7, including the lack of radio and X-ray detections. Alternatively, a young, low-metallicity (~1/200 solar) starburst, modelled including binary stellar pathways, can reproduce the inferred strength of the He II line and simultaneously the strength of the observed [OIII] emission, but only if the gas shows super-solar $\alpha$-element abundances (O/Fe ~ 5 O/Fe solar).

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 new Hubble Space Telescope/Wide Field Camera 3 imaging of 25 extremely luminous (-23.2 < M_ UV < -21.2) Lyman-break galaxies (LBGs) at z ~ 7. The sample was initially selected from 1.65 deg^2 of ground-based imaging in the UltraVISTA/COSMOS and UDS/SXDS fields, and includes the extreme Lyman-alpha emitters, `Himiko' and `CR7'. A deconfusion analysis of the deep Spitzer photometry available suggests that these galaxies exhibit strong rest-frame optical nebular emission lines (EW_0(H_beta + [OIII]) > 600A). We find that irregular, multiple-component morphologies suggestive of clumpy or merging systems are common (f_multi > 0.4) in bright z ~ 7 galaxies, and ubiquitous at the very bright end (M_UV < -22.5). The galaxies have half-light radii in the range r_1/2 ~ 0.5-3 kpc. The size measurements provide the first determination of the size-luminosity relation at z ~ 7 that extends to M_UV ~ -23. We find the relation to be steep with r_1/2 ~ L^1/2. Excluding clumpy, multi-component galaxies however, we find a shallower relation that implies an increased star-formation rate surface density in bright LBGs. Using the new, independent, HST/WFC3 data we confirm that the rest-frame UV luminosity function at z ~ 7 favours a power-law decline at the bright-end, compared to an exponential Schechter function drop-off. Finally, these results have important implications for the Euclid mission, which we predict will detect > 1000 similarly bright galaxies at z ~ 7. Our new HST imaging suggests that the vast majority of these galaxies will be spatially resolved by Euclid, mitigating concerns over dwarf star contamination.

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

We present the results of a new study of the K-band galaxy luminosity function (KLF) at redshifts z<3.75, based on a nested combination of the UltraVISTA, CANDELS and HUDF surveys. The large dynamic range in luminosity spanned by this new dataset (3-4 dex over the full redshift range) is sufficient to clearly demonstrate for the first time that the faint-end slope of the KLF at z>0.25 is relatively steep (-1.3<alpha<-1.5 for a single Schechter function), in good agreement with recent theoretical and phenomenological models. Moreover, based on our new dataset we find that a double Schechter function provides a significantly improved description of the KLF at z<2. At redshifts z>0.25 the evolution of the KLF is remarkably smooth, with little or no evolution evident at faint (M_K>-20.5) or bright magnitudes (M_K<-24.5). Instead, the KLF is seen to evolve rapidly at intermediate magnitudes, with the number density of galaxies at M_K~-23 dropping by a factor of ~5 over the redshift interval 0.25<z<3.75. Motivated by this, we explore a simple description of the evolving KLF based on a double Schechter function with fixed faint-end slopes (alpha_1=-0.5, alpha_2=-1.5) and a shared characteristic magnitude (M_K*). According to this parameterisation, the normalisation of the component which dominates the faint-end of the KLF remains approximately constant, with phi*_2 decreasing by only a factor of ~2 between z~0 and z~3.25. In contrast, the component which dominates the bright end of the KLF at low redshifts evolves dramatically, becoming essentially negligible by z~3. Finally, we note that within this parameterisation, the observed evolution of M_K* between z~0 and z~3.25 is entirely consistent with M_K* corresponding to a constant stellar mass of M*~5x10^10 Msun at all redshifts.

We present the evolution in the number density and stellar mass functions of photometrically selected post-starburst galaxies in the UKIDSS Deep Survey (UDS), with redshifts of 0.5<z<2 and stellar masses logM>10. We find that this transitionary species of galaxy is rare at all redshifts, contributing ~5% of the total population at z~2, to <1% by z~0.5. By comparing the mass functions of quiescent galaxies to post-starburst galaxies at three cosmic epochs, we show that rapid quenching of star formation can account for 100% of quiescent galaxy formation, if the post-starburst spectral features are visible for ~250Myr. The flattening of the low mass end of the quiescent galaxy stellar mass function seen at z~1 can be entirely explained by the addition of rapidly quenched galaxies. Only if a significant fraction of post-starburst galaxies have features that are visible for longer than 250Myr, or they acquire new gas and return to the star-forming sequence, can there be significant growth of the red sequence from a slower quenching route. The shape of the mass function of these transitory post-starburst galaxies resembles that of quiescent galaxies at z~2, with a preferred stellar mass of logM~10.6, but evolves steadily to resemble that of star-forming galaxies at z<1. This leads us to propose a dual origin for post-starburst galaxies: (1) at z>2 they are exclusively massive galaxies that have formed the bulk of their stars during a rapid assembly period, followed by complete quenching of further star formation, (2) at z<1 they are caused by the rapid quenching of gas-rich star-forming galaxies, independent of stellar mass, possibly due to environment and/or gas-rich major mergers.

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

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.

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

We present a new, robust measurement of the evolving rest-frame UV galaxy luminosity function (LF) over the key redshift range z = 2 - 4. Our results are based on the high dynamic range provided by combining the HUDF, CANDELS/GOODS-South, and UltraVISTA/COSMOS surveys. We utilise the unparalleled multi-frequency photometry available in this survey `wedding cake' to compile complete galaxy samples at z ~ 2,3,4 via photometric redshifts (calibrated against the latest spectroscopy) rather than colour-colour selection, and to determine accurate rest-frame UV absolute magnitudes from SED fitting. Our new determinations of the UV LF extend from M_{1500} ~ -22 down to M_{1500} = -14.5, -15.5 and -16 at z ~ 2, 3 and 4 respectively (thus reaching ~ 3 - 4 magnitudes fainter than previous blank-field studies at z ~ 2 - 3). At z ~ 2 - 3 we find a much shallower faint-end slope (alpha = -1.32 +- 0.03) than the steeper values (alpha ~ -1.7) reported by Reddy & Steidel (2009) or by Alavi et al. (2014), and show that this new measurement is robust. By z ~ 4 the faint-end slope has steepened slightly, to alpha = -1.43 +- 0.04, and we show that these measurements are consistent with the overall evolutionary trend from z = 0 to z = 8. Finally, we find that while characteristic number density (phi*) drops from z ~ 2 to z ~ 4, characteristic luminosity (M*) brightens by ~ 1 mag over this redshift range. This, combined with the new flatter faint-end slopes, has the consequence that UV luminosity density (and hence unobscured star-formation density) peaks at z ~ 2.5 - 3, when the Universe was ~ 2.5 Gyr old.

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

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

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

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

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

We present the results of a new search for galaxies at redshift z ~ 9 in the first two Hubble Frontier Fields with completed HST WFC3/IR and ACS imaging. To ensure robust photometric redshift solutions, and to minimize incompleteness, we confine our search to objects with H_{160} < 28.6 (AB mag), consider only image regions with an rms noise sigma_{160} > 30 mag (within a 0.5-arcsec diameter aperture), and insist on detections in both H_{160} and J_{140}. The result is a survey covering an effective area (after accounting for magnification) of 10.9 sq. arcmin, which yields 12 galaxies at 8.4 < z < 9.5. Within the Abell-2744 cluster and parallel fields we confirm the three brightest objects reported by Ishigaki et al. (2014), but recover only one of the four z > 8.4 sources reported by Zheng et al. (2014). In the MACSJ0416.1-240 cluster field we report five objects, and explain why each of these eluded detection or classification as z ~ 9 galaxies in the published searches of the shallower CLASH data. Finally, we uncover four z ~ 9 galaxies from the previously unsearched MACSJ0416.1-240 parallel field. Based on the published magnification maps we find that only one of these 12 galaxies is likely boosted by more than a factor of two by gravitational lensing. Consequently we are able to perform a fairly straightforward reanalysis of the normalization of the z ~ 9 UV galaxy luminosity function as explored previously in the HUDF12 programme. We conclude that the new data strengthen the evidence for a continued smooth decline in UV luminosity density (and hence star-formation rate density) from z ~ 8 to z ~ 9, contrary to recent reports of a marked drop-off at these redshifts. This provides further support for the scenario in which early galaxy evolution is sufficiently extended to explain cosmic reionization.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

Strong gravitational lensing provides a powerful means for studying faint galaxies in the distant universe. By magnifying the apparent brightness of background sources, massive clusters enable the detection of galaxies fainter than the usual sensitivity limit for blank fields. However, this gain in effective sensitivity comes at the cost of a reduced survey volume and, in this {\it Letter}, we demonstrate there is an associated increase in the cosmic variance uncertainty. As an example, we show that the cosmic variance uncertainty of the high redshift population viewed through the Hubble Space Telescope Frontier Field cluster Abell 2744 increases from ~35% at redshift z~7 to >~65% at z~10. Previous studies of high redshift galaxies identified in the Frontier Fields have underestimated the cosmic variance uncertainty that will affect the ultimate constraints on both the faint end slope of the high-redshift luminosity function and the cosmic star formation rate density, key goals of the Frontier Field program.

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

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

We present optical and near-infrared imaging of the field of the z=7.0842 quasar ULAS J112001.48+064124.3 taken with the Hubble Space Telescope. We use these data to search for galaxies that may be physically associated with the quasar, using the Lyman break technique, and find three such objects, although the detection of one in Spitzer Space Telescope imaging strongly suggests it lies at z~2. This is consistent with the field luminosity function and indicates that there is no excess of >L* galaxies within 1 Mpc of the quasar. A detection of the quasar shortward of the Ly-alpha line is consistent with the previously observed evolution of the intergalactic medium at z>5.5.

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.

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

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

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

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.

We present new results on the cosmic star formation history in the SXDS-UDS field out to z=1.6. We compile narrow-band data from the Subaru Telescope and the Visible and Infrared Survey Telescope for Astronomy (VISTA) in conjunction with broad-band data from the SXDS and UDS, to make a selection of 5725 emission-line galaxies in 12 redshift slices, spanning 10 Gyr of cosmic time. We determine photometric redshifts for the sample using 11-band photometry, and use a spectroscopically confirmed subset to fine tune the resultant redshift distribution. We use the maximum-likelihood technique to determine luminosity functions in each redshift slice and model the selection effects inherent in any narrow-band selection statistically, to obviate the retrospective corrections ordinarily required. The deep narrow-band data are sensitive to very low star formation rates (SFRs), and allow an accurate evaluation of the faint end slope of the Schechter function, alpha. We find that alpha is particularly sensitive to the assumed faintest broadband magnitude of a galaxy capable of hosting an emission line, and propose that this limit should be empirically motivated. For this analysis we base our threshold on the limiting observed equivalent widths of emission lines in the local Universe. We compute the characteristic SFR of galaxies in each redshift slice, and the integrated SFR density, rho_{SFR}. We find our results to be in good agreement with the literature and parametrize the evolution of the SFR density as rho_{SFR} proportional to (1 + z)^{4.58} confirming a steep decline in star formation activity since z = 1.6.

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 present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.

We present and analyse the deepest, high-quality Ks-band images ever obtained of luminous quasars at z~4, in an attempt to determine the basic properties of their host galaxies less than 1 Gyr after the first recorded appearance of black holes with Mbh > 10^9 Msol. To maximise the robustness of our results we have carefully selected two SDSS quasars at z~4. These quasars are representative of the most luminous quasars known at this epoch but they also, crucially, lie within 40 arcsec of comparably-bright foreground stars (required for accurate PSF definition), and have redshifts which ensure line-free Ks-band imaging. The data were obtained in excellent seeing (<0.4-arcsec) at the ESO VLT with integration times of ~5.5 hours per source. Via carefully-controlled separation of host-galaxy and nuclear light, we estimate the luminosities and stellar masses of the host galaxies, and set constraints on their half-light radii. The quasar host galaxies have K-band luminosities similar to radio galaxies at comparable redshifts, suggesting that these quasar hosts are also among the most massive galaxies in existence at this epoch. However, the quasar hosts are a factor ~5 smaller than the host galaxies of luminous low-redshift quasars. We estimate the stellar masses of the z~4 host galaxies to lie in the range 2-10x10^11 Msol, and use the CIV emission line in the Sloan spectra to estimate the masses of their black holes. The results imply a black-hole:host-galaxy mass ratio Mbh:Mgal~0.01-0.05. This is an order of magnitude higher than typically seen in the low-redshift Universe, and is consistent with existing evidence for a systematic growth in this mass ratio with increasing redshift, at least for objects selected as powerful AGN.

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.

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.

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

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

We present deep HST/WFPC2, rest-frame U images of 17 ~L* quasars at z=1 and z=2 (V and I bands respectively), designed to explore the host galaxies. We fit the images with simple axisymmetric galaxy models, including a point-source, in order to separate nuclear and host-galaxy emission. We successfully model all of the host galaxies, with luminosities stable to within 0.3 mag. Combining with our earlier NICMOS rest-frame optical study of the same sample, we provide the first rest-frame U-V colours for a sample of quasar host galaxies. While the optical luminosities of their host galaxies indicate that they are drawn purely from the most massive (>~L*) early-type galaxy population, their colours are systematically bluer than those of comparably massive galaxies at the same redshift. The host galaxies of the radio-loud quasars (RLQ) in our sample are more luminous than their radio-quiet quasar (RQQ) counterparts at each epoch, but have indistinguishable colours, confirming that the RLQ's are drawn from only the most massive galaxies (10^{11}-10^{12} M_sun, even at z~2), while the RQQ's are slightly less massive (~10^{11} M_sun). This is consistent with the well-known anticorrelation between radio-loudness and accretion rate. Using simple stellar population "frosting" models we estimate mean star formation rates of ~350 M_sun/yr for the RLQ's and ~100 M_sun/yr for the RQQ's at z~2. By z~1, these rates have fallen to ~150 M_sun/yr for the RLQ's and ~50 M_sun/yr for the RQQ's. We conclude that while the host galaxies are extremely massive, they remain actively star-forming at, or close to, the epoch of the quasar.

We present the results of a search for the most luminous star-forming galaxies at redshifts z~6 based on CFHT Legacy Survey data. We identify a sample of 40 Lyman break galaxies brighter than magnitude z'=25.3 across an area of almost 4 square degrees. Sensitive spectroscopic observations of seven galaxies provide redshifts for four, of which only two have moderate to strong Lyman alpha emission lines. All four have clear continuum breaks in their spectra. Approximately half of the Lyman break galaxies are spatially resolved in 0.7 arcsec seeing images, indicating larger sizes than lower luminosity galaxies discovered with the Hubble Space Telescope, possibly due to on-going mergers. The stacked optical and infrared photometry is consistent with a galaxy model with stellar mass ~ 10^{10} solar masses. There is strong evidence for substantial dust reddening with a best-fit A_V=0.7 and A_V>0.48 at 2 sigma confidence, in contrast to the typical dust-free galaxies of lower luminosity at this epoch. The spatial extent and spectral energy distribution suggest that the most luminous z~6 galaxies are undergoing merger-induced starbursts. The luminosity function of z=5.9 star-forming galaxies is derived. This agrees well with previous work and shows strong evidence for an exponential decline at the bright end, indicating that the feedback processes which govern the shape of the bright end are occurring effectively at this epoch.

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

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

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

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

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 physical properties inferred from the SEDs of z>3 galaxies have been influential in shaping our understanding of early galaxy formation and the role galaxies may play in cosmic reionization. Of particular importance is the stellar mass density at early times which represents the integral of earlier star formation. An important puzzle arising from the measurements so far reported is that the specific star formation rates (sSFR) evolve far less rapidly than expected in most theoretical models. Yet the observations underpinning these results remain very uncertain, owing in part to the possible contamination of rest-optical broadband light from strong nebular emission lines. To quantify the contribution of nebular emission to broad-band fluxes, we investigate the SEDs of 92 spectroscopically-confirmed galaxies in the redshift range 3.8<z<5.0 chosen because the H-alpha line lies within the Spitzer/IRAC 3.6 um filter. We demonstrate that the 3.6 um flux is systematically in excess of that expected from stellar continuum, which we derive by fitting the SED with population synthesis models. No such excess is seen in a control sample at 3.1<z<3.6 in which there is no nebular contamination in the IRAC filters. From the distribution of our 3.6 um flux excesses, we derive an H-alpha equivalent width (EW) distribution. The mean rest-frame H-alpha EW we infer at 3.8<z<5.0 (270 A) indicates that nebular emission contributes at least 30% of the 3.6 um flux. Via our empirically-derived EW distribution we correct the available stellar mass densities and show that the sSFR evolves more rapidly at z>4 than previously thought, supporting up to a 5x increase between z~2 and 7. Such a trend is much closer to theoretical expectations. Given our findings, we discuss the prospects for verifying quantitatively the nebular emission line strengths prior to the launch of the James Webb Space Telescope.

We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8.

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.

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.

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

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.

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.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

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

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

Star forming galaxies represent a valuable tracer of cosmic history. Recent observational progress with Hubble Space Telescope has led to the discovery and study of the earliest-known galaxies corresponding to a period when the Universe was only ~800 million years old. Intense ultraviolet radiation from these early galaxies probably induced a major event in cosmic history: the reionization of intergalactic hydrogen. New techniques are being developed to understand the properties of these most distant galaxies and determine their influence on the evolution of the universe.

We investigate the stellar populations of Lyman alpha emitters (LAEs) at z=5.7 and 6.6 in a 0.65 deg^2 sky of the Subaru/XMM-Newton Deep Survey (SXDS) Field, using deep images taken with Subaru/Suprime-Cam, UKIRT/WFCAM, and Spitzer/IRAC. We produce stacked multiband images at each redshift from 165 (z=5.7) and 91 (z=6.6) IRAC-undetected objects, to derive typical spectral energy distributions (SEDs) of z~6-7 LAEs for the first time. The stacked LAEs have as blue UV continua as the HST/WFC3 z-dropout galaxies of similar Muv, with a spectral slope beta ~ -3, but at the same time they have red UV-to-optical colors with detection in the 3.6um band. Using SED fitting we find that the stacked LAEs have low stellar masses of ~(3-10)*10^7 Msun, very young ages of ~1-3 Myr, negligible dust extinction, and strong nebular emission from the ionized interstellar medium, although the z=6.6 object is fitted similarly well with high-mass models without nebular emission; inclusion of nebular emission reproduces the red UV-to-optical color while keeping the UV color sufficiently blue. We infer that typical LAEs at z~6-7 are building blocks of galaxies seen at lower redshifts. We find a tentative decrease in the Lyman alpha escape fraction from z=5.7 to 6.6, which may imply an increase in the intergalactic medium neutral fraction. From the minimum contribution of nebular emission required to fit the observed SEDs, we place an upper limit on the escape fraction of ionizing photons to be f_esc^ion~0.6 at z=5.7 and ~0.9 at z=6.6. We also compare the stellar populations of our LAEs with that of stacked HST/WFC3 z-dropout galaxies.

We present deep spectroscopic data for a 24-object subsample of our full 41-object z~0.5 radio galaxy sample in order to investigate the evolution of the Fundamental Plane of radio galaxies. We find that the low-luminosity, FRI-type, radio galaxies in our sample are consistent with the local Fundamental Plane of radio galaxies defined by Bettoni et al. when corrected for simple passive evolution of their stellar populations. However, we find that the higher luminosity, FRII-type radio galaxies are inconsistent with the local Fundamental Plane if only passive evolution is considered, and find evidence for a rotation in the Fundamental Plane at z~0.5 when compared with the local relation. We show that neither passive evolution, nor a mass-dependent evolution in the mass-to-light ratio, nor an evolution in the size of the host galaxies can, by themselves, plausibly explain the observed tilt. However, we suggest that some combination of all three effects, with size evolution as the dominant factor, may be sufficient to explain the difference between the planes. We also find evidence for a correlation between host galaxy velocity dispersion and radio luminosity at the 97% significance level within our sub-sample, although further observations are required in order to determine whether this is different for the FRI and FRII radio sources. Assuming that the M_BH - sigma relation still holds at z~0.5, this implies that radio luminosity scales with black hole mass, in agreement with previous studies.

We measure the rest-frame colors (dust-corrected), infrared luminosities, star formation rates, and stellar masses of 92 galaxies in a Spitzer-selected cluster at z=1.62. By fitting spectral energy distributions (SEDs) to 10-band photometry (0.4 micron < lambda(obs) <8 micron) and measuring 24 micron fluxes for the 12 spectroscopically confirmed and 80 photometrically selected members, we discover an exceptionally high level of star formation in the cluster core of ~1700 Msun/yr per Mpc^2. The cluster galaxies define a strong blue sequence in (U-V) color and span a range in color. We identify 17 members with L(IR)>10^(11) Lsun, and these IR luminous members follow the same trend of increasing star formation with stellar mass that is observed in the field at z~2. Using rates derived from both the 24 micron imaging and SED fitting, we find that the relative fraction of star-forming members triples from the lowest to highest galaxy density regions, e.g. the IR luminous fraction increases from ~8% at Sigma~10 gal per Mpc^2 to ~25% at Sigma>100 gal per Mpc^2. The observed increase is a reversal of the well-documented trend at z<1 and signals that we have reached the epoch when massive cluster galaxies are still forming a substantial fraction of their stars.

Wide-field multi-object spectroscopy is a high priority for European astronomy over the next decade. Most 8-10m telescopes have a small field of view, making 4-m class telescopes a particularly attractive option for wide-field instruments. We present a science case and design drivers for a wide-field multi-object spectrograph (MOS) with integral field units for the 4.2-m William Herschel Telescope (WHT) on La Palma. The instrument intends to take advantage of a future prime-focus corrector and atmospheric-dispersion corrector that will deliver a field of view 2 deg in diameter, with good throughput from 370 to 1,000 nm. The science programs cluster into three groups needing three different resolving powers R: (1) high-precision radial-velocities for Gaia-related Milky Way dynamics, cosmological redshift surveys, and galaxy evolution studies (R = 5,000), (2) galaxy disk velocity dispersions (R = 10,000) and (3) high-precision stellar element abundances for Milky Way archaeology (R = 20,000). The multiplex requirements of the different science cases range from a few hundred to a few thousand, and a range of fibre-positioner technologies are considered. Several options for the spectrograph are discussed, building in part on published design studies for E-ELT spectrographs. Indeed, a WHT MOS will not only efficiently deliver data for exploitation of important imaging surveys planned for the coming decade, but will also serve as a test-bed to optimize the design of MOS instruments for the future E-ELT.

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

We present the results of our investigation into the stellar populations of 24 radio galaxies at z~0.5 drawn from four complete, low-frequency selected radio surveys. We use the strength of the 4000A break as an indicator of recent star formation, and compare this with radio luminosity, optical spectral classification and morphological classification. We find evidence of different star formation histories for high- and low-luminosity radio sources; our group of low radio luminosity sources (typically FRI-type sources) has systematically older stellar populations than the higher radio luminosity group. Our sample is also fairly well divided by optical spectral classification. We find that galaxies classified as having low excitation spectra (LEGs) possess older stellar populations than high excitation line objects (HEGs), with the HEGs showing evidence for recent star formation. We also investigate the link between radio morphology, as used by Owen & Laing (1989), and the stellar populations. We find that there is a preference for the "fat-double" sources to have older stellar populations than the "classical double" sources, although this is also linked to these sources lying predominantly in the LEG and HEG categories respectively. These results are consistent with the hypothesis that HEGs are powered by accretion of cold gas, which could be supplied, for example, by recent mergers, secular instabilities, or filamentary cold flows. These processes could also trigger star formation in the host galaxy. The host galaxies of the LEGs do not show evidence for recent star formation and an influx of cold gas, and are consistent with being powered by the accretion of the hot phase of the inter-stellar medium.

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

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 discovery imaging and spectroscopy for nine new z ~ 6 quasars found in the Canada-France High-z Quasar Survey (CFHQS) bringing the total number of CFHQS quasars to 19. By combining the CFHQS with the more luminous SDSS sample we are able to derive the quasar luminosity function from a sample of 40 quasars at redshifts 5.74 < z < 6.42. Our binned luminosity function shows a slightly lower normalisation and flatter slope than found in previous work. The binned data also suggest a break in the luminosity function at M_1450 approx -25. A double power law maximum likelihood fit to the data is consistent with the binned results. The luminosity function is strongly constrained (1 sigma uncertainty < 0.1 dex) over the range -27.5 < M_1450 < -24.7. The best-fit parameters are Phi(M_1450^*) = 1.14 x 10^-8 Mpc^-3 mag^-1, break magnitude M_1450^* = -25.13 and bright end slope beta = -2.81. However the covariance between beta and M_1450^* prevents strong constraints being placed on either parameter. For a break magnitude in the range -26 < M_1450^* < -24 we find -3.8 < beta < -2.3 at 95% confidence. We calculate the z = 6 quasar intergalactic ionizing flux and show it is between 20 and 100 times lower than that necessary for reionization. Finally, we use the luminosity function to predict how many higher redshift quasars may be discovered in future near-IR imaging surveys.

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 investigate the stellar populations of Lyman alpha emitters (LAEs) at z=3.1 and 3.7 in 0.65 deg^2 of the Subaru/XMM-Newton Deep Field, based on rest-frame UV-to-optical photometry obtained from the Subaru/XMM-Newton Deep Survey, the UKIDSS/Ultra Deep Survey, and the Spitzer legacy survey of the UKIDSS/UDS. Among a total of 302 LAEs (224 for z=3.1 and 78 for z=3.7), only 11 are detected in the K band, i.e., brighter than K(3sigma)=24.1 mag. Eight of the 11 K-detected LAEs are spectroscopically confirmed. We find that the K-undetected objects, which should closely represent the LAE population as a whole, have low stellar masses of ~ 10^8 - 10^8.5 Msun, modest SFRs of 1 - 100 Msun yr^-1, and modest dust extinction of E(B-V) < 0.2. The K-detected objects are massive, Mstar ~ 10^9 - 10^10.5 Msun, and have significant dust extinction with a median of E(B-V) ~= 0.3. Four K-detected objects with the reddest spectral energy distributions, two of which are spectroscopically confirmed, are heavily obscured with E(B-V) ~ 0.65, and their continua resemble those of some local ULIRGs. Interestingly, they have large Lyman alpha equivalent widths ~= 70 - 250 A. If these four are excluded, our sample has a weak anti-correlation between EW(Lya) and Mstar. We compare the stellar masses and the specific star formation rates (sSFR) of LAEs with those of Lyman-break galaxies, distant red galaxies, submillimetre galaxies, and I- or K-selected galaxies with z_phot ~ 3. We find that the LAE population is the least massive among all the galaxy populations in question, but with relatively high sSFRs, while NIR-detected LAEs have Mstar and sSFR similar to LBGs. Our reddest four LAEs have very high sSFRs in spite of large Mstar, thus occupying a unique region in the Mstar versus sSFR space. (abridged)

[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 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 present imaging and spectroscopic observations for six quasars at z>5.9 discovered by the Canada-France High-z Quasar Survey (CFHQS). The CFHQS contains sub-surveys with a range of flux and area combinations to sample a wide range of quasar luminosities at z~6. The new quasars have luminosities 10 to 75 times lower than the most luminous SDSS quasars at this redshift. The least luminous quasar, CFHQS J0216-0455 at z=6.01, has absolute magnitude M_1450=-22.21, well below the likely break in the luminosity function. This quasar is not detected in a deep XMM-Newton survey showing that optical selection is still a very efficient tool for finding high redshift quasars.

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 report here on the discovery of a red quasar, J004929.4+351025.7 at a redshift of z = 2.48, situated within a large Lya emission-line halo. The radio spectral energy distribution implies that the radio jets were triggered < 10^4 years prior to the time at which the object is observed, suggesting that the jet triggering of the active galactic nucleus is recent. The loosely biconical structure of the emission-line halo suggests that it is ionised by photons emitted by the central quasar nucleus and that the central nucleus is obscured by a dusty torus with Av ~ 3.0. The large spatial extent of the Lya halo relative to the radio emission means this could only have occurred if the radio jets emerged from an already established highly-accreting black hole. This suggests that the radio-jet triggering is delayed with respect to the onset of accretion activity on to the central supermassive black hole.

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

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.

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

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.

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 basic observational data on the 6C** sample. This is a new sample of radio sources drawn from the 151 MHz 6C survey, which was filtered with radio criteria chosen to optimize the chances of finding radio galaxies at z > 4. The filtering criteria are a steep-spectral index and a small angular size. The final sample consists of 68 sources from a region of sky covering 0.421 sr. We present VLA radio maps, and the results of K-band imaging and optical spectroscopy. Near-infrared counterparts are identified for 66 of the 68 sources, down to a 3-sigma limiting magnitude of K ~ 22 mag in a 3-arcsec aperture. Eight of these identifications are spatially compact, implying an unresolved nuclear source. The K-magnitude distribution peaks at a median K=18.7 mag, and is found to be statistically indistinguishable from that of the similarly selected 6C* sample, implying that the redshift distribution could extend to z > 4. Redshifts determined from spectroscopy are available for 22 (32 per cent) of the sources, over the range of 0.2 < z < 3.3 . We measure 15 of these, whereas the other 7 were previously known. Six sources are at z > 2.5. Four sources show broad emission lines in their spectra and are classified as quasars. Three of these show also an unresolved K-band identification. Eleven sources fail to show any distinctive emission and/or absorption features in their spectra. We suggest that these could be (i) in the so-called `redshift desert' region of 1.2 < z < 1.8, or (ii) at a greater redshift, but feature weak emission line spectra.

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

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

We have performed a deep spectroscopic survey of extremely red galaxies on the GOODS-South field, using GMOS on Gemini South. We present here spectra and redshifts for 16 ERGs at 0.87<z<2.02, to a limit K=20.2. The ERHs are a mixture of spheroidals, mergers and spirals, with one AGN. For at least 10 of these galaxies we observe [OII] emission lines. We perform an age-dating analysis by fitting the spectra and 9-band photometry of the ERGs with models of passively evolving stellar populations combined with a younger star-forming component. The best-fitting ages for the old stellar components range from 0.6 to 4.5 Gyr, with a mean 2.1 Gyr. Masses range from 3 to 20 times 10^10 solar masses. The star-forming component typically forms a few per cent of the total mass, with dust reddening averaging E(B-V)=0.35. Its timescale tends to be short for mergers (<50 Myr) and longer (200-800 Myr) for spiral ERGs.

We investigate the role of the luminous infrared galaxy (LIRG) and ultra-luminous infrared galaxy (ULIRG) phases in the evolution of Ks-selected galaxies and, in particular, Extremely Red Galaxies (ERGs). With this aim, we compare the properties of a sample of 2905 Ks<21.5 (Vega mag) galaxies in the GOODS/CDFS with the sub-sample of those 696 sources which are detected at 24 microns. We find that LIRGs constitute 30% of the galaxies with stellar mass M>1x10^{11} Msun assembled at redshift z=0.5. A minimum of 65% of the galaxies with M>2.5x10^{11} Msun at z~2-3 are ULIRGs at those redshifts. 60% of the ULIRGs in our sample have the characteristic colours of ERGs. Conversely, 40% of the ERGs with stellar mass M>1.3x10^{11} Msun at 1.5<z<2.0 and a minimum of 52% of those with the same mass cut at 2.0<z<3.0 are ULIRGs. The average optical/near-IR properties of the massive ERGs at similar redshifts that are identified with ULIRGs and that are not have basically no difference, suggesting that both populations contain the same kind of objects in different phases of their lives. LIRGs and ULIRGs have an important role in galaxy evolution and mass assembly, and, although they are only able to trace a fraction of the massive (M>1x10^{11} Msun) galaxies present in the Universe at a given time, this fraction becomes very significant (>50%) at redshifts z>~2.

It has long been believed that accretion onto supermassive black holes powers quasars, but there has been relatively few observational constraints on the spins of the black holes. We address this problem by estimating the average radiative efficiencies of a large sample of quasars selected from the Sloan Digital Sky Survey, by combining their luminosity function and their black hole mass function. Over the redshift interval $0.4<z<2.1$, we find that quasars have average radiative efficiencies of $\sim 30% - 35%$, strongly suggesting that their black holes are rotating very fast, with specific angular momentum $a \approx 1$, which stays roughly constant with redshift. The average radiative efficiency could be reduced by a factor of $\sim$2, depending on the adopted zeropoint for the black hole mass scale. The inferred large spins and their lack of significant evolution are in agreement with the predictions of recent semi-analytical models of hierarchical galaxy formation if black holes gain most of their mass through accretion.

In this paper we report new evidence that measurements of the broad-line widths in quasars are dependent on the source orientation, consistent with the idea that the broad-line region is flattened or disc-like. This reinforces the view derived from radio-selected samples, where the radio-core dominance has been used as a measure of orientation. The results presented here show a highly significant (>99.95%) correlation between radio spectral index (which we use as a proxy for source orientation) and broad-line width derived from the H-beta and MgII emission lines. This is the first time that this type of study has used quasars derived from a large optically selected quasar sample, where the radio-loud quasars (RLQs) and radio-quiet quasars (RQQs) have indistinguishable distributions in redshift, bolometric luminosity and colour, and therefore overcomes any biases which may be present in only selecting via radio emission. We find that the mean FWHM for the flat-spectrum radio-loud quasars (FSQs) to be <FWHM(FSQ)>= 4990 +/- 536km/s, which differs significant from the mean FWHM of the steep-spectrum radio-loud quasars (SSQs), where <FWHM(SSQ)>=6464 +/- 506km/s. We also find that the distribution in FWHM for the FSQs is indistinguishable from that of the radio-quiet quasars (RQQs), where <FWHM(RQQ)>= 4831 +/-25km/s. Considering other observational results in the literature we interpret this result in the context of a significant fraction of the FSQs being derived from the underlying RQQ population which have their radio flux Doppler boosted above the RLQ/RQQ divide. Under the assumption of a disc-like broad-line region we find no evidence for a difference in the average line-of-sight angle for RQQs and RLQs, implying that the difference is due to black-hole mass. (abridged)

We present the results of a study which uses the 3CRR sample of radio-loud active galactic nuclei to investigate the evolution of the black-hole:spheroid mass ratio in the most massive early-type galaxies from 0<z<2. Radio-loud unification is exploited to obtain virial (line-width) black-hole mass estimates from the 3CRR quasars, and stellar mass estimates from the 3CRR radio galaxies, thereby providing black-hole and stellar mass estimates for a single population of early-type galaxies. At low redshift (z<1) the 3CRR sample is consistent with a black-hole:spheroid mass ratio of M_bh/M_sph ~0.002, in good agreement with that observed locally for quiescent galaxies of similar stellar mass (M_sph ~5x10^11 M_sun). However, over the redshift interval 0<z<2 the 3CRR black-hole:spheroid mass ratio is found to evolve as M_bh/M_sph \propto (1+z)^{2.07\pm0.76}, reaching M_bh/M_sph ~0.008 by redshift z~2. This evolution is found to be inconsistent with the local black hole:spheroid mass ratio remaining constant at a moderately significant level (98%). If confirmed, the detection of evolution in the 3CRR black-hole:spheroid mass ratio further strengthens the evidence that, at least for massive early-type galaxies, the growth of the central supermassive black hole may be completed before that of the host spheroid.

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 present high-resolution imaging of a sample of 10 SCUBA submillimetre galaxies, obtained using the ACS camera on the Hubble Space Telescope. We find that the majority show compact, disturbed morphologies. Using quantitative morphological classification, we find that at least 6 are classified as major mergers. Simulations suggest that the morphological parameters are unlike local spirals and ellipticals, but similar to those of local ULIRGs. Compared to Lyman-break galaxies, the submillimetre-selected galaxies are on average more asymmetric, but also significantly more concentrated in their light distributions. This is consistent with a higher fraction of the stellar mass being contained within a central spheroid component. Despite their different morphologies, we find that submm galaxies are similar in size to luminous (>~L*) Lyman-break galaxies, with half-light radii in the range 2.8 - 4.6kpc.

We present deep Ks<21.5 (Vega) identifications, redshifts and stellar masses for most of the sources composing the bulk of the 24 micron background in the GOODS/CDFS. Our identified sample consists of 747 Spitzer/MIPS 24 micron objects, and includes ~94% of all the 24 micron sources in the GOODS-South field which have fluxes Snu(24)>83 microJy (the 80% completeness limit of the Spitzer/GTO 24 micron catalog). 36% of our galaxies have spectroscopic redshifts (mostly at z<1.5) and the remaining ones have photometric redshifts of very good quality, with a median of |dz|=|zspec-zphot|/(1+zspec)=0.02. We find that MIPS 24 micron galaxies span the redshift range z~0-4, and that a substantial fraction (28%) lie at high redshifts z>1.5. We determine the existence of a bump in the redshift distribution at z~1.9, indicating the presence of a significant population of galaxies with PAH emission at these redshifts. Massive (M>10^11 Msun) star-forming galaxies at redshifts 2<z<3 are characterized by very high star-formation rates (SFR>500 Msun/yr), and some of them are able to construct a mass of 10^10-10^11 Msun in a single burst lifetime (~0.01-0.1 Gyr). At lower redshifts z<2, massive star-forming galaxies are also present, but appear to be building their stars on long timescales, either quiescently or in multiple modest burst-like episodes. At redshifts z~1-2, the ability of the burst-like mode to produce entire galaxies in a single event is limited to some lower (M<7x10^10 Msun) mass systems, and it is basically negligible at z<1. Our results support a scenario where star-formation activity is differential with assembled stellar mass and redshift, and where the relative importance of the burst-like mode proceeds in a down-sizing way from high to low redshifts. (abridged)

Supermassive black holes in galactic nuclei show a close relationship between the black hole mass, luminosity, and stellar velocity dispersion (sigma_*) of the host galaxy bulge. Probing these relationships at high redshift may shed light on the link between the formation of the galactic bulge and central black hole, but direct measurements of stellar velocity dispersion at high redshift are difficult. We show that [OIII] line widths provide a useful surrogate for sigma_* by comparing sigma_[OIII] with the value of sigma_* predicted by the Faber-Jackson relation for QSOs with measured host galaxy luminosity. Over a wide range of AGN luminosity, sigma_[OIII] tracks sigma_*, albeit with considerable scatter. [OIII] line widths are narrower by ~0.1 dex in radio-loud QSOs than in radio-quiet QSOs of similar host luminosity. In low redshift QSOs, the ratio of star formation rate to black hole growth rate is much smaller than the typical ratio of bulge mass to black hole mass.

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

We have undertaken deep optical imaging observations of three 6.2<z<6.5 quasar fields in the i' and z' filters. These data are used to search for foreground galaxies which are gravitationally lensing the quasars and distant galaxies physically associated with the quasars. Foreground galaxies are found closer than 5 arcsec from the lines-of-sight of two of the three quasars. However, the faintness of these galaxies suggests they have fairly low masses and provide only weak magnifications (mu<1.1). No convincing galaxies physically associated with the quasars are found and the number of i'-band dropouts is consistent with that found in random fields. We consider the expected dark matter halo masses which host these quasars under the assumption that a correlation between black hole mass and dark matter halo mass exists. We show that the steepness of the high-mass tail of the halo mass function at this redshift, combined with realistic amounts of scatter in this correlation, lead to expected halo masses substantially lower than previously believed. This analysis can explain the lack of companion galaxies found here and the low dynamical mass recently published for one of the quasars.

Galaxy bulge luminosity L, black hole mass (Mbh), and stellar velocity dispersion (sigma_*) increase together in a way suggesting a close evolutionary relationship. Measurements of the Mbh-sigma_* relationship as a function of cosmic time may shed light on the origin of this relationship. Direct measurements of sigma_* at high redshift are difficult, and the width of the narrow emission lines of AGN has been proposed as a surrogate. We investigate the utility of using sigma_[OIII] for sigma_* by examining host galaxy magnitudes and [OIII] line widths for low redshift QSOs. For radio-quiet QSOs, sigma_[OIII] is consistent in the mean with the value of sigma_* predicted by the Faber-Jackson relation. For our limited range of host luminosity, scatter obscures the expected increase of sigma_[OIII] with host luminosity. However, for a sample of AGN covering a wide range of measured or inferred sigma_*, there is a clear increase of sigma_[OIII] with sigma_*. Radio-loud QSOs on average have sigma_* smaller by ~0.1 dex than radio-quiet QSOs of similar host luminsotiy, at least for luminosities typical of PG QSOs. Star formation rates in our low redshift QSOs are smaller than required to maintain the typical observed ratio of bulge mass to black hole mass.

We present results of a deep HST/WFPC2 imaging study of 17 quasars at z~0.4, designed to determine the properties of their host galaxies. The sample consists of quasars with absolute magnitudes in the range -24>M_V>-28, allowing us to investigate host galaxy properties across a decade in quasar luminosity, but at a single redshift. We find that the hosts of all the RLQs, and all the RQQs with nuclear luminosities M_V<-24, are massive bulge-dominated galaxies, confirming and extending the trends deduced from our previous studies. From the best-fitting model host galaxies we have estimated spheroid and black-hole masses, and the efficiency (with respect to Eddington luminosity) with which each quasar is radiating. The largest inferred black-hole mass in our sample is \~3.10^9 M_sun, comparable to those at the centres of M87 and Cygnus A. We find no evidence for super-Eddington accretion in even the most luminous objects. We investigate the role of scatter in the black-hole:spheroid mass relation in determining the ratio of quasar to host-galaxy luminosity, by generating simulated populations of quasars lying in hosts with a Schechter mass function. Within the subsample of the highest luminosity quasars, the observed variation in nuclear-host luminosity ratio is consistent with being the result of the scatter in the black-hole:spheroid relation. Quasars with high nuclear-host ratios can be explained by sub-Eddington accretion onto black holes in the high-mass tail of the black-hole:spheroid relation. Our results imply that, owing to the Schechter cutoff, host mass should not continue to increase linearly with quasar luminosity, at the very highest luminosities. Any quasars more luminous than M_V=-27 should be found in massive elliptical hosts which at the present day would have M_V ~ -24.5.

We have obtained deep, multi-band imaging observations around three of the most distant known quasars at redshifts z>6. Standard accretion theory predicts that the supermassive black holes present in these quasars were formed at a very early epoch. If a correlation between black hole mass and dark matter halo mass is present at these early times, then these rare supermassive black holes will be located inside the most massive dark matter halos. These are therefore ideal locations to search for the first clusters of galaxies. We use the Lyman-break technique to identify star-forming galaxies at high redshifts. Our observations show no overdensity of star-forming galaxies in the fields of these quasars. The lack of (dust-free) luminous starburst companions indicates that the quasars may be the only massive galaxies in their vicinity undergoing a period of intense activity.

Using a sample of more than 6000 quasars from the Sloan digital sky survey (SDSS) we compare the black-hole mass distributions of radio-loud and radio-quiet quasars. Based on the virial black-hole mass estimator the radio-loud quasars (RLQs) are found to harbour systematically more massive black holes than radio-quiet quasars (RQQs) with very high significance (>>99.99%), with mean black-hole masses of <\log(M_{bh}/\Msun)>=8.89\pm0.02 and <\log(M_{bh}/\Msun)>=8.69\pm0.01 for the RLQs and RQQs respectively. Crucially, the new RLQ and RQQ samples have indistinguishable distributions on the redshift-optical luminosity plane, excluding the possibility that either parameter is responsible for the observed black-hole mass difference. Moreover, this black-hole mass difference is shown to be in good agreement with the optical luminosity difference observed between RLQ and RQQ host galaxies at low redshift (i.e. \Delta M_{host}=0.4-0.5 magnitudes). Within the SDSS samples, black-hole mass is strongly correlated with both radio luminosity and the radio-loudness $\mathcal{R}$ parameter (>7 sigma significance), although the range in radio luminosity at a given black-hole mass is several orders of magnitude. It is therefore clear that the influence of additional physical parameters or evolution must also be invoked to explain the quasar radio-loudness dichotomy.

We have analysed 5-epoch GOODS HST-ACS B, V, I_775 and z datasets (V1.0 release), in conjunction with existing VLT-ISAAC imaging in the J, H and Ks bands, to derive estimated redshifts for the sample of 198 Extremely Red Galaxies (ERGs) with Ks<22 (Vega) and I_775-Ks>3.92 selected by Roche et al.(2003) from 50.4 sq.arcmin of the GOODS/CDFS field. We find that, at this depth, the ERG population spans the redshift range 0.5<z_phot<4.75 and over two decades in mass (~ 3x10^9 M_sun to ~3x10^11 M_sun). We explore the evolution of the ERG luminosity function (LF) from redshifts <z_phot>=1.0 to <z_phot>=2.5 and compare it with the global Ks-band LF at redshifts 1<z_phot<2. We find that the bright end of the ERG LF does not decrease from redshifts <z_phot>=2.0 to <z_phot>=2.5 and we connect this fact with the presence of progenitors of the local L>L* population at redshifts z_phot>2. We determine lower limits of rho_c=(6.1 +/- 1.9) x 10^{-5} Mpc^{-3} and rho_c=(2.1 +/- 1.1) x 10^{-5}Mpc^{-3} on the comoving densities of progenitors of local massive galaxies already assembled at redshifts <z_phot>=2.5 and <z_phot>=3.5, respectively. We have investigated the existence of high-redshift Lyman break galaxies massive enough to be included in this ERG sample. Out of an initial list of 12 potential very high redshift candidates, we have identified 2 ERGs which have a high probability of lying at z_phot>4. We discuss the advantages of multi-colour to single-colour selection techniques in obtaining reliable lists of very high-redshift candidate sources, and present revised lower redshift estimates for sources previously claimed as potential z> 5 dropouts in recent studies.(abridged)

Virial black-hole mass estimates are presented for 12698 quasars in the redshift interval 0.1<z<2.1, based on modelling of spectra from the Sloan Digital Sky Survey (SDSS) first data release . The black-hole masses of the SDSS quasars are found to lie between $\simeq10^{7}\Msun$ and an upper limit of $\simeq 3\times 10^{9}\Msun$, entirely consistent with the largest black-hole masses found to date in the local Universe. The estimated Eddington ratios of the broad-line quasars (FWHM \geq2000 km s^{-1}) show a clear upper boundary at L_{bol}/L_{Edd}~1, suggesting that the Eddington luminosity is still a relevant physical limit to the accretion rate of luminous broad-line quasars at $z\leq 2$. By combining the black-hole mass distribution of the SDSS quasars with the 2dF quasar luminosity function, the number density of active black holes at $z\simeq 2$ is estimated as a function of mass. By comparing the estimated number density of active black holes at $z\simeq 2$ with the local mass density of dormant black holes, we set lower limits on the quasar lifetimes and find that the majority of black holes with mass $\geq 10^{8.5}\Msun$ are in place by $\simeq 2$.

Virial black-hole mass estimates are presented for 12698 quasars in the redshift interval 0.1<z<2.1, based on modelling of spectra from the Sloan Digital Sky Survey (SDSS) first data release. The black-hole masses of the SDSS quasars are found to lie between $\simeq10^{7}\Msun$ and an upper limit of $\simeq 3\times 10^{9}\Msun$, entirely consistent with the largest black-hole masses found to date in the local Universe. The estimated Eddington ratios of the broad-line quasars (FWHM\geq 2000 kms^{-1}) show a clear upper boundary at $L_{bol}/L_{Edd}\simeq 1$, suggesting that the Eddington luminosity is still a relevant physical limit to the accretion rate of luminous broad-line quasars at $z\leq 2$. By combining the black-hole mass distribution of the SDSS quasars with the 2dF quasar luminosity function, the number density of active black holes at $z\simeq 2$ is estimated as a function of mass. In addition, we independently estimate the local black-hole mass function for early-types using the $M_{bh}-\sigma$ and $M_{bh}-L_{bulge}$ correlations. Based on the SDSS velocity dispersion function and the 2MASS $K-$band luminosity function, both estimates are found to be consistent at the high-mass end ($M_{bh}\geq 10^{8}\Msun$). By comparing the estimated number density of active black holes at $z\simeq 2$ with the local mass density of dormant black holes, we set lower limits on the quasar lifetimes and find that the majority of black holes with mass $\geq 10^{8.5}\Msun$ are in place by $\simeq 2$.

The results of analysis of HST I-band imaging of a sample of 41 radio galaxies spanning three orders of magnitude in radio luminosity at redshift z~0.5 are presented. The full sample is drawn from four complete, low-frequency selected radio samples with progressively fainter flux-density limits (3CRR, 6CE, 7CRS and the new TexOx-1000 sample). Modelling of the HST imaging data shows that the host galaxies have surface-brightness distributions consistent with those expected for classic ellipticals (Sersic parameter, beta~0.25), with beta in the range 0.17<beta<0.30, and a mean of <beta>=0.23\pm0.01. The host-galaxy luminosities are comparable with those of galaxies drawn from the bright end of the local cluster galaxy luminosity function, spanning the range 0.7L^{\star}<L<10L^{\star}, with a mean of 3.2\pm0.3 L^{\star}. In addition, the radio galaxies are shown to follow a Kormendy relation indistinguishable from that of powerful low-redshift radio galaxies. Combining our new results with those in the literature it is found that the scalelengths and Kormendy relations of 3C-class radio galaxies do not vary significantly over the redshift range 0.0<z<0.8, providing no evidence for dynamical evolution of this class of host galaxy within this redshift interval. Converting the host-galaxy luminosities into black-hole mass estimates predicts that the radio galaxies harbour central black holes with masses in the range 10^{8.1}\Msun <M_{bh}< 10^{9.5}\Msun, with a geometric mean of <M_{bh}>=10^{8.87\pm 0.04}\Msun. Finally, a significant (~3 sigma) correlation is found between black-hole mass and 151-MHz radio luminosity for those objects in the sample with either high-excitation nuclear spectra or classical double radio structures. (abridged)

We present the photometric redshift distribution of a sample of 198 Extremely Red Galaxies (ERGs) with Ks<22 and I-Ks>3.92 (Vega), selected by Roche et al. in 50.4 sq. arcmin of the Chandra Deep Field South (CDFS). The sample has been obtained using ISAAC-VLT and ACS-HST GOODS public data. We also show the results of a morphological study of the 72 brightest ERGs in the z band (z<25, AB).

We highlight the importance of the MgII emission-line doublet in probing high-redshift quasars and their supermassive black holes. In the SDSS era, where large scale investigations of quasars across the age of the Universe are possible, this emission-line has the ability to provide accurate systemic redshifts which are important for a variety of follow-up studies, as well as probe the masses of the supermassive black holes that power these phenomena.

We present near-infrared H and K-band spectra of the z=6.41 quasar SDSS J114816.64+525150.3. The spectrum reveals a broad MgII 2799 emission line with a full-width half-maximium of 6000 km/s. From the peak wavelength of this emission line we obtain a more accurate redshift than is possible from the published optical spectrum and determine a redshift of z=6.41+/-0.01. If the true peak of the Lyman alpha emission is at the same redshift, then a large fraction of the flux blueward of the peak is absorbed. The equivalent width of the MgII emission line is similar to that of lower redshift quasars, suggesting that the UV continuum is not dominated by a beamed component. Making basic assumptions about the line-emitting gas we derive an estimate for the central black hole in this quasar of 3x10^9 solar masses. The very high luminosity of the quasar shows that it is accreting at the maximal allowable rate for a black hole of this mass adopting the Eddington limit criterion.

We present the results of deep VLT-ISAAC Ks-band imaging of four z~1.5 flat-spectrum quasars selected from the Parkes half-Jansky flat spectrum sample. We find that the hosts of these flat-spectrum quasars are consistent with lying on the K-z Hubble relation for radio galaxies. This implies that the flat-spectrum quasar hosts fall in line with the expectations from orientation based unified schemes and also that they contain black holes of similar mass. Moreover, the width of the Hbeta broad emission line in these objects tends to be narrower than in their misaligned (low-frequency selected quasar) counterparts, implying that the width of the Hbeta broad emission line depends on source inclination, at least for radio-loud quasars, in line with previous studies.

In this proceedings we summarize our recent letter (Jarvis & McLure 2002) where we suggested that by correcting for the inevitable effects of inclination, the black-hole masses of flat-spectrum quasars (FSQ) with intrinsically powerful radio jets are confined, virtually exclusively, to M_BH >10^{8} Msolar. After considering realistic Doppler boosting factors, many of the FSQ would be more accurately classified as radio-intermediate or radio-quie t quasars. This range in radio luminosity suggests that the FSQ are fully consistent with an upper boundary on radio power of the form L_5GHz proportional to M_BH^2.5.

A reliable method for estimating the black-hole masses of high-redshift quasars would provide crucial new information for understanding the nature and cosmological evolution of quasars. In this proceedings we summarize the results of our recent paper (McLure & Jarvis 2002) which provides a virial black-hole mass estimator based on rest-frame UV observables, thereby allowing reliable black-hole mass estimates to be obtained out to redshifts of $z\sim2.5$ from optical spectra alone.

In this review the correlations between host galaxy properties and black-hole mass determined from nearby quiescent galaxies are briefly discussed, before proceeding to investigate their application to active galactic nuclei (AGN). The recent advances in estimating the black-hole masses of AGN are then reviewed, followed by an examination of the connection between black-hole mass and radio luminosity.

We present the final results from our deep HST imaging study of the hosts of radio-quiet quasars (RQQs), radio-loud quasars (RLQs) and radio galaxies (RGs). We describe new WFPC2 R-band observations for 14 objects and model these images in conjunction with the data already reported in McLure et al (1999). We find that spheroidal hosts become more prevalent with increasing nuclear luminosity such that, for nuclear luminosities M_V < -23.5, the hosts of both radio-loud and radio-quiet AGN are virtually all massive ellipticals. Moreover we demonstrate that the basic properties of these hosts are indistinguishable from those of quiescent, evolved, low-redshift ellipticals of comparable mass. This result kills any lingering notion that radio-loudness is determined by host-galaxy morphology, and also sets severe constraints on evolutionary schemes which attempt to link low-z ULIRGs with RQQs. Instead, we show that our results are as expected given the relationship between black-hole and spheroid mass established for nearby galaxies, and apply this relation to estimate the mass of the black hole in each object. The results agree very well with completely-independent estimates based on nuclear emission-line widths; all the quasars in our sample have M(bh) > 5 x 10^8 solar masses, while the radio-loud objects are confined to M(bh) > 10^9 solar masses. This apparent mass-threshold difference, which provides a natural explanation for why RQQs outnumber RLQs by a factor of 10, appears to reflect the existence of a minimum and maximum level of black-hole radio output which is a strong function of black-hole mass. Finally, we use our results to estimate the fraction of massive spheroids/black-holes which produce quasar-level activity. This fraction is \~0.1% at the present day, rising to > 10% at z = 2-3.

We use a simple model of spheroid formation to explore the relationship between the creation of stars and dust in a massive proto galaxy and the growth of its central black hole. This model predicts that sub-mm luminosity peaks after only ~ 0.2 Gyr. However, without a very massive seed black hole, Eddington-limited growth means that a black hole mass of 10^9 solar masses, and hence very luminous AGN activity, cannot be produced until > 0.5 Gyr after the formation of the first massive stars in the halo. The model thus predicts a time-lag between the peak of sub-mm luminosity and AGN luminosity in a massive proto-elliptical of a few times 10^8 years. For a formation redshift z = 5, this means that powerful AGN activity is delayed until z = 3.5, by which time star formation in the host is 90% complete, and sub-mm luminosity has declined to ~ 25% of its peak value. This provides a natural explanation for why successful sub-mm detections of luminous radio galaxies are largely confined to z > 2.5. Conversely the model also predicts that while all high-redshift luminous sub-mm-selected sources should contain an active (and growing) black hole, the typical luminosity of the AGN in such objects is ~1000 times smaller than that of the most powerful AGN. This is consistent with the almost complete failure to detect sub-mm selected galaxies with existing X-ray surveys. Finally the model yields a black-hole:spheroid mass ratio which evolves rapidly in the first Gyr, but asymptotes to ~ 0.001-0.003 in agreement with results at low redshift. This ratio arises not because the AGN terminates star formation, but because fueling of the massive black hole is linked to the total mass of gas available for star formation in the host.

A new analysis of the connection between black-hole mass and radio luminosity in radio-selected flat-spectrum quasars (FSQ) is presented. In contrast to recent claims in the literature, we find no evidence that the black-hole masses of radio-selected FSQ are systematically lower that those of luminous optically-selected radio-loud quasars. The black-hole masses of the FSQ are estimated via the virial black-hole mass estimator which utilizes the line-width of the H$\beta$ emission line as a tracer of the central gravitational potential. By correcting for the inevitable effects of inclination, incurred due to the FSQ being viewed close to the line of sight, we find that the black-hole masses of the FSQ with intrinsically powerful radio jets are confined, virtually exclusively, to M_bh > 10^8 M_solar. This is in good agreement with previous studies of optically selected FSQ and steep-spectrum radio-loud quasars. Finally, following the application of a realistic Doppler boosting correction, we find that the FSQ occupy a wide range in intrinsic radio luminosity, and that many sources would be more accurately classified as radio-intermediate or radio-quiet quasars. This range in radio luminosity suggests that the FSQ are fully consistent with an upper boundary on radio power of the form L_5GHz \propto M_bh^2.5.

Deep spectroscopy of the two mJy radio galaxies LBDS 53W069 and 53W091 has previously shown them to have old (>3 Gyr) stellar populations at z=1.5. Here we present the results of HST observations of these galaxies. 53W069 has an r^1/4 profile in both the F814W & F110W data with an effective radius of 0.30 arcsec (2.7 kpc). The restframe U-B colour gradient requires a stellar population of super-solar (3Z_sun) metallicity that formed on a very short timescale at high redshift (z>5). 53W091 has a regular r^1/4 profile in F110W with an effective radius of 0.32 arcsec (2.9 kpc). The F814W profile is more extended and is consistent with a blue exponential disk that contributes 20% of the flux within r_e. The restframe U-B colour gradient is significantly larger than that observed in field ellipticals at z<1, implying a stellar population of mixed metallicity (1-3Z_sun) that formed in a high-redshift rapid burst. We have compared these two LBDS radio galaxies with the Kormendy relation of ten 3CR radio galaxies at z=0.8. The LBDS galaxies follow the same relation as the more radio-luminous 3CR galaxies, assuming passive evolution of their stellar populations, although they are smaller than the 3CR galaxies whose mean effective radius is 12 kpc. Their sizes and radio luminosities are consistent with scaling relations applied to the 3CR galaxies, in which both radio power and effective radius scale with galaxy mass. Compared with a sample of z=0.4 cluster ellipticals, 53W069 & 53W091 lie well within the scatter of the Kormendy relation. We conclude that the hosts of these mJy radio sources at z=1.5 are passively-evolving elliptical galaxies that will evolve into ordinary L* ellipticals by the present day. (Abridged)

A new technique is presented for determining the black-hole masses of high-redshift quasars from optical spectroscopy. The new method utilizes the full-width half maximum (FWHM) of the low-ionization MgII emission line and the correlation between broad-line region (BLR) radius and continuum luminosity at 3000 Angstroms. Using archival UV spectra it is found that the correlation between BLR radius and 3000 Angstrom luminosity is tighter than the established correlation with 5100 Angstrom luminosity. Furthermore, it is found that the correlation between BLR radius and 3000 Ang luminosity is consistent with a relation of the form $R_{BLR} \propto lambda L_{lambda}^{0.5}$, as expected for a constant ionization parameter. Using a sample of objects with broad-line radii determined from reverberation mapping it is shown that the FWHM of MgII and H beta are consistent with following an exact one-to-one relation, as expected if both H beta and MgII are emitted at the same radius from the central ionizing source. The resulting virial black-hole mass estimator based on rest-frame UV observables is shown to reproduce black-hole mass measurements based on reverberation mapping to within a factor of 2.5 (1 sigma). Finally, the new UV black-hole mass estimator is shown to produce identical results to the established optical (H beta) estimator when applied to 128 intermediate-redshift (0.3<z<0.9) quasars drawn from the Large Bright Quasar Survey and the radio-selected Molonglo quasar sample. We therefore conclude that the new UV virial black-hole mass estimator can be reliably used to estimate the black-hole masses of quasars from z~0.25 through to the peak epoch of quasar activity at z~2.5 via optical spectroscopy alone.

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

We present high-resolution imaging in the rest-frame optical of the weak radio source LBDS53W091. Previous optical spectroscopy has shown that this object has an evolved stellar population of age >3Gyr at z=1.55, determined from the amplitude of rest-frame UV spectral breaks. We have obtained deep Hubble Space Telescope imaging over 10 orbits with NICMOS camera 2, using the F160W H-band filter (1.6microns) which is a good approximation to the rest-frame R-band. Our observations reveal a radial light profile which is well fit by a de Vaucouleurs' r^{1/4} law, with a scale length of r_e=0.3" (3kpc/h_{50} for Omega_M=0.3 and Omega_{Lambda}=0). The elliptical morphology of the radio galaxy indicates a dynamically-evolved old system, consistent with the spectroscopic results. Some surrounding objects lie on the color:magnitude relation for a cluster at z=1.55, and are likely to be associated. The group of galaxies are somewhat more luminous than the fundamental plane of ellipticals at z=0 subject to (1+z)^4 surface brightness dimming, but are consistent with estimates of the fundamental plane at high redshift subject to passive luminosity evolution from a formation epoch of z>3.

We highlight some of the principal results from our recent Hubble Space Telescope studies of quasars and radio galaxies. The hosts of these powerful AGN are normal massive ellipticals which lie on the region of the fundamental plane populated predominantly by massive ellipticals with boxy isophotes and distinct cores. The hosts of the radio-loud sources are on average 1.5 times brighter than their radio-quiet counterparts and appear to lie above a mass threshold of 4 x 10^11 solar masses. This suggests that black holes more massive than 5 x 10^8 solar masses are required to produce a powerful radio source. However we show that this apparent threshold appears to be a consequence of an upper bound on radio output which is a strong function of black-hole mass, L_5GHz proportional to M_bh^2.5. This steep mass dependence can explain why the hosts of the most powerful radio sources are good standard candles. Such objects were certainly fully assembled by z = 1, and appear to have formed the bulk of their stars prior to z = 3.

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

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

The central black-hole masses of a sample of 30 luminous quasars are estimated using H-beta FWHM measurements from a combination of new and previously-published nuclear spectra. The quasar black-hole mass estimates are combined with reverberation-mapping measurements for a sample of Seyfert galaxies (Wandel 1999) in order to study AGN black-hole masses over a wide range in nuclear luminosity. It is found that black-hole mass and bulge luminosity are well correlated and follow a relation consistent with that expected if black-hole and bulge mass are directly proportional. Contrary to the results of Wandel (1999) no evidence is found that Seyfert galaxies follow a different Mblack-Mbulge relation to quasars. However, the black-hole mass distributions of the radio-loud and radio-quiet quasar sub-samples are found to be significantly different, with the median black-hole mass of the radio-loud quasars a factor of three larger than their radio-quiet counterparts. Finally, utilizing the elliptical galaxy fundamental plane to provide stellar velocity dispersion estimates, a comparison is performed between the virial H-beta black-hole mass estimates and those of the Mblack-sigma correlations of Gebhardt et al. (2000a) and Merritt & Ferrarese (2000). With the disc-like geometry of the broad-line region adopted in this paper, the virial H-beta black-hole masses indicate that the correct normalization of the black-hole vs. bulge mass relation is Mblack=0.0025Mbulge, while the standard assumption of purely random broad-line velocities leads to Mblack=0.0008Mbulge. The normalization of Mblack=0.0025Mbulge provided by the disc model is in remarkably good agreement with that inferred for our quasar sample using the (completely independent) Mblack-sigma correlations.

In this paper, the third and final of a series, we present complete K-band imaging and some complementary I-band imaging of the filtered 6C* sample. We find no systematic differences between the K-z relation of 6C* radio galaxies and those from complete samples, so the near-infrared properties of luminous radio galaxies are not obviously biased by the additional 6C* radio selection criteria (steep spectral index and small angular size). The 6C* K-z data significantly improve delineation of the K-z relation for radio galaxies at high-redshift (z >2). Accounting for non-stellar contamination, and for correlations between radio luminosity and stellar mass, we find little support for previous claims that the underlying scatter in the stellar luminosity of radio galaxies increases significantly at z >2. In a particular spatially-flat universe with a cosmological constant, the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (~5 Lstar), and a low dispersion (sigma ~ 0.5 mag) which formed their stars at epochs corresponding to z >~2.5. This result is in line with recent sub-mm studies of high-redshift radio galaxies and the inferred ages of extremely red objects from faint radio samples.

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

The results of a comprehensive study of the cluster environments and black-hole masses of an optically matched sample of radio-loud and radio-quiet quasars are presented. No evidence is found for a difference in large-scale environments, with both quasar classes found to be located in clusters of Abell class $\sim0$ . Conversely, virial black-hole mass estimates based on H$_{\beta}$ line-widths show a clear difference in the quasar black-hole mass distributions. Our results suggest that a black-hole mass of $\sim10^{9}\Msolar$ is required to produce a powerful radio-loud quasar, and that it is black-hole mass and accretion rate which hold the key to the radio-loudness dichotomy.

We present the results of a deep K-band imaging study which reveals the host galaxies around a sample of luminous radio-quiet quasars. The K-band images, obtained at UKIRT, are of sufficient quality to allow accurate modelling of the underlying host galaxy. Initially, the basic structure of the hosts is revealed using a modified Clean deconvolution routine optimised for this analysis. 2 of the 14 quasars are shown to have host galaxies with violently disturbed morphologies which cannot be modelled by smooth elliptical profiles. For the remainder of our sample, 2D models of the host and nuclear component are fitted to the images using the chi-squared statistic to determine goodness of fit. Host galaxies are detected around all of the quasars. The reliability of the modelling is extensively tested, and we find the host luminosity to be well constrained for 9 quasars. The derived average K-band absolute K-corrected host galaxy magnitude for these luminous radio-quiet quasars is <M_K>=-25.15+/-0.04, slightly more luminous than an L* galaxy. The spread of derived host galaxy luminosities is small, although the spread of nuclear-to-host ratios is not. These host luminosities are shown to be comparable to those derived from samples of quasars of lower total luminosity and we conclude that there is no correlation between host and nuclear luminosity for these quasars. Nuclear-to-host ratios break the lower limit previously suggested from studies of lower nuclear luminosity quasars and Seyfert galaxies. Morphologies are less certain but, on the scales probed by these images, some hosts appear to be dominated by spheroids but others appear to have disk-dominated profiles.

A description is given of the method used to extract quasar host-galaxy parameters from the deep Hubble Space Telescope quasar images presented by McLure et al. (1999) and Dunlop et al. (2000). We then give the results of extensive testing of this technique on a wide range of simulated quasar+host combinations spanning the redshift range of our HST study (0.1 < z < 0.3). These simulations demonstrate that, when applied to our deep HST images, our method of analysis can easily distinguish the morphological type of a given host galaxy, as well as determining its scalelength, luminosity, axial ratio and position angle to within an accuracy of a few percent. We also present new infrared tip-tilt images of 4 of the most luminous quasars in our {\sc hst} sample, along with the outcome of modelling these data in a similar manner. The results provide further confidence in the accuracy of the derived host-galaxy scalelengths, and allow accurate determination of $R-K$ colours for this subset of sources. All 4 of these quasar host galaxies have very similar red colours, $R-K = 2.9 \pm 0.2$, indicative of a well-evolved stellar population.

The spatial clustering amplitude (B_{gq}) is determined for a sample of 44 powerful AGN at z~0.2. No significant difference is detected in the richness of the cluster environments of the radio-loud and radio-quiet sub-samples, both of which typically inhabit environments as rich as Abell Class ~0. Comparison with radio luminosity-matched samples from Hill & Lilly (1991) and Wold et al. (2000a) suggests that there is no epoch-dependent change in environment richness out to at least z>=0.5 for either radio galaxies or radio quasars. Comparison with the APM cluster survey shows that, contrary to current folklore, powerful AGN do not avoid rich clusters, but rather display a spread in cluster environment which is perfectly consistent with being drawn at random from the massive elliptical population. Finally, we argue that virtually all Abell class ~0 clusters contained an active galaxy during the epoch of peak quasar activity at z~2.5.

We present the results of a comprehensive re-analysis of the images of a virtually complete sample of 28 powerful 3CR radio galaxies with redshifts 0.6<z<1.8 from the HST archive. Using a two-dimensional modelling technique we have derived scalelengths and absolute magnitudes for a total of 16 3CR galaxies with a median redshift of z=0.8. Our results confirm the basic conclusions of Best, Longair & R\"{o}ttgering (1997, 1998) in that we also find z=1 3CR galaxies to be massive, well-evolved ellipticals, whose infrared emission is dominated by starlight. However, we in fact find that the scalelength distribution of 3CR galaxies at z \simeq 1 is completely indistinguishable from that derived for their low-redshift counterparts from our own recently-completed HST study of AGN hosts at z \simeq 0.2. There is thus no evidence that 3CR radio galaxies at z \simeq 1 are dynamically different from 3CR galaxies at low redshift. Moreover, for a 10-object sub-sample we have determined the galaxy parameters with sufficient accuracy to demonstrate, for the first time, that the z \simeq 1 3CR galaxies follow a Kormendy relation which is indistinguishable from that displayed by low-redshift ellipticals if one allows for purely passive evolution. The implied rather modest level of passive evolution since z \simeq 1 is consistent with that predicted from spectrophotometric models provided one assumes a high formation redshift (z \ge 4) within a low-density Universe. We conclude that there is no convincing evidence for significant dynamical evolution among 3CR galaxies in the redshift interval 0<z<1, and that simple passive evolution remains an acceptable interpretation of the K-z relation for powerful radio galaxies.

We present the first results from a major HST WFPC2 imaging study aimed at providing the first statistically meaningful comparison of the morphologies, luminosities, scalelengths and colours of the host galaxies of radio-quiet quasars, radio-loud quasars, and radio galaxies. We describe the design of this study and present the images which have been obtained for the first half of our 33-source sample. We find that the hosts of all three classes of luminous AGN are massive elliptical galaxies, with scalelengths ~=10 kpc, and R-K colours consistent with mature stellar populations. Most importantly this is the the first unambiguous evidence that, just like radio-loud quasars, essentially all radio-quiet quasars brighter than M_R = -24 reside in massive ellipticals. This result removes the possibility that radio `loudness' is directly linked to host galaxy morphology, but is however in excellent accord with the black-hole/spheroid mass correlation recently highlighted by Magorrian et al. (1998). We apply the relations given by Magorrian et al. to infer the expected Eddington luminosity of the putative black hole at the centre of each of the spheroidal host galaxies we have uncovered. Comparison with the actual nuclear R-band luminosities suggests that the black holes in most of these galaxies are radiating at a few percent of the Eddington luminosity; the brightest host galaxies in our low-z sample are capable of hosting quasars with M_R = -28, comparable to the most luminous quasars at z = 3. Finally we discuss our host-derived black-hole masses in the context of the radio-luminosity:black-hole mass correlation recently uncovered for nearby galaxies by Franceschini et al. (1998), and the resulting implications for the physical origin of radio loudness.