Measurements of The Hubble-Lemaitre constant from early- and local-universe observations show a significant discrepancy. In an attempt to understand the origin of this mismatch, independent techniques to measure $H_0$ are required. One such technique, strong lensing time delays, is set to become a leading contender amongst the myriad methods due to forthcoming large strong lens samples. It is therefore critical to understand the systematic effects inherent in this method. In this paper, we quantify the influence of additional structures along the line-of-sight by adopting realistic lightcones derived from the \textit{CosmoDC2} semi-analytical extra-galactic catalogue. Using multiple lens plane ray-tracing to create a set of simulated strong lensing systems, we have investigated the impact of line-of-sight structures on time-delay measurements and in turn, on the inferred value of $H_0$. We have also tested the reliability of existing procedures for correcting for line-of-sight effects. We find that if the integrated contribution of the of line-of-sight structures is close to a uniform mass sheet, the bias in $H_0$ can be adequately corrected by including a constant external convergence $\kappa_{\rm ext}$ in the lens model. However, for realistic line-of-sight structures comprising many galaxies at different redshifts, this simple correction over-estimates the bias by a factor of approximately three. We therefore conclude that lens modelling must incorporate multiple lens planes to account for line-of-sight structures for accurate and precise inference of $H_0$.

In this paper we develop a new unsupervised machine learning technique comprised of a feature extractor, a convolutional autoencoder (CAE), and a clustering algorithm consisting of a Bayesian Gaussian mixture model (BGM). We apply this technique to visual band space-based simulated imaging data from the Euclid Space Telescope using data from the Strong Gravitational Lenses Finding Challenge. Our technique promisingly captures a variety of lensing features such as Einstein rings with different radii, distorted arc structures, etc, without using predefined labels. After the clustering process, we obtain several classification clusters separated by different visual features which are seen in the images. Our method successfully picks up $\sim$63\ percent of lensing images from all lenses in the training set. With the assumed probability proposed in this study, this technique reaches an accuracy of $77.25\pm 0.48$\% in binary classification using the training set. Additionally, our unsupervised clustering process can be used as the preliminary classification for future surveys of lenses to efficiently select targets and to speed up the labelling process. As the starting point of the astronomical application using this technique, we not only explore the application to gravitationally lensed systems, but also discuss the limitations and potential future uses of this technique.

Using the IRAM NOrthern Extended Millimeter Array (NOEMA), we conducted a program to measure redshifts for 13 bright galaxies detected in the Herschel Astrophysical Large Area Survey (H-ATLAS) with $S_{500{\mu}\rm m}\ge$80 mJy. We report reliable spectroscopic redshifts for 12 individual sources, which are derived from scans of the 3 and 2 mm bands, covering up to 31 GHz in each band, and are based on the detection of at least two emission lines. The spectroscopic redshifts are in the range $2.08<z<4.05$ with a median value of $z=2.9\pm$0.6. The sources are unresolved or barely resolved on scales of 10 kpc. In one field, two galaxies with different redshifts were detected. In two cases the sources are found to be binary galaxies with projected distances of ~140 kpc. The linewidths of the sources are large, with a mean value for the full width at half maximum of 700$\pm$300 km/s and a median of 800 km/s. We analyse the nature of the sources with currently available ancillary data to determine if they are lensed or hyper-luminous $L_{\rm FIR} > 10^{13}\,L_\odot$ galaxies. We also present a reanalysis of the spectral energy distributions including the continuum flux densities measured at 3 and 2 mm to derive the overall properties of the sources. Future prospects based on these efficient measurements of redshifts of high-z galaxies using NOEMA are outlined, including a comprehensive survey of all the brightest Herschel galaxies.

We present a study of the dust, stars and atomic gas (HI) in an HI-selected sample of local galaxies (z<0.035) in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) fields. This HI-selected sample reveals a population of very high gas fraction (>80 per cent), low stellar mass sources that appear to be in the earliest stages of their evolution. We compare this sample with dust and stellar mass selected samples to study the dust and gas scaling relations over a wide range of gas fraction (proxy for evolutionary state of a galaxy). The most robust scaling relations for gas and dust are those linked to NUV-r (SSFR) and gas fraction, these do not depend on sample selection or environment. At the highest gas fractions, our additional sample shows the dust content is well below expectations from extrapolating scaling relations for more evolved sources, and dust is not a good tracer of the gas content. The specific dust mass for local galaxies peaks at a gas fraction of ~75 per cent. The atomic gas depletion time is also longer for high gas fraction galaxies, opposite to the trend found for molecular gas depletion timescale. We link this trend to the changing efficiency of conversion of HI to H2 as galaxies increase in stellar mass surface density as they evolve. Finally, we show that galaxies start out barely obscured and increase in obscuration as they evolve, yet there is no clear and simple link between obscuration and global galaxy properties.

The largest Herschel extragalactic surveys, H-ATLAS and HerMES, have selected a sample of "ultrared" dusty, star-forming galaxies (DSFGs) with rising SPIRE flux densities ($S_{500} > S_{350} > S_{250}$; so-called "500 $\mu$m-risers") as an efficient way for identifying DSFGs at higher redshift ($z > 4$). In this paper, we present a large Spitzer follow-up program of 300 Herschel ultrared DSFGs. We have obtained high-resolution ALMA, NOEMA, and SMA data for 63 of them, which allow us to securely identify the Spitzer/IRAC counterparts and classify them as gravitationally lensed or unlensed. Within the 63 ultrared sources with high-resolution data, $\sim$65% appear to be unlensed, and $\sim$27% are resolved into multiple components. We focus on analyzing the unlensed sample by directly performing multi-wavelength spectral energy distribution (SED) modeling to derive their physical properties and compare with the more numerous $z \sim 2$ DSFG population. The ultrared sample has a median redshift of 3.3, stellar mass of 3.7 $\times$ 10$^{11}$ $M_{\odot}$, star formation rate (SFR) of 730 $M_{\odot}$yr$^{-1}$, total dust luminosity of 9.0 $\times$ 10$^{12}$ $L_{\odot}$, dust mass of 2.8 $\times$ 10$^9$ $M_{\odot}$, and V-band extinction of 4.0, which are all higher than those of the ALESS DSFGs. Based on the space density, SFR density, and stellar mass density estimates, we conclude that our ultrared sample cannot account for the majority of the star-forming progenitors of the massive, quiescent galaxies found in infrared surveys. Our sample contains the rarer, intrinsically most dusty, luminous and massive galaxies in the early universe that will help us understand the physical drivers of extreme star formation.

We explore the effectiveness of deep learning convolutional neural networks (CNNs) for estimating strong gravitational lens mass model parameters. We have investigated a number of practicalities faced when modelling real image data, such as how network performance depends on the inclusion of lens galaxy light, the addition of colour information and varying signal-to-noise. Our CNN was trained and tested with strong galaxy-galaxy lens images simulated to match the imaging characteristics of the Large Synoptic Survey Telescope (LSST) and Euclid. For images including lens galaxy light, the CNN can recover the lens model parameters with an acceptable accuracy, although a 34 per cent average improvement in accuracy is obtained when lens light is removed. However, the inclusion of colour information can largely compensate for the drop in accuracy resulting from the presence of lens light. While our findings show similar accuracies for single epoch Euclid VIS and LSST r-band datasets, we find a 24 per cent increase in accuracy by adding g- and i-band images to the LSST r-band without lens light and a 20 per cent increase with lens light. The best network performance is obtained when it is trained and tested on images where lens light exactly follows the mass, but when orientation and ellipticity of the light is allowed to differ from those of the mass, the network performs most consistently when trained with a moderate amount of scatter in the difference between the mass and light profiles.

Follow-up observations of (sub-)mm-selected gravitationally-lensed systems have allowed a more detailed study of the dust-enshrouded phase of star-formation up to very early cosmic times. Here, the case of the gravitationally lensed merger in HATLAS J142935.3-002836 (also known as H1429-0028; z_lens=0.218, z_bkg=1.027) is revisited following recent developments in the literature and new APEX observations targeting two carbon monoxide (CO) rotational transitions J_up=3 and 6. We show that the line-profiles comprise three distinct velocity components, where the fainter high-velocity one is less magnified and more compact. The modelling of the observed spectral line energy distribution of CO J_up=2 to 6 and [CI]3P_1-3P_0 assumes a large velocity gradient scenario, where the analysis is based on four statistical approaches. Since the detected gas and dust emission comes exclusively from only one of the two merging components (the one oriented North-South, NS), we are only able to determine upper-limits for the companion. The molecular gas in the NS component in H1429-0028 is found to have a temperature of ~70K, a volume density of log(n[/cm3])~3.7, to be expanding at ~10km/s/pc, and amounts to M_H2=4(-2,+3)*1e9 Msun. The CO to H2 conversion factor is estimated to be alpha_CO=0.4(-0.2,+0.3) Msun/(K.km/s.pc2). The NS galaxy is expected to have a factor of >10x more gas than its companion (M_H2<3e8 Msun). Nevertheless, the total amount of molecular gas in the system comprises only up to 15 per cent (1sigma upper-limit) of the total (dynamical) mass.

Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$\mu$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust continuum and molecular gas emission are resolved into a nearly complete ~1"5 diameter Einstein ring plus a weaker image in the center, which is caused by a special dual deflector lensing configuration. The observed line profiles of the CO, H2O and H2O+ lines are strikingly similar. In the source plane, we reconstruct the dust continuum images and the spectral cubes of the line emission at sub-kpc scales. The reconstructed dust emission in the source plane is dominated by a compact disk with an effective radius of 0.7kpc plus an overlapping extended disk with a radius twice as large. While the average magnification for the dust continuum is $\mu$~10-11, the magnification of the line emission varies 5 to 22 across different velocity components. The emission lines have similar spatial and kinematic distributions. The molecular gas and dust content reveal that G09v1.97 is a gas-rich major merger in its pre-coalescence phase. Both of the merging companions are intrinsically ULIRGs with LIR reaching $\gtrsim 4\times10^{12}L_\odot$, and the total LIR of G09v1.97 is $1.4\times10^{13}L_\odot$. The approaching southern galaxy shows no obvious kinematic structure with a semi-major half-light radius a_s=0.4kpc, while the receding galaxy resembles an a_s=1.2kpc rotating disk. The two galaxies are separated by a projected distance of 1.3kpc, bridged by weak line emission that is co-spatially located with the cold-dust-emission peak, suggesting a large amount of cold ISM in the interacting region. (abridged)

We present the final results from our survey of luminous $z \sim $ 5.5 quasars. This is the first systematic quasar survey focusing on quasars at $z \sim$ 5.5, during the post-reionization epoch. It has been challenging to select quasars at $5.3 < z < 5.7$ using conventional color selections, due to their similar optical colors to those of late-type stars, especially M dwarfs. We developed a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR, and mid-IR photometry, using data from the Sloan Digital Sky Survey (SDSS), PanSTARR1 (PS1), the UKIRT Infrared Deep Sky Surveys - Large Area Survey, the UKIRT Hemisphere Survey, the VISTA Hemisphere Survey, and the Wide Field Infrared Survey Explorer ($WISE$), covering $\sim$ 11000 deg$^2$ of high galactic latitude sky. In this paper, we present the discovery of 15 new quasars at $z\sim 5.5$. Together with results from Yang et al. (2017), our survey provides a complete, flux-limited sample of 31 quasars at $5.3 \le z \le 5.7$. We measure the quasar spatial density at $z \sim 5.5$ and $M_{1450} < -26.2$. Our result is consistent with the rapid decline of quasar spatial density from $z = 5$ to 6, with $k=-0.66 \pm 0.05$ ($\rho(z) \propto 10^{kz}$). In addition, we present a new survey using optical colors only from the full PS1 area for luminous quasars at $z = 5.0 - 5.5$, beyond the SDSS footprint, and report the preliminary results from this survey, including 51 new quasars discovered at $4.61\le z \le5.71$.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

We present Atacama Large Millimeter/submillimeter Array observations of a radio-loud and millimeter-bright galaxy at z=2.6. Gravitational lensing by a foreground galaxy at z~0.2 provides access to physical scales of approximately 360 pc, and we resolve a 2.5 kpc-radius ring of star-forming molecular gas, traced by atomic carbon CI(1-0) and carbon monoxide CO(4-3). We also detect emission from the cyanide radical, CN(4-3). With a velocity width of 680 km/s, this traces dense molecular gas travelling at velocities nearly a factor of two larger than the rotation speed of the molecular ring. While this could indicate the presence of a dynamical and photochemical interaction between the active galactic nucleus and molecular interstellar medium on scales of a few 100 pc, on-going feedback is unlikely to have a significant impact on the assembly of stellar mass in the molecular ring, given the ~10s Myr depletion timescale due to star formation.

The galaxies found in optical surveys fall in two distinct regions of a diagram of optical colour versus absolute magnitude: the red sequence and the blue cloud with the green valley in between. We show that the galaxies found in a submillimetre survey have almost the opposite distribution in this diagram, forming a `green mountain'. We show that these distinctive distributions follow naturally from a single, continuous, curved Galaxy Sequence in a diagram of specific star-formation rate versus stellar mass without there being the need for a separate star-forming galaxy Main Sequence and region of passive galaxies. The cause of the red sequence and the blue cloud is the geometric mapping between stellar mass/specific star-formation rate and absolute magnitude/colour, which distorts a continuous Galaxy Sequence in the diagram of intrinsic properties into a bimodal distribution in the diagram of observed properties. The cause of the green mountain is Malmquist bias in the submillimetre waveband, with submillimetre surveys tending to select galaxies on the curve of the Galaxy Sequence, which have the highest ratios of submillimetre-to-optical luminosity. This effect, working in reverse, causes galaxies on the curve of the Galaxy Sequence to be underrepresented in optical samples, deepening the green valley. The green valley is therefore not evidence (1) for there being two distinct populations of galaxies, (2) for galaxies in this region evolving more quickly than galaxies in the blue cloud and the red sequence, (c) for rapid quenching processes in the galaxy population.

Using high-resolution (sub-kiloparsec scale) submillimeter data obtained by ALMA, we analyze the star formation rate (SFR), gas content and kinematics in SDP 81, a gravitationally-lensed star-forming galaxy at redshift 3. We estimate the SFR surface density ($\Sigma_{\mathrm{SFR}}$) in the brightest clump of this galaxy to be $357^{+135}_{-85}\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$, over an area of $0.07\pm0.02\,\mathrm{kpc}^2$. Using the intensity-weighted velocity of CO$\,$(5-4), we measure the turbulent velocity dispersion in the plane-of-the-sky and find $\sigma_{\mathrm{v,turb}} = 37\pm5\,\mathrm{km\,s}^{-1}$ for the star-forming clump, in good agreement with previous estimates along the line of sight. Our measurements of gas surface density, freefall time and turbulent Mach number reveal that the role of turbulence is vital to explaining the observed SFR in this clump. While the Kennicutt Schmidt (KS) relation predicts a SFR surface density of $\Sigma_{\mathrm{SFR,KS}} = 52\pm17\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$, the single-freefall model by Krumholz, Dekel and McKee (KDM) predicts $\Sigma_{\mathrm{SFR,KDM}} = 106\pm37\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$. In contrast, the multi-freefall (turbulence) model by Salim, Federrath and Kewley (SFK) gives $\Sigma_{\mathrm{SFR,SFK}} = 491^{+139}_{-194}\,\mathrm{M_{\odot}\,yr^{-1}\,kpc^{-2}}$. Although the SFK relation overestimates the SFR in this clump (possibly due to the ignorance of magnetic field), it provides the best prediction among the available models. Finally, we compare the star formation and gas properties of this high-redshift galaxy to local star-forming regions and find that the SFK relation provides the best estimates of SFR in both local and high-redshift galaxies.

This work presents AutoLens, the first entirely automated modeling suite for the analysis of galaxy-scale strong gravitational lenses. AutoLens simultaneously models the lens galaxy's light and mass whilst reconstructing the extended source galaxy on an adaptive pixel-grid. The method's approach to source-plane discretization is amorphous, adapting its clustering and regularization to the intrinsic properties of the lensed source. The lens's light is fitted using a superposition of Sersic functions, allowing AutoLens to cleanly deblend its light from the source. Single component mass models representing the lens's total mass density profile are demonstrated, which in conjunction with light modeling can detect central images using a centrally cored profile. Decomposed mass modeling is also shown, which can fully decouple a lens's light and dark matter and determine whether the two component are geometrically aligned. The complexity of the light and mass models are automatically chosen via Bayesian model comparison. These steps form AutoLens's automated analysis pipeline, such that all results in this work are generated without any user-intervention. This is rigorously tested on a large suite of simulated images, assessing its performance on a broad range of lens profiles, source morphologies and lensing geometries. The method's performance is excellent, with accurate light, mass and source profiles inferred for data sets representative of both existing Hubble imaging and future Euclid wide-field observations.

The {\it Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 660 deg$^2$ with the PACS and SPIRE cameras in five photometric bands: 100, 160, 250, 350 and 500\mic. This is the second of three papers describing the data release for the large fields at the south and north Galactic poles (NGP and SGP). In this paper we describe the catalogues of far-infrared and submillimetre sources for the NGP and SGP, which cover 177 deg$^2$ and 303 deg$^2$, respectively. The catalogues contain 153,367 sources for the NGP field and 193,527 sources for the SGP field detected at more than 4$\sigma$ significance in any of the 250, 350 or 500\mic\ bands. The source detection is based on the 250\mic\ map, and we present photometry in all five bands for each source, including aperture photometry for sources known to be extended. The rms positional accuracy for the faintest sources is about 2.4 arc seconds in both right ascension and declination. We present a statistical analysis of the catalogues and discuss the practical issues -- completeness, reliability, flux boosting, accuracy of positions, accuracy of flux measurements -- necessary to use the catalogues for astronomical projects.

We have modelled high resolution ALMA imaging of six strong gravitationally lensed galaxies detected by the Herschel Space Observatory. Our modelling recovers mass properties of the lensing galaxies and, by determining magnification factors, intrinsic properties of the lensed sub-millimetre sources. We find that the lensed galaxies all have high ratios of star formation rate to dust mass, consistent with or higher than the mean ratio for high redshift sub-millimetre galaxies and low redshift ultra-luminous infra-red galaxies. Source reconstruction reveals that most galaxies exhibit disturbed morphologies. Both the cleaned image plane data and the directly observed interferometric visibilities have been modelled, enabling comparison of both approaches. In the majority of cases, the recovered lens models are consistent between methods, all six having mass density profiles that are close to isothermal. However, one system with poor signal to noise shows mildly significant differences.

With the advent of wide-area submillimeter surveys, a large number of high-redshift gravitationally lensed dusty star-forming galaxies (DSFGs) has been revealed. Due to the simplicity of the selection criteria for candidate lensed sources in such surveys, identified as those with $S_{500\mu m} > 100$ mJy, uncertainties associated with the modelling of the selection function are expunged. The combination of these attributes makes submillimeter surveys ideal for the study of strong lens statistics. We carried out a pilot study of the lensing statistics of submillimetre-selected sources by making observations with the Atacama Large Millimetre Array (ALMA) of a sample of strongly-lensed sources selected from surveys carried out with the Herschel Space Observatory. We attempted to reproduce the distribution of image separations for the lensed sources using a halo mass function taken from a numerical simulation which contains both dark matter and baryons. We used three different density distributions, one based on analytical fits to the halos formed in the EAGLE simulation and two density distributions (Singular Isothermal Sphere (SIS) and SISSA) that have been used before in lensing studies. We found that we could reproduce the observed distribution with all three density distributions, as long as we imposed an upper mass transition of $\sim$$10^{13} M_{\odot}$ for the SIS and SISSA models, above which we assumed that the density distribution could be represented by an NFW profile. We show that we would need a sample of $\sim$500 lensed sources to distinguish between the density distributions, which is practical given the predicted number of lensed sources in the Herschel surveys.

We perform lens modelling and source reconstruction of Submillimeter Array (SMA) data for a sample of 12 strongly lensed galaxies selected at 500$\mu$m in the Herschel Astrophysical Terahertz Large Area Survey H-ATLAS. A previous analysis of the same dataset used a single S\`ersic profile to model the light distribution of each background galaxy. Here we model the source brightness distribution with an adaptive pixel scale scheme, extended to work in the Fourier visibility space of interferometry. We also present new SMA observations for seven other candidate lensed galaxies from the H-ATLAS sample. Our derived lens model parameters are in general consistent with previous findings. However, our estimated magnification factors, ranging from 3 to 10, are lower. The discrepancies are observed in particular where the reconstructed source hints at the presence of multiple knots of emission. We define an effective radius of the reconstructed sources based on the area in the source plane where emission is detected above 5$\sigma$. We also fit the reconstructed source surface brightness with an elliptical Gaussian model. We derive a median value $r_{eff}\,\sim 1.77\,$kpc and a median Gaussian full width at half maximum $\sim1.47\,$kpc. After correction for magnification, our sources have intrinsic star formation rates SFR$\,\sim900-3500\,M_{\odot}yr^{-1}$, resulting in a median star formation rate surface density $\Sigma_{SFR}\sim132\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ (or $\sim 218\,M_{\odot}$ yr$^{-1}$ kpc$^{-2}$ for the Gaussian fit). This is consistent with what observed for other star forming galaxies at similar redshifts, and is significantly below the Eddington limit for a radiation pressure regulated starburst.

By determining the nature of all the Planck compact sources within 808.4 deg^2 of large Herschel surveys, we have identified 27 candidate proto-clusters of dusty star forming galaxies (DSFGs) that are at least 3{\sigma} overdense in either 250, 350 or 500 $\mu$mm sources. We find roughly half of all the Planck compact sources are resolved by Herschel into multiple discrete objects, with the other half remaining unresolved by Herschel. We find a significant difference between versions of the Planck catalogues, with earlier releases hosting a larger fraction of candidate proto-clusters and Galactic Cirrus than later releases, which we ascribe to a difference in the filters used in the creation of the three catalogues. We find a surface density of DSFG candidate proto-clusters of (3.3 $\pm$ 0.7) x 10^(-2) sources deg^(-2), in good agreement with previous similar studies. We find that a Planck colour selection of S_{857}/S_{545} < 2 works well to select candidate proto-clusters, but can miss proto-clusters at z < 2. The Herschel colours of individual candidate proto-cluster members indicate our candidate proto-clusters all likely all lie at z > 1. Our candidate proto-clusters are a factor of 5 times brighter at 353 GHz than expected from simulations, even in the most conservative estimates. Further observations are needed to confirm whether these candidate proto-clusters are physical clusters, multiple proto-clusters along the line of sight, or chance alignments of unassociated sources.

This paper forms part of the second major public data release of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). In this work, we describe the identification of optical and near-infrared counterparts to the submillimetre detected sources in the $177$ deg$^2$ North Galactic Plane (NGP) field. We used the likelihood ratio method to identify counterparts in the Sloan Digital Sky Survey and in the UKIRT Imaging Deep Sky Survey within a search radius of $10$ arcsec of the H-ATLAS sources with a $4\sigma$ detection at $250 \, \mu$m. We obtained reliable ($R \ge 0.8 $) optical counterparts with $r< 22.4$ for 42429 H-ATLAS sources ($37.8$ per cent), with an estimated completeness of $71.7$ per cent and a false identification rate of $4.7$ per cent. We also identified counterparts in the near-infrared using deeper $K$-band data which covers a smaller $\sim25$ deg$^2$. We found reliable near-infrared counterparts to $61.8$ per cent of the $250$-$\mu$m-selected sources within that area. We assessed the performance of the likelihood ratio method to identify optical and near-infrared counterparts taking into account the depth and area of both input catalogues. Using catalogues with the same surface density of objects in the overlapping $\sim25$ deg$^2$ area, we obtained that the reliable fraction in the near-infrared ($54.8$ per cent) is significantly higher than in the optical ($36.4$ per cent). Finally, using deep radio data which covers a small region of the NGP field, we found that $80 - 90$ per cent of our reliable identifications are correct.

We present the largest submillimeter images that have been made of the extragalactic sky. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 660 deg$^2$ with the PACS and SPIRE cameras in five photometric bands: 100, 160, 250, 350, and 500{\mu}m. In this paper we present the images from our two largest fields which account for ~75% of the survey. The first field is 180.1 deg$^2$ in size centered on the North Galactic Pole (NGP) and the second field is 317.6 deg$^2$ in size centered on the South Galactic Pole. The NGP field serendipitously contains the Coma cluster. Over most (~80%) of the images, the pixel noise, including both instrumental noise and confusion noise, is approximately 3.6, and 3.5 mJy/pix at 100 and 160{\mu}m, and 11.0, 11.1 and 12.3 mJy/beam at 250, 350 and 500{\mu}m, respectively, but reaches lower values in some parts of the images. If a matched filter is applied to optimize point-source detection, our total 1{\sigma} map sensitivity is 5.7, 6.0, and 7.3 mJy at 250, 350, and 500{\mu}m, respectively. We describe the results of an investigation of the noise properties of the images. We make the most precise estimate of confusion in SPIRE maps to date finding values of 3.12+/-0.07, 4.13+/-0.02 and 4.45+/-0.04 mJy/beam at 250, 350, and 500{\mu}m in our un-convolved maps. For PACS we find an estimate of the confusion noise in our fast-parallel observations of 4.23 and 4.62 mJy/beam at 100 and 160{\mu}m. Finally, we give recipes for using these images to carry out photometry, both for unresolved and extended sources.

In this work we measure and study the cross-correlation signal between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range $0.2<z<0.8$, and a background sample of H-ATLAS galaxies with photometric redshifts $\gtrsim1.2$. It constitutes a substantial improvement over the cross-correlation measurements made by Gonzalez-Nuevo et al. (2014) with updated catalogues and wider area (with $S/N\gtrsim 5$ below 10' and reaching $S/N\sim 20$ below 30"). The better statistics allow us to split the sample in different redshift bins and to perform a tomographic analysis (with $S/N\gtrsim 3$ below 10 arcmin and reaching $S/N\sim 15$ below 30"). Moreover, we implement a halo model to extract astrophysical information about the background galaxies and the deflectors that are producing the lensing link between the foreground (lenses) and background (sources) samples. In the case of the sources, we find typical mass values in agreement with previous studies: a minimum halo mass to host a central galaxy, $M_{min}\sim 10^{12.26} M_\odot$, and a pivot halo mass to have at least one sub-halo satellite, $M_1\sim 10^{12.84} M_\odot$. However, the lenses are massive galaxies or even galaxy groups/clusters, with minimum mass of $M_{min}^{lens}\sim 10^{13.06} M_\odot$. Above a mass of $M_1^{lens}\sim 10^{14.57} M_\odot$ they contain at least one additional satellite galaxy which contributes to the lensing effect. The tomographic analysis shows that, while $M_1^{lens}$ is almost redshift independent, there is a clear evolution of increase $M_{min}^{lens}$ with redshift in agreement with theoretical estimations. Finally, the halo modeling allows us to identify a strong lensing contribution to the cross-correlation for angular scales below 30". This interpretation is supported by the results of basic but effective simulations.

The Herschel Space Observatory has revealed a very different galaxyscape from that shown by optical surveys which presents a challenge for galaxy-evolution models. The Herschel surveys reveal (1) that there was rapid galaxy evolution in the very recent past and (2) that galaxies lie on a a single Galaxy Sequence (GS) rather than a star-forming `main sequence' and a separate region of `passive' or `red-and-dead' galaxies. The form of the GS is now clearer because far-infrared surveys such as the Herschel ATLAS pick up a population of optically-red star-forming galaxies that would have been classified as passive using most optical criteria. The space-density of this population is at least as high as the traditional star-forming population. By stacking spectra of H-ATLAS galaxies over the redshift range 0.001 < z < 0.4, we show that the galaxies responsible for the rapid low-redshift evolution have high stellar masses, high star-formation rates but, even several billion years in the past, old stellar populations - they are thus likely to be relatively recent ancestors of early-type galaxies in the Universe today. The form of the GS is inconsistent with rapid quenching models and neither the analytic bathtub model nor the hydrodynamical EAGLE simulation can reproduce the rapid cosmic evolution. We propose a new gentler model of galaxy evolution that can explain the new Herschel results and other key properties of the galaxy population.

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

We present the Herschel Bright Sources (HerBS) sample, a sample of bright, high-redshift Herschel sources detected in the 616.4 square degree H-ATLAS survey. The HerBS sample contains 209 galaxies, selected with a 500 {\mu}m flux density greater than 80 mJy and an estimated redshift greater than 2. The sample consists of a combination of HyLIRGs and lensed ULIRGs during the epoch of peak cosmic star formation. In this paper, we present SCUBA-2 observations at 850 ${\mu}$m of 189 galaxies of the HerBS sample, 152 of these sources were detected. We fit a spectral template to the Herschel-SPIRE and 850 ${\mu}$m SCUBA-2 flux densities of 22 sources with spectroscopically determined redshifts, using a two-component modified blackbody spectrum as a template. We find a cold- and hot-dust temperature of 21.29 K and 45.80 K, a cold-to-hot dust mass ratio of 26.62 and a $\beta$ of 1.83. The poor quality of the fit suggests that the sample of galaxies is too diverse to be explained by our simple model. Comparison of our sample to a galaxy evolution model indicates that the fraction of lenses is high. Out of the 152 SCUBA-2 detected galaxies, the model predicts 128.4 $\pm$ 2.1 of those galaxies to be lensed (84.5%). The SPIRE 500 ${\mu}$m flux suggests that out of all 209 HerBS sources, we expect 158.1 $\pm$ 1.7 lensed sources, giving a total lensing fraction of 76 per cent.

We investigate the mean star formation rates (SFRs) in the host galaxies of ~3000 optically selected QSOs from the SDSS survey within the Herschel-ATLAS fields, and a radio-luminous sub-sample, covering the redshift range of z = 0.2-2.5. Using WISE & Herschel photometry (12 - 500{\mu}m) we construct composite SEDs in bins of redshift and AGN luminosity. We perform SED fitting to measure the mean infrared luminosity due to star formation, removing the contamination from AGN emission. We find that the mean SFRs show a weak positive trend with increasing AGN luminosity. However, we demonstrate that the observed trend could be due to an increase in black hole (BH) mass (and a consequent increase of inferred stellar mass) with increasing AGN luminosity. We compare to a sample of X-ray selected AGN and find that the two populations have consistent mean SFRs when matched in AGN luminosity and redshift. On the basis of the available virial BH masses, and the evolving BH mass to stellar mass relationship, we find that the mean SFRs of our QSO sample are consistent with those of main sequence star-forming galaxies. Similarly, the radio-luminous QSOs have mean SFRs that are consistent with both the overall QSO sample and with star-forming galaxies on the main sequence. In conclusion, on average QSOs reside on the main sequence of star-forming galaxies, and the observed positive trend between the mean SFRs and AGN luminosity can be attributed to BH mass and redshift dependencies.

Combining large-area optical quasar surveys with the new far-infrared Herschel-ATLAS Data Release 1, we search for an observational signature associated with the minority of quasars possessing bright far-infrared (FIR) luminosities. We find that FIR-bright quasars show broad CIV emission line blueshifts in excess of that expected from the optical luminosity alone, indicating particularly powerful nuclear outflows. The quasars show no signs of having redder optical colours than the general ensemble of optically-selected quasars, ruling out differences in line-of-sight dust within the host galaxies. We postulate that these objects may be caught in a special evolutionary phase, with unobscured, high black hole accretion rates and correspondingly strong nuclear outflows. The high FIR emission found in these objects is then either a result of star formation related to the outflow, or is due to dust within the host galaxy illuminated by the quasar. We are thus directly witnessing coincident small-scale nuclear processes and galaxy-wide activity, commonly invoked in galaxy simulations which rely on feedback from quasars to influence galaxy evolution.

We present an extragalactic survey using observations from the Atacama Large Millimeter/submillimeter Array (ALMA) to characterise galaxy populations up to $z=0.35$: the Valpara\'iso ALMA Line Emission Survey (VALES). We use ALMA Band-3 CO(1--0) observations to study the molecular gas content in a sample of 67 dusty normal star-forming galaxies selected from the $Herschel$ Astrophysical Terahertz Large Area Survey ($H$-ATLAS). We have spectrally detected 49 galaxies at $>5\sigma$ significance and 12 others are seen at low significance in stacked spectra. CO luminosities are in the range of $(0.03-1.31)\times10^{10}$ K km s$^{-1}$ pc$^2$, equivalent to $\log({\rm M_{gas}/M_{\odot}}) =8.9-10.9$ assuming an $\alpha_{\rm CO}$=4.6(K km s$^{-1}$ pc$^{2}$)$^{-1}$, which perfectly complements the parameter space previously explored with local and high-z normal galaxies. We compute the optical to CO size ratio for 21 galaxies resolved by ALMA at $\sim 3$."$5$ resolution (6.5 kpc), finding that the molecular gas is on average $\sim$ 0.6 times more compact than the stellar component. We obtain a global Schmidt-Kennicutt relation, given by $\log [\Sigma_{\rm SFR}/({\rm M_{\odot} yr^{-1}kpc^{-2}})]=(1.26 \pm 0.02) \times \log [\Sigma_{\rm M_{H2}}/({\rm M_{\odot}\,pc^{-2}})]-(3.6 \pm 0.2)$. We find a significant fraction of galaxies lying at `intermediate efficiencies' between a long-standing mode of star-formation activity and a starburst, specially at $\rm L_{IR}=10^{11-12} L_{\odot}$. Combining our observations with data taken from the literature, we propose that star formation efficiencies can be parameterised by $\log [{\rm SFR/M_{H2}}]=0.19 \times {\rm (\log {L_{IR}}-11.45)}-8.26-0.41 \times \arctan[-4.84 (\log {\rm L_{IR}}-11.45) ]$. Within the redshift range we explore ($z<0.35$), we identify a rapid increase of the gas content as a function of redshift.

We present the first discoveries from a survey of $z\gtrsim6$ quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg$^2$ of sky down to $z_{\rm AB}\sim23.0$, and UKIDSS and UHS, which will map the northern sky at $0<DEC<+60^{\circ}$, reaching $J_{\rm VEGA}\sim19.6$ (5-$\sigma$). The combination of these datasets allows us to discover quasars at redshift $z\gtrsim7$ and to conduct a complete census of the faint quasar population at $z\gtrsim6$. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint $z\gtrsim6$ quasars and one new $z=6.63$ quasar in our pilot spectroscopic observations. The two new $z\sim6$ quasars are at $z=6.07$ and $z=6.17$ with absolute magnitudes at rest-frame wavelength 1450 \AA\ being $M_{1450}=-25.83$ and $M_{1450}=-25.76$, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at $z\gtrsim6$. The new $z=6.63$ quasar has an absolute magnitude of $M_{1450}=-25.95$. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS datasets to find $z\gtrsim7$ quasars. Extrapolating from previous QLF measurements, we predict that these combined datasets will yield $\sim200$ $z\sim6$ quasars to $z_{\rm AB} < 21.5$, $\sim1{,}000$ $z\sim6$ quasars to $z_{\rm AB}<23$, and $\sim 30$ quasars at $z>6.5$ to $J_{\rm VEGA}<19.5$.

We present initial results from the first systematic survey of luminous $z\sim 5.5$ quasars. Quasars at $z \sim$ 5.5, the post-reionization epoch, are crucial tools to explore the evolution of intergalactic medium, quasar evolution and the early super-massive black hole growth. However, it has been very challenging to select quasars at redshifts 5.3 $\le z \le$ 5.7 using conventional color selections, due to their similar optical colors to late-type stars, especially M dwarfs, resulting in a glaring redshift gap in quasar redshift distributions. We develop a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR and mid-IR photometric data from Sloan Digital Sky Survey (SDSS), UKIRT InfraRed Deep Sky Surveys - Large Area Survey (ULAS), VISTA Hemisphere Survey (VHS) and Wide field Infrared Survey Explorer (WISE). From our pilot observations in SDSS-ULAS/VHS area, we have discovered 15 new quasars at 5.3 $\le z \le$ 5.7 and 6 new lower redshift quasars, with SDSS z band magnitude brighter than 20.5. Including other two $z \sim$ 5.5 quasars already published in our previous work, we now construct an uniform quasar sample at 5.3 $\le z \le$ 5.7 with 17 quasars in a $\sim$ 4800 square degree survey area. For further application in a larger survey area, we apply our selection pipeline to do a test selection by using the new wide field J band photometric data from a preliminary version of the UKIRT Hemisphere Survey (UHS). We successfully discover the first UHS selected $z \sim$ 5.5 quasar.

We present the calibration between the dust continuum luminosity and interstellar gas content obtained from the Valpara\'{i}so ALMA Line Emission Survey (VALES) sample of 67 main-sequence star-forming galaxies at 0.02<$z$<0.35. We use CO(1-0) observations from the Atacama Large Millimetre/submillimetre Array (ALMA) to trace the molecular gas mass, $M_{\mathrm{H}_{2}}$, and estimate the rest-frame monochromatic luminosity at 850 $\mu$m, $L_{\nu_{850}}$, by extrapolating the dust continuum from MAGPHYS modelling of the far-ultraviolet to submillimetre spectral energy distribution sampled by the Galaxy And Mass Assembly (GAMA) survey. Adopting $\alpha_{\rm CO}$ = 6.5 (K km s$^{-1}$ pc$^{2}$)$^{-1}$, the average ratio of $L_{\nu_{850}}/M_{\mathrm{H}_{2}}$ = (6.4$\pm$1.4)$\times10^{19}$ erg s$^{-1}$ Hz$^{-1}$ $\mathrm{M}_{\odot}^{-1}$, in excellent agreement with literature values. We obtain a linear fit of $\log_{10}$ ($M_{\mathrm{H}_{2}}/\mathrm{M}_{\odot}$) = (0.92$\pm$0.02) $\log_{10}$ ($L_{\nu_{850}}$/erg s$^{-1}$ Hz$^{-1}$)-(17.31$\pm$0.59). We provide relations between $L_{\nu_{850}}$, $M_{\mathrm{H}_{2}}$ and $M_{\mathrm{ISM}}$ when combining the VALES and literature samples, and adopting a Galactic $\alpha_{\rm CO}$ value.

Submillimeter galaxies (SMGs) at $z\gtrsim1$ are luminous in the far-infrared and have star-formation rates, SFR, of hundreds to thousands of solar masses per year. However, it is unclear whether they are true analogs of local ULIRGs or whether the mode of their star formation is more similar to that in local disk galaxies. We target these questions by using Herschel-PACS to examine the conditions in the interstellar medium (ISM) in far-infrared luminous SMGs at z~1-4. We present 70-160 micron photometry and spectroscopy of the [OIV]26 micron, [FeII]26 micron, [SIII]33 micron, [SiII]34 micron, [OIII]52 micron, [NIII]57 micron, and [OI]63 micron fine-structure lines and the S(0) and S(1) hydrogen rotational lines in 13 lensed SMGs identified by their brightness in early Herschel data. Most of the 13 targets are not individually spectroscopically detected and we instead focus on stacking these spectra with observations of an additional 32 SMGs from the \herschel\ archive -- representing a complete compilation of PACS spectroscopy of SMGs. We detect [OI]63 micron, [SiII]34 micron, and [NIII]57 micron at >3sigma in the stacked spectra, determining that the average strengths of these lines relative to the far-IR continuum are $(0.36\pm0.12)\times10^{-3}$, $(0.84\pm0.17)\times10^{-3}$, and $(0.27\pm0.10)\times10^{-3}$, respectively. Using the [OIII]52/[NIII]57 emission line ratio we show that SMGs have average gas-phase metallicities $\gtrsim Z_{\rm sun}$. By using PDR modelling and combining the new spectral measurements with integrated far-infrared fluxes and existing [CII]158 micron data we show that SMGs have average gas densities, n, of $\sim10^{1-3}{\rm cm^{-3}}$ and FUV field strengths, $G_0\sim10^{2.2-4.5}$ (in Habing units: $1.6\times10^{-3}{\rm erg~cm^{-2}~s^{-1}}$), consistent with both local ULIRGs and lower luminosity star-forming galaxies.

We use new Band-3 CO(1-0) observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to study the physical conditions in the interstellar gas of a sample of 27 dusty main-sequence star-forming galaxies at 0.03<$z$<0.2 present in the Valpara\'iso ALMA Line Emission Survey (VALES). The sample is drawn from far-IR bright galaxies over $\sim$160 deg$^{2}$ in the Herschel Astrophysical Terahertz Large Area Survey (HATLAS), which is covered by Herschel [CII] 158 $\mu$m spectroscopy and far-infrared (FIR) photometry. The [CII] and CO lines are both detected at >5$\sigma$ in 26 sources. We find an average [CII] to CO(1-0) luminosity ratio of 3500$\pm$1200 for our sample that is consistent with previous studies. Using the [CII], CO and FIR measurements as diagnostics of the physical conditions of the interstellar medium, we compare these observations to the predictions of a photodissociation region (PDR) model to determine the gas density, surface temperature, pressure, and the strength of the incident far-ultraviolet (FUV) radiation field, $G_{0}$, normalised to the Habing Field. The majority of our sample exhibit hydrogen densities of 4 < $\log n/\mathrm{cm}^{3}$ < 5.5 and experience an incident FUV radiation field with strengths of 2 < $\log G_0$ < 3 when adopting standard adjustments. A comparison to galaxy samples at different redshifts indicates that the average strength of the FUV radiation field appears constant up to redshift $z\sim$6.4, yet the neutral gas density increases with redshift by a factor of $\sim$100, that persists regardless of various adjustments to our observable quantities. This evolution could provide an explanation for the observed evolution of the star formation rate density with cosmic time, yet could arise from a combination of observational biases when using different suites of emission lines as diagnostic tracers of PDR gas.

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

We present a source-plane reconstruction of a ${\it Herschel}$ and ${\it Planck}$-detected gravitationally-lensed dusty star-forming galaxy (DSFG) at $z=1.68$ using {\it Hubble}, Sub-millimeter Array (SMA), and Keck observations. The background sub-millimeter galaxy (SMG) is strongly lensed by a foreground galaxy cluster at $z=0.997$ and appears as an arc of length $\sim 15^{\prime \prime}$ in the optical images. The continuum dust emission, as seen by SMA, is limited to a single knot within this arc. We present a lens model with source plane reconstructions at several wavelengths to show the difference in magnification between the stars and dust, and highlight the importance of a multi-wavelength lens models for studies involving lensed DSFGs. We estimate the physical properties of the galaxy by fitting the flux densities to model SEDs leading to a magnification-corrected star formation rate of $390 \pm 60$ M$_{\odot}$ yr$^{-1}$ and a stellar mass of $1.1 \pm 0.4\times 10^{11}$ M$_{\odot}$. These values are consistent with high-redshift massive galaxies that have formed most of their stars already. The estimated gas-to-baryon fraction, molecular gas surface density, and SFR surface density have values of $0.43 \pm 0.13$, $350 \pm 200$ M$_{\odot}$ pc$^{-2}$, and $\sim 12 \pm 7~$M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, respectively. The ratio of star formation rate surface density to molecular gas surface density puts this among the most star-forming systems, similar to other measured SMGs and local ULIRGS.

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $\mu L_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/\Delta V \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $\mu$m continuum opacity $\tau_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.

We compare common star-formation rate (SFR) indicators in the local Universe in the GAMA equatorial fields (around 160 sq. deg.), using ultraviolet (UV) photometry from GALEX, far-infrared (FIR) and sub-millimetre (sub-mm) photometry from H-ATLAS, and Halpha spectroscopy from the GAMA survey. With a high-quality sample of 745 galaxies (median redshift 0.08), we consider three SFR tracers: UV luminosity corrected for dust attenuation using the UV spectral slope beta (SFRUV,corr), Halpha line luminosity corrected for dust using the Balmer decrement (BD) (SFRHalpha,corr), and the combination of UV and IR emission (SFRUV+IR). We demonstrate that SFRUV,corr can be reconciled with the other two tracers after applying attenuation corrections by calibrating IRX (i.e. the IR to UV luminosity ratio) and attenuation in the Halpha (derived from BD) against beta. However, beta on its own is very unlikely to be a reliable attenuation indicator. We find that attenuation correction factors depend on parameters such as stellar mass, z and dust temperature (Tdust), but not on Halpha equivalent width (EW) or Sersic index. Due to the large scatter in the IRX vs beta correlation, when compared to SFRUV+IR, the beta-corrected SFRUV,corr exhibits systematic deviations as a function of IRX, BD and Tdust.

We present the first major data release of the largest single key-project in area carried out in open time with the Herschel Space Observatory. The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) is a survey of 600 deg^2 in five photometric bands - 100, 160, 250, 350 and 500 um - with the PACS and SPIRE cameras. In this paper and a companion paper (Bourne et al. 2016) we present the survey of three fields on the celestial equator, covering a total area of 161.6 deg^2 and previously observed in the Galaxy and Mass Assembly (GAMA) spectroscopic survey. This paper describes the Herschel images and catalogues of the sources detected on the SPIRE 250 um images. The 1-sigma noise for source detection, including both confusion and instrumental noise, is 7.4, 9.4 and 10.2 mJy at 250, 350 and 500 um. Our catalogue includes 120230 sources in total, with 113995, 46209 and 11011 sources detected at >4-sigma at 250, 350 and 500 um. The catalogue contains detections at >3-sigma at 100 and 160 um for 4650 and 5685 sources, and the typical noise at these wavelengths is 44 and 49 mJy. We include estimates of the completeness of the survey and of the effects of flux bias and also describe a novel method for determining the true source counts. The H-ATLAS source counts are very similar to the source counts from the deeper HerMES survey at 250 and 350 um, with a small difference at 500 um. Appendix A provides a quick start in using the released datasets, including instructions and cautions on how to use them.

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

Aims. We aim to study the 250 micron luminosity function (LF) down to much fainter luminosities than achieved by previous efforts. Methods. We developed a modified stacking method to reconstruct the 250 micron LF using optically selected galaxies from the SDSS survey and Herschel maps of the GAMA equatorial fields and Stripe 82. Our stacking method not only recovers the mean 250 micron luminosities of galaxies that are too faint to be individually detected, but also their underlying distribution functions. Results. We find very good agreement with previous measurements in the overlapping luminosity range. More importantly, we are able to derive the LF down to much fainter luminosities (around 25 times fainter) than achieved by previous studies. We find strong positive luminosity evolution \propto (1 + z)^4.89\pm1.07 and moderate negative density evolution \propto (1 + z)^-1.02\pm0.54 over the redshift range z=[0.02, 0.5].

We describe a far infrared survey of the Coma cluster and the galaxy filament it resides within. Our survey covers an area of $\sim$150 deg$^2$ observed by $Herschel$ H-ATLAS in five bands at 100, 160, 250, 350 and 500$\mu$m. The SDSS spectroscopic survey ($m_{r} \le 17.8)$ is used to define an area (within the Virial radius) and redshift selected ($4268 < v < 9700$ km s$^{-1}$) sample of 744 Coma cluster galaxies - the Coma Cluster Catalogue (CCC). For comparison we also define a sample of 951 galaxies in the connecting filament - the Coma Filament Catalogue (CFC). The optical positions and parameters are used to define appropriate apertures to measure each galaxy's far infrared emission. We have detected 99 of 744 (13\%) and 422 of 951 (44\%) of the cluster and filament galaxies in the SPIRE 250$\mu$m band. We consider the relative detection rates of galaxies of different morphological types finding that it is only the S0/Sa population that shows clear differences between the cluster and filament. We find no differences between the dust masses and temperatures of cluster and filament galaxies with the exception of early type galaxy dust temperatures, which are significantly hotter in the cluster than in the filament (X-ray heating?). From a chemical evolution model we find no evidence for different evolutionary processes (gas loss or infall) between galaxies in the cluster and filament.

Using results from the Herschel Astrophysical Terrahertz Large-Area Survey and the Galaxy and Mass Assembly project, we show that, for galaxy masses above approximately 1.0e8 solar masses, 51% of the stellar mass-density in the local Universe is in early-type galaxies (ETGs: Sersic n > 2.5) while 89% of the rate of production of stellar mass-density is occurring in late-type galaxies (LTGs: Sersic n < 2.5). From this zero-redshift benchmark, we have used a calorimetric technique to quantify the importance of the morphological transformation of galaxies over the history of the Universe. The extragalactic background radiation contains all the energy generated by nuclear fusion in stars since the Big Bang. By resolving this background radiation into individual galaxies using the deepest far-infrared survey with the Herschel Space Observatory and a deep near-infrared/optical survey with the Hubble Space Telescope (HST), and using measurements of the Sersic index of these galaxies derived from the HST images, we estimate that approximately 83% of the stellar mass-density formed over the history of the Universe occurred in LTGs. The difference between this and the fraction of the stellar mass-density that is in LTGs today implies there must have been a major transformation of LTGs into ETGs after the formation of most of the stars.

We present a new pixelized method for the inversion of gravitationally lensed extended source images which we term adaptive semi-linear inversion (SLI). At the heart of the method is an h-means clustering algorithm which is used to derive a source plane pixelization that adapts to the lens model magnification. The distinguishing feature of adaptive SLI is that every pixelization is derived from a random initialization, ensuring that data discretization is performed in a completely different and unique way for every lens model parameter set. We compare standard SLI on a fixed source pixel grid with the new method and demonstrate the shortcomings of the former when modeling singular power law ellipsoid (SPLE) lens profiles. In particular, we demonstrate the superior reliability and efficiency of adaptive SLI which, by design, fixes the number of degrees of freedom (NDOF) of the optimization and thereby removes biases present with other methods that allow the NDOF to vary. In addition, we highlight the importance of data discretization in pixel-based inversion methods, showing that adaptive SLI averages over significant systematics that are present when a fixed source pixel grid is used. In the case of the SPLE lens profile, we show how the method successfully samples its highly degenerate posterior probability distribution function with a single non-linear search. The robustness of adaptive SLI provides a firm foundation for the development of a strong lens modeling pipeline, which will become necessary in the short-term future to cope with the increasing rate of discovery of new strong lens systems.

We have modelled Atacama Large Millimeter/sub-millimeter Array (ALMA) long baseline imaging of the strong gravitational lens system H-ATLAS J090311.6+003906 (SDP.81). We have reconstructed the distribution of band 6 and 7 continuum emission in the z=3.042 source and we have determined its kinematic properties by reconstructing CO(5-4) and CO(8-7) line emission in bands 4 and 6. The continuum imaging reveals a highly non-uniform distribution of dust with clumps on scales of 200 pc. In contrast, the CO line emission shows a relatively smooth, disk-like velocity field which is well fit by a rotating disk model with an inclination angle of 40+/-5 degrees and an asymptotic rotation velocity of 320 km/s. The inferred dynamical mass within 1.5kpc is (3.5+/-0.5)x10^{10} M_sol which is comparable to the total molecular gas masses of (2.7+/-0.5)x10^{10} M_sol and (3.5+/-0.6)x10^{10} M_sol from the dust continuum emission and CO emission respectively. Our new reconstruction of the lensed HST near-infrared emission shows two objects which appear to be interacting, with the rotating disk of gas and dust revealed by ALMA distinctly offset from the near-infrared emission. The clumpy nature of the dust and a low value of the Toomre parameter of Q=0.3 suggest that the disk is in a state of collapse. We estimate a star formation rate in the disk of 470+/-80 M_sol/yr with an efficiency 65 times greater than typical low-redshift galaxies. Our findings add to the growing body of evidence that the most infra-red luminous, dust obscured galaxies in the high redshift Universe represent a population of merger induced starbursts.

We present the properties of the first 250 $\mu$m blind sample of nearby galaxies (15 < D < 46 Mpc) containing 42 objects from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Herschel's sensitivity probes the faint end of the dust luminosity function for the first time, spanning a range of stellar mass (7.4 < log$_{10}$ M$_{\star}$ < 11.3 M$_{\odot}$), star formation activity (-11.8 < log$_{10}$ SSFR < -8.9 yr$^{-1}$), gas fraction (3-96 per cent), and colour (0.6 < FUV-Ks < 7.0 mag). The median cold dust temperature is 14.6 K, colder than in the Herschel Reference Survey (18.5 K) and Planck Early Release Compact Source Catalogue 17.7 K. The mean dust-to-stellar mass ratio in our sample is higher than these surveys by factors of 3.7 and 1.8, with a dust mass volume density of (3.7 $\pm$ 0.7) x 10$^{5}$ M$_{\odot}$ Mpc$^{-3}$. Counter-intuitively, we find that the more dust rich a galaxy, the lower its UV attenuation. Over half of our dust-selected sample are very blue in FUV-Ks colour, with irregular and/or highly flocculent morphology, these galaxies account for only 6 per cent of the sample's stellar mass but contain over 35 per cent of the dust mass. They are the most actively star forming galaxies in the sample, with the highest gas fractions and lowest UV attenuation. They also appear to be in an early stage of converting their gas into stars, providing valuable insights into the chemical evolution of young galaxies.

The Herschel Space Observatory has had a tremendous impact on the study of extragalactic dust. Specifically, early-type galaxies (ETG) have been the focus of several studies. In this paper we combine results from two Herschel studies - a Virgo cluster study HeViCS and a broader, low-redshift H-ATLAS/GAMA study - and contrast the dust and associated properties for similar mass galaxies. This comparison is motivated by differences in results exhibited between multiple Herschel studies of early-type galaxies. A comparison between consistent modified blackbody derived dust mass is carried out, revealing strong differences between the two samples in both dust mass and dust-to-stellar mass ratio. In particular, the HeViCS sample lacks massive ETG with as high a specific dust content as found in H-ATLAS. This is most likely connected with the difference in environment for the two samples. We calculate nearest neighbour environment densities in a consistent way, showing that H-ATLAS ETG occupy sparser regions of the local Universe, whereas HeViCS ETG occupy dense regions. This is also true for ETG that are not Herschel-detected but are in the Virgo and GAMA parent samples. Spectral energy distributions are fit to the panchromatic data. From these we find that in H-ATLAS the specific star formation rate anticorrelates with stellar mass and reaches values as high as in our Galaxy. On the other hand HeViCS ETG appear to have little star formation. Based on the trends found here, H-ATLAS ETG are thought to have more extended star formation histories and a younger stellar population than HeViCS ETG.

We exploit long-baseline ALMA sub-mm observations of the lensed star-forming galaxy SDP 81 at z=3.042 to investigate the properties of inter-stellar medium on scales of 50-100pc. The kinematics of the CO gas within this system are well described by a rotationally-supported disk with an inclination-corrected rotation speed, v=320+/-20km/s and a dynamical mass of M=(3.5+/-1.0)x10^10Mo within a radius of 1.5 kpc. The disk is gas rich and unstable, with a Toomre parameter, Q=0.30+/-0.10 and so should collapse in to star-forming regions with Jeans length L_J~130pc. We identify five star-forming regions within the ISM on these scales and show that their scaling relations between luminosity, line-widths and sizes are significantly offset from those typical of molecular clouds in local Galaxies (Larson's relations). These offsets are likely to be caused by the high external hydrostatic pressure for the interstellar medium (ISM), P/kB=(40+/-20)x10^7K/cm3, which is ~10,000x higher than the typical ISM pressure in the Milky Way. The physical conditions of the star-forming ISM and giant molecular clouds appears to be similar to the those found in the densest environments in the local Universe, such as those in the Galactic center.

Interstellar dust in galaxies can be traced either through its extinction effects on the star light, or through its thermal emission at infrared wavelengths. Recent radiative transfer studies of several nearby edge-on galaxies have found an apparent inconsistency in the dust energy balance: the radiative transfer models that successfully explain the optical extinction underestimate the observed fluxes by an average factor of three. We investigate the dust energy balance for IC4225 and NGC5166, two edge-on spiral galaxies observed by the Herschel Space Observatory in the frame of the H-ATLAS survey. We start from models which were constrained from optical data and extend them to construct the entire spectral energy distribution of our galaxies. These predicted values are subsequently compared to the observed far-infrared fluxes. We find that including a young stellar population in the modelling is necessary as it plays a non-negligible part in the heating of the dust grains. While the modelling approach for both galaxies is nearly identical, we find two very different results. As is often seen in other edge-on spiral galaxies, the far-infrared emission of our radiative transfer model of IC4225 underestimates the observed fluxes by a factor of about three. For NGC5166 on the other hand, we find that both the predicted spectral energy distribution as well as the simulated images match the observations particularly well. We explore possible reasons for this difference and conclude that it is unlikely that one single mechanism is the cause of the dust energy balance problem in spiral galaxies. We discuss the different approaches that can be considered in order to get a conclusive answer on the origin this discrepancy.

We present the first measurement of the correlation between the map of the cosmic microwave background (CMB) lensing potential derived from the \emph{Planck} nominal mission data and $z\gtrsim 1.5$ galaxies detected by the \emph{Herschel}-ATLAS (H-ATLAS) survey covering about $600\,\hbox{deg}^2$, i.e. about 1.4\% of the sky. We reject the hypothesis that there is no correlation between CMB lensing and galaxy detection at a $20\,\sigma$ significance, checking the result by performing a number of null tests. The significance of the detection of the theoretically expected cross-correlation signal is found to be $10\,\sigma$. The galaxy bias parameter, $b$, derived from a joint analysis of the cross-power spectrum and of the auto-power spectrum of the galaxy density contrast is found to be $b=2.80^{+0.12}_{-0.11}$, consistent with earlier estimates for H-ATLAS galaxies at similar redshifts. On the other hand, the amplitude of the cross-correlation is found to be a factor $1.62 \pm 0.16$ higher than expected from the standard model and also found by cross-correlation analyses with other tracers of the large-scale structure. The enhancement due to lensing magnification can account for only a fraction of the excess cross-correlation signal. We suggest that part of it may be due to an incomplete removal of the contamination of the CIB, that includes the H-ATLAS sources we are cross-correlating with. In any case, the highly significant detection reported here using a catalog covering only 1.4\% of the sky demonstrates the potential of CMB lensing correlations with submillimeter surveys.

We present Hubble Space Telescope (HST) WFC3 imaging and grism spectroscopy observations of the {\it Herschel}-selected gravitationally-lensed starburst galaxy HATLASJ1429-0028. The lensing system consists of an edge-on foreground disk galaxy at $z=0.218$ with a nearly complete Einstein ring of the infrared luminous galaxy at $z=1.027$. The WFC3 spectroscopy with G102 and G141 grisms, covering the wavelength range of 0.8 to 1.7 $\mu$m, resulted in detections of H$\alpha$+[NII], H$\beta$, [SII], and [OIII] for the background galaxy from which we measure line fluxes and ratios. The Balmer line ratio H$\alpha$/H$\beta$ of $7.5 \pm 4.4$, when corrected for [NII], results in an extinction for the starburst galaxy of $E(B-V)=0.8 \pm 0.5$. The $H\alpha$ based star-formation rate, when corrected for extinction, is $60 \pm 50$ M$_{\odot}$ yr$^{-1}$, lower than the instantaneous star-formation rate of 390 $\pm$ 90 M$_{\odot}$ yr$^{-1}$ from the total IR luminosity. We also compare the nebular line ratios of HATLASJ1429-0028 with other star-forming and sub-mm bright galaxies. The nebular line ratios are consistent with an intrinsic ultra-luminous infrared galaxy with no evidence for excitation by an active galactic nucleus (AGN). We estimate the metallicity, $12 + log(O/H)$, of HATLASJ1429-0028 to be 8.49 $\pm$ 0.16. Such a low value is below the average relations for stellar mass vs. metallicity of galaxies at $z \sim 1$ for a galaxy with stellar mass of $\sim 2 \times 10^{11}$ M$_{\odot}$. The combination of high stellar mass, lack of AGN indicators, low metallicity, and the high star-formation rate of HATLASJ1429-0028 suggest that this galaxy is currently undergoing a rapid formation.

We explore the behaviour of [CII]-157.74um forbidden fine-structure line observed in a sample of 28 galaxies selected from ~50deg^2 of the H-ATLAS survey. The sample is restricted to galaxies with flux densities higher than S_160um>150mJy and optical spectra from the GAMA survey at 0.02<z<0.2. Far-IR spectra centred on this redshifted line were taken with the PACS instrument on-board the Herschel Space Observatory. The galaxies span 10<log(L_IR/Lo)<12 (where L_IR=L_IR[8-1000um]) and 7.3<log(L_[CII]/Lo)<9.3, covering a variety of optical galaxy morphologies. The sample exhibits the so-called [CII] deficit at high IR luminosities, i.e. L_[CII]/L_IR (hereafter [CII]/IR) decreases at high L_IR. We find significant differences between those galaxies presenting [CII]/IR>2.5x10^-3 with respect to those showing lower ratios. In particular, those with high ratios tend to have: (1) L_IR<10^11Lo; (2) cold dust temperatures, T_d<30K; (3) disk-like morphologies in r-band images; (4) a WISE colour 0.5<S_12um/S_22um<1.0; (5) low surface brightness Sigma_IR~10^8-9 Lo kpc^-2, (6) and specific star-formation rates of sSFR~0.05-3 Gyr^-1. We suggest that the strength of the far-UV radiation fields (<G_O>) is main parameter responsible for controlling the [CII]/IR ratio. It is possible that relatively high <G_O> creates a positively charged dust grain distribution, impeding an efficient photo-electric extraction of electrons from these grains to then collisionally excite carbon atoms. Within the brighter IR population, 11<log(L_IR/Lo)<12, the low [CII]/IR ratio is unlikely to be modified by [CII] self absorption or controlled by the presence of a moderately luminous AGN (identified via the BPT diagram).

Gamma-ray bursts (GRBs) are the most energetic phenomena in the Universe; believed to result from the collapse and subsequent explosion of massive stars. Even though it has profound consequences for our understanding of their nature and selection biases, little is known about the dust properties of the galaxies hosting GRBs. We present analysis of the far-infrared properties of an unbiased sample of 20 \textit{BeppoSAX} and \textit{Swift} GRB host galaxies (at an average redshift of $z\,=\,3.1$) located in the {\it Herschel} Astrophysical Terahertz Large Area Survey, the {\it Herschel} Virgo Cluster Survey, the {\it Herschel} Fornax Cluster Survey, the {\it Herschel} Stripe 82 Survey and the {\it Herschel} Multi-tiered Extragalactic Survey, totalling $880$ deg$^2$, or $\sim 3$\% of the sky in total. Our sample selection is serendipitous, based only on whether the X-ray position of a GRB lies within a large-scale {\it Herschel} survey -- therefore our sample can be considered completely unbiased. Using deep data at wavelengths of 100\,--\,500$\,\mu$m, we tentatively detected 1 out of 20 GRB hosts located in these fields. We constrain their dust masses and star formation rates (SFRs), and discuss these in the context of recent measurements of submillimetre galaxies and ultraluminous infrared galaxies. The average far-infrared flux of our sample gives an upper limit on SFR of $<114\,{\rm M}\odot\,\mbox{yr}^{-1}$. The detection rate of GRB hosts is consistent with that predicted assuming that GRBs trace the cosmic SFR density in an unbiased way, i.e. that the fraction of GRB hosts with $\mbox{SFR}>500\,{\rm M}\odot\,\mbox{yr}^{-1}$ is consistent with the contribution of such luminous galaxies to the cosmic star formation density.

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

We report an unexpected variation in the positional offset distributions between Herschel-ATLAS sub-millimetre (submm) sources and their optical associations, depending on both 250-{\mu}m signal-to-noise ratio and 250/350-{\mu}m colour. We show that redder and brighter submm sources have optical associations with a broader distribution of positional offsets than would be expected if these offsets were due to random positional errors in the source extraction. The observation can be explained by two possible effects: either red submm sources trace a more clustered population than blue ones, and their positional errors are increased by confusion; or red submm sources are generally at high redshifts and are frequently associated with low-redshift lensing structures which are identified as false counterparts. We perform various analyses of the data, including the multiplicity of optical associations, the redshift and magnitude distributions in H-ATLAS in comparison to HerMES, and simulations of weak lensing, and we conclude that the effects are most likely to be explained by widespread weak lensing of Herschel-SPIRE sources by foreground structures. This has important consequences for counterpart identification and derived redshift distributions and luminosity functions of submm surveys.

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

[Abridged] Aims: This work focuses on one lensed system, HATLAS J142935.3-002836 (H1429-0028), selected in the Herschel-ATLAS field. Gathering a rich, multi-wavelength dataset, we aim to confirm the lensing hypothesis and model the background source's morphology and dynamics, as well as to provide a full physical characterisation. Methods: Multi-wavelength high-resolution data is utilised to assess the nature of the system. A lensing-analysis algorithm which simultaneously fits different wavebands is adopted to characterise the lens. The background galaxy dynamical information is studied by reconstructing the 3-D source-plane of the ALMA CO(J:4-3) transition. Near-IR imaging from HST and Keck-AO allows to constrain rest-frame optical photometry independently for the foreground and background systems. Physical parameters (such as stellar and dust masses) are estimated via modelling of the spectral energy distribution taking into account source blending, foreground obscuration, and differential magnification. Results: The system comprises a foreground edge-on disk galaxy (at z_sp=0.218) with an almost complete Einstein ring around it. The background source (at z_sp=1.027) is magnified by a factor of ~8-10 depending on wavelength. It is comprised of two components and a tens of kpc long tidal tail resembling the Antennae merger. As a whole, the system is a massive stellar system (1.32[-0.41,+0.63] x1E11 Mo) forming stars at a rate of 394+-90 Mo/yr, and has a significant gas reservoir M_ISM = 4.6+-1.7 x1E10 Mo. Its depletion time due to star formation alone is thus expected to be tau_SF=M_ISM/SFR=117+-51 Myr. The dynamical mass of one of the components is estimated to be 5.8+-1.7 x1E10 Mo, and, together with the photometric total mass estimate, it implies that H1429-0028 is a major merger system (1:2.8[-1.5,+1.8]).

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

We present a comparison of the physical properties of a rest-frame $250\mu$m selected sample of massive, dusty galaxies from $0<z<5.3$. Our sample comprises 29 high-redshift submillimetre galaxies (SMGs) from the literature, and 843 dusty galaxies at $z<0.5$ from the Herschel-ATLAS, selected to have a similar stellar mass to the SMGs. The $z>1$ SMGs have an average SFR of $390^{+80}_{-70}\,$M$_\odot$yr$^{-1}$ which is 120 times that of the low-redshift sample matched in stellar mass to the SMGs (SFR$=3.3\pm{0.2}$ M$_\odot$yr$^{-1}$). The SMGs harbour a substantial mass of dust ($1.2^{+0.3}_{-0.2}\times{10}^9\,$M$_\odot$), compared to $(1.6\pm0.1)\times{10}^8\,$M$_\odot$ for low-redshift dusty galaxies. At low redshifts the dust luminosity is dominated by the diffuse ISM, whereas a large fraction of the dust luminosity in SMGs originates from star-forming regions. At the same dust mass SMGs are offset towards a higher SFR compared to the low-redshift H-ATLAS galaxies. This is not only due to the higher gas fraction in SMGs but also because they are undergoing a more efficient mode of star formation, which is consistent with their bursty star-formation histories. The offset in SFR between SMGs and low-redshift galaxies is similar to that found in CO studies, suggesting that dust mass is as good a tracer of molecular gas as CO.

We apply a chemical evolution model to investigate the sources and evolution of dust in a sample of 26 high-redshift ($z>1$) submillimetre galaxies (SMGs) from the literature, with complete photometry from ultraviolet to the submillimetre. We show that dust produced only by low-intermediate mass stars falls a factor 240 short of the observed dust masses of SMGs, the well-known `dust-budget crisis'. Adding an extra source of dust from supernovae can account for the dust mass in 19 per cent of the SMG sample. Even after accounting for dust produced by supernovae the remaining deficit in the dust mass budget provides support for higher supernova yields, substantial grain growth in the interstellar medium or a top-heavy IMF. Including efficient destruction of dust by supernova shocks increases the tension between our model and observed SMG dust masses. The models which best reproduce the physical properties of SMGs have a rapid build-up of dust from both stellar and interstellar sources and minimal dust destruction. Alternatively, invoking a top-heavy IMF or significant changes in the dust grain properties can solve the dust budget crisis only if dust is produced by both low mass stars and supernovae and is not efficiently destroyed by supernova shocks.

We have determined the mass-density radial profiles of the first five strong gravitational lens systems discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We present an enhancement of the semi-linear lens inversion method of Warren & Dye which allows simultaneous reconstruction of several different wavebands and apply this to dual-band imaging of the lenses acquired with the Hubble Space Telescope. The five systems analysed here have lens redshifts which span a range, 0.22<z<0.94. Our findings are consistent with other studies by concluding that: 1) the logarithmic slope of the total mass density profile steepens with decreasing redshift; 2) the slope is positively correlated with the average total projected mass density of the lens contained within half the effective radius and negatively correlated with the effective radius; 3) the fraction of dark matter contained within half the effective radius increases with increasing effective radius and increases with redshift.

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}

We report the redshift of HATLAS J132427.0+284452 (hereafter HATLAS J132427), a gravitationally lensed starburst galaxy, the first determined 'blind' by the Herschel Space Observatory. This is achieved via the detection of [C II] consistent with z = 1.68 in a far-infrared spectrum taken with the SPIRE Fourier Transform Spectrometer. We demonstrate that the [C II] redshift is secure via detections of CO J = 2 - 1 and 3 - 2 using the Combined Array for Research in Millimeter-wave Astronomy and the Institut de Radioastronomie Millimetrique's Plateau de Bure Interferometer. The intrinsic properties appear typical of high-redshift starbursts despite the high lensing-amplified fluxes, proving the ability of the FTS to probe this population with the aid of lensing. The blind detection of [C II] demonstrates the potential of the SAFARI imaging spectrometer, proposed for the much more sensitive SPICA mission, to determine redshifts of multiple dusty galaxies simultaneously without the benefit of lensing.

We present an analysis of CO molecular gas tracers in a sample of 500{\mu}m-selected Herschel-ATLAS galaxies at z<0.05 (cz<14990km/s). Using 22-500{\mu}m photometry from WISE, IRAS and Herschel, with HI data from the literature, we investigate correlations between warm and cold dust, and tracers of the gas in different phases. The correlation between global CO(3-2) line fluxes and FIR-submillimetre fluxes weakens with increasing IR wavelength ({\lambda}>60{\mu}m), as a result of colder dust being less strongly associated with dense gas. Conversely, CO(2-1) and HI line fluxes both appear to be better correlated with longer wavelengths, suggesting that cold dust is more strongly associated with diffuse atomic and molecular gas phases, consistent with it being at least partially heated by radiation from old stellar populations. The increased scatter at long wavelengths implies that submillimetre fluxes are a poorer tracer of SFR. Fluxes at 22 and 60{\mu}m are also better correlated with diffuse gas tracers than dense CO(3-2), probably due to very-small-grain emission in the diffuse interstellar medium, which is not correlated with SFR. The FIR/CO luminosity ratio and the dust mass/CO luminosity ratio both decrease with increasing luminosity, as a result of either correlations between mass and metallicity (changing CO/H2) or between CO luminosity and excitation [changing CO(3-2)/CO(1-0)].

We present the results from nine years of optically monitoring the gravitationally lensed z=0.658 quasar RX J1131-1231. The R-band light curves of the four individual images of the quasar were obtained using deconvolution photometry for a total of 707 epochs. Several sharp quasar variability features strongly constrain the time delays between the quasar images. Using three different numerical techniques, we measure these delays for all possible pairs of quasar images while always processing the four light curves simultaneously. For all three methods, the delays between the three close images A, B, and C are compatible with being 0, while we measure the delay of image D to be 91 days, with a fractional uncertainty of 1.5% (1 sigma), including systematic errors. Our analysis of random and systematic errors accounts in a realistic way for the observed quasar variability, fluctuating microlensing magnification over a broad range of temporal scales, noise properties, and seasonal gaps. Finally, we find that our time-delay measurement methods yield compatible results when applied to subsets of the data.

Using IRAM PdBI we report the detection of H2O in six new lensed ultra-luminous starburst galaxies at high redshift, discovered in the Herschel H-ATLAS survey. The sources are detected either in the 2_{02}-1_{11} or 2_{11}-2_{02} H_2O emission lines with integrated line fluxes ranging from 1.8 to 14 Jy.km/s. The corresponding apparent luminosities are mu x L_H2O ~ 3-12 x 10^8 Lo, where mu is the lensing magnification factor (3 < mu < 12). These results confirm that H2O lines are among the strongest molecular lines in such galaxies, with intensities almost comparable to those of the high-J CO lines, and same profiles and line widths (200-900 km/s) as the latter. With the current sensitivity of PdBI, H2O can therefore easily be detected in high-z lensed galaxies (with F(500um) > 100 mJy) discovered in the Herschel surveys. Correcting the luminosities for lensing amplification, L_H2O is found to have a strong dependence on the IR luminosity, varying as ~L_IR^{1.2}. This relation which needs to be confirmed with better statistics, may indicate a role of radiative (IR) excitation of the H2O lines, and implies that high-z galaxies with L_IR >~ 10^13 Lo tend to be very strong emitters in H2O, that have no equivalent in the local universe.

We have constructed a sample of radio-loud objects with optical spectroscopy from the Galaxy and Mass Assembly (GAMA) project over the Herschel-ATLAS Phase 1 fields. Classifying the radio sources in terms of their optical spectra, we find that strong-emission-line sources (`high-excitation radio galaxies') have, on average, a factor ~4 higher 250-micron Herschel luminosity than weak-line (`low-excitation') radio galaxies and are also more luminous than magnitude-matched radio-quiet galaxies at the same redshift. Using all five H-ATLAS bands, we show that this difference in luminosity between the emission-line classes arises mostly from a difference in the average dust temperature; strong-emission-line sources tend to have comparable dust masses to, but higher dust temperatures than, radio galaxies with weak emission lines. We interpret this as showing that radio galaxies with strong nuclear emission lines are much more likely to be associated with star formation in their host galaxy, although there is certainly not a one-to-one relationship between star formation and strong-line AGN activity. The strong-line sources are estimated to have star-formation rates at least a factor 3-4 higher than those in the weak-line objects. Our conclusion is consistent with earlier work, generally carried out using much smaller samples, and reinforces the general picture of high-excitation radio galaxies as being located in lower-mass, less evolved host galaxies than their low-excitation counterparts.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

Chemical and physical Earth models agree little as to the radioactive power of the planet. Each predicts a range of radioactive powers, overlapping slightly with the other at about 24 TW, and together spanning 14-46 TW. Approximately 20 % of this radioactive power (3-8 TW) escapes to space in the form of geo-neutrinos. The remaining 11-38 TW heats the planet with significant geo-dynamical consequences, appearing as the radiogenic component of the 43-49 TW surface heat flow. The non-radiogenic component of the surface heat flow (5-38 TW) is presumably primordial, a legacy of the formation and early evolution of the planet. A constraining measurement of radiogenic heating provides insights to the thermal history of the Earth and potentially discriminates chemical and physical Earth models. Radiogenic heating in the planet primarily springs from unstable nuclides of uranium, thorium, and potassium. The paths to their stable daughter nuclides include nuclear beta decays, producing geo-neutrinos. Large sub-surface detectors efficiently record the energy but not the direction of the infrequent interactions of the highest energy geo-neutrinos, originating only from uranium and thorium. The measured energy spectrum of the interactions estimates the relative amounts of these heat-producing elements, while the intensity estimates planetary radiogenic power. Recent geo-neutrino observations in Japan and Italy find consistent values of radiogenic heating. The combined result mildly excludes the lowest model values of radiogenic heating and, assuming whole mantle convection, identifies primordial heat loss. Future observations have the potential to measure radiogenic heating with better precision, further constraining geological models and the thermal evolution of the Earth.

We have used data from the Herschel-ATLAS at 250, 350 and 500 \mu m to determine the far-infrared (FIR) properties of 50 Broad Absorption Line Quasars (BAL QSOs). Our sample contains 49 high-ionization BAL QSOs (HiBALs) and 1 low-ionization BAL QSO (LoBAL) which are compared against a sample of 329 non-BAL QSOs. These samples are matched over the redshift range 1.5 \leq z < 2.3 and in absolute i-band magnitude over the range -28 \leq M_{i} \leq -24. Of these, 3 BAL QSOs (HiBALs) and 27 non-BAL QSOs are detected at the > 5 sigma level. We calculate star-formation rates (SFR) for our individually detected HiBAL QSOs and the non-detected LoBAL QSO as well as average SFRs for the BAL and non-BAL QSO samples based on stacking the Herschel data. We find no difference between the HiBAL and non-BAL QSO samples in the FIR, even when separated based on differing BAL QSO classifications. Using Mrk 231 as a template, the weighted mean SFR is estimated to be \approx240\pm21 M_{\odot} yr^{-1} for the full sample, although this figure should be treated as an upper limit if AGN-heated dust makes a contribution to the FIR emission. Despite tentative claims in the literature, we do not find a dependence of {\sc C\,iv} equivalent width on FIR emission, suggesting that the strength of any outflow in these objects is not linked to their FIR output. These results strongly suggest that BAL QSOs (more specifically HiBALs) can be accommodated within a simple AGN unified scheme in which our line-of-sight to the nucleus intersects outflowing material. Models in which HiBALs are caught towards the end of a period of enhanced spheroid and black-hole growth, during which a wind terminates the star-formation activity, are not supported by the observed FIR properties.

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 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 use spitzer-IRAC data to identify near-infrared counterparts to submillimeter galaxies detected with Herschel-SPIRE at 250um in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). Using a likelihood ratio analysis we identify 146 reliable IRAC counterparts to 123 SPIRE sources out of the 159. We find that, compared to the field population, the SPIRE counterparts occupy a distinct region of 3.6 and 4.5um color-magnitude space, and we use this property to identify a further 23 counterparts to 13 SPIRE sources. The IRAC identification rate of 86% is significantly higher than those that have been demonstrated with wide-field ground-based optical and near-IR imaging of Herschel fields. We estimate a false identification rate of 3.6%, corresponding to 4 to 5 sources. Among the 73 counterparts that are undetected in SDSS, 57 have both 3.6 and 4.5um coverage. Of these 43 have [3.6] - [4.5]> 0 indicating that they are likely to be at z > 1.4. Thus, ~ 40% of identified SPIRE galaxies are likely to be high redshift (z > 1.4) sources. We discuss the statistical properties of the IRAC-identified SPIRE galaxy sample including far-IR luminosities, dust temperatures, star-formation rates, and stellar masses. The majority of our detected galaxies have 10^10 to 10^11 L_sun total IR luminosities and are not intense starbursting galaxies as those found at z ~ 2, but they have a factor of 2 to 3 above average specific star-formation rates compared to near-IR selected galaxy samples.

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

The Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) provides an unprecedented opportunity to search for blazars at sub-mm wavelengths. We cross-matched the FIRST radio source catalogue with the 11655 sources brighter than 35 mJy at 500{\mu}m in the \sim 135 square degrees of the sky covered by the H-ATLAS equatorial fields at 9 h and 15 h, plus half of the field at 12 h. We found that 379 of the H-ATLAS sources have a FIRST counterpart within 10 arcsec, including 8 catalogued blazars (plus one known blazar that was found at the edge of one the H-ATLAS maps). To search for additional blazar candidates we have devised new diagnostic diagrams and found that known blazars occupy a region of the log(S500{\mu}m/S350{\mu}m) vs. log(S500{\mu}m/S1.4GHz) plane separated from that of the other sub-mm sources with radio counterparts. Using this diagnostic we have selected 12 further candidates that turn out to be scattered in the (r-z) vs. (u-r) plane or in the WISE colour-colour diagram proposed by Massaro et al. (2012), where known blazars are concentrated in well defined strips. This suggests that the majority of them either are not blazars or have spectral energy distributions contaminated by their host galaxies. A significant fraction of true blazars are found to be hosted by star-forming galaxies. This finding, supported by an analysis of blazars detected in Planck 545 and 857 GHz bands, is at odds with the notion that blazar hosts are passive ellipticals and indicates that the sub-mm selection is providing a novel prospect on blazar properties. Based on an inspection of the available photometric data, including the WISE all-sky survey, the unpublished VIKING survey and new radio observations, we tentatively estimate that there are 11 blazars with synchrotron flux density S500{\mu}m > 35mJy over the considered area. This result already allows us to constrain blazar evolution models.

We present high-resolution maps of stars, dust, and molecular gas in a strongly lensed submillimeter galaxy (SMG) at z = 3.259. HATLAS12--00 is selected from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) as a strong lens candidate mainly based on its unusually high 500um flux density (~300 mJy). It is the only high-redshift Planck detection in the 130 deg^2 H-ATLAS Phase 1 area. Keck Adaptive Optics images reveal a quadruply imaged galaxy in the K-band while the Submillimeter Array and the Extended Very Large Array show doubly imaged 880um and CO(1-0) sources, indicating differentiated distributions of the various components in the galaxy. In the source plane, the stars reside in three major kpc-scale clumps extended over ~1.6 kpc, the dust in a compact (~1 kpc) region ~3 kpc north of the stars, and the cold molecular gas in an extended (~7 kpc) disk ~5 kpc northeast of the stars. The emission from the stars, dust, and gas are magnified by ~17, 8, and 7 times, respectively, by four lensing galaxies at z ~ 1. Intrinsically, the lensed galaxy is a warm (T_dust ~ 40-65 K), hyper-luminous (L_IR ~ 1.7e13 Lsun; SFR ~ 2000 Msun/yr), gas-rich (M_gas/M_baryon ~ 70%), young (M_stellar/SFR ~ 20 Myr), and short-lived (M_gas/SFR ~ 40 Myr) starburst, without a significant active galactic nucleus. With physical properties similar to unlensed z > 2 SMGs, HATLAS12--00 offers a detailed view of a typical SMG through a powerful cosmic microscope.

We report measurements of the carbon monoxide ground state rotational transition (12C16O J = 1--0) with the Zpectrometer ultra-wideband spectrometer on the 100-m diameter Green Bank Telescope. The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 24 targets identified by Herschel-ATLAS photometric colors from the SPIRE instrument. Nine of the CO measurements are new redshift determinations, substantially adding to the number of detections of galaxies with rest-frame peak submillimeter emission near 100um. The CO detections confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies (DSFGs). The CO redshift distribution of the 350um-selected galaxies is strikingly similar to the optical redshifts of 850um-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by foreground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding approximate gravitational lens magnifications mu from spectroscopic data alone, yielding mu ~ 3--20. Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 10^13 L_sun.

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.

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

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

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.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

We present the dust properties and star-formation histories of local submillimetre-selected galaxies in Herschel-ATLAS, classified by optical morphology. The early-type galaxies (ETGs) that are detected contain as much dust as typical spirals, and form a unique sample that has been blindly selected at submillimetre wavelengths. Comparing H-ATLAS galaxies to a control sample of optically selected galaxies, we find 5.5% of luminous ETGs are detected in H-ATLAS. The H-ATLAS ETGs contain a significant mass of cold dust: the mean dust mass is 5.5x10^7 Msun, with individual galaxies ranging from 9x10^5-4x10^8 Msun. This is comparable to that of spirals in our sample, and is an order of magnitude more dust than that found for the control ETGs, which have a median dust mass inferred from stacking of (0.8-4.0)x10^6 Msun. The ETGs detected in H-ATLAS have bluer NUV-r colours, higher specific star-formation rates and younger stellar populations than ETGs which are optically selected, and may be transitioning from the blue cloud to the red sequence. We also find that H-ATLAS and control ETGs inhabit similar low-density environments. We conclude that the dust in H-ATLAS and control ETGs cannot be solely from stellar sources, and a large contribution from dust formed in the ISM or external sources is required. Alternatively, dust destruction may not be as efficient as predicted. We also explore the properties of the most passive spiral galaxies in our sample with SSFR<10^-11/yr. We find these passive spirals have lower dust-to-stellar mass ratios, higher stellar masses and older stellar population ages than normal spirals. The passive spirals inhabit low density environments similar to those of the normal spiral galaxies in our sample. This shows that the processes which turn spirals passive do not occur solely in the intermediate density environments of group and cluster outskirts. (Abridged)

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.

The Herschel survey, H-ATLAS, with its large areal coverage, has recently discovered a number of bright, strongly lensed high-z submillimeter galaxies. The strong magnification makes it possible to study molecular species other than CO, which are otherwise difficult to observe in high-z galaxies. Among the lensed galaxies already identified by H-ATLAS, the source J090302.9-014127B (SDP.17b) at z = 2.305 is remarkable due to its excitation conditions and a tentative detection of the H2O 202-111 emission line (Lupu et al. 2010). We report observations of this line in SDP.17b using the IRAM interferometer equipped with its new 277- 371GHz receivers. The H2O line is detected at a redshift of z = 2.3049+/-0.0006, with a flux of 7.8+/-0.5 Jy km s-1 and a FWHM of 250+/-60 km s-1. The new flux is 2.4 times weaker than the previous tentative detection, although both remain marginally consistent within 1.6-sigma. The intrinsic line luminosity and ratio of H2O(202-111)/CO8-7 seem comparable with those of the nearby starburst/enshrouded-AGN Mrk 231, suggesting that SDP.17b could also host a luminous AGN. The detection of a strong H2O 202-111 line in SDP.17b implies an efficient excitation mechanism of the water levels that must occur in very dense and warm interstellar gas.

We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.

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 report ground-based follow-up observations of the exceptional source, ID141, one the brightest sources detected so far in the H-ATLAS cosmological survey. ID141 was observed using the IRAM 30-meter telescope and Plateau de Bure interferometer (PdBI), the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX) submillimeter telescope to measure the dust continuum and emission lines of the main isotope of carbon monoxide and carbon ([C I] and [C II]). The detection of strong CO emission lines with the PdBI confirms that ID141 is at high redshift (z=4.243 +/- 0.001). The strength of the continuum and emission lines suggests that ID141 is gravitationally lensed. The width (Delta V (FWHM) ~ 800 km/s}) and asymmetric profiles of the CO and carbon lines indicate orbital motion in a disc or a merger. The properties derived for ID141 are compatible with a ultraluminous (L_FIR ~ 8.5 +/- 0.3 x 10^13/mu_L Lsun, where mu_L is the amplification factor, dense (n ~ 10^4 cm^-3) and warm (T_kin ~ 40K) starburst galaxy, with an estimated star-formation rate of (0.7 to 1.7) x 10^4/mu_L Msun/yr. The carbon emission lines indicate a dense (n ~ 10^4 cm^-3) Photo-Dominated Region, illuminated by a far-UV radiation field a few thousand times more intense than that in our Galaxy. In conclusion, the physical properties of the high-z galaxy, ID141, are remarkably similar to those of local ultraluminous infrared galaxies.

We investigate the ultraviolet and optical properties and environment of low redshift galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration data. We use the Sloan Digital Sky Survey seventh release and the Galaxy And Mass Assembly database to select galaxies with r_Petro < 19.0 mag in the redshift range 0.02 < z < 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor contribution from red systems which are highly obscured by dust. In addition we find that the colour of a galaxy rather than the local density of its environment determines whether it is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously detected by both PACS and SPIRE is 25K \pm 4K, independent of environment. This analysis provides a glimpse of the potential of the H-ATLAS data to investigate the submillimeter properties of galaxies in the local universe.

The intergalactic medium was not completely reionized until approximately a billion years after the Big Bang, as revealed by observations of quasars with redshifts of less than 6.5. It has been difficult to probe to higher redshifts, however, because quasars have historically been identified in optical surveys, which are insensitive to sources at redshifts exceeding 6.5. Here we report observations of a quasar (ULAS J112001.48+064124.3) at a redshift of 7.085, which is 0.77 billion years after the Big Bang. ULAS J1120+0461 had a luminosity of 6.3x10^13 L_Sun and hosted a black hole with a mass of 2x10^9 M_Sun (where L_Sun and M_Sun are the luminosity and mass of the Sun). The measured radius of the ionized near zone around ULAS J1120+0641 is 1.9 megaparsecs, a factor of three smaller than typical for quasars at redshifts between 6.0 and 6.4. The near zone transmission profile is consistent with a Ly alpha damping wing, suggesting that the neutral fraction of the intergalactic medium in front of ULAS J1120+0641 exceeded 0.1.

We report the detection of a significant excess in the surface density of far-infrared sources from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) within ~1 Mpc of the centres of 66 optically-selected clusters of galaxies in the SDSS with <z>~0.25. From the analysis of the multiwavelength properties of their counterparts we conclude that the far-infrared emission is associated with dust-obscured star formation and/or active galactic nuclei within galaxies in the clusters themselves. The excess reaches a maximum at a radius of ~0.8 Mpc, where we find 1.0\pm0.3 S_250um>34 mJy sources on average per cluster above what would be expected for random field locations. If the far-infrared emission is dominated by star formation (as opposed to AGN) then this corresponds to an average star formation rate of ~7 M_sun/yr per cluster in sources with L_IR>5d10 L_sun. Although lensed sources make a negligible contribution to the excess signal, a fraction of the sources around the clusters could be gravitationally lensed, and we have identified a sample of potential cases of cluster-lensed Herschel sources that could be targeted in follow-up studies.

The Herschel-ATLAS is a survey of 550 square degrees with the Herschel Space Observatory in five far--infrared and submillimetre bands. The first data for the survey, observations of a field 4x4 sq. degrees in size, were taken during the Science Demonstration Phase, and reach a 5 sigma noise level of 33 mJy/beam at 250 microns. This paper describes the source extraction methods used to create the corresponding Science Demonstration Phase catalogue, which contains 6876 sources, selected at 250 microns, within ~14 sq. degrees. SPIRE sources are extracted using a new method specifically developed for Herschel data; PACS counterparts of these sources are identified using circular apertures placed at the SPIRE positions. Aperture flux densities are measured for sources identified as extended after matching to optical wavelengths. The reliability of this catalogue is also discussed, using full simulated maps at the three SPIRE bands. These show that a significant number of sources at 350 and 500 microns have undergone flux density enhancements of up to a factor of ~2, due mainly to source confusion. Correction factors are determined for these effects. The SDP dataset and corresponding catalogue will be available from http://www.h-atlas.org/.

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 use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

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.

We carry out a multi-wavelength study of individual galaxies detected by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) and identified at other wavelengths, using data spanning the radio to the ultraviolet (UV). We develop a Monte Carlo method to account for flux boosting, source blending, and correlations among bands, which we use to derive deboosted far-infrared (FIR) luminosities for our sample. We estimate total star-formation rates for BLAST counterparts with z < 0.9 by combining their FIR and UV luminosities. Star formation is heavily obscured at L_FIR > 10^11 L_sun, z > 0.5, but the contribution from unobscured starlight cannot be neglected at L_FIR < 10^11 L_sun, z < 0.25. We assess that about 20% of the galaxies in our sample show indication of a type-1 active galactic nucleus (AGN), but their submillimeter emission is mainly due to star formation in the host galaxy. We compute stellar masses for a subset of 92 BLAST counterparts; these are relatively massive objects, with a median mass of ~10^11 M_sun, which seem to link the 24um and SCUBA populations, in terms of both stellar mass and star-formation activity. The bulk of the BLAST counterparts at z<1 appear to be run-of-the-mill star-forming galaxies, typically spiral in shape, with intermediate stellar masses and practically constant specific star-formation rates. On the other hand, the high-z tail of the BLAST counterparts significantly overlaps with the SCUBA population, in terms of both star-formation rates and stellar masses, with observed trends of specific star-formation rate that support strong evolution and downsizing.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

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 study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $\beta$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = \^a H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

We have reduced the data taken with the Spectral and Photometric Imaging Receiver (SPIRE) photometer on board the Herschel Space Observatory in the Science Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We describe the data reduction, which poses specific challenges, both because of the sheer size of the data, and because only two scans are made for each region. We implement effective solutions to process the bolometric timelines into maps, and show that correlations among detectors are negligible, and that the photometer is stable on time scales up to 250 s. This is longer than the time the telescope takes to cross the observed sky region, and it allows us to use naive binning methods for an optimal reconstruction of the sky emission. The maps have equal contribution of confusion and white instrumental noise, and the latter is estimated to 5.3, 6.4, and 6.7 mJy/beam (1-{\sigma}), at 250, 350, and 500 \mu{m}, respectively. This pipeline is used to reduce other H-ATLAS observations, as they became available, and we discuss how it can be used with the optimal map maker implemented in the Herschel Interactive Processing Environment (HIPE), to improve computational efficiency and stability. The SDP dataset is available from http://www.h-atlas.org/.

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.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 2008 on the 3.9-m Anglo-Australian Telescope using the AAOmega fibre-fed spectrograph facility to acquire spectra with a resolution of R~1300 for 120,862 SDSS selected galaxies. The target catalogue constitutes three contiguous equatorial regions centred at 9h (G09), 12h (G12) and 14.5h (G15) each of 12 x 4 sq.deg to limiting fluxes of r < 19.4, r < 19.8, and r < 19.4 mag respectively (and additional limits at other wavelengths). Spectra and reliable redshifts have been acquired for over 98 per cent of the galaxies within these limits. Here we present the survey footprint, progression, data reduction, redshifting, re-redshifting, an assessment of data quality after 3 years, additional image analysis products (including ugrizYJHK photometry, Sersic profiles and photometric redshifts), observing mask, and construction of our core survey catalogue (GamaCore). From this we create three science ready catalogues: GamaCoreDR1 for public release, which includes data acquired during year 1 of operations within specified magnitude limits (February 2008 to April 2008); GamaCoreMainSurvey containing all data above our survey limits for use by the GAMA team and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched to Herschel-ATLAS Science Demonstration data. These catalogues along with the associated spectra, stamps and profiles can be accessed via the GAMA website: http://www.gama-survey.org/

In order to generate credible 0.1-2 {\mu}m SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate S\'ersic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, {\alpha} +/- 0.014, {\Phi}* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or S\'ersic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We have measured star formation histories (SFHs) and stellar masses of galaxies detected by the Balloon-borne Large Aperture Sub-millimetre Telescope (BLAST) over approximately 9 square degrees centred on the Chandra Deep Field South. We have applied the recently developed SFH reconstruction method of Dye et al. to optical, near-infrared and mid-infrared photometry of 92 BLAST galaxies. We find significant differences between the SFHs of low mass (<10^11 M_sol) and high mass (>10^11 M_sol) galaxies. On average, low mass galaxies exhibit a dominant late burst of star formation which creates a large fraction of their stellar mass. Conversely, high mass systems tend to have a significant amount of stellar mass that formed much earlier. We also find that the high mass SFHs evolve more strongly than the low mass SFHs. These findings are consistent with the phenomenon of downsizing observed in optically selected samples of galaxies.

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.

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.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

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.

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

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 have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

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 new efficient technique for measuring evolution of the galaxy luminosity function. The method reconstructs the evolution over the luminosity-redshift plane using any combination of three input dataset types: 1) number counts, 2) galaxy redshifts, 3) integrated background flux measurements. The evolution is reconstructed in adaptively sized regions of the plane according to the input data as determined by a Bayesian formalism. We demonstrate the performance of the method using a range of different synthetic input datasets. We also make predictions of the accuracy with which forthcoming surveys conducted with SCUBA2 and the Herschel Space Satellite will be able to measure evolution of the sub-millimetre luminosity function using the method.

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

We present a comparison of the SCUBA Half Degree Extragalactic Survey (SHADES) at 450, 850 and 1100 microns with deep guaranteed time 15 microns AKARI FU-HYU survey data and Spitzer guaranteed time data at 3.6-24 microns in the Lockman Hole East. The AKARI data was analysed using bespoke software based in part on the drizzling and minimum-variance matched filtering developed for SHADES, and was cross calibrated against Infrared Space Observatory (ISO) fluxes. Our stacking analyses find AKARI 15um galaxies with >~200 microJy contribute >10% of the 450 micron background, but only <4% of the 1100 micron background, suggesting that different populations contribute at mm-wavelengths. We confirm our earlier result that the ultra-deep 450 micron SCUBA-2 Cosmology Survey will be dominated by populations already detected by AKARI and Spitzer mid-infrared surveys. The superb mid-infrared wavelength coverage afforded by combining Spitzer and AKARI photometry is an excellent diagnostic of AGN contributions, and we find that (23-52)% of submm-selected galaxies have AGN bolometric fractions f_AGN>0.3.

We describe the spectroscopic target selection for the Galaxy And Mass Assembly (GAMA) survey. The input catalogue is drawn from the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS). The aim is to measure redshifts for galaxies in three 4x12 degree regions at 9h, 12h and 14.5h, on the celestial equator, with magnitude selections r < 19.4, z < 18.2 and K(AB) < 17.6 over all three regions, and r < 19.8 in the 12-h region. The target density is 1080 per sq. deg. in the 12-h region and 720 per sq. deg. in the other regions. The average GAMA target density and area are compared with completed and ongoing galaxy redshift surveys. The GAMA survey implements a highly complete star-galaxy separation that jointly uses an intensity-profile separator (delta(sg) = r-band psf mag - model mag) as per the SDSS) and a colour separator. The colour separator is defined as delta(sg,jk) = J - K - f(g-i), where f(g-i) is a quadratic fit to the J-K colour of the stellar locus over the range 0.3 < g-i < 2.3. All galaxy populations investigated are well separated with delta(sg,jk) > 0.2. From two years out of a three-year AAOmega program on the Anglo-Australian Telescope, we have obtained 79599 unique galaxy redshifts. Previously known redshifts in the GAMA region bring the total up to 98497. The median galaxy redshift is 0.2 with 99% at z < 0.5. We present some of the global statistical properties of the survey, including K-band galaxy counts, colour-redshift relations and preliminary n(z).

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has recently surveyed ~8.7 deg^2 centered on GOODS-South at 250, 350, and 500 microns. In Dye et al. (2009) we presented the catalogue of sources detected at 5-sigma in at least one band in this field and the probable counterparts to these sources in other wavebands. In this paper, we present the results of a redshift survey in which we succeeded in measuring redshifts for 82 of these counterparts. The spectra show that the BLAST counterparts are mostly star-forming galaxies but not extreme ones when compared to those found in the Sloan Digital Sky Survey. Roughly one quarter of the BLAST counterparts contain an active nucleus. We have used the spectroscopic redshifts to carry out a test of the ability of photometric redshift methods to estimate the redshifts of dusty galaxies, showing that the standard methods work well even when a galaxy contains a large amount of dust. We have also investigated the cases where there are two possible counterparts to the BLAST source, finding that in at least half of these there is evidence that the two galaxies are physically associated, either because they are interacting or because they are in the same large-scale structure. Finally, we have made the first direct measurements of the luminosity function in the three BLAST bands. We find strong evolution out to z=1, in the sense that there is a large increase in the space-density of the most luminous galaxies. We have also investigated the evolution of the dust-mass function, finding similar strong evolution in the space-density of the galaxies with the largest dust masses, showing that the luminosity evolution seen in many wavebands is associated with an increase in the reservoir of interstellar matter in galaxies.

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

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.

[abridged] We present constraints on the nature of the first galaxies selected at 350 microns. The sample includes galaxies discovered in the deepest blank-field survey at 350 microns (in the Bootes Deep Field) and also later serendipitous detections in the Lockman Hole. Spectral energy distribution templates are fit to identified counterparts, and the sample is found to comprise IR-luminous galaxies at 1<z<3 predominantly powered by star formation. The first spectrum of a 350-micron-selected galaxy provides an additional confirmation, showing prominent dust grain features typically associated with star-forming galaxies. Compared to submillimeter galaxies selected at 850 and 1100 microns, galaxies selected at 350 microns have a similar range of far-infrared color temperatures. However, no 350-micron-selected sources are reliably detected at 850 or 1100 microns. Galaxies in our sample with redshifts 1<z<2 show a tight correlation between the far- and mid-infrared flux densities, but galaxies at higher redshifts show a large dispersion in their mid- to far-infrared colors. The 350-micron flux densities (15<S(350)<40 mJy) place these objects near the Herschel/SPIRE 350-micron confusion threshold, with the lower limit on the star formation rate density suggesting the bulk of the 350-micron contribution will come from less luminous infrared sources and normal galaxies.

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 directly measure redshift evolution in the mean physical properties (far-infrared luminosity, temperature, and mass) of the galaxies that produce the cosmic infrared background (CIB), using measurements from the Balloon-borne Large Aperture Sub-millimeter Telescope (BLAST), and Spitzer which constrain the CIB emission peak. This sample is known to produce a surface brightness in the BLAST bands consistent with the full CIB, and photometric redshifts are identified for all of the objects. We find that most of the 70 micron background is generated at z <~ 1 and the 500 micron background generated at z >~ 1. A significant growth is observed in the mean luminosity from ~ 10^9 - 10^12 L_sun, and in the mean temperature by 10 K, from redshifts 0< z < 3. However, there is only weak positive evolution in the comoving dust mass in these galaxies across the same redshift range. We also measure the evolution of the far-infrared luminosity density, and the star-formation rate history for these objects, finding good agreement with other infrared studies up to z ~1, exceeding the contribution attributed to optically-selected galaxies.

Optical and near-infrared (NIR) spectra are presented for ULAS J131911.29+095951.4 (hereafter ULAS J1319+0950), a new redshift z = 6.127 +/- 0.004 quasar discovered in the Third Data Release (DR3) of the UKIRT Infrared Deep Sky Survey (UKIDSS). The source has Y_Vega = 19.10 +/- 0.03, corresponding to M_1450,AB = -27.12, which is comparable to the absolute magnitudes of the z ~= 6 quasars discovered in the SDSS. ULAS J1319+0950 was, in fact, registered by SDSS as a faint source with z_AB = 20.1 +/- 0.1, just below the signal--to--noise ratio limit of the SDSS high-redshift quasar survey. The faint z-band magnitude is a consequence of the weak Ly alpha/N V emission line, which has a rest-frame equivalent width of only ~20 A and provides only a small boost to the z-band flux. Nevertheless, there is no evidence from this UKIDSS-based search for a significant population of high-redshift quasars with weak emission lines. The Ly alpha optical depth to ULAS J1319+0950 is consistent with that measured towards similarly distant SDSS quasars, implying that results from optical- and NIR-selected quasars may be combined. Also presented is a new NIR-spectrum of the previously discovered UKIDSS quasar ULAS J020332.38+001229.2, which reveals the object to be a broad absorption line quasar. The new spectrum shows that the emission line previously identified as Ly alpha is actually N V, leading to a revised redshift of z = 5.72, rather than z = 5.86 as previously estimated.

We have identified radio and/or mid-infrared counterparts to 198 out of 350 sources detected at >=5 sigma over ~ 9 square degrees centered on the Chandra Deep Field South (CDFS) by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) at 250, 350 and 500 um. We have matched 114 of these counterparts to optical sources with previously derived photometric redshifts and fitted SEDs to the BLAST fluxes and fluxes at 70 and 160 um acquired with the Spitzer Space Telescope. In this way, we have constrained dust temperatures, total far-infrared/sub-millimeter luminosities and star formation rates for each source. Our findings show that on average, the BLAST sources lie at significantly lower redshifts and have significantly lower rest-frame dust temperatures compared to submm sources detected in surveys conducted at 850 um. We demonstrate that an apparent increase in dust temperature with redshift in our sample arises as a result of selection effects. Finally, we provide the full multi-wavelength catalog of >= 5 sigma BLAST sources contained within the complete ~ 9 square degree survey area.

We present new optical and infrared photometry for a statistically complete sample of seven 1.1 mm selected sources with accurate Submillimetre Array coordinates. We determine photometric redshifts for four of the seven sources of 4.47, 4.50, 1.49 and 0.64. Of the other three sources two are undetected at optical wavelengths down to the limits of very deep Subaru and Canada-France-Hawaii Telescope images ($\sim$27 mag AB, i band) and the photometry of the remaining source is corrupted by a bright nearby galaxy. The sources with the highest redshifts are at higher redshifts than all but one of the $\sim$200 sources taken from the largest recent 850 $\mu$m surveys, which may indicate that 1.1 mm surveys are more efficient at finding sources at very high redshifts than 850 $\mu$m surveys. We investigate the evolution of the number density with redshift of our sample using a banded $V_{e}/V_{a}$ analysis and find no evidence for a redshift cutoff, although the number of sources is very small. We also perform the same analysis on a statistically complete sample of 38 galaxies selected at 850$\mu$m from the GOODS-N field and find evidence for a drop-off in the number density beyond $z\sim1$ and 2 for hot and cold dust dominated SMGs respectively, confirming the earlier conclusion of Wall, Pope & Scott. We also find strong evidence for the existence of two differently evolving sub-populations separated in luminosity, with a higher relative density of the high luminosity galaxies at higher redshifts.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

If Type-II supernovae - the evolutionary end points of short-lived, massive stars - produce a significant quantity of dust (>0.1 M_sun) then they can explain the rest-frame far-infrared emission seen in galaxies and quasars in the first Gyr of the Universe. Submillimetre observations of the Galactic supernova remnant, Cas A, provided the first observational evidence for the formation of significant quantities of dust in Type-II supernovae. In this paper we present new data which show that the submm emission from Cas A is polarised at a level significantly higher than that of its synchrotron emission. The orientation is consistent with that of the magnetic field in Cas A, implying that the polarised submm emission is associated with the remnant. No known mechanism would vary the synchrotron polarisation in this way and so we attribute the excess polarised submm flux to cold dust within the remnant, providing fresh evidence that cosmic dust can form rapidly. This is supported by the presence of both polarised and unpolarised dust emission in the north of the remnant, where there is no contamination from foreground molecular clouds. The inferred dust polarisation fraction is unprecedented (f_pol ~ 30%) which, coupled with the brief timescale available for grain alignment (<300 yr), suggests that supernova dust differs from that seen in other Galactic sources (where f_pol=2-7%), or that a highly efficient grain alignment process must operate in the environment of a supernova remnant.

Recent studies of galaxies ~2-3 Gyr after the Big Bang have revealed large, turbulent rotating systems. The existence of well-ordered rotation in galaxies during this peak epoch of cosmic star formation may suggest that gas accretion through cold streams is likely to be the dominant mode by which most star-forming galaxies at high redshift since major mergers can completely disrupt the observed velocity fields. However poor spatial resolution and sensitivity have hampered this interpretation, limiting the study to the largest and most luminous galaxies, which may have fundamentally different modes of assembly than more typical star forming galaxies. Here we report observations of a typical star forming galaxy at z=3.07 with a linear resolution of ~100 parsec. This spatial sampling is made possible by the combination of gravitational lensing and laser guide star adaptive optics. We find a well-ordered compact source in which molecular gas is being converted efficiently into stars, likely assembling a spheroidal bulge and disk similar to those seen in spiral galaxies at the present day.

The UKIDSS Galactic Plane Survey (GPS) is one of the five near infrared Public Legacy Surveys that are being undertaken by the UKIDSS consortium, using the Wide Field Camera on the United Kingdom Infrared Telescope. It is surveying 1868 sq.deg. of the northern and equatorial Galactic plane at Galactic latitudes -5<b<5 in the J, H and K filters and a ~200 sq.deg. area of the Taurus-Auriga-Perseus molecular cloud complex in these three filters and the 2.12 um (1-0) H_2 filter. It will provide data on ~2 billion sources. Here we describe the properties of the dataset and provide a user's guide for its exploitation. We also present brief Demonstration Science results from DR2 and from the Science Verification programme. These results illustrate how GPS data will frequently be combined with data taken in other wavebands to produce scientific results. The Demonstration Science includes studies of: (i) the star formation region G28.983-0.603, cross matching with Spitzer-GLIMPSE data to identify YSOs; (ii) the M17 nebula; (iii) H_2 emission in the rho Ophiuchi dark cloud; (iv) X-ray sources in the Galactic Centre; (v) external galaxies in the Zone of Avoidance; (vi) IPHAS-GPS optical-infrared spectrophotometric typing. (abridged).

A new method of determining galaxy star-formation histories (SFHs) is presented. Using the method, the feasibility of recovering SFHs with multi-band photometry is investigated. The method divides a galaxy's history into discrete time intervals and reconstructs the average rate of star formation in each interval. This directly gives the total stellar mass. A simple linear inversion solves the problem of finding the most likely discretised SFH for a given set of galaxy parameters. It is shown how formulating the method within a Bayesian framework lets the data simultaneously select the optimal regularisation strength and the most appropriate number of discrete time intervals for the reconstructed SFH. The method is demonstrated by applying it to mono-metallic synthetic photometric catalogues created with different input SFHs, assessing how the accuracy of the recovered SFHs and stellar masses depend on the photometric passband set, signal-to-noise and redshift. The results show that reconstruction of SFHs using multi-band photometry is possible, being able to distinguish an early burst of star formation from a late one, provided an appropriate passband set is used. Although the resolution of the recovered SFHs is on average inferior compared to what can be achieved with spectroscopic data, the multi-band approach can process a significantly larger number of galaxies per unit exposure time.

We model the extremely massive and luminous lens galaxy in the Cosmic Horseshoe Einstein ring system, recently discovered in the Sloan Digital Sky Survey. We use the semi-linear method of Warren & Dye (2003), which pixelises the source surface brightness distribution, to invert the Einstein ring for sets of parameterised lens models. Here, the method is refined by exploiting Bayesian inference to optimise adaptive pixelisation of the source plane and to choose between three differently parameterised models: a singular isothermal ellipsoid, a power law model and a NFW profile. The most probable lens model is the power law with a volume mass density that scales as r^(-1.96+/-0.02) and an axis ratio of ~0.8. The mass within the Einstein ring (i.e., within a cylinder with projected distance of ~30 kpc from the centre of the lens galaxy) is (5.02+/-0.09)*10^12 M_solar, and the mass-to-light ratio is ~30. Even though the lens lies in a group of galaxies, the preferred value of the external shear is almost zero. This makes the Cosmic Horseshoe unique amongst large separation lenses, as almost all the deflection comes from a single, very massive galaxy with little boost from the environment.

We present the discovery of fifteen new T2.5-T7.5 dwarfs (with estimated distances between ~24-93pc, identified in the first three main data releases of the UKIRT Infrared Deep Sky Survey. This brings the total number of T dwarfs discovered in the Large Area Survey (to date) to 28. These discoveries are confirmed by near infrared spectroscopy, from which we derive spectral types on the unified scheme of Burgasser et al. (2006). Seven of the new T dwarfs have spectral types of T2.5-T4.5, five have spectral types of T5-T5.5, one is a T6.5p, and two are T7-7.5. We assess spectral morphology and colours to identify T dwarfs in our sample that may have non-typical physical properties (by comparison to solar neighbourhood populations). The colours of the full sample of LAS T dwarfs show a possible trend to bluer Y-J with decreasing effective temperature beyond T8. By accounting for the main sources of incompleteness (selection, follow-up and spatial) as well as the effects of unresolved binarity and Malmquist bias, we estimate that there are 17+-4 >=T4 dwarfs in the J<=19 volume of the LAS second data release. Comparing this to theoretical predictions is most consistent with a sub-stellar mass function exponent alpha between -1.0 and 0. This is consistent with the latest 2MASS/SDSS constraint (which is based on lower number statistics), and is significantly lower than the alpha~1.0 suggested by L dwarf field populations, possibly a result of the lower mass range probed by the T dwarf class.

We present high-resolution interferometric imaging of LH850.02, the brightest 850- and 1200-micron submillimetre (submm) galaxy in the Lockman Hole. Our observations were made at 890 micron with the Submillimetre Array (SMA). Our high-resolution submm imaging detects LH850.02 at >6-sigma as a single compact (size < 1 arcsec or < 8 kpc) point source and yields its absolute position to ~0.2-arcsec accuracy. LH850.02 has two alternative radio counterparts within the SCUBA beam (LH850.02N & S), both of which are statistically very unlikely to be so close to the SCUBA source position by chance. However, the precise astrometry from the SMA shows that the submm emission arises entirely from LH850.02N, and is not associated with LH850.02S (by far the brighter of the two alternative identifications at 24-micron). Fits to the optical-infrared multi-colour photometry of LH850.02N & S indicate that both lie at z~3.3, and are therefore likely to be physically associated. At these redshifts, the 24 micron--to--submm flux density ratios suggest that LH850.02N has an Arp220-type starburst-dominated far-IR SED, while LH850.02S is more similar to Mrk231, with less dust-enshrouded star-formation activity, but a significant contribution at 24-micron (rest-frame ~5-6 micron) from an active nucleus. This complex mix of star-formation and AGN activity in multi-component sources may be common in the high redshift ultraluminous galaxy population, and highlights the need for precise astrometry from high resolution interferometric imaging for a more complete understanding.

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

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

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

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 microns has been carried out to map the regions around several 850 micron-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared luminosities and hence star-formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350 micron flux densities for 24 850 micron sources, seven of which are detected at >2.5-sigma within a 10 arcsec search radius of the 850 micron positions. When results from the literature are included the total number of 350 micron flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-infrared (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich (Mdust~9x10^8 Msun), luminous (Lfir~2x10^12 Lsun), star-forming galaxies with intrinsic dust temperatures of ~35 K and star-formation rates of ~400 Msun/yr. We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350 micron source counts, N350(>25mJy)~200-500 deg^-2.

We describe the WFCAM Science Archive (WSA), which is the primary point of access for users of data from the wide-field infrared camera WFCAM on the United Kingdom Infrared Telescope (UKIRT), especially science catalogue products from the UKIRT Infrared Deep Sky Survey (UKIDSS). We describe the database design with emphasis on those aspects of the system that enable users to fully exploit the survey datasets in a variety of different ways. We give details of the database-driven curation applications that take data from the standard nightly pipeline-processed and calibrated files for the production of science-ready survey datasets. We describe the fundamentals of querying relational databases with a set of astronomy usage examples, and illustrate the results.

We describe the goals, design, implementation, and initial progress of the UKIRT Infrared Deep Sky Survey (UKIDSS), a seven year sky survey which began in May 2005, using the UKIRT Wide Field Camera. It is a portfolio of five survey components covering various combinations of the filter set ZYJHK and H_2. The Large Area Survey, the Galactic Clusters Survey, and the Galactic Plane Survey cover approximately 7000 square degrees to a depth of K~18; the Deep Extragalactic Survey covers 35 square degrees to K~21, and the Ultra Deep Survey covers 0.77 square degrees to K~23. Summed together UKIDSS is 12 times larger in effective volume than the 2MASS survey. The prime aim of UKIDSS is to provide a long term astronomical legacy database; the design is however driven by a series of specific goals -- for example to find the nearest and faintest sub-stellar objects; to discover Population II brown dwarfs, if they exist; to determine the substellar mass function; to break the z=7 quasar barrier; to determine the epoch of re-ionisation; to measure the growth of structure from z=3 to the present day; to determine the epoch of spheroid formation; and to map the Milky Way through the dust, to several kpc. The survey data are being uniformly processed, and released in stages through the WFCAM Science Archive (WSA : http://surveys.roe.ac.uk/wsa). Before the formal survey began, UKIRT and the UKIDSS consortium collaborated in obtaining and analysing a series of small science verification (SV) projects to complete the commissioning of the camera. We show some results from these SV projects in order to demonstrate the likely power of the eventual complete survey. Finally, using the data from the First Data Release we assess how well UKIDSS is meeting its design targets so far.

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

We give a non-exhaustive review of the use of strong gravitational lensing in placing constraints on the quantity of dark and visible mass in galaxies. We discuss development of the methodology and summarise some recent results.

(abridged) We report the discovery of a very cool brown dwarf, ULAS J003402.77-005206.7 (ULAS J0034-00), identified in UKIDSS DR1. We provide optical, near-infrared, and mid-infrared photometry of the source, and two near-infrared spectra. Comparing the spectral energy distribution of ULAS J0034-00 to that of the T8 brown dwarf 2MASS J0415-09, the latest-type and coolest well-studied brown dwarf to date, with Teff~750 K, we find evidence that ULAS J0034-00 is significantly cooler. First, the measured values of the near-infrared absorption spectral indices imply a later classification, of T8.5. Second, the H-[4.49] colour provides an empirical estimate of the temperature of 540<Teff<660 K (+/-2sig range). Third, the J- and H-band peaks are somewhat narrower in ULAS J0034-00, and detailed comparison against spectral models calibrated to 2MASS J0415-09 yields an estimated temperature lower by 60-120 K relative to 2MASS J0415-09 i.e. 630<Teff<690 K (+/-2sig), and lower gravity or higher metallicity according to the degenerate combination -0.5<delta(log g-2[m/H])<-0.25 (+/-2sig). Combining these estimates, and considering systematics, it is likely the temperature lies in the range 600<Teff<700 K. Despite the low inferred Teff we find no evidence for strong absorption by NH3 over the wavelength range 1.51-1.56 um. Evolutionary models imply that the mass and age are in the ranges 15-36 M(Jup) and 0.5-8 Gyr, respectively. The measured proper motion, of (0.37+/-0.07)arcsec/yr, combined with the photometrically estimated distance of 14-22 pc, implies a tangential velocity of ~30 km/s. ULAS J0034-00 is significantly bluer than 2MASS J0415-09 in Y-J, so future searches should allow for the possibility that cooler T dwarfs are bluer still.

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.

The UKIRT Infrared Deep Sky Survey (UKIDSS) is a set of five large near-infrared surveys, covering a complementary range of areas, depths, and Galactic latitudes. The UKIDSS Second Data Release (DR2) includes the First Data Release (DR1), with minor improvements, plus new data for the LAS, GPS, GCS, and DXS, from observations made over 2006 May through July (when the UDS was unobservable). DR2 was staged in two parts. The first part excluded the GPS, and took place on 2007 March 1. The GPS was released on 2007 April 12. DR2 includes 282 sq. degs of multicolour data to (Vega) K=18, complete in the full YJHK set for the LAS, 57 sq. degs in the ZYJHK set for the GCS, and 236 sq. degs in the JHK set for the GPS. DR2 includes nearly 7 sq. degs of deep JK data (DXS, UDS) to an average depth K=21. In addition the release includes a comparable quantity of data where coverage of the filter set for any survey is incomplete. We document changes that have occurred since DR1 to the pipeline, calibration, and archive procedures. The two most noteworthy changes are presentation of the data in a single database (compared to two previously), and provision of additional error flags for detected sources, flagging potentially spurious artifacts, corrupted data and suspected cross-talk sources. We summarise the contents of each of the surveys in terms of filters, areas, and depths.

We analyse the first publicly released deep field of the UKIDSS Deep eXtragalactic Survey (DXS) to identify candidate galaxy over-densities at z~1 across ~1 sq. degree in the ELAIS-N1 field. Using I-K, J-K and K-3.6um colours we identify and spectroscopically follow-up five candidate structures with Gemini/GMOS and confirm they are all true over-densities with between five and nineteen members each. Surprisingly, all five structures lie in a narrow redshift range at z=0.89+/-0.01, although they are spread across 30Mpc on the sky. We also find a more distant over-density at z=1.09 in one of the spectroscopic survey regions. These five over-dense regions lying in a narrow redshift range indicate the presence of a supercluster in this field and by comparing with mock cluster catalogs from N-body simulations we discuss the likely properties of this structure. Overall, we show that the properties of this supercluster are similar to the well-studied Shapley and Hercules superclusters at lower redshift.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

We model the mass distribution in the recently discovered Einstein ring LBG J213512.73-010143 (the `Cosmic Eye') using archival Hubble Space Telescope imaging. We reconstruct the mass density profile of the z=0.73 lens and the surface brightness distribution of the z=3.07 source and find that the observed ring is best fit with a dual-component lens model consisting of a baryonic Sersic component nested within a dark matter halo. The dark matter halo has an inner slope of 1.42+/-0.23, consistent with CDM simulations after allowing for baryon contraction. The baryonic component has a B-band mass-to-light ratio of 1.71+0.28-0.38 (solar units) which when evolved to the present day is in agreement with local ellipticals. Within the Einstein radius of 0.77'' (5.6 kpc), the baryons account for (46+/-11)% of the total lens mass. External shear from a nearby foreground cluster is accurately predicted by the model. The reconstructed surface brightness distribution in the source plane clearly shows two peaks. Through a generalisation of our lens inversion method, we conclude that the redshifts of both peaks are consistent with each other, suggesting that we are seeing structure within a single galaxy.

We report the detection of CO(3-2) emission from a bright, gravitationally lensed Lyman Break Galaxy, LBGJ213512.73-010143 (the 'Cosmic Eye'), at z=3.07 using the Plateau de Bure Interferometer. This is only the second detection of molecular gas emission from an LBG and yields an intrinsic molecular gas mass of (2.4+/-0.4)x10^9 Mo. The lens reconstruction of the UV morphology of the LBG indicates that it comprises two components separated by ~2 kpc. The CO emission is unresolved, and appears to be centered on the intrinsically fainter (and also less highly magnified) of the two UV components. The width of the CO line indicates a dynamical mass of (8+/-2)x10^9csc(i)^2 Mo within the central 2 kpc. Employing mid-infrared observations from Spitzer we derive a stellar mass of ~(6+/-2)x10^9 Mo and a star-formation rate of ~60 Mo/yr, indicating that the molecular gas will be consumed in ~40 Myr. The gas fractions, star-formation efficiencies and line widths suggests that LBGJ213512 is a high-redshift, gas-rich analog of a local luminous infrared galaxy. This galaxy has a similar gas-to-dynamical mass fraction as observed in the submillimeter-selected population, although the gas surface density and star-formation efficiency is a factor of 3x less, suggesting less vigorous activity. We discuss the uncertainties in our conclusions arising from adopting a CO-to-H2 conversion factor appropriate for either the Milky Way or local luminous infrared galaxies. These observations demonstrate that current facilities, when aided by fortuitous gravitational magnification, can study 'ordinary' galaxies at high-redshift and so act as pathfinders for ALMA.

We report on the first ultracool dwarf discoveries from the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey Early Data Release (LAS EDR), in particular the discovery of T dwarfs which are fainter and more distant than those found using the 2MASS and SDSS surveys. We aim to show that our methodologies for searching the ~27 sq degs of the LAS EDR are successful for finding both L and T dwarfs $via$ cross-correlation with the Sloan Digital Sky Survey (SDSS) DR4 release. While the area searched so far is small, the numbers of objects found shows great promise for near-future releases of the LAS and great potential for finding large numbers of such dwarfs. Ultracool dwarfs are selected by combinations of their YJH(K) UKIDSS colours and SDSS DR4 z-J and i-z colours, or, lower limits on these red optical/infrared colours in the case of DR4 dropouts. After passing visual inspection tests, candidates have been followed up by methane imaging and spectroscopy at 4m and 8m-class facilities. Our main result is the discovery following CH4 imaging and spectroscopy of a T4.5 dwarf, ULASJ 1452+0655, lying ~80pc distant. A further T dwarf candidate, ULASJ 1301+0023, has very similar CH4 colours but has not yet been confirmed spectroscopically. We also report on the identification of a brighter L0 dwarf, and on the selection of a list of LAS objects designed to probe for T-like dwarfs to the survey J-band limit. Our findings indicate that the combination of the UKIDSS LAS and SDSS surveys provide an excellent tool for identifying L and T dwarfs down to much fainter limits than previously possible. Our discovery of one confirmed and one probable T dwarf in the EDR is consistent with expectations from the previously measured T dwarf density on the sky.

We present the initial results from our search for high redshift, z > 6, quasars using near infrared data from the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey (LAS). Our analysis of 106 deg^2 of sky from Data Release 1 (DR1) has resulted in the discovery of ULAS J020332.38+001229.2, a luminous (J_AB=20.0, J_Vega=19.1, M_1450=-26.2) quasar at z=5.86. Following candidate selection from the combined IR and optical catalogue data and stacking of multiple epoch Sloan Digital Sky Survey (SDSS) data, we have obtained optical spectroscopy for the only two z > 6 quasar candidates. The VLT FORS2 spectrum of ULAS J020332.38+001229.2 shows broad Ly alpha + NV1240 emission at ~8350 A and an abrupt continuum break due to absorption by the Ly alpha forest. The quasar is not present in the SDSS DR5 catalogue and the continuum spectral index of alpha=-1.4 (F_nu nu^alpha) is redder than a composite of SDSS quasars at similar redshifts (alpha=-0.5). The discovery of one z~6 quasar in ~100 deg^2 in a complete sample within our selection criteria down to a median depth of Y_AB=20.4 (7 sigma) is consistent with existing SDSS results. We describe our survey methodology, including the use of optical data from the SDSS and the highly effective procedures developed to isolate the very small surface density of high-probability quasar candidates.

We have measured the contribution of submillimeter and mid-infrared sources to the extragalactic background radiation at 70 and 160um. Specifically, we have stacked flux in 70 and 160um Spitzer Space Telescope (Spitzer) observations of the Canada-UK Deep Sub-millimeter Survey 14h field at the positions of 850um sources detected by SCUBA and also 8 and 24um sources detected by Spitzer. We find that per source, the SCUBA galaxies are the strongest and the 8um sources the weakest contributors to the background flux at both 70 and 160um. Our estimate of the contribution of the SCUBA sources is higher than previous estimates. However, expressed as a total contribution, the full 8um source catalogue accounts for twice the total 24um source contribution and ~10 times the total SCUBA source contribution. The 8um sources account for the majority of the background radiation at 160um with a flux of 0.87+/-0.16 MJy/sr and at least a third at 70um with a flux of 0.103+/-0.019 MJy/sr. These measurements are consistent with current lower limits on the background at 70 and 160um. Finally, we have investigated the 70 and 160um emission from the 8 and 24um sources as a function of redshift. We find that the average 70um flux per 24um source and the average 160um flux per 8 and 24um source is constant over all redshifts, up to z~4. In contrast, the low-redshift half (z<1) of the of 8um sample contributes approximately four times the total 70um flux of the high-redshift half. These trends can be explained by a single non-evolving SED.

We report the discovery using Hubble Space Telescope imaging and Keck spectroscopy of a very bright, highly magnified (~30*) Lyman Break Galaxy (LBG) at z=3.07 in the field of the massive z=0.33 cluster MACSJ2135.2-0102. The system comprises two high-surface brightness arcs with a maximum extent of 3", bracketing a central object which we identify as a massive early-type galaxy at z=0.73. We construct a lens model which reproduces the main features of the system using a combination of a galaxy-scale lens and the foreground cluster. We show that the morphological, spectral and photometric properties of the arcs are consistent with them arising from the lensing of a single ~L* LBG. The most important feature of this system is that the lensing magnification results in an apparent magnitude of r=20.3, making this one of the brightest LBGs known. Such a high magnification provides the opportunity of obtaining very high signal to noise (and potentially spatially resolved) spectroscopy of a high redshift galaxy to study its physical properties. We present initial imaging and spectroscopy demonstrating the basic properties of the system and discuss the opportunities for future observations.

This paper defines the UKIRT Infrared Deep Sky Survey (UKIDSS) Early Data Release (EDR). UKIDSS is a set of five large near-infra-red surveys defined by Lawrence et al. (2006), being undertaken with the UK Infra-red Telescope (UKIRT) Wide Field Camera (WFCAM). The programme began in May 2005 and has an expected duration of seven years. Each survey uses some or all of the broadband filter complement ZYJHK. The EDR is the first public release of data to the European Southern Observatory (ESO) community. All worldwide releases occur after a delay of 18 months from the ESO release. The EDR provides a small sample dataset, ~50 sq.deg (about 1% of the whole of UKIDSS), that is a lower limit to the expected quality of future survey data releases. In addition, an EDR+ dataset contains all EDR data plus extra data of similar quality, but for areas not observed in all of the required filters (amounting to ~220 sq.deg). The first large data release, DR1, will occur in mid-2006. We provide details of the observational implementation, the data reduction, the astrometric and photometric calibration, and the quality control procedures. We summarise the data coverage and quality (seeing, ellipticity, photometricity, depth) for each survey and give a brief guide to accessing the images and catalogues from the WFCAM Science Archive.

We identify 17 possible 8.0 micron-selected counterparts to the submillimeter galaxies in the CUDSS 14-hour field, derived from deep imaging carried out with the IRAC and MIPS instruments aboard the Spitzer Space Telescope. Ten of the 17 counterparts are not the same as those previously identified at shorter wavelengths. We argue that 8.0 micron selection offers a better means for identifying counterparts to submillimeter galaxies than near-infrared or optical selection. Based on the panchromatic SEDs, most counterparts appear to be powered by ongoing star formation. Power-law fits to the SEDs suggest that five objects in the 8.0 micron-selected counterpart sample harbor dominant AGNs; a sixth object is identified as a possible AGN. The 3.6 to 8.0 micron colors of the infrared-selected counterparts are significantly redder than the general IRAC galaxy population in the CUDSS 14-hour field.

We present an analysis of stellar mass estimates for a sample of 25000 galaxies from the COMBO-17 survey over the interval 0.2<z<1.0. We have developed, implemented, and tested a new method of estimating stellar mass-to-light ratios, which relies on redshift and spectral energy distribution (SED) classification from 5 broadband and 12 medium band filters. We find that the majority (>60%) of massive galaxies with M_* > 10^{11} solar masses at all z<1 are non-star-forming; blue star-forming galaxies dominate at lower masses. We have used these mass estimates to explore the evolution of the stellar mass function since z=1. We find that the total stellar mass density of the universe has roughly doubled since z~1. Our measurements are consistent with other measurements of the growth of stellar mass with cosmic time and with estimates of the time evolution of the cosmic star formation rate. Intriguingly, the integrated stellar mass of blue galaxies with young stars has not significantly changed since z~1, even though these galaxies host the majority of the star formation: instead, the growth of the total stellar mass density is dominated by the growth of the total mass in the largely passive galaxies on the red sequence.

Bright galaxy-galaxy strong lenses are much more powerful than lensed quasars for measuring the mass profiles of galaxies, but until this year only a handful have been known. Here we present five new examples, identified via the optimal line-of-sight gravitational lens search strategy applied to luminous red galaxies in the Sloan Digital Sky Survey (SDSS). Our survey largely complements a similar survey by Bolton et al., who recently presented several new lenses. The lensed background galaxies are selected from the SDSS spectra via the presence of narrow emission line signatures, including the [OII] 3726,3729, Hb and [OIII] 4960,5008 lines, superposed on the spectra of the bright, intervening, deflector galaxies. Our five confirmed new systems include deflector galaxies with redshifts z=0.17-0.28 and lensed galaxies with redshifts z=0.47-1.18. Simulations of moderately deep (few orbits) HST-ACS imaging of systems such as these, where the lensed source is brighter than r~23, are presented. These demonstrate the feasibility of measuring accurately the inner slope of the dark matter halo to within an uncertainty sigma(gamma)~0.1, the dark matter fraction within the Einstein radius, and the mass-to-light ratio of the stars alone, independently of dynamical measurements. The high success rate of our search so far, >60%, and the relatively modest observational resources necessary to confirm the gravitational lens nature of the candidates, demonstrate that compilation of a sample of ~100 galaxy-galaxy lenses from the SDSS is readily achievable, opening up a rich new field in dark matter studies.

(Abridged) We present our VLT/FORS1 deep spectroscopic observations of the gravitationally lensed quasar SDSS J0924+0219, as well as archival HST/NICMOS and ACS images of the same object. The two-epoch spectra, obtained in the Multi Object Spectroscopy (MOS) mode, allow for very accurate flux calibration, spatial deconvolution of the data, and provide the redshift of the lensing galaxy z=0.394 +/- 0.001. These spectra, taken 15 days apart, show only slight continuum variations, while the broad emission lines display obvious changes in the red wing of the Mg II line, in the Fe II bands, and in the central part of the C III] line. Even though variations in the line profiles are present, we do not see any significant differences between the continuum and emission line flux ratios of images A and B of the quasar. Spatial deconvolution of the HST images reveals a double Einstein ring. One ring corresponds to the lensed quasar host galaxy at z=1.524 and a second bluer one, is the image either of a star-forming region in the host galaxy, or of another unrelated lower redshift object. We find that a broad range of lens models gives a satisfactory fit to the data. However, they predict very different time delays, making SDSS J0924+0219 an object of particular interest for photometric monitoring. In addition, the lens models reconstructed using exclusively the constraints from the Einstein rings, or using exclusively the astrometry of the quasar images, are not compatible. This suggests that substructures play an important role in SDSS J0924+0219.

Aims: To predict time delays for a sample of gravitationally lensed quasars and to evaluate the accuracy that can be realistically achieved on the value of H_0. Methods: We consider 14 lensed quasars that are candidates for time-delay monitoring and model them in detail using pixelized lens models. For each system, we provide a mass map, arrival-time surface and the distribution of predicted time-delays in a concordance cosmology, assuming H_0^{-1}=14 Gyr (H_0=70 in local units). Based on the predicted time-delays and on the observational circumstances, we rate each lens as `excellent' or `good' or `unpromising' for time-delay monitoring. Finally, we analyze simulated time delays for the 11 lens rated excellent or good, and show that H_0 can be recovered to a precision of 5%. Results: In combination with COSMOGRAIL paper I on the temporal sampling of lensed quasar light curves, the present work will help design monitoring campaigns of lensed quasars.

We investigate the redshift dependence of the contribution to the extragalactic far-infrared/sub-millimeter background from galaxies detected by the Spitzer Space Telescope at 8um and 24um. Using seven-band optical to mid-infrared photometry, we estimate photometric redshifts for the Spitzer sources which appear to be mostly L* galaxies at a median redshift of z=1.0. These sources, extracted from deep 8um and 24um mosaics of the CUDSS 14-hour field with 5sigma limits of 5.8uJy and 70uJy respectively, exhibit significant 850um and 450um emission as observed by SCUBA. At 850um, after removing >=4sigma sources and those securely identified in our previous cross-matching paper, we measure stacked flux at the significance level of 4.4sigma and 2.9sigma from the full 8um and 24um galaxy catalogue respectively. At 450um, flux is detected from all 8um galaxies at the level of 3.5sigma, while there is no significant emission from the 24um galaxies. We find that the 850um flux is emitted almost exclusively at z>~1.3 from the Spitzer sources with 0.44mJy (4.7sigma) per 8um source and 0.51mJy (2.8sigma) per 24um source. This corresponds to a contribution of (16+/-3)% toward the 850um extra-galactic background from the 8um sources and (5.0+/-1.8)% from the 24um sources. At 450um, only the 8um sources within the redshift interval 1<z<2 exhibit significant emission with an average flux per source of 3.35mJy (3.0sigma). This is a contribution of (37+/-12)% to the 450um background.

We measure the mass density profile of the lens galaxy in the Einstein ring system 0047-2808 using our semi-linear inversion method developed in an earlier paper. By introducing an adaptively gridded source plane, we are able to eliminate the need for regularisation of the inversion. This removes the problem of a poorly defined number of degrees of freedom, encountered by inversion methods that employ regularisation, and so allows a proper statistical comparison between models. We confirm the results of Wayth et al. (2004), that the source is double, and that a power-law model gives a significantly better fit that the singular isothermal ellipsoid model. We measure a slope alpha=2.11+/-0.04. We find, further, that a dual-component constant M/L baryonic + dark halo model gives a significantly better fit than the power-law model, at the 99.7% confidence level. The inner logarithmic slope of the dark halo profile is found to be 0.87^{+0.69}_{-0.61} (95% CL), consistent with the predictions of CDM simulations of structure formation. We determine an unevolved B-band mass to light ratio for the baryons (only) of 3.05^{+0.53}_{-0.90} h_65 M_sol/L_Bsol (95% CL). This is the first measurement of the baryonic M/L of a single galaxy by purely gravitational lens methods. The baryons account for 65^{+10}_{-18}% (95% CL) of the total projected mass, or, assuming spherical symmetry, 84^{+12}_{-24}% (95% CL) of the total three-dimensional mass within the mean radius of 1.16'' (7.5 h_65^{-1} kpc) traced by the ring. Finally, at the level of >3sigma, we find that the halo mass is rounder than the baryonic distribution and that the two components are offset in orientation from one another.

We use numerical simulations to test a broad range of plausible observational strategies designed to measure the time delay between the images of gravitationally lensed quasars. Artificial quasar light curves are created along with Monte-Carlo simulations in order to determine the best temporal sampling to adopt when monitoring the photometric variations of systems with time delays between 5 and 120 days, i.e., always shorter than the visibility window across the year. Few and realistic assumptions are necessary on the quasar photometric variations (peak-to-peak amplitude and time-scale of the variations) and on the accuracy of the individual photometric points. The output of the simulations is the (statistical) relative error made on the time delay measurement, as a function of 1- the object visibility over the year, 2- the temporal sampling of the light curves and 3- the time delay. Also investigated is the effect of long term microlensing variations which must be below the 5 % level (either intrinsically or by subtraction) if the goal is to measure time delays with an accuracy of 1-2 %. However, while microlensing increases the random error on the time delay, it does not significantly increase the systematic error, which is always a factor 5 to 10 smaller than the random error. Finally, it is shown that, when the time delay is comparable to the visibility window of the object, a logarithmic sampling can significantly improve the time delay determination. All results are presented in the form of compact plots to be used to optimize the observational strategy of future monitoring programs.

We investigate how galaxy-galaxy lensing measurements depend on the knowledge of redshifts for lens and source galaxies. Galaxy-galaxy lensing allows one to study dark matter halos of galaxies statistically using weak gravitational lensing. Redshift information is required to reliably distinguish foreground lens galaxies from background source galaxies and to convert the measured shear into constraints on the lens model. Without spectroscopy or multi-colour information, redshifts can be drawn from independently estimated probability distributions. The COMBO-17 survey provides redshifts for both lens and source galaxies. It thus offers the unique possibility to do this investigation with observational data. We find that it is of great importance to know the redshifts of individual lens galaxies in order to constrain the properties of their dark matter halos. Whether the redshifts are derived from $UBVRI$ or the larger number of filters available in COMBO-17 is not very important. In contrast, knowledge of individual source redshifts improves the measurements only very little over the use of statistical source redshift distributions as long as the source redshift distribution is known accurately.

We present a measurement of mass estimates for dark matter halos around galaxies from the COMBO-17 survey using weak gravitational lensing. COMBO-17 is particularly useful for this kind of investigation because it covers observations in 17 optical filters from which accurate photometric redshifts and spectral classification for objects with $R<24$ are derived. This allows us to select lens and source galaxies from their redshifts and to thus avoid any uncertainties from estimates of the source redshift distribution. We study galaxy lenses at redshifts $z_\mathrm{d}=0.2-0.7$ by fitting the normalization of either singular isothermal spheres (SIS) or Navarro-Frenk-White (NFW) profiles to the whole lens sample; we then consider halos around blue and red subsamples separately. We also constrain the scaling of halo mass with light. For the NFW model, we find virial masses $M_\mathrm{vir}^*=3.9^{+3.3}_{-2.4}\times 10^{11}h^{-1}M_{\sun}$ (1-$\sigma$) for blue and $M_\mathrm{vir}^*=7.1^{+7.1}_{-3.8}\times 10^{11}h^{-1}M_{\sun}$ for red galaxies of $L_\star=10^{10}h^{-2}L_{\sun}$, respectively. The derived mass-to-light scaling relations suggest that the mass-to-light ratio might decrease with increasing luminosity for blue galaxies but increase with increasing luminosity for red galaxies. However, these differences between blue and red galaxies are only marginally significant and both subsamples are consistent with having the same mass-to-light ratio at all luminosities. Finally, we compare our results to those obtained from the Red-Sequence Cluster Survey (RCS) and the Sloan Digital Sky Survey (SDSS). Taking differences in the actual modelling into account, we find very good agreement with these surveys.

We present the COMBO-17 object catalogue of the Chandra Deep Field South for public use, covering a field which is 31.5' x 30' in size. This catalogue lists astrometry, photometry in 17 passbands from 350 to 930 nm, and ground-based morphological data for 63,501 objects. The catalogue also contains multi-colour classification into the categories 'Star', 'Galaxy' and 'Quasar' as well as photometric redshifts. We include restframe luminosities in Johnson, SDSS and Bessell passbands and estimated errors. The redshifts are most reliable at R<24, where the sample contains approximately 100 quasars, 1000 stars and 10000 galaxies. We use nearly 1000 spectroscopically identified objects in conjunction with detailed simulations to characterize the performance of COMBO-17. We show that the selection of quasars, more generally type-1 AGN, is nearly complete and minimally contaminated at z=[0.5,5] for luminosities above M_B=-21.7. Their photometric redshifts are accurate to roughly 5000 km/sec. Galaxy redshifts are accurate to 1% in dz/(1+z) at R<21. They degrade in quality for progressively fainter galaxies, reaching accuracies of 2% for galaxies with R~222 and of 10% for galaxies with R>24. The selection of stars is complete to R~23, and deeper for M stars. We also present an updated discussion of our classification technique with maps of survey completeness, and discuss possible failures of the statistical classification in the faint regime at R>24.

We present a direct detection of the growth of large-scale structure, using weak gravitational lensing and photometric redshift data from the COMBO-17 survey. We use deep R-band imaging of two 0.25 square degree fields, affording shear estimates for over 52000 galaxies; we combine these with photometric redshift estimates from our 17 band survey, in order to obtain a 3-D shear field. We find theoretical models for evolving matter power spectra and correlation functions, and fit the corresponding shear correlation functions to the data as a function of redshift. We detect the evolution of the power at the 7.7 sigma level given minimal priors, and measure the rate of evolution for 0<z<1. We also fit correlation functions to our 3-D data as a function of cosmological parameters sigma_8 and Omega_Lambda. We find joint constraints on Omega_Lambda and sigma_8, demonstrating an improvement in accuracy by a factor of 2 over that available from 2D weak lensing for the same area.

We present the rest-frame colors and luminosities of ~25000 m_R<24 galaxies in the redshift range 0.2<z<1.1, drawn from 0.78 square degrees of the COMBO-17 survey. We find that the rest-frame color distribution of these galaxies is bimodal at all redshifts out to z~1. This bimodality permits a model-independent definition of red, early-type galaxies and blue, late-type galaxies at any given redshift. The colors of the blue peak become redder towards the present day, and the number density of blue luminous galaxies has dropped strongly since z~1. Focusing on the red galaxies, we find that they populate a color-magnitude relation. Such red sequences have been identified in galaxy cluster environments, but our data show that such a sequence exists over this redshift range even when averaging over all environments. The mean color of the red galaxy sequence evolves with redshift in a way that is consistent with the aging of an ancient stellar population. The rest-frame B-band luminosity density in red galaxies evolves only mildly with redshift in a Lambda-dominated cold dark matter universe. Accounting for the change in stellar mass-to-light ratio implied by the redshift evolution in red galaxy colors, the COMBO-17 data indicate an increase in stellar mass on the red sequence by a factor of two since z~1. The largest source of uncertainty is large-scale structure, implying that considerably larger surveys are necessary to further refine this result. We explore mechanisms that may drive this evolution in the red galaxy population, finding that both galaxy merging and truncation of star formation in some fraction of the blue, star-forming population are required to fully explain the properties of these galaxies.

We present a 3-dimensional lensing analysis of the z=0.16 supercluster A901/2, resulting in a 3-D map of the dark matter distribution within a 3 X 10^{5} [Mpc]^3 volume from the COMBO-17 survey. We perform a chi^2-fit of isothermal spheres to the tangential shear pattern around each cluster as a function of redshift to estimate the 3-D positions and masses of the main clusters in the supercluster from lensing alone. We then present the first 3-D map of the dark matter gravitational potential field, Phi, using the Kaiser-Squires (1993) and Taylor (2001) inversion methods. These maps clearly show the potential wells of the main supercluster components, including a new cluster behind A902, and demonstrates the applicability of 3-D dark matter mapping and projection free-mass-selected cluster finding to current data. Finally, we develop the halo model of dark matter and galaxy clustering and compare this with the auto-and cross-correlation functions of the 3-D gravitational potential, galaxy number densities and galaxy luminosity densities measured in the A901/2 field. We find significant anti-correlations between the gravitational potential field and the galaxy number density and luminosities, as expected due to baryonic infall into dark matter concentrations. We find good agreement with the halo model for the number densities and luminosity correlation functions.

We investigate correlations between the location of galaxies in dense environments and their degree of star-formation activity. Using photometric redshifts and spectral classifications from the unique 17-band COMBO-17 survey we are able to precisely isolate galaxies from the Abell 901/902 supercluster within a thin redshift slice around z=0.16. We compare the detailed photometric properties of the supercluster galaxies with the underlying dark matter density field as revealed by weak gravitational lensing. We find strong evidence for segregation by type, with the highest density regions populated almost exclusively by galaxies classified according to their rest-frame U-V colours as quiescent. We also observe a threshold surface mass density from lensing, \kappa\sim 0.05 (corresponding to a physical density \Sigma = 2.5x10^(14)h M_\sun Mpc^(-2)), above which star-formation activity is rapidly suppressed. This abrupt transformation affects primarily the faint end of the star-forming galaxy population and occurs at a local surface number density corresponding to roughly 400h^2 Mpc^(-2) to a limit of M*_V+6. When only galaxies brighter than M*+1 are considered the trends with environment remain, but are more gradual and extend beyond 2h^(-1) Mpc radius.

Using our semi-linear inversion method, we measure the mass profile of the lens galaxy in the Einstein ring system 0047-2808. The lens is modelled as a baryonic component following the observed light, embedded in a dark matter generalised NFW halo. The semi-linear method makes full use of the information content in the ring image. We determine an unevolved B-band mass to light ratio for the baryons of 3.05 +0.45 -0.86 h_65 M_sol/L_sol (95% CL), accounting for 65% of the total projected mass within the radius of 1.16'' traced by the ring. This result is obtained without need of dynamical measurements. The inner logarithmic slope of the halo is found to be 0.87 +0.69 -0.61 (95% CL). We find that the halo is fairly well aligned with the light but has only half the ellipticity.

We present a determination of the optical/UV AGN luminosity function and its evolution, based on a large sample of faint (R<24) QSOs identified in the COMBO-17 survey. Using multi-band photometry in 17 filters within 350nm < lambda_obs < 930nm, we could simultaneously determine photometric redshifts with an accuracy of sigma_z<0.03 and obtain spectral energy distributions. The redshift range covered by the sample is 1.2<z<4.8, which implies that even at z~3, the sample reaches below luminosities corresponding to M_B = -23, conventionally employed to distinguish between Seyfert galaxies and quasars. We clearly detect a broad plateau-like maximum of quasar activity around z~2 and map out the smooth turnover between z~1 and z~4. The shape of the LF is characterised by some mild curvature, but no sharp `break' is present within the range of luminosities covered. Using only the COMBO-17 data, the evolving LF can be adequately described by either a pure density evolution (PDE) or a pure luminosity evolution (PLE) model. However, the absence of a strong L*-like feature in the shape of the LF inhibits a robust distinction between these modes. We present a robust estimate for the integrated UV luminosity generation by AGN as a function of redshift. We find that the LF continues to rise even at the lowest luminosities probed by our survey, but that the slope is sufficiently shallow that the contribution of low-luminosity AGN to the UV luminosity density is negligible. Although our sample reaches much fainter flux levels than previous data sets, our results on space densities and LF slopes are completely consistent with extrapolations from recent major surveys such as SDSS and 2QZ.

We perform a cosmic shear analysis of the COMBO-17 survey -- a unique dataset with shear quality R-band imaging and accurate photometric redshift estimates (dz=0.05) for ~90% of galaxies to R=24.0. We undertake a full maximum likelihood analysis to measure the weak lensing power spectra, Cl^kk, Cl^bb & Cl^kb from l=400 to l=10^4. We find a strong measurement of the convergence power over five fields. The b-field has a much lower significance, indicating our data is free of major systematics, while the cross-correlation of k & b is consistent with zero. We have also calculated the shear correlation functions and variance over a range of scales between 0.5 and 20 arcmin. In addition, we have used our results to measure cosmological parameters, constraining the normalisation of the matter power spectrum to be sigma_8=(0.72 \pm 0.09)(Omega_m/0.3)^-0.49, where the errors quoted are 1-sigma due to the intrinsic dispersion in galaxy ellipticities, cosmic and sampling variance. We have significantly reduced the usual additional uncertainty in the median redshift (z_m) of the survey by estimating z_m directly from the data. To demonstrate the power of accurate redshift information, we have also measured parameters from a shear analysis of only those galaxies with accurate redshifts. In this case, we have eliminated the uncertainty in the redshift distribution of sources and we show that the uncertainty in the resulting parameter constraints are reduced by more than a factor of 2 compared to the typical uncertainties found in cosmic shear surveys to date. Finally, we combine our parameter measurements with constraints from the 2dFGRS and from the CMB. With these additional constraints, we measure sigma_8=0.73 +0.05/-0.03 and Omega_m=0.27 +/- 0.02.

We describe a new method for analyzing gravitational lens images, for the case where the source light distribution is pixelized. The method is suitable for high resolution, high S/N data of a multiply-imaged extended source. For a given mass distribution, we show that the step of inverting the image to obtain the deconvolved pixelized source light distribution, and the uncertainties, is a linear one. This means that the only parameters of the non-linear problem are those required to model the mass distribution. This greatly simplifies the search for a min.-chi^2 fit to the data and speeds up the inversion. The method is extended in a straightforward way to include linear regularization. We apply the method to simulated Einstein ring images and demonstrate the effectiveness of the inversion for both the unregularized and regularized cases.

We present a detailed empirical assessment of how the galaxy luminosity function and stellar luminosity density evolves over the last half of the universe's age (0.2<z<1.2) for galaxies of different spectral energy distributions (SED). The results are based on ~25,000 galaxies (R<24) with redshift measurements (sigma_z~0.03) and SEDs across 350..930 nm, derived from medium-band photometry in 17 filters, observed as part of the COMBO-17 survey (``Classifying Objects by Medium-Band Observations in 17 Filters'') over three disjoint fields with a total area of 0.78 square degrees. Luminosity functions (LF), binned in redshift and SED-type, are presented in the restframe passbands of the SDSS r-band, the Johnson B-band and a synthetic UV continuum band at 280 nm. We find that the luminosity function depends strongly on SED-type at all redshifts covered. The shape of the LF, i.e. the faint-end power-law slope, does depend on SED type, but not on redshift. However, the redshift evolution of the characteristic luminosity M* and density phi* depends strongly on SED-type: (1) Early-type galaxies, defined as redder than a present-day reference Sa spectrum, become drastically more abundant towards low redshift, by a factor of 10 in the number density phi* from z=1.1 to now, and by a factor of 4 in their contribution to the co-moving r-band luminosity density, j_r. (2) Galaxies resembling present-day Sa- to Sbc-colours show a co-moving number density and contribution to j_r that does not vary much with redshift. (3) Galaxies with blue spectra reflecting strong star formation decrease towards low redshift both in luminosity and density, and by a factor of 4 in their j_r contribution. (abridged)

[ABRIDGED] We estimate the total mass distribution of the galaxy cluster CL0024+1654 from the measured source depletion due to lens magnification in the R band. Within a radius of 0.54Mpc/h, a total projected mass of (8.1+/-3.2)*10^14 M_sol/h (EdS) is measured, which corresponds to a mass- to-light ratio of M/L(B)=470+/-180. We compute the luminosity function of CL0024+1654 in order to estimate contamination of the background source counts from cluster galaxies. Three different magnification-based reconstruction methods are employed using both local and non-local techniques. We have modified the standard single power-law slope number count theory to incorporate a break and applied this to our observations. Fitting analytical magnification profiles of different cluster models to the observed number counts, we find that the cluster is best described either by a NFW model with scale radius r_s=334+/-191 kpc/h and normalisation kappa_s=0.23+/-0.08 or a power-law profile with slope xi=0.61+/-0.11, central surface mass density kappa_0=1.52+/-0.20 and assuming a core radius of r_core=35 kpc/h. The NFW model predicts that the cumulative projected mass contained within a radius R scales as M(<R)=2.9*10^14*(R/1')^[1.3-0.5lg (R/1')] M_sol/h. Finally, we have exploited the fact that flux magnification effectively enables us to probe deeper than the physical limiting magnitude of our observations in searching for a change of slope in the U band number counts. We rule out both a total flattening of the counts with a break up to U_AB<=26.6 and a change of slope, reported by some studies, from dlog N/dm=0.4->0.15 up to U_AB<=26.4 with 95% confidence.

We measure the alignment of galaxy ellipticities in the local universe over a range of scales using digitized photographic data from the SuperCOSMOS Sky Survey. We find for a magnitude cut of b_J < 20.5, corresponding to a median galaxy redshift of z = 0.1, and 2x10^6 galaxies, that the galaxy ellipticities exhibit a non-zero correlation over a range of scales between 1 and 100 arcminutes. In particular, we measure the variance of mean galaxy ellipticities, sg^2(theta), in square angular cells on the sky as a function of cell size and find it lies in the range, 2 x 10^{-4} > sg^2(theta) > 1 x 10^{-5} for cell side lengths between 15 < theta < 100 arcminutes. Considering the low median redshift of the galaxies in the sample and hence the relatively low effective cross-section for lensing of these galaxies by the large-scale structure of the Universe, we propose that we have detected an intrinsic alignment of galaxy ellipticities. We compare our results to recent analytical and numerical predictions made for the intrinsic galaxy alignment and find good agreement. We discuss the importance of these results for measuring cosmic shear from upcoming shallow surveys (e.g. Sloan Digital Sky Survey) and we outline how these measurements could possibly be used to constrain models of galaxy formation and/or measure the mass distribution in the local universe.

We present a weak shear analysis of the Abell 901/902 supercluster, composed of three rich clusters at z=0.16. Using a deep R-band image from the 0.5 x 0.5 degree MPG/ESO Wide Field Imager together with supplementary B-band observations, we build up a comprehensive picture of the light and mass distributions in this region. We find that, on average, the light from the early-type galaxies traces the dark matter fairly well, although one cluster is a notable exception to this rule. The clusters themselves exhibit a range of mass-to-light (M/L) ratios, X-ray properties, and galaxy populations. We attempt to model the relation between the total mass and the light from the early-type galaxies with a simple scale-independent linear biasing model. We find M/L_B=130h for the early type galaxies with zero stochasticity, which, if taken at face value, would imply Omega_m < 0.1. However, this linear relation breaks down on small scales and on scales equivalent to the average cluster separation (approximately 1 Mpc), demonstrating that a single M/L ratio is not adequate to fully describe the mass-light relation in the supercluster. Rather, the scatter in M/L ratios observed for the clusters supports a model incorporating non-linear biasing or stochastic processes. Finally, there is a clear detection of filamentary structure connecting two of the clusters, seen in both the galaxy and dark matter distributions, and we discuss the effects of cluster-cluster and cluster-filament interactions as a means to reconcile the disparate descriptions of the supercluster.

We report on deep multi-color imaging (R_5-sigma = 26) of the Chandra Deep Field South, obtained with the Wide Field Imager (WFI) at the MPG/ESO 2.2 m telescope on La Silla as part of the COMBO-17 survey. As a result we present a catalog of 63501 objects in a field measuring 31.5' x 30' with astrometry and BVR photometry. A sample of 37 variable objects is selected from two-epoch photometry. We try to give interpretations based on color and variation amplitude.

Gravitational lensing magnifies the luminosity of galaxies behind the lens. We use this effect to constrain the total mass in the cluster Abell 1689 by comparing the lensed luminosities of background galaxies with the luminosity function of an undistorted field. Since galaxies are assumed to be a random sampling of luminosity space, this method is not limited by clustering noise. We use photometric redshift information to estimate galaxy distance and intrinsic luminosity. Knowing the redshift distribution of the background population allows us to lift the mass/background degeneracy common to lensing analysis. In this paper we use 9 filters observed over 12 hours with the Calar Alto 3.5m telescope to determine the redshifts of 1000 galaxies in the field of Abell 1689. Using a complete sample of 151 background galaxies we measure the cluster mass profile. We find that the total projected mass interior to 0.25h^(-1)Mpc is (0.48 +/- 0.16) * 10^(15)h^(-1) solar masses, where our error budget includes uncertainties from the photometric redshift determination, the uncertainty in the off-set calibration and finite sampling. This result is in good agreement with that found by number count and shear-based methods and provides a new and independent method to determine cluster masses.

We present a new method for directly determining accurate, self-consistent cluster lens mass and shear maps in the strong lensing regime from the magnification bias of background galaxies. The method relies upon pixellisation of the surface mass density distribution which allows us to write down a simple, solvable set of equations. We also show how pixellisation can be applied to methods of mass determination from measurements of shear and present a simplified method of application. The method is demonstrated with cluster models and applied to magnification data from the lensing cluster Abell 1689.

We present the first application of lens magnification to measure the absolute mass of a galaxy cluster; Abell 1689. The absolute mass of a galaxy cluster can be measured by the gravitational lens magnification of a background galaxy population by the cluster potential. The lensing signal is complicated by the variation in number counts due to galaxy clustering and shot-noise, and by additional uncertainties in relating magnification to mass in the strong lensing regime. Clustering and shot-noise can be dealt with using maximum likelihood methods. Local approximations can then be used to estimate the mass from magnification. Alternatively if the lens is axially symmetric we show that the amplification equation can be solved nonlocally for the surface mass density and the tangential shear. In this paper we present the first maps of the total mass distribution in Abell 1689, measured from the deficit of lensed red galaxies behind the cluster. Although noisier, these reproduce the main features of mass maps made using the shear distortion of background galaxies but have the correct normalisation, finally breaking the ``sheet-mass'' degeneracy that has plagued lensing methods based on shear. We derive the cluster mass profile in the inner 4' (0.48 Mpc/h). These show a profile with a near isothermal surface mass density \kappa = (0.5+/-0.1)(\theta/1')^{-1} out to a radius of 2.4' (0.28Mpc/h), followed by a sudden drop into noise. We find that the projected mass interior to 0.24 h^{-1}$Mpc is M(<0.24 Mpc/h)=(0.50+/- 0.09) \times 10^{15} Msol/h. We compare our results with masses estimated from X-ray temperatures and line-of-sight velocity dispersions, as well as weak shear and lensing arclets and find all are in fair agreement for Abell 1698.

We report the results of a search for neutrino-induced particle cascades using a deep ocean water Cherenkov detector. The effective mass of the detector, a string of seven 40 cm diameter photomultipliers at 5.2 m spacing, is found through simulation analysis to be surprisingly large: greater than 1 megaton of water at incident neutrino energies of 1 PeV. We find no evidence for neutrino-induced cascades in 18.6 hours of observation. Although the limit implied by this observation is the strongest yet for predictions of active galatic nuclei (AGN) neutrinos at energies above 100 TeV, perhaps the more intriguing result is that the power of these techniques can be exploited to test these AGN models in a relatively short time.