We measure star-formation rates (SFRs) and specific SFRs (SSFRs) of Ks-selected galaxies from the VIDEO survey by stacking 1.4-GHz Very Large Array data. We split the sample, which spans 0 < z < 3 and stellar masses 10**8.0 < Mstellar/Msol < 10**11.5, into elliptical, irregular or starburst galaxies based on their spectral-energy distributions. We find that SSFR falls with stellar mass, in agreement with the `downsizing' paradigm. We consider the dependence of the SSFR-mass slope on redshift: for our full and elliptical samples the slope flattens, but for the irregular and starburst samples the slope is independent of redshift. The rate of SSFR evolution reduces slightly with stellar mass for ellipticals, but irregulars and starbursts co-evolve across stellar masses. Our results for SSFR as a function of stellar mass and redshift are in agreement with those derived from other radio-stacking measurements of mass-selected passive and star-forming galaxies, but inconsistent with those generated from semi-analytic models, which tend to underestimate SFRs and SSFRs. There is a need for deeper high-resolution radio surveys such as those from telescopes like the next-generation MeerKAT in order to probe lower masses at earlier times and to permit direct detections, i.e. to study individual galaxies in detail.

There is an ongoing debate concerning the driver of nuclear activity in galaxies, with AGN either being triggered by major or minor galactic mergers or, alternatively, through secular processes like cooling gas accretion and/or formation of bars. We investigate the close environment of active galaxies selected in the X-ray, the radio, and the mid-IR. We utilise the first data release of the new near-IR VISTA Deep Extragalactic Observations (VIDEO) survey of the XMM-Large Scale Structure (LSS) field. We use two measures of environment density, namely counts within a given aperture and a finite redshift slice (pseudo-3D density) and closest neighbour density measures ${\Sigma}_{2}$ and ${\Sigma}_{5}$. We select both AGN and control samples, matching them in redshift and apparent Ks-band magnitude. We find that AGN are found in a range of environments, with a subset of the AGN samples residing in over-dense environments. Seyfert-like X-ray AGN and flat-spectrum radio-AGN are found to inhabit significantly over-dense environments compared to their control sample. The relation between over-densities and AGN luminosity does not however reveal any positive correlation. Given the absence of an environment density-AGN luminosity relation,we find no support for a scheme where high luminosity AGN are preferentially triggered by mergers. On the contrary, we find that AGN likely trace over dense environments at high redshift due to the fact that they inhabit the most massive galaxies, rather than being an AGN.

We present a stellar mass-metallicity relation at z~1.4 with an unprecedentedly large sample of ~340 star-forming galaxies obtained with FMOS on the Subaru Telescope. We observed K-band selected galaxies at 1.2 < z_{ph} < 1.6 in the SXDS/UDS fields with M_{*} > 10^{9.5} M_{\sun}, and expected F(H\alpha) > 5 \times 10^{-17} erg s^{-1} cm^{-2}. Among the observed ~1200 targets, 343 objects show significant H\alpha emission lines. The gas-phase metallicity is obtained from [NII]\lambda 6584/H\alpha line ratio, after excluding possible active galactic nuclei (AGNs). Due to the faintness of the [NII]\lambda 6584 lines, we apply the stacking analysis and derive the mass-metallicity relation at z~1.4. Our results are compared to past results at different redshifts in the literature. The mass-metallicity relation at z~1.4 is located between those at z~0.8 and z~2.2; it is found that the metallicity increases with decreasing redshift from z~3 to z~0 at fixed stellar mass. Thanks to the large size of the sample, we can study the dependence of the mass-metallicity relation on various galaxy physical properties. The average metallicity from the stacked spectra is close to the local FMR in the higher metallicity part but >0.1 dex higher in metallicity than the FMR in the lower metallicity part. We find that galaxies with larger E(B-V), B-R, and R-H colours tend to show higher metallicity by ~0.05 dex at fixed stellar mass. We also find relatively clearer size dependence that objects with smaller half light radius tend to show higher metallicity by ~0.1 dex at fixed stellar mass, especially in the low mass part.

We use galaxies from the Herschel-ATLAS survey, and a suite of ancillary simulations based on an isothermal dust model, to study our ability to determine the effective dust temperature, luminosity and emissivity index of 250um selected galaxies in the local Universe (z < 0.5). As well as simple far-infrared SED fitting of individual galaxies based on chi^2 minimisation, we attempt to derive the best global isothermal properties of 13,826 galaxies with reliable optical counterparts and spectroscopic redshifts. Using our simulations, we highlight the fact that applying traditional SED fitting techniques to noisy observational data in the Herschel Space Observatory bands introduces artificial anti-correlation between derived values of dust temperature and emissivity index. This is true even for galaxies with the most robust detections in our sample, making the results hard to interpret. We apply a method to determine the best-fit global values of isothermal effective temperature and emissivity index for z < 0.5 galaxies in H-ATLAS, deriving T = 22.3 +/- 0.1K and Beta = 1.98 +/- 0.02 (or T = 23.5 +/- 0.1K and Beta = 1.82 +/- 0.02 if we attempt to correct for bias by assuming that T and Beta are independent and normally distributed). We use our technique to test for an evolving emissivity index, finding only weak evidence. The median dust luminosity of our sample is log(Ldust/Lsolar) = 10.72 +/- 0.05 which (unlike T) shows little dependence on the choice of Beta used in our analysis, including whether it is variable or fixed. We use a further suite of simulations to emphasise the importance of the H-ATLAS PACS data for deriving dust temperatures at these redshifts, even though they are less sensitive than the SPIRE data. The majority of galaxies detected by H-ATLAS are normal star-forming galaxies, though a substantial minority (~31%) fall in the Luminous Infrared Galaxy category. (Abridged)

It has been speculated that low luminosity radio-loud AGN have the potential to serve as an important source of AGN feedback, and may be responsible for suppressing star-formation activity in massive elliptical galaxies at late times. As such the cosmic evolution of these sources is vitally important to understanding the significance of such AGN feedback processes and their influence on the global star-formation history of the universe. In this paper we present a new investigation of the evolution of faint radio sources out to $z{\sim}2.5$. We combine a 1 square degree VLA radio survey, complete to a depth of 100 $\mu$Jy, with accurate 10 band photometric redshifts from the VIDEO and CFHTLS surveys. The results indicate that the radio population experiences mild positive evolution out to $z{\sim}1.2$ increasing their space density by a factor of $\sim$3, consistent with results of several previous studies. Beyond $z$=1.2 there is evidence of a slowing down of this evolution. Star-forming galaxies drive the more rapid evolution at low redshifts, $z{<}$1.2, while more slowly evolving AGN populations dominate at higher redshifts resulting in a decline in the evolution of the radio luminosity function at $z{>}$1.2. The evolution is best fit by pure luminosity evolution with star-forming galaxies evolving as $(1+z)^{2.47\pm0.12}$ and AGN as $(1+z)^{1.18\pm0.21}$.

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 compare the environmental and star formation properties of far-infrared detected and non--far-infrared detected galaxies out to $z \sim0.5$. Using optical spectroscopy and photometry from the Galaxy And Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS), with far-infrared observations from the {\em Herschel}-ATLAS Science Demonstration Phase (SDP), we apply the technique of Voronoi Tessellations to analyse the environmental densities of individual galaxies. Applying statistical analyses to colour, $r-$band magnitude and redshift-matched samples, we show there to be a significant difference at the 3.5$\sigma$ level between the normalized environmental densities of these two populations. This is such that infrared emission (a tracer of star formation activity) favours underdense regions compared to those inhabited by exclusively optically observed galaxies selected to be of the same $r-$band magnitude, colour and redshift. Thus more highly star-forming galaxies are found to reside in the most underdense environments, confirming previous studies that have proposed such a correlation. However, the degeneracy between redshift and far-infrared luminosity in our flux-density limited sample means that we are unable to make a stronger statement in this respect. We then apply our method to synthetic light cones generated from semi-analytic models, finding that over the whole redshift distribution the same correlations between star-formation rate and environmental density are found.

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

In order to investigate the growth of super-massive black holes (SMBHs), we construct the black hole mass function (BHMF) and Eddington ratio distribution function (ERDF) of X-ray-selected broad-line AGNs at z~1.4 in the Subaru XMM-Newton Deep Survey field. In this redshift range, a significant part of the accretion growth of SMBHs is thought to be taking place. Black hole masses of X-ray-selected broad-line AGNs are estimated using the width of the broad MgII line and the 3000A monochromatic luminosity. We supplement the MgII FWHM values with the Ha FWHM obtained from our NIR spectroscopic survey. Using the black hole masses of broad-line AGNs at redshifts between 1.18 and 1.68, the binned broad-line AGN BHMF and ERDF are calculated using the Vmax method. To properly account for selection effects that impact the binned estimates, we derive the corrected broad-line AGN BHMF and ERDF by applying the Maximum Likelihood method, assuming that the ERDF is constant regardless of the black hole mass. We do not correct for the non-negligible uncertainties in virial BH mass estimates. If we compare the corrected broad-line AGN BHMF with that in the local Universe, the corrected BHMF at z~1.4 has a higher number density above 10^8 Msolar but a lower number density below that mass range. The evolution may be indicative of a down-sizing trend of accretion activity among the SMBH population. The evolution of broad-line AGN ERDF from z=1.4 to 0 indicates that the fraction of broad-line AGNs with accretion rate close to the Eddington-limit is higher at higher redshifts.

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)

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

We investigate the [O II] emission line properties of 18,508 quasars at z<1.6 drawn from the Sloan Digital Sky Survey (SDSS) quasar sample. The quasar sample has been separated into 1,692 radio-loud and 16,816 radio-quiet quasars (RLQs and RQQs hereafter) matched in both redshift and i'-band absolute magnitude. We use the [O II]\lambda3726+3729 line as an indicator of star formation. Based on these measurements we find evidence that star-formation activity is higher in the RLQ population. The mean equivalent widths (EW) for [O II] are EW([O II])_RL=7.80\pm0.30 \AA, and EW([O II])_RQ=4.77\pm0.06 \AA, for the RLQ and RQQ samples respectively. The mean [O II] luminosities are \log[L([O II])_RL/W]=34.31\pm0.01 and \log[L([O II])_RQ/W]=34.192\pm0.004 for the samples of RLQs and RQQs respectively. Finally, to overcome possible biases in the EW measurements due to the continuum emission below the [O II] line being contaminated by young stars in the host galaxy, we use the ratio of the [O II] luminosity to rest-frame i'-band luminosity, in this case, we find for the RLQs \log[L([O II])_RL/L_opt]=-3.89\pm0.01 and \log[L([O II])_RQ/L_opt]=-4.011\pm0.004 for RQQs. However the results depend upon the optical luminosity of the quasar. RLQs and RQQs with the same high optical luminosity \log(L_opt/W)>38.6, tend to have the same level of [O II] emission. On the other hand, at lower optical luminosities \log(L_opt/W)<38.6, there is a clear [O II] emission excess for the RLQs. As an additional check of our results we use the [O III] emission line as a tracer of the bolometric accretion luminosity, instead of the i'-band absolute magnitude, and we obtain similar results. Radio jets appear to be the main reason for the [O II] emission excess in the case of RLQs. In contrast, we suggest AGN feedback ensures that the two populations acquire the same [O II] emission at higher optical luminosities.

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

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

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

We present high-spatial resolution imaging obtained with the Submillimeter Array (SMA) at 880um and the Keck Adaptive Optics (AO) system at Ks-band of a gravitationally lensed sub-millimeter galaxy (SMG) at z=4.243 discovered in the Herschel-Astrophysical Terahertz Large Area Survey. The SMA data (angular resolution ~0.6") resolve the dust emission into multiple lensed images, while the Keck AO Ks-band data (angular resolution ~0.1") resolve the lens into a pair of galaxies separated by 0.3". We present an optical spectrum of the foreground lens obtained with the Gemini-South telescope that provides a lens redshift of z_lens = 0.595 +/- 0.005. We develop and apply a new lens modeling technique in the visibility plane that shows that the SMG is magnified by a factor of mu = 4.1 +/- 0.2 and has an intrinsic infrared (IR) luminosity of L_IR = (2.1 +/- 0.2) x 10^13 Lsun. We measure a half-light radius of the background source of r_s = 4.4 +/- 0.5 kpc which implies an IR luminosity surface density of Sigma_IR = (3.4 +/- 0.9) x 10^11 Lsun kpc^-2, a value that is typical of z > 2 SMGs but significantly lower than IR luminous galaxies at z~0. The two lens galaxies are compact (r_lens ~ 0.9 kpc) early-types with Einstein radii of theta_E1 = 0.57 +/- 0.01 and theta_E2 = 0.40 +/- 0.01 that imply masses of M_lens1 = (7.4 +/- 0.5) x 10^10 Msun and M_lens2 = (3.7 +/- 0.3) x 10^10 Msun. The two lensing galaxies are likely about to undergo a dissipationless merger, and the mass and size of the resultant system should be similar to other early-type galaxies at z~0.6. This work highlights the importance of high spatial resolution imaging in developing models of strongly lensed galaxies discovered by Herschel.

We present IRAM Plateau de Bure Interferometer observations of the 12CO(3-2) emission from two far-infrared luminous QSOs at z ~ 2.5 selected from the Herschel-ATLAS survey. These far-infrared bright QSOs were selected to have supermassive black holes (SMBHs) with masses similar to those thought to reside in sub-millimetre galaxies (SMGs) at z ~ 2.5; making them ideal candidates as systems in transition from an ultraluminous infrared galaxy phase to a sub-mm faint, unobscured, QSO. We detect 12CO(3-2) emission from both QSOs and we compare their baryonic, dynamical and SMBH masses to those of SMGs at the same epoch. We find that these far-infrared bright QSOs have similar dynamical but lower gas masses than SMGs. In particular we find that far-infrared bright QSOs have ~50+-23% less warm/dense gas than SMGs, which combined with previous results showing the QSOs lack the extended, cool reservoir of gas seen in SMGs, suggests that they are at a different evolutionary stage. This is consistent with the hypothesis that far-infrared bright QSOs represent a short (~1Myr) but ubiquitous phase in the transformation of dust obscured, gas-rich, starburst-dominated SMGs into unobscured, gas-poor, QSOs.

We identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

In this paper we investigate the performance of the likelihood ratio method as a tool for identifying optical and infrared counterparts to proposed radio continuum surveys with SKA precursor and pathfinder telescopes. We present a comparison of the infrared counterparts identified by the likelihood ratio in the VISTA Deep Extragalactic Observations (VIDEO) survey to radio observations with 6, 10 and 15 arcsec resolution. We cross-match a deep radio catalogue consisting of radio sources with peak flux density $>$ 60 $\mu$Jy with deep near-infrared data limited to $K_{\mathrm{s}}\lesssim$ 22.6. Comparing the infrared counterparts from this procedure to those obtained when cross-matching a set of simulated lower resolution radio catalogues indicates that degrading the resolution from 6 arcsec to 10 and 15 arcsec decreases the completeness of the cross-matched catalogue by approximately 3 and 7 percent respectively. When matching against shallower infrared data, comparable to that achieved by the VISTA Hemisphere Survey, the fraction of radio sources with reliably identified counterparts drops from $\sim$89%, at $K_{\mathrm{s}}\lesssim$22.6, to 47% with $K_{\mathrm{s}}\lesssim$20.0. Decreasing the resolution at this shallower infrared limit does not result in any further decrease in the completeness produced by the likelihood ratio matching procedure. However, we note that radio continuum surveys with the MeerKAT and eventually the SKA, will require long baselines in order to ensure that the resulting maps are not limited by instrumental confusion noise.

PDS 144 is a pair of Herbig Ae stars that are separated by 5.35" on the sky. It has previously been shown to have an A2Ve Herbig Ae star viewed at 83\circ inclination as its northern member and an A5Ve Herbig Ae star as its southern member. Direct imagery revealed a disk occulting PDS 144 N - the first edge-on disk observed around a Herbig Ae star. The lack of an obvious disk in direct imagery suggested PDS 144 S might be viewed face-on or not physically associated with PDS 144 N. Multi-epoch Hubble Space Telescope imagery of PDS 144 with a 5 year baseline demonstrates PDS 144 N & S are comoving and have a common proper motion with TYC 6782-878-1. TYC 6782-878-1 has previously been identified as a member of Upper Sco sub-association A at d = 145 \pm 2 pc with an age of 5-10 Myr. Ground-based imagery reveals jets and a string of Herbig-Haro knots extending 13' (possibly further) which are aligned to within 7\circ \pm 6\circ on the sky. By combining proper motion data and the absence of a dark mid-plane with radial velocity data, we measure the inclination of PDS 144 S to be i = 73\circ \pm 7\circ. The radial velocity of the jets from PDS 144 N & S indicates they, and therefore their disks, are misaligned by 25\circ \pm 9\circ. This degree of misalignment is similar to that seen in T Tauri wide binaries.

We present near-infrared spectroscopic observations of star-forming galaxies at z~1.4 with FMOS on the Subaru Telescope. We observed K-band selected galaxies in the SXDS/UDS fields with K<23.9 mag, 1.2<z_ph<1.6, M*>10^{9.5} Msun, and expected F(Halpha)>10^{-16} erg s^{-1} cm^{-2}. 71 objects in the sample have significant detections of Halpha. For these objects, excluding possible AGNs identified from the BPT diagram, gas-phase metallicities are obtained from [NII]/Halpha line ratio. The sample is split into three stellar mass bins, and the spectra are stacked in each stellar mass bin. The mass-metallicity relation obtained at z~1.4 is located between those at z~0.8 and z~2.2. We constrain an intrinsic scatter to be ~0.1 dex or larger in the mass-metallicity relation at z~1.4; the scatter may be larger at higher redshifts. We found trends that the deviation from the mass-metallicity relation depends on the SFR and the half light radius: Galaxies with higher SFR and larger half light radii show lower metallicities at a given stellar mass. One possible scenario for the trends is the infall of pristine gas accreted from IGM or through merger events. Our data points show larger scatter than the fundamental metallicity relation (FMR) at z~0.1 and the average metallicities slightly deviate from the FMR. The compilation of the mass-metallicity relations at z~3 to z~0.1 shows that they evolve smoothly from z~3 to z~0 without changing the shape so much except for the massive part at z~0.

The European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) is a 4-m class survey telescope for wide-field near-infrared imaging. VISTA is currently running a suite of six public surveys, which will shortly deliver their first Europe wide public data releases to ESO. The VISTA Kilo-degree Infrared Galaxy Survey (VIKING) forms a natural intermediate between current wide shallow, and deeper more concentrated surveys, by targeting two patches totalling 1500 sq.deg in the northern and southern hemispheres with measured 5-sigma limiting depths of Z ~ 22.4, Y ~ 21.4, J ~ 20.9, H ~ 19.9 and Ks ~19.3 (Vega). This architecture forms an ideal working parameter space for the discovery of a significant sample of 6.5 <= z <= 7.5 quasars. In the first data release priority has been placed on small areas encompassing a number of fields well sampled at many wavelengths, thereby optimising science gains and synergy whilst ensuring a timely release of the first products. For rare object searches e.g. high-z quasars, this policy is not ideal since photometric selection strategies generally evolve considerably with the acquisition of data. Without a reasonably representative data set sampling many directions on the sky it is not clear how a rare object search can be conducted in a highly complete and efficient manner. In this paper, we alleviate this problem by supplementing initial data with a realistic model of the spatial, luminosity and colour distributions of sources known to heavily contaminate photometric quasar selection spaces, namely dwarf stars of spectral type M, L and T. We use this model along with a subset of available data to investigate contamination of quasar selection space by cool stars and galaxies and lay down a set of benchmark selection constraints that limit contamination to reasonable levels whilst maintaining high completeness...

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

We present the first direct and unbiased measurement of the evolution of the dust mass function of galaxies over the past 5 billion years of cosmic history using data from the Science Demonstration Phase of the Herschel-ATLAS. The sample consists of galaxies selected at 250{\mu}m which have reliable counterparts from SDSS at z < 0.5, and contains 1867 sources. Dust masses are calculated using both a single temperature grey-body model for the spectral energy distribution and also using a model with multiple temperature components. The dust temperature for either model shows no trend with redshift. Splitting the sample into bins of redshift reveals a strong evolution in the dust properties of the most massive galaxies. At z = 0.4 - 0.5, massive galaxies had dust masses about five times larger than in the local Universe. At the same time, the dust-to-stellar mass ratio was about 3-4 times larger, and the optical depth derived from fitting the UV-sub-mm data with an energy balance model was also higher. This increase in the dust content of massive galaxies at high redshift is difficult to explain using standard dust evolution models and requires a rapid gas consumption timescale together with either a more top-heavy IMF, efficient mantle growth, less dust destruction or combinations of all three. This evolution in dust mass is likely to be associated with a change in overall ISM mass, and points to an enhanced supply of fuel for star formation at earlier cosmic epochs.

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

This paper presents a study of the environments of SDSS Quasi-Stellar-Objects (QSOs) in the Spitzer Extragalactic Representative Volume Survey (SERVS). We concentrate on the high-redshift QSOs as these have not been studied in large numbers with data of this depth before. We use the IRAC 3.6-4.5{\mu}m colour of objects and ancillary r-band data to filter out as much foreground contamination as possible. This technique allows us to find a significant (> 4-{\sigma}) over-density of galaxies around QSOs in a redshift bin centred on z ~ 2.0 and a (> 2-{\sigma}) over-density of galaxies around QSOs in a redshift bin centred on z ~ 3.3. We compare our findings to the predictions of a semi-analytic galaxy formation model, based on the {\Lambda}CDM millennium simulation, and find for both redshift bins that the model predictions match well the source-density we have measured from the SERVS data.

We use the science demonstration field data of the Herschel-ATLAS to study how star formation, traced by the far-infrared Herschel data, is related to both the accretion luminosity and redshift of quasars selected from the Sloan Digital Sky Survey and the 2SLAQ survey. By developing a maximum likelihood estimator to investigate the presence of correlations between the far-infrared and optical luminosities we find evidence that the star-formation in quasar hosts is correlated with both redshift and quasar accretion luminosity. Assuming a relationship of the form L_IR \propto L_QSO^{\theta} (1 + z)^{\zeta}, we find {\theta} = 0.22 +/- 0.08 and {\zeta} = 1.6 +/- 0.4, although there is substantial additional uncertainty in {\zeta} of order +/- 1, due to uncertainties in the host galaxy dust temperature. We find evidence for a large intrinsic dispersion in the redshift dependence, but no evidence for intrinsic dispersion in the correlation between L_QSO and L_IR, suggesting that the latter may be due to a direct physical connection between star formation and black hole accretion. This is consistent with the idea that both the quasar activity and star formation are dependent on the same reservoir of cold gas, so that they are both affected by the influx of cold gas during mergers or heating of gas via feedback processes.

We present physical properties of two submillimeter selected gravitationally lensed sources, identified in the Herschel Astrophysical Terahertz Large Area Survey. These submillimeter galaxies (SMGs) have flux densities > 100 mJy at 500 um, but are not visible in existing optical imaging. We fit light profiles to each component of the lensing systems in Spitzer IRAC 3.6 and 4.5 um data and successfully disentangle the foreground lens from the background source in each case, providing important constraints on the spectral energy distributions (SEDs) of the background SMG at rest-frame optical-near-infrared wavelengths. The SED fits show that these two SMGs have high dust obscuration with Av ~4 to 5 and star formation rates of ~100 M_sun/yr. They have low gas fractions and low dynamical masses compared to 850 um selected galaxies.

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

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

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

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

Fibre Multi-Object Spectrograph (FMOS) is the first near-infrared instrument with a wide field of view capable of acquiring spectra simultaneously from up to 400 objects. It has been developed as a common-use instrument for the F/2 prime-focus of the Subaru Telescope. The field coverage of 30' diameter is achieved using a new 3-element corrector optimized in the near-infrared (0.9-1.8um) wavelength range. Due to limited space at the prime-focus, we have had to develop a novel fibre positioner called "Echidna" together with two OH-airglow suppressed spectrographs. FMOS consists of three subsystems: the prime focus unit for IR, the fibre positioning system/connector units, and the two spectrographs. After full systems integration, FMOS was installed on the telescope in late 2007. Many aspects of performance were checked through various test and engineering observations. In this paper, we present the optical and mechanical components of FMOS and show the results of our on-sky engineering observations to date.

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

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

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

We present a comparison between Gaussian processes (GPs) and artificial neural networks (ANNs) as methods for determining photometric redshifts for galaxies, given training set data. In particular, we compare their degradation in performance as the training set size is degraded in ways which might be caused by the observational limitations of spectroscopy. We find that performance with large, complete training sets is very similar, although the ANN achieves slightly smaller root mean square errors. If the size of the training set is reduced by random sampling, the RMS errors of both methods increase, but they do so to a lesser extent and in a much smoother manner for the case of GP regression. When training objects are removed at redshifts 1.3<z<1.7, to simulate the effects of the "redshift desert" of optical spectroscopy, the GP regression is successful at interpolating across the redshift gap, while the ANN suffers from strong bias for test objects in this redshift range. Overall, GP regression has attractive properties for photometric redshift estimation, particularly for deep, high-redshift surveys where it is difficult to obtain a large, complete training set. At present, unlike the ANN code, public GP regression codes do not take account of inhomogeneous measurement errors on the photometric data, and thus cannot estimate reliable uncertainties on the predicted redshifts. However, a better treatment of errors is in principle possible, and the promising results in this paper suggest that such improved GP algorithms should be pursued. (abridged)

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.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We present the results of our search for high-redshift Lyman-break galaxies over the GOODS-South field. We use HST-ACS data in B, V, i' & z', VLT-ISAAC J and Ks, Spitzer-IRAC 3.6, 4.5, 5.8 and 8.0 micron data in conjunction with the new HAWK-I Y-band science verification data to search for dropout galaxies in the redshift range 6<z<9. We survey ~119 arcmin^2 to Y_AB=25.7 (5sigma), of which 37.5 arcmin^2 reaches Y_AB=25.9. Candidate z' and Y drop-outs were selected on the basis of a colour cut of (Y-J)_{AB}>0.75 mag and (z'-Y)_{AB}>1.0 mag respectively. We find no robust Y-drops (z~9) brighter than J_{AB}<25.4. In our search for z'-band dropouts (z~6.5-7.5), we identify four possible candidates, two with z'-drop colours and clear Spitzer-IRAC detections and two less likely candidates. We also identify two previously-known Galactic T-dwarf stellar contaminants with these colours, and two likely transient objects seen in the Y-band data. The implications if all or none of our candidates are real on the Ultra-Violet galaxy luminosity functions at z>6.5 are explored. We find our number of z'-drop candidates to be insufficient based on the expected number of z' drops in a simple no-evolution scenario from the z=3 Lyman-break galaxy luminosity function but we are consistent with the observed luminosity function at z~6 (if all our candidates are real). However, if one or both of our best z'-drop candidates are not z>6.5 galaxies, this would demand evolution of the luminosity function at early epochs. We show that the future surveys to be conducted with the ESO VISTA telescope over the next five years will be able to measure the bulk of the luminosity function for both z' and Y drop-outs.

In this letter we report the discovery of a z=4.88 radio galaxy discovered with a new technique which does not rely on pre-selection of a sample based on radio properties such as steep-spectral index or small angular size. This radio galaxy was discovered in the Elais-N2 field and has a spectral index of alpha = 0.75, i.e. not ultra-steep spectrum. It also has a luminosity consistent with being drawn from the break of the radio luminosity function and can therefore be considered as a typical radio galaxy. Using the Spitzer-SWIRE data over this field we find that the host galaxy is consistent with being similarly massive to the lower redshift powerful radio galaxies (~1-3L*). We note however, that at z=4.88 the H-alpha line is redshifted into the IRAC 3.6micron filter and some of the flux in this band may be due to this rather than stellar continuum emission. The discovery of such a distant radio source from our initial spectroscopic observations demonstrate the promise of our survey for finding the most distant radio sources.

We detect the optical afterglow and host galaxy of GRB 070714B. Our observations of the afterglow show an initial plateau in the lightcurve for approximately the first 5 to 25 minutes, then steepening to a powerlaw decay with index alpha= 0.86 +/- 0.10 for the period between 1 to 24 hours post burst. This is consistent with the X-ray light-curve which shows an initial plateau followed by a similar subsequent decay. At late time, we detect a host galaxy at the location of the optical transient. Gemini Nod & Shuffle spectroscopic observations of the host show a single emission line at 7167 angstroms which, based on a grizJHK photometric redshift, we conclude is the 3727 angstrom [O II] line. We therefore find a redshift of z=0.923. This redshift, as well as a subsequent probable spectroscopic redshift determination of GRB 070429B at z=0.904 by two other groups, significantly exceeds the previous highest spectroscopically confirmed short burst redshift of z=0.546 for GRB 051221. This dramatically moves back the time at which we know short bursts were being formed, and suggests that the present evidence for an old progenitor population may be observationally biased.

We discuss the optical properties, X-ray detections, and Active Galactic Nucleus (AGN) populations of four clusters at z ~ 1 in the Subaru-XMM Deep Field (SXDF). The velocity distribution and plausible extended X-ray detections are examined, as well as the number of X-ray point sources and radio sources associated with the clusters. We find that the two clusters that appear virialised and have an extended X-ray detection contain few, if any, AGN, whereas the two pre-virialised clusters have a large AGN population. This constitutes evidence that the AGN fraction in clusters is linked to the clusters' evolutionary stage. The number of X-ray AGN in the pre-virialised clusters is consistent with an overdensity of factor ~ 200; the radio AGN appear to be clustered with a factor of three to six higher. The median K-band luminosities of L_K = 1.7 +/- 0.7 L* for the X-ray sources and L_K = 2.3 +/- 0.1 L* for the radio sources support the theory that these AGN are triggered by galaxy interaction and merging events in sub-groups with low internal velocity distributions, which make up the cluster environment in a pre-virialisation evolutionary stage.

We study the 37 brightest radio sources in the Subaru/XMM-Newton Deep Field (SXDF). We have spectroscopic redshifts for 24 of 37 objects and photometric redshifts for the remainder, yielding a median redshift z_med for the whole sample of z_med ~= 1.1 and a median radio luminosity close to the `FRI/FRII' luminosity divide. Using mid-IR (Spitzer MIPS 24 um) data we expect to trace nuclear accretion activity, even if it is obscured at optical wavelengths, unless the obscuring column is extreme. Our results suggest that above the FRI/FRII radio luminosity break most of the radio sources are associated with objects that have excess mid-IR emission, only some of which are broad-line objects, although there is one clear low-accretion-rate object with an FRI radio structure. For extended steep-spectrum radio sources, the fraction of objects with mid-IR excess drops dramatically below the FRI/FRII luminosity break, although there exists at least one high-accretion-rate `radio-quiet' QSO. We have therefore shown that the strong link between radio luminosity (or radio structure) and accretion properties, well known at z ~ 0.1, persists to z ~ 1. Investigation of mid-IR and blue excesses shows that they are correlated as predicted by a model in which, when significant accretion exists, a torus of dust absorbs ~30% of the light, and the dust above and below the torus scatters >~1% of the light.

We have obtained deep optical spectroscopic data of the highest-redshift cluster candidate (z ~ 1.4, CVB13) selected by Van Breukelen et al. (2006) in a photometric optical/infrared catalogue of the Subaru XMM-Newton Deep Field. The data, which comprise 104 targeted galaxies, were taken with the DEep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck 2 telescope and yielded 31 secure redshifts in the range 1.25 < z < 1.54 within a 7' x 4' field centred on CVB13. Instead of one massive cluster at z = 1.4, we find evidence for three projected structures at z = 1.40, z = 1.45, and z = 1.48. The most statistically robust of these structures, at z = 1.454, has six spectroscopically confirmed galaxies. Its total mass is estimated at > 10^14 M_sun and it may therefore be termed a cluster. There is an X-ray source at the cluster position which is marginally spatially resolved but whose X-ray spectrum is too hard to be thermal cluster emission. Its origin could be the summed X-ray emission from active galaxies in, and projected onto, the cluster. Serendipitously we have discovered a cluster at z = 1.28 with a mass of > 10^14 M_sun at the same position on the sky, comprising six spectroscopically confirmed cluster galaxies and at least one additional radio source. The selection of CVB13 for the cluster catalogue was evidently aided by the superposition of other, presumably lower-mass, structures, whereas the single cluster at z = 1.28 contained too few galaxies to be isolated by the same algorithm. Given the complicated nature of such structures, caution must be employed when measuring the mass function of putative high-redshift clusters with photometric techniques alone.

Using mid-infrared and radio selection criteria, we pre-select a sample of candidate high-redshift type-2 quasars in the Subaru XMM-Newton Deep Field (SXDF). To filter out starburst contaminants, we use a bayesian method to fit the spectral energy distributions (SEDs) between 24-microns and B-band, obtain photometric redshifts, and identify the best candidates for high-z type-2 quasars. This leaves us with 12 z_phot >= 1.7 type-2 quasar candidates in an area ~0.8 deg^2, of which only two have secure X-ray detections. The two detected sources have estimated column densities N_H~2 & 3x10^27 m^-2, i.e. heavily obscured but Compton-thin quasars. Given the large bolometric luminosities and redshifts of the undetected objects, the lack of X-ray detections suggests extreme absorbing columns N_H >= 10^28 m^-2 are typical. We have found evidence for a population of ``Compton-thick'' high-redshift type-2 quasars, at least comparable to, and probably larger than the type-1 quasar population, although spectroscopic confirmation of their AGN nature is important.

We present a new cluster-finding algorithm based on a combination of the Voronoi Tessellation and Friends-Of-Friends methods. The algorithm utilises probability distribution functions derived from a photometric redshift analysis and is tested on simulated cluster-catalogues. We use a 9 band photometric catalogue over 0.5 square degrees in the Subaru XMM-Newton Deep Field. The photometry is comprised of 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. The cluster catalogue contains 13 clusters at redshifts 0.61 <= z <= 1.39 with luminosities 10L* < L_tot < 50 L*.

We present a new cluster-finding algorithm based on a combination of the Voronoi Tessellation and Friends-Of-Friends methods. The algorithm utilises probability distribution functions derived from a photometric redshift analysis. We test our algorithm on a set of simulated cluster-catalogues and have published elsewhere its employment on 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. This pilot study has detected clusters over 0.5 square degrees in the Subaru XMM-Newton Deep Field. The resulting cluster catalogue contains 13 clusters at redshifts 0.61 <= z <= 1.39 with luminosities 10 L* <~ L_tot <~ 50 L*.

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