Matching the number counts of high-$z$ sub-millimetre-selected galaxies (SMGs) has been a long standing problem for galaxy formation models. In this paper, we use 3D dust radiative transfer to model the sub-mm emission from galaxies in the SIMBA cosmological hydrodynamic simulations, and compare predictions to the latest single-dish observational constraints on the abundance of 850$\mathrm{\mu m}$-selected sources. We find unprecedented agreement with the integrated 850$\mathrm{\mu m}$ luminosity function, along with good agreement in the redshift distribution of bright SMGs. The excellent agreement is driven primarily by SIMBA's good match to infrared measures of the star formation rate (SFR) function between $z = 2-4$ at high SFRs. Also important is the self-consistent on-the-fly dust model in SIMBA, which predicts, on average, higher dust masses (by up to a factor of 7) compared to using a fixed dust-to-metals ratio of 0.3. We construct a lightcone to investigate the effect of far-field blending, and find minimal contribution to the shape and normalisation of the luminosity function. We provide new fits to the 850$\mathrm{\mu m}$ luminosity as a function of SFR and dust mass. Our results demonstrate that exotic solutions to the discrepancy between sub-mm counts in simulations and observations, such as a top-heavy IMF, are unnecessary, and that sub-millimetre-bright phases are a natural consequence of massive galaxy evolution.

High-redshift Lyman-alpha blobs (LABs) are an enigmatic class of objects that have been the subject of numerous observational and theoretical investigations. It is of particular interest to determine the dominant power sources for the copious luminosity, as direct emission from HII regions, cooling gas, and fluorescence due to the presence of active galactic nuclei (AGN) can all contribute significantly. In this paper, we present the first theoretical model to consider all of these physical processes in an attempt to develop an evolutionary model for the origin of high-z LABs. This is achieved by combining a series of high-resolution cosmological zoom-in simulations with ionization and Lyman-alpha (Lya) radiative transfer models. We find that massive galaxies display a range of Lya luminosities and spatial extents (which strongly depend on the limiting surface brightness used) over the course of their lives, though regularly exhibit luminosities and sizes consistent with observed LABs. The model LABs are typically powered from a combination of recombination in star-forming galaxies, as well as cooling emission from gas associated with accretion. When AGN are included in the model, the fluorescence caused by AGN-driven ionization can be a significant contributor to the total Lya luminosity as well. We propose that the presence of an AGN may be predicted from the Gini coefficient of the blob's surface brightness. Within our modeled mass range, there are no obvious threshold physical properties that predict appearance of LABs, and only weak correlations of the luminosity with the physical properties of the host galaxy. This is because the emergent Lya luminosity from a system is a complex function of the gas temperature, ionization state, and Lya escape fraction.

We demonstrate how galaxy morphologies can be represented by weighted sums of 'eigengalaxies' and how eigengalaxies can be used in a probabilistic framework to enable principled and simplified approaches in a variety of applications. Eigengalaxies can be derived from a Principal Component Analysis (PCA) of sets of single- or multi-band images. They encode the image space equivalent of basis vectors that can be combined to describe the structural properties of large samples of galaxies in a massively reduced manner. As an illustration, we show how a sample of 10,243 galaxies in the Hubble Space Telescope CANDELS survey, which have visual classifications from Galaxy Zoo, can be represented by just 12 eigengalaxies whilst retaining 96% explained variance. We show how the emphasis of certain eigengalaxy components correspond to visual features (e.g. disc, point source, etc.), and explore the correspondence between combinations of eigengalaxies and features defined in the Galaxy Zoo-CANDELS catalogue. We also describe a probabilistic extension to PCA (PPCA) which enables the eigengalaxy framework to assign probabilities to galaxies and characterise a whole population as a generative distribution. We present four practical applications of the probabilistic eigengalaxy framework that are particularly relevant for the next generation of large imaging surveys: we (i) show how low probability galaxies make for natural candidates for outlier detection (ii) demonstrate how missing data can be predicted (iii) show how a similarity search can be performed on exemplars (iv) demonstrate how unsupervised clustering of objects can be implemented.

Lyman-$\alpha$ blobs (LABs) are spatially extended nebulae of emission in the Ly$\alpha$ line of hydrogen, seen at high redshifts$^{1,2}$, and most commonly found in the dense environment of star-forming galaxies$^{3,4}$. The origin of Ly$\alpha$ emission in the LABs is still unclear and under debate$^{5}$. Proposed powering sources generally fall into two categories: (1) photoionization, galactic super-winds/outflows, resonant scattering of Ly$\alpha$ photons from starbursts or active galactic nuclei (AGNs)$^{6,7,8,9,10}$ and (2) cooling radiation from cold streams of gas accreting onto galaxies$^{12}$. Here we analyze the gas kinematics within a LAB providing rare observational evidence for infalling gas. This is consistent with the release of gravitational accretion energy as cold streams radiate Ly$\alpha$ photons. It also provides direct evidence for possible cold streams feeding the central galaxies. The infalling gas is not important by mass but hints at more than one mechanism to explain the origin of the extended Ly$\alpha$ emission around young galaxies. It is also possible that the infalling gas may represent material falling back to the galaxy from where it originated, forming a galactic fountain.

We analyse 870um Atacama Large Millimetre Array (ALMA) dust continuum detections of 41 canonically-selected z~3 Lyman-break galaxies (LBGs), as well as 209 ALMA-undetected LBGs, in follow-up of SCUBA-2 mapping of the UKIDSS Ultra Deep Survey (UDS) field. We find that our ALMA-bright LBGs lie significantly off the locally calibrated IRX-beta relation and tend to have relatively bluer rest-frame UV slopes (as parametrised by beta), given their high values of the 'infrared excess' (IRX=L_IR/L_UV), relative to the average 'local' IRX-beta relation. We attribute this finding in part to the young ages of the underlying stellar populations but we find that the main reason behind the unusually blue UV slopes are the relatively shallow slopes of the corresponding dust attenuation curves. We show that, when stellar masses are being established via SED fitting, it is absolutely crucial to allow the attenuation curves to vary (rather than fixing it on Calzetti-like law), where we find that the inappropriate curves may underestimate the resulting stellar masses by a factor of ~2-3x on average. In addition, we find these LBGs to have relatively high specific star-formation rates (sSFRs), dominated by the dust component, as quantified via the fraction of obscured star formation ( f_obs = SFR_IR/SFR_(UV+IR)). We conclude that the ALMA-bright LBGs are, by selection, massive galaxies undergoing a burst of a star formation (large sSFRs, driven, for example, by secular or merger processes), with a likely geometrical disconnection of the dust and stars, responsible for producing shallow dust attenuation curves.

We present sensitive 850$\mu$m imaging of the COSMOS field using 640hr of new and archival observations taken with SCUBA-2 at the East Asian Observatory's James Clerk Maxwell Telescope. The SCUBA-2 COSMOS survey (S2COSMOS) achieves a median noise level of $\sigma_{850\mu{\mathrm{m}}}$=1.2mJy/beam over an area of 1.6 sq. degree (MAIN; HST/ACS footprint), and $\sigma_{850\mu{\mathrm{m}}}$=1.7mJy/beam over an additional 1 sq. degree of supplementary (SUPP) coverage. We present a catalogue of 1020 and 127 sources detected at a significance level of >4$\sigma$ and >4.3$\sigma$ in the MAIN and SUPP regions, respectively, corresponding to a uniform 2% false-detection rate. We construct the single-dish 850$\mu$m number counts at $S_{850}$>2mJy and show that these S2COSMOS counts are in agreement with previous single-dish surveys, demonstrating that degree-scale fields are sufficient to overcome the effects of cosmic variance in the $S_{850}$=2-10mJy population. To investigate the properties of the galaxies identified by S2COSMOS sources we measure the surface density of near-infrared-selected galaxies around their positions and identify an average excess of 2.0$\pm$0.2 galaxies within a 13$''$ radius (~100kpc at $z$~2). The bulk of these galaxies represent near-infrared-selected SMGs and/or spatially-correlated sources and lie at a median photometric redshift of $z$=2.0$\pm$0.1. Finally, we perform a stacking analysis at sub-millimeter and far-infrared wavelengths of stellar-mass-selected galaxies ($M_{\star}$=10$^{10}$-10$^{12}{\rm M_{\odot}}$) from $z$=0-4, obtaining high-significance detections at 850um in all subsets (SNR=4-30), and investigate the relation between far-infrared luminosity, stellar mass, and the peak wavelength of the dust SED. The publication of this survey adds a new deep, uniform sub-millimeter layer to the wavelength coverage of this well-studied COSMOS field.

Ninety per cent of baryons are located outside galaxies, either in the circumgalactic or intergalactic medium. Theory points to galactic winds as the primary source of the enriched and massive circumgalactic medium. Winds from compact starbursts have been observed to flow to distances somewhat greater than ten kiloparsecs, but the circumgalactic medium typically extends beyond a hundred kiloparsecs. Here we report optical integral field observations of the massive but compact galaxy SDSS J211824.06+001729.4. The oxygen [O II] lines at wavelengths of 3726 and 3729 angstroms reveal an ionized outflow spanning 80 by 100 square kiloparsecs, depositing metal-enriched gas at 10,000 kelvin through an hourglass-shaped nebula that resembles an evacuated and limb-brightened bipolar bubble. We also observe neutral gas phases at temperatures of less than 10,000 kelvin reaching distances of 20 kiloparsecs and velocities of around 1,500 kilometres per second. This multi-phase outflow is probably driven by bursts of star formation, consistent with theory.

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

Galaxy morphology is a fundamental quantity, that is essential not only for the full spectrum of galaxy-evolution studies, but also for a plethora of science in observational cosmology. While a rich literature exists on morphological-classification techniques, the unprecedented data volumes, coupled, in some cases, with the short cadences of forthcoming 'Big-Data' surveys (e.g. from the LSST), present novel challenges for this field. Large data volumes make such datasets intractable for visual inspection (even via massively-distributed platforms like Galaxy Zoo), while short cadences make it difficult to employ techniques like supervised machine-learning, since it may be impractical to repeatedly produce training sets on short timescales. Unsupervised machine learning, which does not require training sets, is ideally suited to the morphological analysis of new and forthcoming surveys. Here, we employ an algorithm that performs clustering of graph representations, in order to group image patches with similar visual properties and objects constructed from those patches, like galaxies. We implement the algorithm on the Hyper-Suprime-Cam Subaru-Strategic-Program Ultra-Deep survey, to autonomously reduce the galaxy population to a small number (160) of 'morphological clusters', populated by galaxies with similar morphologies, which are then benchmarked using visual inspection. The morphological classifications (which we release publicly) exhibit a high level of purity, and reproduce known trends in key galaxy properties as a function of morphological type at z<1 (e.g. stellar-mass functions, rest-frame colours and the position of galaxies on the star-formation main sequence). Our study demonstrates the power of unsupervised machine learning in performing accurate morphological analysis, which will become indispensable in this new era of deep-wide surveys.

We analyse the physical properties of a large, homogeneously selected sample of ALMA-located sub-mm galaxies (SMGs) detected in the SCUBA-2 CLS 850-$\mu$m map of the UKIDSS/UDS field. This survey, AS2UDS, identified 707 SMGs across the ~1 sq.deg. field, including ~17 per cent which are undetected in the optical/near-infrared to $K$>~25.7 mag. We interpret the UV-to-radio data using a physically motivated model, MAGPHYS and determine a median photometric redshift of z=2.61+-0.08, with a 68th percentile range of z=1.8-3.4 and just ~6 per cent at z>4. The redshift distribution is well fit by a model combining evolution of the gas fraction in halos with the growth of halo mass past a threshold of ~4x10$^{12}$M$_\odot$, thus SMGs may represent the highly efficient collapse of gas-rich massive halos. Our survey provides a sample of the most massive, dusty galaxies at z>~1, with median dust and stellar masses of $M_d$=(6.8+-0.3)x10$^{8}$M$_\odot$ (thus, gas masses of ~10$^{11}$M$_\odot$) and $M_\ast=$(1.26+-0.05)x10$^{11}$M$_\odot$. These galaxies have gas fractions of $f_{gas}=$0.41+-0.02 with depletion timescales of ~150Myr. The gas mass function evolution at high masses is consistent with constraints at lower masses from blind CO-surveys, with an increase to z~2-3 and then a decline at higher redshifts. The space density and masses of SMGs suggests that almost all galaxies with $M_\ast$>~2x10$^{11}$M$_\odot$ have passed through an SMG-like phase. We find no evolution in dust temperature at a constant far-infrared luminosity across z~1.5-4. We show that SMGs appear to behave as simple homologous systems in the far-infrared, having properties consistent with a centrally illuminated starburst. Our study provides strong support for an evolutionary link between the active, gas-rich SMG population at z>1 and the formation of massive, bulge-dominated galaxies across the history of the Universe.

We identify multi-wavelength counterparts to 1,147 submillimeter sources from the S2COSMOS SCUBA-2 survey of the COSMOS field by employing a recently developed radio$+$machine-learning method trained on a large sample of ALMA-identified submillimeter galaxies (SMGs), including 260 SMGs identified in the AS2COSMOS pilot survey. In total, we identify 1,222 optical/near-infrared(NIR)/radio counterparts to the 897 S2COSMOS submillimeter sources with S$_{850}$>1.6mJy, yielding an overall identification rate of ($78\pm9$)%. We find that ($22\pm5$)% of S2COSMOS sources have multiple identified counterparts. We estimate that roughly 27% of these multiple counterparts within the same SCUBA-2 error circles very likely arise from physically associated galaxies rather than line-of-sight projections by chance. The photometric redshift of our radio$+$machine-learning identified SMGs ranges from z=0.2 to 5.7 and peaks at $z=2.3\pm0.1$. The AGN fraction of our sample is ($19\pm4$)%, which is consistent with that of ALMA SMGs in the literature. Comparing with radio/NIR-detected field galaxy population in the COSMOS field, our radio+machine-learning identified counterparts of SMGs have the highest star-formation rates and stellar masses. These characteristics suggest that our identified counterparts of S2COSMOS sources are a representative sample of SMGs at z<3. We employ our machine-learning technique to the whole COSMOS field and identified 6,877 potential SMGs, most of which are expected to have submillimeter emission fainter than the confusion limit of our S2COSMOS surveys (S$_{850}$<1.5mJy). We study the clustering properties of SMGs based on this statistically large sample, finding that they reside in high-mass dark matter halos ($(1.2\pm0.3)\times10^{13}\,h^{-1}\,\rm M_{\odot}$), which suggests that SMGs may be the progenitors of massive ellipticals we see in the local Universe.

The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.

We present James Clerk Maxwell Telescope Submillimetre Common-User Bolometer Array 2 (SCUBA-2) 850 & 450 $\mu$m observations ($\sigma_{850}\sim0.5$ mJy, $\sigma_{450}\sim5$ mJy) of the HS1549+19 and HS1700+64 survey fields containing two of the largest known galaxy over-densities at $z=2.85$ and $2.30$, respectively. We detect 56 sub-millimetre galaxies (SMGs) with SNR $> 4$ over $\sim 50$ arcmin$^2$ at 850 $\mu$m with flux densities of 3 - 17 mJy. The number counts indicate over-densities in the $3$-arcmin diameter core region ($\sim 1.5$ Mpc at $z=2.5$) of $6^{+4}_{-2}\times$ (HS1549) and $4^{+6}_{-2}\times$ (HS1700) compared to blank field surveys. Within these core regions, we spectroscopically confirm that approximately one third of the SMGs lie at the protocluster redshifts for both HS1549 and HS1700. We use statistical identifications of other SMGs in the wider fields to constrain an additional four candidate protocluster members in each system. We combine multi wavelength estimates of the star-formation rates (SFRs) from Lyman-break dropout- and narrowband-selected galaxies, and the SCUBA-2 SMGs, to estimate total SFRs of 12,$500\pm2800$ M$_\odot$ yr$^{-1}$ ($4900\pm1200$ M$_\odot$ yr$^{-1}$) in HS1549 (HS1700), and SFR densities (SFRDs) within the central 1.5-Mpc diameter of each protocluster to be $3000\pm900$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ ($1300\pm400$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$) in the HS1549 (HS1700) protocluster, $\sim10^4\times$ larger than the global SFRDs found at their respective epochs, due to the concentration of star-forming galaxies in the small volume of the dense cluster cores. Our results suggest centrally concentrated starbursts within protoclusters may be a relatively common scenario for the build up of mass in rich clusters assembling at $z\gtrsim2$.

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

We present adaptive optics assisted integral field spectroscopy of 34 star-forming galaxies at $z$ = 0.8-3.3 selected from the HiZELS narrow-band survey. We measure the kinematics of the ionised interstellar medium on $\sim$1 kpc scales, and show that the galaxies are turbulent, with a median ratio of rotational to dispersion support of $v$/$\sigma$=0.82$\pm$0.13. We combine the dynamics with high-resolution rest-frame optical imaging and extract emission line rotation curves. We show that high-redshift star-forming galaxies follow a similar power-law trend in specific angular momentum with stellar mass as that of local late type galaxies. We exploit the high resolution of our data and examine the radial distribution of angular momentum within each galaxy by constructing total angular momentum profiles. Although the stellar mass of a typical star-forming galaxy is expected to grow by a factor $\sim$8 in the $\sim$5 Gyrs between $z$$\sim$3.3 and $z$$\sim$0.8, we show that the internal distribution of angular momentum becomes less centrally concentrated in this period i.e the angular momentum grows outwards. To interpret our observations, we exploit the EAGLE simulation and trace the angular momentum evolution of star forming galaxies from $z$$\sim$3 to $z$$\sim$0, identifying a similar trend of decreasing angular momentum concentration. This change is attributed to a combination of gas accretion in the outer disk, and feedback that preferentially arises from the central regions of the galaxy. We discuss how the combination of the growing bulge and angular momentum stabilises the disk and gives rise to the Hubble sequence.

Generative Adversarial Networks (GANs) are a class of artificial neural network that can produce realistic, but artificial, images that resemble those in a training set. In typical GAN architectures these images are small, but a variant known as Spatial-GANs (SGANs) can generate arbitrarily large images, provided training images exhibit some level of periodicity. Deep extragalactic imaging surveys meet this criteria due to the cosmological tenet of isotropy. Here we train an SGAN to generate images resembling the iconic Hubble Space Telescope eXtreme Deep Field (XDF). We show that the properties of 'galaxies' in generated images have a high level of fidelity with galaxies in the real XDF in terms of abundance, morphology, magnitude distributions and colours. As a demonstration we have generated a 7.6-billion pixel 'generative deep field' spanning 1.45 degrees. The technique can be generalised to any appropriate imaging training set, offering a new purely data-driven approach for producing realistic mock surveys and synthetic data at scale, in astrophysics and beyond.

Feedback processes from active galactic nuclei (AGN) are thought to play a crucial role in regulating star formation in massive galaxies. Previous studies using \textit{Herschel} have resulted in conflicting conclusions as to whether star formation is quenched, enhanced, or not affected by AGN feedback. We use new deep 850 $\mu$m observations from the SCUBA-2 Cosmology Legacy survey (S2CLS) to investigate star formation in a sample of X-ray selected AGN, probing galaxies up to $L_{0.5-7~\rm keV} = 10^{46}$ erg s$^{-1}$. Here we present the results of our analysis on a sample of 1957 galaxies at $ 1 < z < 3 $, using both S2CLS and ancilliary data at seven additional wavelengths (24--500 \mcm) from \textit{Herschel} and \textit{Spitzer}. We perform a stacking analysis, binning our sample by redshift and X-ray luminosity. By fitting analytical spectral energy distributions (SEDs) to decompose contributions from cold and warm dust, we estimate star-formation rates for each `average' source. We find that the average AGN in our sample resides in a star-forming host galaxy, with SFRs ranging from 80 - 600 $M_{\odot}$ year$^{-1}$. Within each redshift bin, we see no trend of SFR with X-ray luminosity, instead finding a flat distribution of SFR across $\sim$3 orders of magnitude of AGN luminosity. By studying instantaneous X-ray luminosities and SFRs, we find no evidence that AGN activity affects star formation in host galaxies.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

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

We measure the average deflection of cosmic microwave background photons by quasars at $\langle z \rangle =1.7$. Our sample is selected from the Sloan Digital Sky Survey to cover the redshift range $0.9\leq z\leq2.2$ with absolute i-band magnitudes of $M_i\leq-24$ (K-corrected to z=2). A stack of nearly 200,000 targets reveals an 8$\sigma$ detection of Planck's estimate of the lensing convergence towards the quasars. We fit the signal with a model comprising a Navarro-Frenk-White density profile and a 2-halo term accounting for correlated large scale structure, which dominates the observed signal. The best-fitting model is described by an average halo mass $\log_{10}(M_{\rm h}/h^{-1}M_\odot)=12.6\pm0.2$ and linear bias $b=2.7\pm0.3$ at $z=1.7$, in excellent agreement with clustering studies. We also report of a hint, at a 90% confidence level, of a correlation between the convergence amplitude and luminosity, indicating that quasars brighter than $M_i\lesssim -26$ reside in halos of typical mass ${M_{\rm h}\approx 10^{13}\,h^{-1}M_\odot}$, scaling roughly as ${M_{\rm h}\propto L_{\rm opt}^{3/4}}$ at ${M_i\lesssim-24}$, in good agreement with physically-motivated quasar demography models. Although we acknowledge this luminosity dependence is a marginal result, the observed $M_{\rm h}$-$L_{\rm opt}$ relationship could be interpreted as a reflection of the cutoff in the distribution of black hole accretion rates towards high Eddington ratios: the weak trend of $M_{\rm h}$ with $L_{\rm opt}$ observed at low luminosity becomes stronger for the most powerful quasars, which tend to be accreting close to the Eddington limit.

We present a new measurement of the linear bias of radio loud active galactic nuclei (RLAGN) at $z\approx0.3$ and $L_{\rm 1.4GHz}>10^{23}\,{\rm W\,Hz^{-1}}$ selected from the Best & Heckman (2012) sample, made by cross-correlating the RLAGN surface density with a map of the convergence of the weak lensing field of the cosmic microwave background from Planck. We detect the cross-power signal at a significance of $3\sigma$ and use the amplitude of the cross-power spectrum to estimate the linear bias of RLAGN, $b=2.5 \pm 0.8$, corresponding to a typical dark matter halo mass of $\log_{10}(M_{\rm h} /h^{-1} M_\odot)=14.0^{+0.3}_{-0.5}$. When RLAGN associated with optically-selected clusters are removed we measure a lower bias corresponding to $\log_{10}(M_{\rm h} /h^{-1} M_\odot)=13.7^{+0.4}_{-1.0}$. These observations support the view that powerful RLAGN typically inhabit rich group and cluster environments.

Future radio surveys will generate catalogues of tens of millions of radio sources, for which redshift estimates will be essential to achieve many of the science goals. However, spectroscopic data will be available for only a small fraction of these sources, and in most cases even the optical and infrared photometry will be of limited quality. Furthermore, radio sources tend to be at higher redshift than most optical sources and so a significant fraction of radio sources hosts differ from those for which most photometric redshift templates are designed. We therefore need to develop new techniques for estimating the redshifts of radio sources. As a starting point in this process, we evaluate a number of machine-learning techniques for estimating redshift, together with a conventional template-fitting technique. We pay special attention to how the performance is affected by the incompleteness of the training sample and by sparseness of the parameter space or by limited availability of ancillary multi-wavelength data. As expected, we find that the quality of the photometric-redshift degrades as the quality of the photometry decreases, but that even with the limited quality of photometry available for all sky-surveys, useful redshift information is available for the majority of sources, particularly at low redshift. We find that a template-fitting technique performs best with high-quality and almost complete multi-band photometry, especially if radio sources that are also X-ray emitting are treated separately. When we reduced the quality of photometry to match that available for the EMU all-sky radio survey, the quality of the template-fitting degraded and became comparable to some of the machine learning methods. Machine learning techniques currently perform better at low redshift than at high redshift, because of incompleteness of the currently available training data at high redshifts.

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.

We present Atacama Large Millimeter/submillimeter Array observations of a massive (M_stars~10^11 M_Sun) compact (r_e,UV~100 pc) merger remnant at z=0.66 that is driving a 1000 km/s outflow of cool gas, with no observational trace of an active galactic nucleus (AGN). We resolve molecular gas on scales of approximately 1-2 kpc, and our main finding is the discovery of a wing of blueshifted CO(2-1) emission out to -1000 km/s relative to the stars. We argue that this is the molecular component of a multiphase outflow, expelled from the central starburst within the past 5 Myr through stellar feedback, although we cannot rule out previous AGN activity as a launching mechanism. If the latter is true, then this is an example of a relic multiphase AGN outflow. We estimate a molecular mass outflow rate of approximately 300 M_Sun/yr, or about one third of the 10 Myr-averaged star formation rate. This system epitomizes the multiphase 'blowout' episode following a dissipational major merger - a process that has violently quenched central star formation and supermassive black hole growth.

Measuring molecular gas mass is vital for understanding the evolution of galaxies at high redshifts (z$\geq$1). Most measurements rely on CO as a tracer, but dependences on metallicity, dynamics and surface density lead to systematic uncertainties in high-z galaxies, where these physical properties are difficult to observe, and where the physical environments can differ systematically from those at z=0. Dust continuum emission provides a potential alternative assuming a known dust/gas ratio, but this must be calibrated on a direct gas tracer at z$\geq$1. In this paper we consider the [CI] 492-GHz emission line, which has been shown to trace molecular gas closely throughout Galactic clouds and has the advantages of being optically thin in typical conditions (unlike CO), and being observable at accessible frequencies at high redshifts (in contrast to the low-excitation lines of CO). We use the Atacama Large Millimetre/submillimetre Array (ALMA) to measure [CI], CO(4-3) and dust emission in a representative sample of star-forming galaxies at z=1, and combine these data with multi-wavelength spectral energy distributions to study relationships between dust and gas components of galaxies. We uncover a strong [CI]-dust correlation, suggesting that both trace similar phases of the gas. By incorporating other samples from the literature, we show that this correlation persists over a wide range of luminosities and redshifts up to z$\sim$4. Finally we explore the implications of our results as an independent test of literature calibrations for dust as a tracer of gas mass, and for predicting the CI abundance.

We present the detection of a bright 3mm continuum source in the Brightest Cluster Galaxy (BCG) in Abell 851 (z=0.411) with the NOrthern Extended Millimeter Array (NOEMA). When this detection is compared to other multi-frequency observations across 21cm-100um, including new Arcminute Microkelvin Imager 15GHz observations, we find evidence for a relatively flat, variable core source associated with the BCG. The radio power and amplitude of variability observed in this galaxy is consistent with the cores in lower redshift BCGs in X-ray-selected clusters, and the flat mm-cm spectrum is suggestive of the BCG being a low luminosity AGN archetype. The discovery of this system could provide a basis for a long-term study of the role of low luminosity radio mode 'regulatory' feedback in massive clusters.

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

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

We have used the Submillimeter Array at 860$\,\mu$m to observe the brightest SCUBA-2 sources in 4$\,$deg$^{2}$ of the Cosmology Legacy Survey. We have targeted 75 of the brightest single-dish SCUBA-2 850$\,\mu$m sources down to $S_{850}\,{\approx}\,8\,$mJy, achieving an average synthesized beam of 2.4$^{\prime\prime}$ and an average rms of $\sigma_{860}\,{=}\,1.5\,$mJy in our primary beam-corrected maps. We searched our maps for $4\sigma$ peaks, corresponding to $S_{860}\,{\gtrsim}\,6\,$mJy sources, and detected 59 single galaxies and three pairs of galaxies. We include in our study 28 archival observations, bringing our sample size to 103 bright single-dish submillimetre sources with interferometric follow-up. We compute the cumulative and differential number counts of our sample, finding them to overlap with previous single-dish survey number counts within the uncertainties, although our cumulative number count is systematically lower than the parent SCUBA-2 cumulative number count by $24\,{\pm}\,6$ per cent between 11 and 15$\,$mJy. We estimate the probability that a ${\gtrsim}\,10\,$mJy single-dish submillimetre source resolves into two or more galaxies with similar flux densities, causing a significant change in the number counts, to be about 15 per cent. Assuming the remaining 85 per cent of the targets are ultra-luminous starburst galaxies between $z\,{=}\,2$-3, we find a likely volume density of ${\gtrsim}\,400\,$M$_{\odot}\,$yr$^{-1}$ sources to be ${\sim}\,3^{+0.7}_{-0.6}\,{\times}\,10^{-7}\,$Mpc$^{-3}$. We show that the descendants of these galaxies could be ${\gtrsim}\,4\,{\times}\,10^{11}\,$M$_{\odot}$ local quiescent galaxies, and that about 10 per cent of their total stellar mass would have formed during these short bursts of star-formation.

We present deep high resolution (0.03", 200pc) ALMA Band 7 observations covering the dust continuum and [CII] $\lambda157.7\mu$m emission in four $z\sim4.4-4.8$ sub-millimeter galaxies (SMGs) selected from the ALESS and AS2UDS surveys. The data show that the rest-frame 160$\mu$m (observed 345 GHz) dust emission is consistent with smooth morphologies on kpc scales for three of the sources. One source, UDS47.0, displays apparent substructure but this is also consistent with a smooth morphology, as indicated by simulations showing that smooth exponential disks can appear clumpy when observed at high angular resolution (0.03") and depth of these observations ($\sigma_{345\text{GHz}} \sim27-47\mu$Jy beam$^{-1}$). The four SMGs are bright [CII] emitters, and we extract [CII] spectra from the high resolution data, and recover $\sim20-100$% of the [CII] flux and $\sim40-80$% of the dust continuum emission, compared to the previous lower resolution observations. When tapered to 0.2" resolution our maps recover $\sim80-100$% of the continuum emission, indicating that $\sim60$% of the emission is resolved out on $\sim200$pc scales. We find that the [CII] emission in high-redshift galaxies is more spatially extended than the rest-frame 160$\mu$m dust continuum by a factor of $1.6\pm0.4$. By considering the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of the star-formation rate surface density ($\Sigma_{\text{SFR}}$) we revisit the [CII] deficit, and suggest that the decline in the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of $\Sigma_{\text{SFR}}$ is consistent with local processes. We also explore the physical drivers that may be responsible for these trends and can give rise to the properties found in the densest regions of SMGs.

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

Interferometric observations have demonstrated that a significant fraction of single-dish submillimetre (submm) sources are blends of multiple submm galaxies (SMGs), but the nature of this multiplicity, i.e. whether the galaxies are physically associated or chance projections, has not been determined. We performed spectroscopy of 11 SMGs in six multi-component submm sources, obtaining spectroscopic redshifts for nine of them. For an additional two component SMGs, we detected continuum emission but no obvious features. We supplement our observed sources with four sources from the literature. This sample allows us to statistically constrain the physical nature of single-dish submm source multiplicity for the first time. In three [3/7, or 43 (-33/+39) per cent at 95% confidence] of the single-dish sources for which the nature of the blending is unambiguous, the components for which spectroscopic redshifts are available are physically associated, whereas 4/7 [57 (-39/+33) per cent] have at least one unassociated component. When components whose spectra exhibit continuum but no features and for which the photometric redshift is significantly different from the spectroscopic redshift of the other component are also considered, 6/9 [67 (-37/+26) per cent] of the single-dish sources are comprised of at least one unassociated component. The nature of the multiplicity of one single-dish source is ambiguous. We conclude that physically associated systems and chance projections both contribute to the multi-component single-dish submm source population. This result contradicts the conventional wisdom that bright submm sources are solely a result of merger-induced starbursts, as blending of unassociated galaxies is also important.

High-redshift, luminous, dusty star forming galaxies (DSFGs) constrain the extremity of galaxy formation theories. The most extreme are discovered through follow-up on candidates in large area surveys. Here we present 850 $\mu$m SCUBA-2 follow-up observations of 188 red DSFG candidates from the \textit{Herschel} Multi-tiered Extragalactic Survey (HerMES) Large Mode Survey, covering 274 deg$^2$. We detected 87 per cent with a signal-to-noise ratio $>$ 3 at 850~$\mu$m. We introduce a new method for incorporating the confusion noise in our spectral energy distribution fitting by sampling correlated flux density fluctuations from a confusion limited map. The new 850~$\mu$m data provide a better constraint on the photometric redshifts of the candidates, with photometric redshift errors decreasing from $\sigma_z/(1+z)\approx0.21$ to $0.15$. Comparison spectroscopic redshifts also found little bias ($\langle (z-z_{\rm spec})/(1+z_{\rm spec})\rangle = 0.08 $). The mean photometric redshift is found to be 3.6 with a dispersion of $0.4$ and we identify 21 DSFGs with a high probability of lying at $z > 4$. After simulating our selection effects we find number counts are consistent with phenomenological galaxy evolution models. There is a statistically significant excess of WISE-1 and SDSS sources near our red galaxies, giving a strong indication that lensing may explain some of the apparently extreme objects. Nevertheless, our sample should include examples of galaxies with the highest star formation rates in the Universe ($\gg10^3$ M$_\odot$yr$^{-1}$).

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

We study the heating mechanisms and Ly{\alpha} escape fractions of 35 Ly{\alpha} blobs (LABs) at z = 3.1 in the SSA22 field. Dust continuum sources have been identified in 11 of the 35 LABs, all with star formation rates (SFRs) above 100 Msun/yr. Likely radio counterparts are detected in 9 out of 29 investigated LABs. The detection of submm dust emission is more linked to the physical size of the Ly{\alpha} emission than to the Ly{\alpha} luminosities of the LABs. A radio excess in the submm/radio detected LABs is common, hinting at the presence of active galactic nuclei. Most radio sources without X-ray counterparts are located at the centers of the LABs. However, all X-ray counterparts avoid the central regions. This may be explained by absorption due to exceptionally large column densities along the line-of-sight or by LAB morphologies, which are highly orientation dependent. The median Ly{\alpha} escape fraction is about 3\% among the submm-detected LABs, which is lower than a lower limit of 11\% for the submm-undetected LABs. We suspect that the large difference is due to the high dust attenuation supported by the large SFRs, the dense large-scale environment as well as large uncertainties in the extinction corrections required to apply when interpreting optical data.

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

We present an unsupervised machine learning technique that automatically segments and labels galaxies in astronomical imaging surveys using only pixel data. Distinct from previous unsupervised machine learning approaches used in astronomy we use no pre-selection or pre-filtering of target galaxy type to identify galaxies that are similar. We demonstrate the technique on the HST Frontier Fields. By training the algorithm using galaxies from one field (Abell 2744) and applying the result to another (MACS0416.1-2403), we show how the algorithm can cleanly separate early and late type galaxies without any form of pre-directed training for what an 'early' or 'late' type galaxy is. We then apply the technique to the HST CANDELS fields, creating a catalogue of approximately 60,000 classifications. We show how the automatic classification groups galaxies of similar morphological (and photometric) type, and make the classifications public via a catalogue, a visual catalogue and galaxy similarity search. We compare the CANDELS machine-based classifications to human-based classifications from the Galaxy Zoo: CANDELS project. Although there is not a direct mapping between Galaxy Zoo and our hierarchical labelling, we demonstrate a good level of concordance between human and machine classifications. Finally, we show how the technique can be used to identify rarer objects and present new lensed galaxy candidates from the CANDELS imaging.

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

Identifying galaxy clusters through overdensities of galaxies in photometric surveys is the oldest and arguably the most economic and mass-sensitive detection method, compared to X-ray and Sunyaev-Zel'dovich Effect surveys that detect the hot intracluster medium. However, a perennial problem has been the mapping of optical 'richness' measurements on to total cluster mass. Emitted at a conformal distance of 14 Gpc, the cosmic microwave background acts as a backlight to all intervening mass in the Universe, and therefore has been gravitationally lensed. Here we present a calibration of cluster optical richness at the 10 per cent level by measuring the average cosmic microwave background lensing convergence measured by Planck towards the positions of large numbers of optically-selected clusters, detecting the deflection of photons by haloes of total mass of the order 10**14 solar masses. Although mainly aimed at the study of larger-scale structures, the Planck lensing reconstruction can yield nearly unbiased results for stacked clusters on arcminute scales. The lensing convergence only depends on the redshift integral of the fractional overdensity of matter, so this approach offers a clean measure of cluster mass over most of cosmic history, largely independent of baryon physics.

We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dusty interstellar medium. Our main goals are to understand the origin of, and scatter in the IRX-beta relation; to assess the efficacy of simplified stellar population synthesis screen models in capturing the essential physics in the IRX-beta relation; and to understand systematic deviations from the canonical local IRX-beta relations in particular populations of high-redshift galaxies. Our main results follow. Galaxies that have young stellar populations with relatively cospatial UV and IR emitting regions and a Milky Way-like extinction curve fall on or near the standard Meurer relation. This behaviour is well captured by simplified screen models. Scatter in the IRX-beta relation is dominated by three major effects: (i) older stellar populations drive galaxies below the relations defined for local starbursts due to a reddening of their intrinsic UV SEDs; (ii) complex geometries in high-z heavily star forming galaxies drive galaxies toward blue UV slopes owing to optically thin UV sightlines; (iii) shallow extinction curves drive galaxies downward in the IRX-beta plane due to lowered NUV/FUV extinction ratios. We use these features of the UV slopes of galaxies to derive a fitting relation that reasonably collapses the scatter back toward the canonical local relation. Finally, we use these results to develop an understanding for the location of two particularly enigmatic populations of galaxies in the IRX-beta plane: z~2-4 dusty star forming galaxies, and z>5 star forming galaxies.

We present the VISTA-CFHT Stripe 82 (VICS82) survey: a near-infrared (J+Ks) survey covering 150 square degrees of the Sloan Digital Sky Survey (SDSS) equatorial Stripe 82 to an average depth of J=21.9 AB mag and Ks=21.4 AB mag (80% completeness limits; 5-sigma point source depths are approximately 0.5 mag brighter). VICS82 contributes to the growing legacy of multi-wavelength data in the Stripe 82 footprint. The addition of near-infrared photometry to the existing SDSS Stripe 82 coadd ugriz photometry reduces the scatter in stellar mass estimates to delta log(M_stellar)~0.3 dex for galaxies with M_stellar>10^9M_sun at z~0.5, and offers improvement compared to optical-only estimates out to z~1, with stellar masses constrained within a factor of approximately 2.5. When combined with other multi-wavelength imaging of the Stripe, including moderate-to-deep ultraviolet (GALEX), optical and mid-infrared (Spitzer IRAC) coverage, as well as tens of thousands of spectroscopic redshifts, VICS82 gives access to approximately 0.5 Gpc^3 of comoving volume. Some of the main science drivers of VICS82 include (a) measuring the stellar mass function of L^star galaxies out to z~1; (b) detecting intermediate redshift quasars at 2<z<3.5; (c) measuring the stellar mass function and baryon census of clusters of galaxies, and (d) performing optical/near-infrared-cosmic microwave background lensing cross-correlation experiments linking stellar mass to large-scale dark matter structure. Here we define and describe the survey, highlight some early science results and present the first public data release, which includes an SDSS-matched catalogue as well as the calibrated pixel data itself.

We present clustering analyses of identically-selected star-forming galaxies in 3 narrow redshift slices (at z=0.8, z=1.47 and z=2.23), from HiZELS, a deep, near-infrared narrow-band survey. The HiZELS samples span the peak in the cosmic star-formation rate density, identifying typical star-forming galaxies at each epoch. Narrow-band samples have well-defined redshift distributions and are therefore ideal for clustering analyses. We quantify the clustering of the three samples, and of H-alpha luminosity-selected subsamples, initially using simple power law fits to the two-point correlation function. We extend this work to link the evolution of star-forming galaxies and their host dark matter halos over cosmic time using sophisticated dark matter halo models. We find that the clustering strength, r0, and the bias of galaxy populations relative to the clustering of dark matter increase linearly with H-alpha luminosity (and, by implication, star-formation rate) at all three redshifts, as do the host dark matter halo masses of the HiZELS galaxies. The typical galaxies in our samples are star-forming centrals, residing in halos of mass M_halo ~ a few times 10^12M_solar. We find a remarkably tight redshift-independent relation between the H-alpha luminosity scaled by the characteristic luminosity, L(H-alpha)/L(H-alpha)*(z), and the minimum host dark matter halo mass of central galaxies. This reveals that the dark matter halo environment is a strong driver of galaxy star-formation rate and therefore of the evolution of the star-formation rate density in the Universe.

The identification of high-redshift massive galaxies with old stellar populations may pose challenges to some models of galaxy formation. However, to securely classify a galaxy as quiescent, it is necessary to exclude significant ongoing star formation, something that can be challenging to achieve at high redshift. In this letter, we analyse deep ALMA/870um and SCUBA-2/450um imaging of the claimed "post-starburst" galaxy ZF-20115 at z=3.717 that exhibits a strong Balmer break and absorption lines. The far-infrared imaging reveals a luminous starburst located 0.4+/-0.1 arcsec (~3kpc in projection) from the position of the rest-frame ultra-violet/optical emission, with an obscured star-formation rate of 100 Mo/yr. This star-forming component is undetected in the rest-frame ultraviolet but contributes significantly to the lower angular resolution photometry at restframe wavelengths >3500A, significantly complicating the determination of a reliable stellar mass. Importantly, in the presence of dust obscuration, strong Balmer features are not a unique signature of a post-starburst galaxy and are indeed frequently observed in infrared-luminous galaxies. We conclude that the ZF20015 system does not pose a challenge to current models of galaxy formation and that deep sub-/millimeter observations are a prerequisite for any claims of quiescence. The multi-wavelength observations of ZF20115 unveil a complex system with an intricate and spatially-varying star-formation history. ZF20115 demonstrates that understanding high-redshift obscured starbursts will only be possible with multi-wavelength studies that include high-resolution observations, available with the JWST, at mid-infrared wavelengths.

We present the results of a Spitzer/Herschel infrared photometric analysis of the largest (716) and highest-redshift (z=1.8) sample of Brightest Cluster Galaxies (BCGs), those from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). Given the tension that exists between model predictions and recent observations of BCGs at z<2, we aim to uncover the dominant physical mechanism(s) guiding the stellar-mass buildup of this special class of galaxies, the most massive in the Universe uniquely residing at the centres of galaxy clusters. Through a comparison of their stacked, broadband, infrared spectral energy distributions (SEDs) to a variety of SED model templates in the literature, we identify the major sources of their infrared energy output, in multiple redshift bins between 0 < z < 1.8. We derive estimates of various BCG physical parameters from the stacked {\nu}L{\nu} SEDs, from which we infer a star-forming, as opposed to a 'red and dead' population of galaxies, producing tens to hundreds of solar masses per year down to z=0.5. This discovery challenges the accepted belief that BCGs should only passively evolve through a series of gas-poor, minor mergers since z~4 (De Lucia & Blaizot 2007), but agrees with the improved semi-analytic model of hierarchical structure formation of Tonini et al. (2012), which predicts star-forming BCGs throughout the epoch considered. We attribute the star formation inferred from the stacked infrared SEDs to both major and minor 'wet' (gas-rich) mergers, based on a lack of key signatures (to date) of the cluster cooling flows to which BCG star formation is typically attributed, as well as a number of observational and simulation-based studies that support this scenario.

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

We investigate the far-infrared properties of galaxies selected via deep, narrow-band imaging of the H$\alpha$ emission line in four redshift slices from $z=0.40$--$2.23$ over $\sim 1$deg$^2$ as part of the High-redshift Emission Line Survey (HiZELS). We use a stacking approach in the Herschel PACS/SPIRE bands, along with $850\,\mu$m imaging from SCUBA-2 to study the evolution of the dust properties of H$\alpha$-emitters selected above an evolving characteristic luminosity threshold, $0.2L^\star_{{\rm H}\alpha}(z)$. We investigate the relationship between the dust temperatures and the far-infrared luminosities of our stacked samples, finding that H$\alpha$-selection identifies cold, low-$L_{\rm IR}$ galaxies ($T_{\rm dust}\sim 14$k; $\log[L_{\rm IR}/{\rm L}_\odot]\sim 9.9$) at $z=0.40$, and more luminous, warmer systems ($T_{\rm dust}\sim 34$k; $\log[L_{\rm IR}/{\rm L}_\odot]\sim 11.5$) at $z=2.23$. Using a modified greybody model, we estimate "characteristic sizes" for the dust-emitting regions of HiZELS galaxies of $\sim 0.5$kpc, nearly an order of magnitude smaller than their stellar continuum sizes, which may provide indirect evidence of clumpy ISM structure. Lastly, we measure the dust masses from our far-IR SEDs along with metallicity-dependent gas-to-dust ratios ($\delta_{\rm GDR}$) to measure typical molecular gas masses of $\sim 10^{10}$M$_\odot$ for these bright H$\alpha$-emitters. The gas depletion timescales are shorter than the Hubble time at each redshift, suggesting probable replenishment of their gas reservoirs from the intergalactic medium. Based on the number density of H$\alpha$-selected galaxies, we find that typical star-forming galaxies brighter than $0.2L^{\star}_{{\rm H}\alpha}(z)$ host a significant fraction ($35\pm10$%) of the total gas content of the Universe, consistent with the predictions of the latest cosmological simulations.

We present a new exploration of the cosmic star-formation history and dust obscuration in massive galaxies at redshifts $0.5< z<6$. We utilize the deepest 450 and 850$\mu$m imaging from SCUBA-2 CLS, covering 230arcmin$^2$ in the AEGIS, COSMOS and UDS fields, together with 100-250$\mu$m imaging from Herschel. We demonstrate the capability of the T-PHOT deconfusion code to reach below the confusion limit, using multi-wavelength prior catalogues from CANDELS/3D-HST. By combining IR and UV data, we measure the relationship between total star-formation rate (SFR) and stellar mass up to $z\sim5$, indicating that UV-derived dust corrections underestimate the SFR in massive galaxies. We investigate the relationship between obscuration and the UV slope (the IRX-$\beta$ relation) in our sample, which is similar to that of low-redshift starburst galaxies, although it deviates at high stellar masses. Our data provide new measurements of the total SFR density (SFRD) in $M_\ast>10^{10}M_\odot$ galaxies at $0.5<z<6$. This is dominated by obscured star formation by a factor of $>10$. One third of this is accounted for by 450$\mu$m-detected sources, while one fifth is attributed to UV-luminous sources (brighter than $L^\ast_{UV}$), although even these are largely obscured. By extrapolating our results to include all stellar masses, we estimate a total SFRD that is in good agreement with previous results from IR and UV data at $z\lesssim3$, and from UV-only data at $z\sim5$. The cosmic star-formation history undergoes a transition at $z\sim3-4$, as predominantly unobscured growth in the early Universe is overtaken by obscured star formation, driven by the build-up of the most massive galaxies during the peak of cosmic assembly.

We report the result from observations conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) to detect [CII] 158 um fine structure line emission from galaxies embedded in one of the most spectacular Lyman-alpha blobs (LABs) at z=3.1, SSA22-LAB1. Of three dusty star-forming galaxies previously discovered by ALMA 860 um dust continuum survey toward SSA22-LAB1, we detected the [CII] line from one, LAB1-ALMA3 at z=3.0993+/-0.0004. No line emission was detected, associated with the other ALMA continuum sources or from three rest-frame UV/optical selected z_spec~3.1 galaxies within the field of view. For LAB1-ALMA3, we find relatively bright [CII] emission compared to the infrared luminosity (L_[CII]/L_[CII]) and an extremely high [CII] 158 um and [NII] 205 um emission line ratio (L_[CII]/L_[NII]>55). The relatively strong [CII] emission may be caused by abundant photodissociation regions and sub-solar metallicity, or by shock heating. The origin of the unusually strong [CII] emission could be causally related to the location within the giant LAB, although the relationship between extended Lyman-alpha emission and ISM conditions of associated galaxies is yet to be understand.

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

In these proceedings we highlight the primary scientific goals and design of the WEAVE-LOFAR survey, which will use the new WEAVE spectrograph on the 4.2m William Herschel Telescope to provide the primary source of spectroscopic information for the LOFAR Surveys Key Science Project. Beginning in 2018, WEAVE-LOFAR will generate more than 10$^6$ R=5000 365-960 nm spectra of low-frequency selected radio sources, across three tiers designed to efficiently sample the redshift-luminosity plane, and produce a data set of enormous legacy value. The radio frequency selection, combined with the high multiplex and throughput of the WEAVE spectrograph, make obtaining redshifts in this way very efficient, and we expect that the redshift success rate will approach 100 per cent at $z < 1$. This unprecedented spectroscopic sample - which will be complemented by an integral field component - will be transformational in key areas, including studying the star formation history of the Universe, the role of accretion and AGN-driven feedback, properties of the epoch of reionisation, cosmology, cluster haloes and relics, as well as the nature of radio galaxies and protoclusters. Each topic will be addressed in unprecedented detail, and with the most reliable source classifications and redshift information in existence.

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

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $\sigma_{450}=1.2$ mJy/beam and $\sigma_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

Recent measurements of the dark matter halo masses of infrared-selected obscured quasars are in tension --- some indicate that obscured quasars have higher halo mass compared to their unobscured counterparts, while others find no difference. The former result is inconsistent with the simplest models of quasar unification that rely solely on viewing angle, while the latter may support such models. Here, using empirical relationships between dark matter halo and supermassive black hole masses, we provide a simple evolutionary picture that naturally explains these findings and is motivated by more sophisticated merger-driven quasar fueling models. The model tracks the growth rate of haloes, with the black hole growing in spurts of quasar activity in order to "catch-up" with the $M_{\textrm{BH}}$ - $M_{\textrm{star}}$ - $M_{\textrm{halo}}$ relationship. The first part of the quasar phase is obscured and is followed by an unobscured phase. Depending on the luminosity limit of the sample, driven by observational selection effects, a difference in halo masses may or may not be significant. For high luminosity samples, the difference can be large (a few to 10 times higher masses in obscured quasars), while for lower luminosity samples the halo mass difference is very small, much smaller than current observational constraints. Such a simple model provides a qualitative explanation for the higher mass haloes of obscured quasars, as well as rough quantitative agreement with seemingly disparate results.

Submillimetre galaxies (SMGs) are among the most luminous dusty galaxies in the Universe, but their true nature remains unclear; are SMGs the progenitors of the massive elliptical galaxies we see in the local Universe, or are they just a short-lived phase among more typical star-forming galaxies? To explore this problem further, we investigate the clustering of SMGs identified in the SCUBA-2 Cosmology Legacy Survey. We use a catalogue of submillimetre ($850\mu$m) source identifications derived using a combination of radio counterparts and colour/IR selection to analyse a sample of 914 SMGs in the UKIDSS Ultra Deep Survey (UDS), making this the largest high redshift sample of these galaxies to date. Using angular cross-correlation techniques, we estimate the halo masses for this large sample of SMGs and compare them with passive and star-forming galaxies selected in the same field. We find that SMGs, on average, occupy high-mass dark matter halos (M$_{\text{halo}} >10^{13}$M$_{\odot}$) at redshifts $z > 2.5$, consistent with being the progenitors of massive quiescent galaxies in present-day galaxy clusters. We also find evidence of downsizing, in which SMG activity shifts to lower mass halos at lower redshifts. In terms of their clustering and halo masses, SMGs appear to be consistent with other star-forming galaxies at a given redshift.

We have observed the dust continuum of ten z=3.1 Lyman Break Galaxies with the Atacama Large Millimeter/Submillimeter Array at ~450 mas resolution in Band 7. We detect and resolve the 870um emission in one of the targets with an integrated flux density of S(870)=(192+/-57) uJy, and measure a stacked 3-sigma signal of S(870)=(67+/-23) uJy for the remaining nine. The total infrared luminosities estimated from full spectral energy distribution fits are L(8-1000um)=(8.4+/-2.3)x10^10 Lsun for the detection and L(8-1000um)=(2.9+/-0.9)x10^10 Lsun for the stack. With HST ACS I-band imaging we map the rest-frame UV emission on the same scale as the dust, effectively resolving the 'infrared excess' (IRX=L_FIR/L_UV) in a normal galaxy at z=3. Integrated over the galaxy we measure IRX=0.56+/-0.15, and the galaxy-averaged UV slope is beta=-1.25+/-0.03. This puts the galaxy a factor of ~10 below the IRX-beta relation for local starburst nuclei of Meurer et al. (1999). However, IRX varies by more than a factor of 3 across the galaxy, and we conclude that the complex relative morphology of the dust relative to UV emission is largely responsible for the scatter in the IRX-beta relation at high-z. A naive application of a Meurer-like dust correction based on the UV slope would dramatically over-estimate the total star formation rate, and our results support growing evidence that when integrated over the galaxy, the typical conditions in high-z star-forming galaxies are not analogous to those in the local starburst nuclei used to establish the Meurer relation.

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

We investigate the properties of the galaxies selected from the deepest 850-micron survey undertaken to date with SCUBA-2 on the JCMT. This deep 850-micron imaging was taken in parallel with deep 450-micron imaging in the very best observing conditions as part of the SCUBA-2 Cosmology Legacy Survey. A total of 106 sources were uncovered at 850 microns from ~150, sq. arcmin in the centre of the COSMOS/UltraVISTA/CANDELS field, imaged to a typical rms depth of ~0.25 mJy. We utilise the wealth of available deep multi-frequency data to establish the complete redshift distribution for this sample, yielding <z> = 2.38 +- 0.09, a mean redshift comparable with that derived for all but the brightest previous sub-mm samples. We have also been able to establish the stellar masses of the majority of the galaxy identifications, enabling us to explore their location on the star-formation-rate:stellar-mass (SFR:M*) plane. Crucially, our new deep sample reaches flux densities equivalent to SFR ~ 100 Msun/yr, enabling us to confirm that sub-mm galaxies form the high-mass end of the `main sequence' (MS) of star-forming galaxies at z > 1.5 (with a mean specific SFR of sSFR = 2.25 +- 0.19 /Gyr at z ~ 2.5). Our results are consistent with no significant flattening of the MS towards high masses at these redshifts, suggesting that reports of such flattening possibly arise from under-estimates of dust-enshrouded star-formation activity in massive star-forming galaxies. However, our findings add to the growing evidence that average sSFR rises only slowly at high redshift, resulting in log(sSFR) being an apparently simple linear function of the age of the Universe.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

We exploit ALMA 870um observations to measure the star-formation rates (SFRs) of eight X-ray detected Active Galactic Nuclei (AGNs) in a z~3.1 protocluster, four of which reside in extended Ly-alpha haloes (often termed Ly-alpha blobs: LABs). Three of the AGNs are detected by ALMA and have implied SFRs of ~220-410~M_sun/yr; the non detection of the other five AGNs places SFR upper limits of <210 M_sun/yr. The mean SFR of the protocluster AGNs (~110-210 M_sun/yr) is consistent (within a factor of ~0.7-2.3) with that found for co-eval AGNs in the field, implying that galaxy growth is not significantly accelerated in these systems. However, when also considering ALMA data from the literature, we find evidence for elevated mean SFRs (up-to a factor of ~5.9 over the field) for AGNs at the protocluster core, indicating that galaxy growth is significantly accelerated in the central regions of the protocluster. We also show that all of the four protocluster LABs are associated with an ALMA counterpart within the extent of their Ly-alpha emission. The SFRs of the ALMA sources within the LABs (~150-410 M_sun/yr) are consistent with those expected for co-eval massive star-forming galaxies in the field. Furthermore, the two giant LABs (with physical extents of >100 kpc) do not host more luminous star formation than the smaller LABs, despite being an order of magnitude brighter in Ly-alpha emission. We use these results to discuss star formation as the power source of LABs.

We present a Herschel/SPIRE survey of three protoclusters at z=2-3 (2QZCluster, HS1700, SSA22). Based on the SPIRE colours (S350/S250 and S500/S350) of 250 $\mu$m sources, we selected high redshift dusty star-forming galaxies potentially associated with the protoclusters. In the 2QZCluster field, we found a 4-sigma overdensity of six SPIRE sources around 4.5' (~2.2 Mpc) from a density peak of H$\alpha$ emitters at z=2.2. In the HS1700 field, we found a 5-sigma overdensity of eight SPIRE sources around 2.1' (~1.0 Mpc) from a density peak of LBGs at z=2.3. We did not find any significant overdensities in SSA22 field, but we found three 500 $\mu$m sources are concentrated 3' (~1.4 Mpc) east to the LAEs overdensity. If all the SPIRE sources in these three overdensities are associated with protoclusters, the inferred star-formation rate densities are 10$^3$-10$^4$ times higher than the average value at the same redshifts. This suggests that dusty star-formation activity could be very strongly enhanced in z~2-3 protoclusters. Further observations are needed to confirm the redshifts of the SPIRE sources and to investigate what processes enhance the dusty star-formation activity in z~2-3 protoclusters.

Ly-alpha blobs (LABs) offer insight into the complex interface between galaxies and their circumgalactic medium. Whilst some LABs have been found to contain luminous star-forming galaxies and active galactic nuclei that could potentially power the Ly-alpha emission, others appear not to be associated with obvious luminous galaxy counterparts. It has been speculated that LABs may be powered by cold gas streaming on to a central galaxy, providing an opportunity to directly observe the `cold accretion' mode of galaxy growth. Star-forming galaxies in LABs could be dust obscured and therefore detectable only at longer wavelengths. We stack deep SCUBA-2 observations of the SSA22 field to determine the average 850um flux density of 34 LABs. We measure S_850 = 0.6 +/- 0.2mJy for all LABs, but stacking the LABs by size indicates that only the largest third (area > 1794 kpc^2) have a mean detection, at 4.5 sigma, with S_850 = 1.4 +/- 0.3mJy. Only two LABs (1 and 18) have individual SCUBA-2 > 3.5 sigma detections at a depth of 1.1mJy/beam. We consider two possible mechanisms for powering the LABs and find that central star formation is likely to dominate the emission of Ly-alpha, with cold accretion playing a secondary role.

To investigate what drives the reversal of the morphology-density relation at intermediate/high redshift, we present a multi-wavelength analysis of 27 dusty starburst galaxies in the massive cluster Cl 0024+17 at z = 0.4. We combine H-alpha dynamical maps from the VLT/FLAMES multi-IFU system with far-infrared imaging using Herschel SPIRE and millimetre spectroscopy from IRAM/NOEMA, in order to measure the dynamics, star formation rates and gas masses of this sample. Most galaxies appear to be rotationally supported, with a median ratio of rotational support to line-of-sight velocity dispersion v/sigma ~ 5 +/- 2, and specific angular momentum lambda_R = 0.83 +/- 0.06 - comparable to field spirals of a similar mass at this redshift. The star formation rates of 3 - 26 M_solar/yr and average 12 CO derived gas mass of 1 x 10^10 M_solar suggest gas depletion timescales of ~ 1Gyr (~ 0.25 of the cluster crossing time). We derive characteristic dust temperatures (mean T_dust = 26 +/- 1 K) consistent with local galaxies of similar far-infrared luminosity, suggesting that the low density gas is yet to be stripped. Taken together, these results suggest that these starbursts have only recently accreted from the field, with star formation rates likely enhanced due to the effects of ram pressure. In order to make the transition to cluster S0s these galaxies must lose ~ 40% of their specific angular momentum. We suggest this must occur > 1 Gyr later, after the molecular gas has been depleted and/or stripped, via multiple tidal interactions with other cluster members.

Recent studies have found that obscured quasars cluster more strongly and are thus hosted by dark matter haloes of larger mass than their unobscured counterparts. These results pose a challenge for the simplest unification models, in which obscured objects are intrinsically the same as unobscured sources but seen through a dusty line of sight. There is general consensus that a structure like a "dusty torus" exists, meaning that this intrinsic similarity is likely the case for at least some subset of obscured quasars. However, the larger host halo masses of obscured quasars implies that there is a second obscured population that has an even higher clustering amplitude and typical halo mass. Here, we use simple assumptions about the host halo mass distributions of quasars, along with analytical methods and cosmological $N$-body simulations to isolate the signal from this population. We provide values for the bias and halo mass as a function of the fraction of the "non-torus obscured" population. Adopting a reasonable value for this fraction of $\sim$25% implies a non-torus obscured quasar bias that is much higher than the observed obscured quasar bias, because a large fraction of the obscured population shares the same clustering strength as the unobscured objects. For this non-torus obscured population, we derive a bias of $\sim$3, and typical halo masses of $\sim3\times10^{13}$ M$_{\odot}/h$ at $z=1$. These massive haloes are likely the descendants of high-mass unobscured quasars at high redshift, and will evolve into members of galaxy groups at $z=0$.

The future space missions Euclid and WFIRST-AFTA will use the Halpha emission line to measure the redshifts of tens of millions of galaxies. The Halpha luminosity function at z>0.7 is one of the major sources of uncertainty in forecasting cosmological constraints from these missions. We construct unified empirical models of the Halpha luminosity function spanning the range of redshifts and line luminosities relevant to the redshift surveys proposed with Euclid and WFIRST-AFTA. By fitting to observed luminosity functions from Halpha surveys, we build three models for its evolution. Different fitting methodologies, functional forms for the luminosity function, subsets of the empirical input data, and treatment of systematic errors are considered to explore the robustness of the results. Functional forms and model parameters are provided for all three models, along with the counts and redshift distributions up to z~2.5 for a range of limiting fluxes (F_Halpha>0.5 - 3 x 10^-16 erg cm^-2 s^-1) that are relevant for future space missions. For instance, in the redshift range 0.90<z<1.8, our models predict an available galaxy density in the range 7700--13300 and 2000--4800 deg^-2 respectively at fluxes above F_Halpha>1 and 2 x 10^-16 erg cm^-2 s^-1, and 32000--48000 for F_Halpha>0.5 x 10^-16 erg cm^-2 s^-1 in the extended redshift range 0.40<z<1.8. We also consider the implications of our empirical models for the total Halpha luminosity density of the Universe, and the closely related cosmic star formation history.

We report the discovery of a gravitationally lensed hyperluminous infrared galaxy (L_IR~10^13 L_sun) with strong radio emission (L_1.4GHz~10^25 W/Hz) at z=2.553. The source was identified in the citizen science project SpaceWarps through the visual inspection of tens of thousands of iJKs colour composite images of Luminous Red Galaxies (LRGs), groups and clusters of galaxies and quasars. Appearing as a partial Einstein ring (r_e~3") around an LRG at z=0.2, the galaxy is extremely bright in the sub-millimetre for a cosmological source, with the thermal dust emission approaching 1 Jy at peak. The redshift of the lensed galaxy is determined through the detection of the CO(3-2) molecular emission line with the Large Millimetre Telescope's Redshift Search Receiver and through [OIII] and H-alpha line detections in the near-infrared from Subaru/IRCS. We have resolved the radio emission with high resolution (300-400 mas) eMERLIN L-band and JVLA C-band imaging. These observations are used in combination with the near-infrared imaging to construct a lens model, which indicates a lensing magnification of ~10x. The source reconstruction appears to support a radio morphology comprised of a compact (<250 pc) core and more extended component, perhaps indicative of an active nucleus and jet or lobe.

We present multiwavelength identifications for the counterparts of 1088 submillimeter sources detected at 850$\mu$m in the SCUBA-2 Cosmology Legacy Survey study of the UKIDSS-UDS field. By utilising an ALMA pilot study on a subset of our bright SCUBA-2 sample as a training set, along with the deep optical-near-infrared data available in this field, we develop a novel technique, Optical-IR Triple Color (OIRTC), using $z-K$, $K-[3.6]$, $[3.6]-[4.5]$ colors to select the candidate submillimeter galaxy (SMG) counterparts. By combining radio identification and the OIRTC technique, we find counterpart candidates for 80% of the Class = 1 $\geq4\,\sigma$ SCUBA-2 sample, defined as those that are covered by both radio and OIR imaging and the base sample for our scientific analyses. Based on the ALMA training set, we expect the accuracy of these identifications to be $82\pm20$%, with a completeness of $69\pm16$%, essentially as accurate as the traditional $p$-value technique but with higher completeness. We find that the fraction of SCUBA-2 sources having candidate counterparts is lower for fainter 850$\mu$m sources, and we argue that for follow-up observations sensitive to SMGs with $S_{850}\gtrsim 1$ mJy across the whole ALMA beam, the fraction with multiple counterparts is likely to be $>40$% for SCUBA-2 sources at $S_{850} \gtrsim 4$ mJy. We find that the photometric redshift distribution for the SMGs is well fit by a lognormal distribution, with a median redshift of $z=2.3\pm0.1$. After accounting for the sources without any radio and/or OIRTC counterpart, we estimate the median redshift to be $z=2.6\pm0.1$ for SMGs with $S_{850} >1$ mJy. We also use this new large sample to study the clustering of SMGs and the the far-infrared properties of the unidentified submillimeter sources by stacking their Herschel SPIRE far-infrared emission.

It is thought that a class of broad absorption line (BAL) QSOs, characterised by Fe absorption features in their UV spectra (called `FeLoBALs'), could mark a transition stage between the end of an obscured starburst event and a youthful QSO beginning to shed its dust cocoon, where Fe has been injected into the interstellar medium by the starburst. To test this hypothesis we have undertaken deep SCUBA-2 850 $\mu$m observations of a sample of 17 FeLoBAL QSOs with 0.89 $\leq$ z $\leq$ 2.78 and -23.31 $\leq$ M$_{B}$ $\leq$-28.50 to directly detect an excess in the thermal emission of the dust which would probe enhanced star-formation activity. We find that FeLoBALs are not luminous sources in the submillimetre, none of them are individually detected at 850 $\mu$m, nor as a population through stacking ($F_{s}=1.14\pm0.58$ mJy). Statistical and survival analyses reveal that FeLoBALs have sub-mm properties consistent with BAL and non-BAL QSOs with matched redshifts and magnitudes. An SED fitting analysis shows that the FIR emission is dominated by AGN activity, and a starburst component is required only in 6/17 sources of our sample; moreover the integrated total luminosity of 16/17 sources is L$\geq$10$^{12}$L$_{\odot}$, high enough to classify FeLoBALs as infrared luminous. In conclusion, we do not find any evidence in support of FeLoBAL QSOs being a transition population between a ULIRG and an unobscured QSO; in particular, FeLoBALs are not characterized by a cold starburst which would support this hypothesis.

We present millimetre (SMA) and sub-millimetre (SCUBA-2) continuum observations of the peculiar star KIC 8462852 which displayed several deep and aperiodic dips in brightness during the Kepler mission. Our observations are approximately confusion-limited at 850 $\mu$m and are the deepest millimetre and sub-millimetre photometry of the star that has yet been carried out. No significant emission is detected towards KIC 8462852. We determine upper limits for dust between a few 10$^{-6}$ M$_{\oplus}$ and 10$^{-3}$ M$_{\oplus}$ for regions identified as the most likely to host occluding dust clumps and a total overall dust budget of $<$7.7 M$_{\oplus}$ within a radius of 200 AU. Such low limits for the inner system make the catastrophic planetary disruption hypothesis unlikely. Integrating over the Kepler lightcurve we determine that at least 10$^{-9}$ M$_{\oplus}$ of dust is required to cause the observed Q16 dip. This is consistent with the currently most favoured cometary breakup hypothesis, but nevertheless implies the complete breakup of $\sim$ 30 Comet 1/P Halley type objects. Finally, in the wide SCUBA-2 field-of-view we identify another candidate debris disc system that is potentially the largest yet discovered.

Radio emission from radio-quiet quasars may be due to star formation in the quasar host galaxy, to a jet launched by the supermassive black hole, or to relativistic particles accelerated in a wide-angle radiatively-driven outflow. In this paper we examine whether radio emission from radio-quiet quasars is a byproduct of star formation in their hosts. To this end we use infrared spectroscopy and photometry from Spitzer and Herschel to estimate or place upper limits on star formation rates in hosts of ~300 obscured and unobscured quasars at z<1. We find that low-ionization forbidden emission lines such as [NeII] and [NeIII] are likely dominated by quasar ionization and do not provide reliable star formation diagnostics in quasar hosts, while PAH emission features may be suppressed due to the destruction of PAH molecules by the quasar radiation field. While the bolometric luminosities of our sources are dominated by the quasars, the 160 micron fluxes are likely dominated by star formation, but they too should be used with caution. We estimate median star formation rates to be 6-29 Msun/year, with obscured quasars at the high end of this range. This star formation rate is insufficient to explain the observed radio emission from quasars by an order of magnitude, with log(L_radio, observed/L_radio, SF)=0.6-1.3 depending on quasar type and star formation estimator. Although radio-quiet quasars in our sample lie close to the 8-1000 micron infrared / radio correlation characteristic of the star-forming galaxies, both their infrared emission and their radio emission are dominated by the quasar activity, not by the host galaxy.

The over-dense environments of protoclusters of galaxies in the early Universe (z>2) are expected to accelerate the evolution of galaxies, with an increased rate of stellar mass assembly and black hole accretion compared to co-eval galaxies in the average density `field'. These galaxies are destined to form the passive population of massive systems that dominate the cores of rich clusters today. While signatures of accelerated growth of galaxies in the SSA22 protocluster (z=3.1) have been observed, the mechanism driving this remains unclear. Here we show an enhanced rate of galaxy-galaxy mergers could be responsible. We morphologically classify Lyman-break Galaxies (LBGs) in the SSA22 protocluster and compare these to those of galaxies in the field at z=3.1 as either active mergers or non-merging using Hubble Space Telescope ACS/F814W imaging, probing the rest-frame ultraviolet stellar light. We measure a merger fraction of 48+/-10 per cent for LBGs in the protocluster compared to 30+/-6 per cent for the field. Although the excess is marginal, an enhanced rate of mergers in SSA22 hints that galaxy-galaxy mergers are one of the key channels driving accelerated star formation and AGN growth in protocluster environments.

We report the results of $1^{\prime}.5 \times3^{\prime}$ mapping at 1.1~mm with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the central region of the $z=3.09$ SSA22 protocluster. By combining our source catalog with archival spectroscopic redshifts, we find that eight submillimeter galaxies (SMGs) with flux densities, $S_{\rm 1.1~mm}=0.7-6.4$~mJy ($L_{\rm IR}\sim10^{12.1}-10^{13.1}L_\odot$) are at $z=3.08-3.10$. Not only are these SMGs members of the protocluster but they in fact reside within the node at the junction of the 50 Mpc-scale filamentary three-dimensional structure traced by Lyman-$\alpha$ emitters (LAEs) in this field. The eight SMGs account for a star formation rate density (SFRD) $\sim$10 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$ in the node, which is two orders of magnitudes higher than the global SFRD at this redshift. We find that four of the eight SMGs host a X-ray luminous active galactic nuclei (AGN). Our results suggest that the vigorous star formation activity and the growth of super massive black holes (SMBHs) occurred simultaneously in the densest regions at $z\sim3$, which may correspond to the most active historical phase of the massive galaxy population found in the core of the clusters in the present universe. Two SMGs are associated with Lyman-$\alpha$ blobs (LABs), implying that the two populations coexist in high density environments for a few cases.

We present VLT VIMOS, Keck DEIMOS and Keck LRIS multi-object spectra of 367 sources in the field of the z ~ 3.09 protocluster SSA22. Sources are spectroscopically classified via template matching, allowing new identifications for 206 extragalactic sources, including 36 z > 2 Lyman-break galaxies (LBGs) and Lyman-\alpha\ emitters (LAEs), 8 protocluster members, and 94 X-ray sources from the ~ 400 ks Chandra deep survey of SSA22. Additionally, in the area covered by our study, we have increased by ~ 4, 13, and 6 times the number of reliable redshifts of sources at 1.0 < z < 2.0, at z > 3.4, and with X-Ray emission, respectively. We compare our results with past spectroscopic surveys of SSA22 to investigate the completeness of the LBGs and the X-Ray properties of the new spectroscopically-classified sources in the SSA22 field.

We study the 850um emission in X-ray selected AGN in the 2 sq-deg COSMOS field using new data from the SCUBA-2 Cosmology Legacy Survey. We find 19 850um bright X-ray AGN in a high-sensitivity region covering 0.89 sq-deg with flux densities of S850=4-10 mJy. The 19 AGN span the full range in redshift and hard X-ray luminosity covered by the sample - 0.7<z<3.5 and 43.2<log10(LX) <45. We report a highly significant stacked 850um detection of a hard X-ray flux-limited population of 699 z>1 X-ray AGN - S850=0.71+/-0.08mJy. We explore trends in the stacked 850um flux densities with redshift, finding no evolution in the average cold dust emission over the redshift range probed. For Type 1 AGN, there is no significant correlation between the stacked 850um flux and hard X-ray luminosity. However, in Type 2 AGN the stacked submm flux is a factor of 2 higher at high luminosities. When averaging over all X-ray luminosities, no significant differences are found in the stacked submm fluxes of Type 1 and Type 2 AGN as well as AGN separated on the basis of X-ray hardness ratios and optical-to-infrared colours. However, at log10(LX) >44.4, dependences in average submm flux on the optical-to-infrared colours become more pronounced. We argue that these high luminosity AGN represent a transition from a secular to a merger-driven evolutionary phase where the star formation rates and accretion luminosities are more tightly coupled. Stacked AGN 850um fluxes are compared to the stacked fluxes of a mass-matched sample of K-band selected non-AGN galaxies. We find that at 10.5<log10(M*/M0)<11.5, the non-AGN 850um fluxes are 1.5-2x higher than in Type 2 AGN of equivalent mass. We suggest these differences are due to the presence of massive dusty, red starburst galaxies in the K-band selected non-AGN sample, which are not present in optically selected catalogues covering a smaller area.

We present the results of a MIPS-24um study of the Brightest Cluster Galaxies (BCGs) of 535 high-redshift galaxy clusters. The clusters are drawn from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS), which effectively provides a sample selected on total stellar mass, over 0.2 < z < 1.8 within the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Survey fields. 20%, or 106 clusters have spectroscopically confirmed redshifts, and the rest have redshifts estimated from the color of their red sequence. A comparison with the public SWIRE images detects 125 individual BCGs at 24um > 100uJy, or 23%. The luminosity-limited detection rate of BCGs in similar richness clusters (Ngal> 12) increases rapidly with redshift. Above z ~ 1, an average of ~20\% of the sample have 24um-inferred infrared luminosities of LIR > 10^12 Lsun, while the fraction below z ~ 1 exhibiting such luminosities is < 1 \%. The Spitzer-IRAC colors indicate the bulk of the 24um-detected population is predominantly powered by star formation, with only 7/125 galaxies lying within the color region inhabited by Active Galactic Nuclei (AGN). Simple arguments limit the star-formation activity to several hundred million years and this may therefore be indicative of the timescale for AGN feedback to halt the star formation. Below redshift z ~ 1 there is not enough star formation to significantly contribute to the overall stellar mass of the BCG population, and therefore BCG growth is likely dominated by dry-mergers. Above z~ 1, however, the inferred star formation would double the stellar mass of the BCGs and is comparable to the mass assembly predicted by simulations through dry mergers. We cannot yet constrain the process driving the star formation for the overall sample, though a single object studied in detail is consistent with a gas-rich merger.

We have discovered an optically rich galaxy cluster at z=1.7089 with star formation occurring in close proximity to the central galaxy. The system, SpARCS104922.6+564032.5, was detected within the Spitzer Adaptation of the red-sequence Cluster Survey, (SpARCS), and confirmed through Keck-MOSFIRE spectroscopy. The rest-frame optical richness of Ngal(500kpc) = 30+/-8 implies a total halo mass, within 500kpc, of ~3.8+/-1.2 x 10^14 Msun, comparable to other clusters at or above this redshift. There is a wealth of ancillary data available, including Canada-France-Hawaii Telescope optical, UKIRT-K, Spitzer-IRAC/MIPS, and Herschel-SPIRE. This work adds submillimeter imaging with the SCUBA2 camera on the James Clerk Maxwell Telescope and near-infrared imaging with the Hubble Space Telescope (HST). The mid/far-infrared (M/FIR) data detect an Ultra-luminous Infrared Galaxy spatially coincident with the central galaxy, with LIR = 6.2+/-0.9 x 10^12 Lsun. The detection of polycyclic aromatic hydrocarbons (PAHs) at z=1.7 in a Spitzer-IRS spectrum of the source implies the FIR luminosity is dominated by star formation (an Active Galactic Nucleus contribution of 20%) with a rate of ~860+/-30 Msun/yr. The optical source corresponding to the IR emission is likely a chain of of > 10 individual clumps arranged as "beads on a string" over a linear scale of 66 kpc. Its morphology and proximity to the Brightest Cluster Galaxy imply a gas-rich interaction at the center of the cluster triggered the star formation. This system indicates that wet mergers may be an important process in forming the stellar mass of BCGs at early times.

We present Early Science observations with the Large Millimeter Telescope, AzTEC 1.1 mm continuum images and wide bandwidth spectra (73-111 GHz) acquired with the Redshift Search Receiver, towards four bright lensed submillimetre galaxies identified through the Herschel Lensing Survey-snapshot and the SCUBA-2 Cluster Snapshot Survey. This pilot project studies the star formation history and the physical properties of the molecular gas and dust content of the highest redshift galaxies identified through the benefits of gravitational magnification. We robustly detect dust continuum emission for the full sample and CO emission lines for three of the targets. We find that one source shows spectroscopic multiplicity and is a blend of three galaxies at different redshifts (z=2.040, 3.252 and 4.680), reminiscent of previous high-resolution imaging follow-up of unlensed submillimetre galaxies, but with a completely different search method, that confirm recent theoretical predictions of physically unassociated blended galaxies. Identifying the detected lines as 12CO (J_up=2-5) we derive spectroscopic redshifts, molecular gas masses, and dust masses from the continuum emission. The mean H_2 gas mass of the full sample is (2.0 +- 0.2) x 10^11 M_sun/\mu, and the mean dust mass is (2.0+-0.2) x 10^9 M_sun/\mu, where \mu=2-5 is the expected lens amplification. Using these independent estimations we infer a gas-to-dust ratio of \delta_GDR=55-75, in agreement with other measurements of submillimetre galaxies. Our magnified high-luminosity galaxies fall on the same locus as other high-redshift submillimetre galaxies, extending the L'_CO - L_FIR correlation observed for local luminous and ultraluminous infrared galaxies to higher FIR and CO luminosities.

We consider the high radio frequency (15 GHz - 353 GHz) properties and variability of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing us to relate our results to the general population. We find that >6.0% of our parent sample (>15.1% if only cool-core clusters are considered) contain a radio-source at 150 GHz of at least 3mJy (~1x10^23 W/Hz at our median redshift of z~0.13). Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component (Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increasing to >8.5% if only cool-core clusters are considered. We see little evidence for strong variability at 15 GHz on short (week-month) time-scales although we see variations greater than 20% at 150 GHz over 6-month times-frames for 4 of the 23 sources with multi-epoch observations. Much more prevalent is long-term (year-decade time-scale) variability, with average annual amplitude variations greater than 1% at 15 GHz being commonplace. There is a weak trend towards higher variability as the peak of the GPS-like component occurs at higher frequency. We demonstrate the complexity that is seen in the radio spectra of BCGs and discuss the potentially significant implications of these high-peaking components for Sunyaev-Zel'dovich cluster searches.

We report the redshift of an unlensed, highly obscured submillimetre galaxy (SMG), HS1700.850.1, the brightest SMG (S850um =19.1 mJy) detected in the JCMT/SCUBA-2 Baryonic Structure Survey, based on the detection of its CO line emission. Using the IRAM PdBI-WIDEX with 3.6GHz band width, we serendipitously detect an emission line at 150.6 GHz. From a search over 14.5 GHz in the 3-mm and 2-mm atmospheric windows, we confirm the identification of this line as CO(5-4) at z = 2.816, meaning that it does not reside in the z~2.30 proto-cluster in this field. Measurement of the 870um source size (<0.85") from the Sub-Millimeter Array (SMA) confirms a compact emission in a S870um =14.5mJy, LIR~10^13 Lsun component, suggesting an Eddington-limited starburst. We use the double-peaked CO line profile measurements along with the SMA size constraints to study the gas dynamics of a HyLIRG, estimating the gas and dynamical masses of HS1700.850.1. While HS1700.850.1 is one of the most extreme galaxies known in the Universe, we find that it occupies a relative void in the Lyman-Break Galaxy distribution in this field. Comparison with other extreme objects at similar epochs (HyLIRG Quasars), and cosmological simulations, suggests such an anti-bias of bright SMGs could be relatively common, with the brightest SMGs rarely occupying the most overdense regions at z=2-4.

We present a novel unsupervised learning approach to automatically segment and label images in astronomical surveys. Automation of this procedure will be essential as next-generation surveys enter the petabyte scale: data volumes will exceed the capability of even large crowd-sourced analyses. We demonstrate how a growing neural gas (GNG) can be used to encode the feature space of imaging data. When coupled with a technique called hierarchical clustering, imaging data can be automatically segmented and labelled by organising nodes in the GNG. The key distinction of unsupervised learning is that these labels need not be known prior to training, rather they are determined by the algorithm itself. Importantly, after training a network can be be presented with images it has never 'seen' before and provide consistent categorisation of features. As a proof-of-concept we demonstrate application on data from the Hubble Space Telescope Frontier Fields: images of clusters of galaxies containing a mixture of galaxy types that would easily be recognised and classified by a human inspector. By training the algorithm using one field (Abell 2744) and applying the result to another (MACS0416.1-2403), we show how the algorithm can cleanly separate image features that a human would associate with early and late type galaxies. We suggest that the algorithm has potential as a tool in the automatic analysis and data mining of next-generation imaging and spectral surveys, and could also find application beyond astronomy.

We present high-resolution 870-um ALMA continuum maps of 30 bright sub-millimeter sources in the UKIDSS UDS field. These sources are selected from deep, 1-square degrees 850-um maps from the SCUBA--2 Cosmology Legacy Survey, and are representative of the brightest sources in the field (median SCUBA2 flux S_850=8.7+/-0.4 mJy). We detect 52 sub-millimeter galaxies (SMGs) at >4-sigma significance in our 30 ALMA maps. In 61+/-17% of the ALMA maps the single-dish source comprises a blend of >=2 SMGs, where the secondary SMGs are Ultra--Luminous Infrared Galaxies (ULIRGs) with L_IR>10^12 Lo. The brightest SMG contributes on average 80+/-4% of the single-dish flux density, and in the ALMA maps containing >=2 SMGs the secondary SMG contributes 25+/-3% of the integrated ALMA flux. We construct source counts and show that multiplicity boosts the apparent single-dish cumulative counts by 20% at S_870>7.5mJy, and by 60% at S_870>12mJy. We combine our sample with previous ALMA studies of fainter SMGs and show that the counts are well-described by a double power-law with a break at 8.5+/-0.6mJy. The break corresponds to a luminosity of ~6x10^12Lsol or a star-formation rate of ~1000Mo/yr. For the typical sizes of these SMGs, which are resolved in our ALMA data with r=1.2+/-0.1kpc, this yields a limiting SFR density of ~100Msol/yr/kpc2. Finally, the number density of S_870>2mJy SMGs is 80+/-30 times higher than that derived from blank-field counts. An over-abundance of faint SMGs is inconsistent with line-of-sight projections dominating multiplicity in the brightest SMGs, and indicates that a significant proportion of these high-redshift ULIRGs must be physically associated.

We report the blind detection of 12CO emission from a Distant Red Galaxy, HS1700.DRG55. We have used the IRAM PdBI-WIDEX, with its 3.6GHz of instantaneous dual-polarization bandwidth, to target 12CO(3--2) from galaxies lying in the proto-cluster at z=2.300 in the field HS1700+64. If indeed this line in DRG55 is 12CO(3--2), it's detection at 104.9GHz indicates a z_CO=2.296. None of the other eight known z~2.30 proto-cluster galaxies lying within the primary beam (PB) are detected in 12CO, although the limits are ~2x worse towards the edge of the PB where several lie. The optical/near-IR magnitudes of DRG55 (R_AB>27, K_AB=22.3) mean that optical spectroscopic redshifts are difficult with 10m-class telescopes, but near-IR redshifts would be feasible. The 24um-implied SFR (210 M_odot yr-1), stellar mass (~10^11 M-odot) and 12CO line luminosity (3.6x10^10 K km s-1 pc^2) are comparable to other normal 12CO-detected star forming galaxies in the literature, although the galaxy is some ~2 mag (~6x) fainter in the rest-frame UV than 12CO-detected galaxies at z>2. The detection of DRG55 in 12CO complements three other 12CO detected UV-bright galaxies in this proto-cluster from previous studies, and suggests that many optically faint galaxies in the proto-cluster may host substantial molecular gas reservoirs, and a full blind census of 12CO in this overdense environment is warranted.

The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.

Recent observations have revealed that starburst galaxies can drive molecular gas outflows through stellar radiation pressure. Molecular gas is the phase of the interstellar medium from which stars form, so these outflows curtail stellar mass growth in galaxies. Previously known outflows, however, involve small fractions of the total molecular gas content and are restricted to sub-kiloparsec scales. It is also apparent that input from active galactic nuclei is in at least some cases dynamically important, so pure stellar feedback has been considered incapable of aggressively terminating star formation on galactic scales. Extraplanar molecular gas has been detected in the archetype starburst galaxy M82, but so far there has been no evidence that starbursts can propel significant quantities of cold molecular gas to the same galactocentric radius (~10 kpc) as the warmer gas traced by metal absorbers. Here we report observations of molecular gas in a compact (effective radius 100 pc) massive starburst galaxy at z=0.7, which is known to drive a fast outflow of ionized gas. We find that 35 per cent of the total molecular gas is spatially extended on a scale of approximately 10 kpc, and one third of this has a velocity of up to 1000 km/s. The kinetic energy associated with this high-velocity component is consistent with the momentum flux available from stellar radiation pressure. This result demonstrates that nuclear bursts of star formation are capable of ejecting large amounts of cold gas from the central regions of galaxies, thereby strongly affecting their evolution.

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

We cross-correlate a cosmic microwave background (CMB) lensing map with the projected space densities of quasars to measure the bias and halo masses of a quasar sample split into obscured and unobscured populations, the first application of this method to distinct quasar subclasses. Several recent studies of the angular clustering of obscured quasars have shown that these objects likely reside in higher-mass halos compared to their unobscured counterparts. This has important implications for models of the structure and geometry of quasars, their role in growing supermassive black holes, and mutual quasar/host galaxy evolution. However, the magnitude and significance of this difference has varied from study to study. Using data from \planck, \wise, and SDSS, we follow up on these results using the independent method of CMB lensing cross-correlations. The region and sample are identical to that used for recent angular clustering measurements, allowing for a direct comparison of the CMB-lensing and angular clustering methods. At $z \sim 1$, we find that the bias of obscured quasars is $b_q = 2.57 \pm 0.24$, while that of unobscured quasars is $b_q = 1.89 \pm 0.19$. This corresponds to halo masses of $\log (M_h / M_{\odot} h^{-1}) = 13.24_{-0.15}^{+0.14}$ (obscured) and $\log (M_h / M_{\odot} h^{-1}) = 12.71_{-0.13}^{+0.15}$ (unobscured). These results agree well with with those from angular clustering (well within $1\sigma$), and confirm that obscured quasars reside in host halos $\sim$3 times as massive as halos hosting unobscured quasars. This implies that quasars spend a significant portion of their lifetime in an obscured state, possibly more than one half of the entire active phase.

We describe a search for submillimeter emission in the vicinity of one of the most distant, luminous galaxies known, HerMES FLS3 at z=6.34, exploiting it as a signpost to a potentially biased region of the early Universe, as might be expected in hierarchical structure formation models. Imaging to the confusion limit with the innovative, wide-field submillimeter bolometer camera, SCUBA-2, we are sensitive to colder and/or less luminous galaxies in the surroundings of HFLS3. We use the Millennium Simulation to illustrate that HFLS3 may be expected to have companions if it is as massive as claimed, but find no significant evidence from the surface density of SCUBA-2 galaxies in its vicinity, or their colors, that HFLS3 marks an over-density of dusty, star-forming galaxies. We cannot rule out the presence of dusty neighbours with confidence, but deeper 450-um imaging has the potential to more tightly constrain the redshifts of nearby galaxies, at least one of which likely lies at z>~5. If associations with HFLS3 can be ruled out, this could be taken as evidence that HFLS3 is less biased than a simple extrapolation of the Millennium Simulation may imply. This could suggest either that it represents a rare short-lived, but highly luminous, phase in the evolution of an otherwise typical galaxy, or that this system has suffered amplification due to a foreground gravitational lens and so is not as intrinsically luminous as claimed.

We present statistically significant detections at 850um of the Lyman Break Galaxy (LBG) population at z=3, 4, and 5 using data from the Submillimetre Common User Bolometer Array 2 (SCUBA-2) Cosmology Legacy Survey (S2CLS) in the United Kingdom Infrared Deep Sky Survey Ultra Deep Survey (UKIDSS-UDS) field. We employ a stacking technique to probe beneath the survey limit to measure the average 850um flux density of LBGs at z=3, 4, and 5 with typical ultraviolet luminosities of L(1700A)~10^29 erg/s/Hz. We measure 850um flux densities of (0.25 +/- 0.03, (0.41 +/- 0.06), and (0.88 +/- 0.23) mJy respectively, and find that they contribute at most 20 per cent to the cosmic far-infrared background at 850um. Fitting an appropriate range of spectral energy distributions to the z=3, 4, and 5 LBG stacked 24-850um fluxes, we derive infrared (IR) luminosities of L(8-1000um)~3.2, 5.5, and 11.0x10^11 Lsun (corresponding to star formation rates of ~50-200 Msun/yr) respectively. We find that the evolution in the IR luminosity density of LBGs is broadly consistent with model predictions for the expected contribution of luminous IR galaxy (LIRG) to ultraluminous IR galaxy (ULIRG) type systems at these epochs. We also see a strong positive correlation between stellar mass and IR luminosity. Our data are consistent with the main sequence of star formation showing little or no evolution from z=3 to 5. We have also confirmed that, for a fixed mass, the reddest LBGs (UV slope Beta -> 0) are indeed redder due to dust extinction, with SFR(IR)/SFR(UV) increasing by approximately an order of magnitude over -2<Beta<0 such that SFR(IR)/SFR(UV)~20 for the reddest LBGs. Furthermore, the most massive LBGs also tend to have higher obscured-to-unobscured ratio, hinting at a variation in the obscuration properties across the mass range.

We have observed a sample of typical z=1 star forming galaxies, selected from the HiZELS survey, with the new KMOS near-infrared, multi-IFU instrument on the VLT, in order to obtain their dynamics and metallicity gradients. The majority of our galaxies have a metallicity gradient consistent with being flat or negative (i.e. higher metallicity cores than outskirts). Intriguingly, we find a trend between metallicity gradient and specific star formation rate (sSFR), such that galaxies with a high sSFR tend to have relatively metal-poor centres, a result which is strengthened when combined with datasets from the literature. This result appears to explain the discrepancies reported between different high redshift studies and varying claims for evolution. From a galaxy evolution perspective, the trend we see would mean that a galaxy's sSFR is governed by the amount of metal poor gas that can be funnelled into its core, triggered either by merging or through efficient accretion. In fact merging may play a significant role as it is the starburst galaxies at all epochs, which have the more positive metallicity gradients. Our results may help to explain the origin of the fundamental metallicity relation, in which galaxies at a fixed mass are observed to have lower metallicities at higher star formation rates, especially if the metallicity is measured in an aperture encompassing only the central regions of the galaxy. Finally, we note that this study demonstrates the power of KMOS as an efficient instrument for large scale resolved galaxy surveys.

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

Recent studies of luminous infrared-selected active galactic nuclei (AGN) suggest that the reddest, most obscured objects display a higher angular clustering amplitude, and thus reside in higher-mass dark matter halos. This is a direct contradiction to the prediction of the simplest unification-by-orientation models of AGN and quasars. However, clustering measurements depend strongly on the "mask" that removes low-quality data and describes the sky and selection function. We find that applying a robust, conservative mask to WISE-selected quasars yields a weaker but still significant difference in the bias between obscured and unobscured quasars. These findings are consistent with results from previous Spitzer surveys, and removes any scale dependence of the bias. For obscured quasars with $\langle z \rangle = 0.99$ we measure a bias of $b_q = 2.67 \pm 0.16$, corresponding to a halo mass of $\log (M_h / M_{\odot} h^{-1}) = 13.3 \pm 0.1$, while for unobscured sources with $\langle z \rangle = 1.04$ we find $b_q = 2.04 \pm 0.17$ with a halo mass $\log (M_h / M_{\odot} h^{-1} )= 12.8 \pm 0.1$. This improved measurement indicates that WISE-selected obscured quasars reside in halos only a few times more massive than the halos of their unobscured counterparts, a reduction in the factor of $\sim$10 larger halo mass as has been previously reported using WISE-selected samples. Additionally, an abundance matching analysis yields lifetimes for both obscured and unobscured quasar phases on the order of a few 100 Myr ($\sim$ 1\% of the Hubble time) --- however, the obscured phase lasts roughly twice as long, in tension with many model predictions.

We report Plateau de Bure Interferometer (PdBI) 1.1 mm continuum imaging towards two extremely red H-[4.5]>4 (AB) galaxies at z>3, which we have previously discovered making use of Spitzer SEDS and Hubble Space Telescope CANDELS ultra-deep images of the Ultra Deep Survey field. One of our objects is detected on the PdBI map with a 4.3 sigma significance, corresponding to Snu(1.1mm)=(0.78 +/- 0.18) mJy. By combining this detection with the Spitzer 8 and 24 micron photometry for this source, and SCUBA2 flux density upper limits, we infer that this galaxy is a composite active galactic nucleus/star-forming system. The infrared (IR)-derived star formation rate is SFR~(200 +/- 100) Msun/yr, which implies that this galaxy is a higher-redshift analogue of the ordinary ultra-luminous infrared galaxies (ULIRGs) more commonly found at z~2-3. In the field of the other target, we find a tentative 3.1 sigma detection on the PdBI 1.1 mm map, but 3.7 arcsec away of our target position, so it likely corresponds to a different object. In spite of the lower significance, the PdBI detection is supported by a close SCUBA2 3.3 sigma detection. No counterpart is found on either the deep SEDS or CANDELS maps, so, if real, the PdBI source could be similar in nature to the sub-millimetre source GN10. We conclude that the analysis of ultra-deep near- and mid-IR images offers an efficient, alternative route to discover new sites of powerful star formation activity at high redshifts.

We investigate the process of rapid star formation quenching in a sample of 12 massive galaxies at intermediate redshift (z~0.6) that host high-velocity ionized gas outflows (v>1000 km/s). We conclude that these fast outflows are most likely driven by feedback from star formation rather than active galactic nuclei (AGN). We use multiwavelength survey and targeted observations of the galaxies to assess their star formation, AGN activity, and morphology. Common attributes include diffuse tidal features indicative of recent mergers accompanied by bright, unresolved cores with effective radii less than a few hundred parsecs. The galaxies are extraordinarily compact for their stellar mass, even when compared with galaxies at z~2-3. For 9/12 galaxies, we rule out an AGN contribution to the nuclear light and hypothesize that the unresolved core comes from a compact central starburst triggered by the dissipative collapse of very gas-rich progenitor merging disks. We find evidence of AGN activity in half the sample but we argue that it accounts for only a small fraction (<10%) of the total bolometric luminosity. We find no correlation between AGN activity and outflow velocity and we conclude that the fast outflows in our galaxies are not powered by on-going AGN activity, but rather by recent, extremely compact starbursts.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

We analyse new SCUBA-2 submillimeter and archival SPIRE far-infrared imaging of a z=1.62 cluster, Cl0218.3-0510, which lies in the UKIDSS/UDS field of the SCUBA-2 Cosmology Legacy Survey. Combining these tracers of obscured star formation activity with the extensive photometric and spectroscopic information available for this field, we identify 31 far-infrared/submillimeter-detected probable cluster members with bolometric luminosities >1e12 Lo and show that by virtue of their dust content and activity, these represent some of the reddest and brightest galaxies in this structure. We exploit Cycle-1 ALMA submillimeter continuum imaging which covers one of these sources to confirm the identification of a SCUBA-2-detected ultraluminous star-forming galaxy in this structure. Integrating the total star-formation activity in the central region of the structure, we estimate that it is an order of magnitude higher (in a mass-normalised sense) than clusters at z~0.5-1. However, we also find that the most active cluster members do not reside in the densest regions of the structure, which instead host a population of passive and massive, red galaxies. We suggest that while the passive and active populations have comparable near-infrared luminosities at z=1.6, M(H)~-23, the subsequent stronger fading of the more active galaxies means that they will evolve into passive systems at the present-day which are less luminous than the descendants of those galaxies which were already passive at z~1.6 (M(H)~-20.5 and M(H)~-21.5 respectively at z~0). We conclude that the massive galaxy population in the dense cores of present-day clusters were already in place at z=1.6 and that in Cl0218.3-0510 we are seeing continuing infall of less extreme, but still ultraluminous, star-forming galaxies onto a pre-existing structure.

We present the spatially resolved H-alpha (Ha) dynamics of sixteen star-forming galaxies at z~0.81 using the new KMOS multi-object integral field spectrograph on the ESO VLT. These galaxies were selected using 1.18 um narrow-band imaging from the 10 deg^2 CFHT-HiZELS survey of the SA22hr field, are found in a ~4Mpc over-density of Ha emitters and likely reside in a group/intermediate environment, but not a cluster. We confirm and identify a rich group of star-forming galaxies at z=0.813+-0.003, with thirteen galaxies within 1000 km/s of each other, and 7 within a diameter of 3Mpc. All our galaxies are "typical" star-forming galaxies at their redshift, 0.8+-0.4 SFR*(z=0.8), spanning a range of specific star formation rate of sSFR=0.2-1.1 Gyr^-1 and have a median metallicity very close to solar of 12+log(O/H)=8.62+-0.06. We measure the spatially resolved Ha dynamics of the galaxies in our sample and show that thirteen out of sixteen galaxies can be described by rotating disks and use the data to derive inclination corrected rotation speeds of 50-275 km/s. The fraction of disks within our sample is 75+-8, consistent with previous results based on HST morphologies of Ha selected galaxies at z~1 and confirming that disks dominate the star formation rate density at z~1. Our Ha galaxies are well fitted by the z~1-2 Tully-Fisher relation, confirming the evolution seen in the zero-point. Apart from having, on average, higher stellar masses and lower sSFRs, our group galaxies at z=0.813 present the same mass-metallicity and TF relation as z~1 field galaxies, and are all disk galaxies.

We measure the cross-power spectrum of the projected mass density as traced by the convergence of the cosmic microwave background lensing field from the South Pole Telescope (SPT) and a sample of Type 1 and 2 (unobscured and obscured) quasars at z~1 selected with the Wide-field Infrared Survey Explorer, over 2500 deg^2. The cross-power spectrum is detected at ~7-sigma, and we measure a linear bias b=1.67+/-0.24, consistent with clustering analyses. Using an independent lensing map, derived from Planck observations, to measure the cross-spectrum, we find excellent agreement with the SPT analysis. The bias of the combined sample of Type 1 and 2 quasars determined in this work is similar to that previously determined for Type 1 quasars alone; we conclude that that obscured and unobscured quasars must trace the matter field in a similar way. This result has implications for our understanding of quasar unification and evolution schemes.

We obtained Subaru FMOS observations of Halpha emitting galaxies selected from the HiZELS narrow-band survey, to investigate the relationship between stellar mass, metallicity and star-formation rate at z = 0.84 - 1.47, for comparison with the Fundamental Metallicity Relation seen at low redshift. Our findings demonstrate, for the first time with a homogeneously selected sample, that a relationship exists for typical star-forming galaxies at z = 1 - 1.5 and that it is surprisingly similar to that seen locally. Therefore, star-forming galaxies at z = 1 - 1.5 are no less metal abundant than galaxies of similar mass and star formation rate (SFR) at z = 0.1, contrary to claims from some earlier studies. We conclude that the bulk of the metal enrichment for this star-forming galaxy population takes place in the 4 Gyr before z = 1.5. We fit a new mass-metallicity-SFR plane to our data which is consistent with other high redshift studies. However, there is some evidence that the mass-metallicity component of this high redshift plane is flattened, at all SFR, compared with z = 0.1, suggesting that processes such as star-formation driven winds, thought to remove enriched gas from low mass halos, are yet to have as large an impact at this early epoch. The negative slope of the SFR-metallicity relation from this new plane is consistent with the picture that the elevation in the SFR of typical galaxies at z > 1 is fuelled by the inflow of metal-poor gas and not major merging.

We investigate the multi-wavelength properties of a sample of 450-\mu m selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4\sigma\ in deep SCUBA-2 450-\mu m observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of \sigma 450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450\mu m-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <\beta_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between \beta_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 \mu m sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

We present a wide-field (30' diameter) 850um SCUBA-2 map of the spectacular three-component merging supercluster, RCS 231953+00, at z=0.9. The brightest submillimetre galaxy (SMG) in the field (S_850=12mJy) is within 30" of one of the cluster cores (RCS 2319-C), and is likely to be a more distant, lensed galaxy. Interestingly, the wider field around RCS 2319-C reveals a local overdensity of SMGs, exceeding the average source density by a factor of 4.5, with a <1 per cent chance of being found in a random field. Utilizing Herschel-SPIRE observations, we find three of these SMGs have similar submillimetre colours. We fit their observed 250-850um spectral energy distributions to estimate their redshift, yielding 2.5<z<3.5, and calculate prodigious star formation rates (SFRs) ranging from 500-2500 solar masses per year. We speculate that these galaxies are either lensed SMGs, or signpost a physical structure at z~3: a 'protocluster' inhabited by young galaxies in a rapid phase of growth, destined to form the core of a massive galaxy cluster by z=0.

This paper discusses the evolution of the correlation between galaxy star formation rates (SFRs) and stellar mass (M*) over the last ~10 Gyrs, particularly focusing on its environmental dependence. We first present the mid-infrared (MIR) properties of the H-alpha-selected galaxies in a rich cluster Cl0939+4713 at z=0.4. With a MIR stacking analysis, we find that the median SFR of H-alpha emitters is higher in higher-density environment at z=0.4. We also find that star-forming galaxies in high-density environment tend to have higher specific SFR (SSFR), although the trend is much less significant compared to that of SFR. This increase of SSFR in high-density environment is not visible when we consider the SFR derived from H-alpha alone, suggesting that the dust attenuation in galaxies depends on environment; galaxies in high-density environment tend to be dustier (by up to ~0.5 mag). We then discuss the environmental dependence of the SFR-M* relation for star-forming galaxies since z~2, by compiling our comparable, narrow-band selected, large H-alpha emitter samples in both distant cluster environments (from MAHALO-Subaru) and field environments (from HiZELS). We find that the SSFR of H-alpha-selected galaxies (at the fixed mass of log(M*/Msun)=10) rapidly evolves as (1+z)^3, but the SFR-M* relation is independent of the environment since z~2, as far as we rely on the H-alpha-based SFRs (with M*-dependent extinction correction). Even if we consider the possible environmental variation in the dust attenuation, we conclude that the difference in the SFR-M* relation between cluster and field star-forming galaxies is always small (~<0.2dex level) at any time in the history of the Universe since z~2. [abridged]

Panchromatic observations of the best candidate HyLIRG from the widest Herschel extragalactic imaging survey have led to the discovery of at least four intrinsically luminous z=2.41 galaxies across a ~100-kpc region - a cluster of starbursting proto-ellipticals. Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densities. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO J=4-3, and double that in CO J=1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to large rotational velocities in two counter-rotating gas disks -- a scenario predicted to lead to the most intense starbursts, which will therefore come in pairs. The disks have M(dyn) of several x 10^11 Msun, and gas fractions of ~40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close to the Eddington limit. Their specific star-formation rates place them ~>5x above the main sequence, which supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile with a simple volumetric star-formation law. N-body and dark matter simulations suggest this system is the progenitor of a B(inary)-type ~10^14.6-Msun cluster.

We present an IR study of high-redshift galaxy clusters with the MIPS camera on the Spitzer Space Telescope. Employing a sample of 42 clusters from the RCS-1 over the redshift range 0.3 < z < 1.0 and spanning an approximate range in mass of 10^{14-15} Msun, we show the number of IR-luminous galaxies in clusters above a fixed IR luminosity of 2x10^{11} Msun per unit cluster mass evolves as (1+z)^{5.1+/-1.9}. These results assume a single star forming galaxy template; the presence of AGN, and an evolution in their relative contribution to the mid-IR galaxy emission, will alter the overall number counts per cluster and their rate of evolution. We infer the total SFR per unit cluster mass and find T_SFR/M_c ~ (1+z)^{5.4+/-1.9}. This evolution can be attributed entirely to the change in the in-falling field galaxy population. The T_SFR/M_c (binned over all redshift) decreases with increasing cluster mass with a slope (T_SFR/M_c ~ M_c^{-1.5+/-0.4}) consistent with the dependence of the stellar-to-total mass per unit cluster mass seen locally. The inferred star formation seen here could produce ~5-10% of the total stellar mass in massive clusters at z = 0. Finally, we show a clear decrease in the number of IR-bright galaxies per unit optical galaxy in the cluster cores, confirming star formation continues to avoid the highest density regions of the universe at z ~ 0.75 (the average redshift of the high-redshift clusters). While several previous studies appear to show enhanced star formation in high-redshift clusters relative to the field we note that these papers have not accounted for the overall increase in galaxy or dark matter density at the location of clusters. Once this is done, clusters at z ~ 0.75 have the same or less star formation per unit mass or galaxy as the field.

We report observations of the CO(2-1) emission of SDSSJ1506+54, a compact (r_e~135pc) starburst galaxy at z=0.6. SDSSJ1506+54 appears to be forming stars close to the limit allowed by stellar radiation pressure feedback models: the measured L_IR/L'_CO 1500 is one of the highest measured for any galaxy. With its compact optical morphology but extended low surface brightness envelope, post-starburst spectral features, high infrared luminosity (L_IR>10^12.5 L_Sun), low gas fraction (M_H2/M_stars~15%), and short gas depletion time (tens of Myr), we speculate that this is a feedback- limited central starburst episode at the conclusion of a major merger. Taken as such, SDSSJ1504+54 epitomizes the brief closing stage of a classic model of galaxy growth: we are witnessing a key component of spheroid formation during what we term a 'redline' starburst.

The RCS 2319+00 supercluster is a massive supercluster at z=0.9 comprising three optically selected, spectroscopically confirmed clusters separated by <3 Mpc on the plane of the sky. This supercluster is one of a few known examples of the progenitors of present-day massive clusters (10^{15} Msun by z~0.5). We present an extensive spectroscopic campaign carried out on the supercluster field resulting, in conjunction with previously published data, in 1961 high confidence galaxy redshifts. We find 302 structure members spanning three distinct redshift walls separated from one another by ~65 Mpc. The component clusters have spectroscopic redshifts of 0.901, 0.905 and 0.905. The velocity dispersions are consistent with those predicted from X-ray data, giving estimated cluster masses of ~10^{14.5} - 10^{14.9} Msun. The Dressler-Shectman test finds evidence of substructure in the supercluster field and a friends-of-friends analysis identified 5 groups in the supercluster, including a filamentary structure stretching between two cluster cores previously identified in the infrared by Coppin et al. (2012). The galaxy colors further show this filamentary structure to be a unique region of activity within the supercluster, comprised mainly of blue galaxies compared to the ~43-77% red-sequence galaxies present in the other groups and cluster cores. Richness estimates from stacked luminosity function fits results in average group mass estimates consistent with ~10^{13} Msun halos. Currently, 22% of our confirmed members reside in >~10^{13} Msun groups/clusters destined to merge onto the most massive cluster, in agreement with the massive halo galaxy fractions important in cluster galaxy pre-processing in N-body simulation merger tree studies.

We present results from a ~100 ks Chandra observation of the 2QZ Cluster 1004+00 structure at z = 2.23 (hereafter, 2QZ Clus). 2QZ Clus was originally identified as an overdensity of four optically-selected QSOs at z = 2.23 within a 15x15 arcmin^2 region. Narrow-band imaging in the near-IR revealed that the structure contains an additional overdensity of 22 z = 2.23 Halpha-emitting galaxies (HAEs), resulting in 23 unique z = 2.23 HAEs/QSOs. Our Chandra observations reveal that 3 HAEs in addition to the 4 QSOs harbor powerfully accreting supermassive black holes (SMBHs), with 2-10 keV luminosities of ~(8-60) x 10^43 ergs/s and X-ray spectral slopes consistent with unobscured AGN. Using a large comparison sample of 210 HAEs in Chandra-COSMOS (C-COSMOS), we find suggestive evidence that the AGN fraction increases with local HAE galaxy density. The 2QZ Clus HAEs reside in a moderately overdense environment (a factor of ~2 times over the field), and after excluding optically-selected QSOs, we find the AGN fraction is a factor of ~3.5^+3.8_-2.2 times higher than C-COSMOS HAEs in similar environments. Using stacking analyses of the Chandra data and Herschel SPIRE observations at 250 um, we respectively estimate mean SMBH accretion rates (Mdot) and star-formation rates (SFRs) for the 2QZ Clus and C-COSMOS samples. We find that the mean 2QZ Clus HAE stacked 2-10 keV luminosity is QSO-like (~6-10 x 10^43 ergs/s), and the implied Mdot/SFR \sim (1.6-3.2) x10^-3 is broadly consistent with the local MBH/M* relation and z ~ 2 X-ray selected AGN. The C-COSMOS HAEs are on average an order of magnitude less X-ray luminous and have Mdot/SFR \sim (0.2-0.4) x10^-3, comparable to z ~ 1-2 star-forming galaxies with similar mean X-ray luminosities. We estimate that a periodic QSO phase with duty cycle ~2-8% would be sufficient to bring the star-forming galaxies onto the local MBH/M* relation.

We use the HiZELS narrow-band H-alpha survey in combination with CANDELS, UKIDSS and WIRDS near-infrared imaging, to investigate the morphologies, merger rates and sizes of a sample of H-alpha emitting galaxies in the redshift range z=0.40 - 2.23, an epoch encompassing the rise to the peak of the star formation rate density. Merger rates are estimated from space- and ground-based imaging using the M20 coefficient. To account for the increase in the specific star-formation rate (sSFR) of the star forming `main-sequence' with redshift, we normalise the star-formation rate of galaxies at each epoch to the typical value derived from the H-alpha luminosity function. Once this trend in sSFR is removed we see no evidence for an increase in the number density of star-forming galaxies or the merger rate with redshift. We thus conclude that neither is the main driver of the enhanced star-formation rate density at z=1-2, with secular processes such as instabilities within efficiently fuelled, gas-rich discs or multiple minor mergers the most likely alternatives. However, we find that 40-50% of starburst galaxies, those with enhanced specific star formation at their epoch, are major mergers and this fraction is redshift independent. Finally, we find the surprising result that the typical size of a star-forming galaxy of a given mass does not evolve across the redshift range considered, suggesting a universal size-mass relation. Taken in combination, these results indicate a star-forming galaxy population that is statistically similar in physical size, merger rate and mass over the 6 Gyr covered in this study, despite the increase in typical sSFR.

The clustering of galaxies observed in future redshift surveys will provide a wealth of cosmological information. Matching the signal at different redshifts constrains the dark energy driving the acceleration of the expansion of the Universe. In tandem with these geometrical constraints, redshift-space distortions (RSD) depend on the build up of large-scale structure. As pointed out by many authors measurements of these effects are intrinsically coupled. We investigate this link, and argue that it strongly depends on the cosmological assumptions adopted when analysing data. Using representative assumptions for the parameters of the Euclid survey in order to provide a baseline future experiment, we show how the derived constraints change due to different model assumptions. We argue that even the assumption of a Friedman-Robertson-Walker (FRW) space-time is sufficient to reduce the importance of the coupling to a significant degree. Taking this idea further, we consider how the data would actually be analysed and argue that we should not expect to be able to simultaneously constrain multiple deviations from the standard $\Lambda$CDM model. We therefore consider different possible ways in which the Universe could deviate from the $\Lambda$CDM model, and show how the coupling between geometrical constraints and structure growth affects the measurement of such deviations.

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

We present adaptive optics assisted integral field spectroscopy of nine Halpha-selected galaxies at z=0.84--2.23 selected from the HiZELS narrow-band survey. Our observations map the star-formation and kinematics of these representative star-forming galaxies on ~kpc-scales. We demonstrate that within the ISM of these galaxies, the velocity dispersion of the star-forming gas (\sigma) follows a scaling relation \sigma\propto\Sigma_SFR^(1/n)+constant (where \Sigma_SFR is the star formation surface density and the constant includes the stellar surface density). Assuming the disks are marginally stable (Toomre Q=1), we show that this follows from the Kennicutt-Schmidt relation (\Sigma_SFR = A\Sigma_gas^n), and we use the data to derive best fit parameters of n=1.34+/-0.15 and A=3.4_(-1.6)^(+2.5)x10^(-4)Mo/yr/kpc^2, consistent with the local relation and implying cold molecular gas masses of M_gas=10^(9-10)Mo and molecular gas fractions M_gas/(M_gas+Mstars)=0.3+/-0.1, with a range of 10-75%. These values confirm the high gas fractions for high-redshift star-forming galaxies, independent of CO-H_2 conversion factor. We also identify eleven ~kpc-scale star-forming regions (clumps) within our sample and show that their sizes are comparable to the wavelength of the fastest growing unstable mode. The luminosities and velocity dispersions of these clumps follow the same scaling relations as local HII regions, although their star formation densities are a factor 15+/-5x higher than typically found locally. We discuss how the clump properties are related to the disk, and show that their high masses and luminosities are a consequence of the high disk surface density.

We compare cosmic microwave background lensing convergence maps derived from South Pole Telescope (SPT) data with galaxy survey data from the Blanco Cosmology Survey, the Wide-field Infrared Survey Explorer, and a new large Spitzer/IRAC field designed to overlap with the SPT survey. Using optical and infrared catalogs covering between 17 and 68 square degrees of sky, we detect correlation between the SPT convergence maps and each of the galaxy density maps at >4 sigma, with zero cross-correlation robustly ruled out in all cases. The amplitude and shape of the cross-power spectra are in good agreement with theoretical expectations and the measured galaxy bias is consistent with previous work. The detections reported here utilize a small fraction of the full 2500 square degree SPT survey data and serve as both a proof of principle of the technique and an illustration of the potential of this emerging cosmological probe.

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

We present the results of a search for molecular gas emission from a star-forming galaxy at z = 4.9. The galaxy benefits from magnification of 22 +/- 5x due to strong gravitational lensing by the foreground cluster MS1358+62. We target the CO(5-4) emission at a known position and redshift from existing Hubble Space Telescope/ACS imaging and Gemini/NIFS [OII]3727 imaging spectroscopy, and obtain a tentative detection at the 4.3sigma level with a flux of 0.104 +/- 0.024Jkm/s. From the CO line luminosity and assuming a CO-to-H2 conversion factor alpha=2, we derive a gas mass M_gas ~ 1^{+1}_{-0.6} x 10^9 M_sun. Combined with the existing data, we derive a gas fraction Mgas/(Mgas + M*) = 0.59^{+0.11}_{-0.06}. The faint line flux of this galaxy highlights the difficulty of observing molecular gas in representative galaxies at this epoch, and suggests that routine detections of similar galaxies in the absence of gravitational lensing will remain challenging even with ALMA in full science operations.

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

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

We present adaptive optics assisted, spatially resolved spectroscopy of a sample of nine H-alpha-selected galaxies at z=0.84--2.23 drawn from the HiZELS narrow-band survey. These galaxies have star-formation rates of 1-27Mo/yr and are therefore representative of the typical high-redshift star-forming population. Our ~kpc-scale resolution observations show that approximately half of the sample have dynamics suggesting that the ionised gas is in large, rotating disks. We model their velocity fields to infer the inclination-corrected, asymptotic rotational velocities. We use the absolute B-band magnitudes and stellar masses to investigate the evolution of the B-band and stellar mass Tully-Fisher relationships. By combining our sample with a number of similar measurements from the literature, we show that, at fixed circular velocity, the stellar mass of star-forming galaxies has increased by a factor 2.5 between z=2 and z=0, whilst the rest-frame B-band luminosity has decreased by a factor ~6 over the same period. Together, these demonstrate a change in mass-to-light ratio in the B-band of Delta(M/L_B)/(M/L_B)_(z=0) \sim 3.5 between z=1.5 and z=0, with most of the evolution occurring below z=1. We also use the spatial variation of [NII]/Halpha to show that the metallicity of the ionised gas in these galaxies declines monotonically with galacto-centric radius, with an average Delta(log O/H)/DeltaR=-0.027+/-0.005dex/kpc. This gradient is consistent with predictions for high-redshift disk galaxies from cosmologically based hydrodynamic simulations.

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.

We make use of our 'minimal' cold interstellar medium (ISM) emission line model that predicts the molecular and atomic line emission per unit dense, star-forming gas mass (Geach & Papadopoulos 2012; Paper I) to examine the utility of key line ratios in surveys of the so-called star formation 'mode' as traced by xi_SF = M_dense(H_2)/M_total(H_2). We argue that xi_SF and its proxies provide very sensitive, extinction-free discriminators of rapid starburst/merger-driven versus secular quiescent/disk-like stellar mass assembly, with the most promising diagnostic to be applied in the near-future being CO(4-3)/[CI](1-0). These lines are accessible across nearly the full range 0<z<2 (thus covering the bulk of galaxy evolution) with the Atacama Large Millimeter Array. In addition to their diagnostic power, another advantage of this combination is the similar observed frequencies (Delta nu_0 ~ 30 GHz) of the lines, resulting in nearly spatially-matched beams for a fixed aperture, thus mitigating the effects of resolution/morphology bias in the interpretation of galaxy-averaged line ratios. Finally we discuss the capability of deep blind redshift surveys with the high frequency component of the Square Kilometer Arrray (SKA) in discovering H_2-rich galaxies with very low xi_SF values. These could be the progenitors of starburst galaxies seen prior to the onset of star formation; such galaxies could be a class of extreme outliers from local (gas surface density)-(star formation rate) scaling laws, which would exclude them from current star formatation or stellar mass selected samples. Our conservative model suggests that SKA could detect such systems residing at z~3 at a rate of 20-200 per hour.

We predict the space density of molecular gas reservoirs in the Universe, and place a lower limit on the number counts of carbon monoxide (CO), hydrogen cyanide (HCN) molecular and [CII] atomic emission lines in blind redshift surveys in the submillimeter-centimeter spectral regime. Our model uses: (a) recently available HCN Spectral Line Energy Distributions (SLEDs) of local Luminous Infrared Galaxies (LIRGs, L_IR>10^11 L_sun), (b) a value for epsilon=SFR/M_dense(H_2) provided by new developments in the study of star formation feedback on the interstellar medium and (c) a model for the evolution of the infrared luminosity density. Minimal 'emergent' CO SLEDs from the dense gas reservoirs expected in all star-forming systems in the Universe are then computed from the HCN SLEDs since warm, HCN-bright gas will necessarily be CO-bright, with the dense star-forming gas phase setting an obvious minimum to the total molecular gas mass of any star-forming galaxy. We include [CII] as the most important of the far-infrared cooling lines. Optimal blind surveys with the Atacama Large Millimeter Array (ALMA) could potentially detect very distant (z~10-12) [CII] emitters in the >ULIRG galaxy class at a rate of ~0.1-1 per hour (although this prediction is strongly dependent on the star formation and enrichment history at this early epoch), whereas the (high-frequency) Square Kilometer Array (SKA) will be capable of blindly detecting z>3 low-J CO emitters at a rate of ~40-70 per hour. The [CII] line holds special promise for the detection of metal-poor systems with extensive reservoirs of CO-dark molecular gas where detection rates with ALMA can reach up to 2-7 per hour in Bands 4-6.

(abridged) We present a clustering analysis of 370 high-confidence H-alpha emitters (HAEs) at z=2.23. The HAEs are detected in the Hi-Z Emission Line Survey (HiZELS), a large-area blank field 2.121um narrowband survey using the United Kingdom Infrared Telescope (UKIRT) Wide Field Camera (WFCAM). Averaging the two-point correlation function of HAEs in two ~1 degree scale fields (United Kingdom Infrared Deep Sky Survey/Ultra Deep Survey [UDS] and Cosmological Evolution Survey [COSMOS] fields) we find a clustering amplitude equivalent to a correlation length of r_0=3.7+/-0.3 Mpc/h for galaxies with star formation rates of >7 M_sun/yr. The data are well-fitted by the expected correlation function of Cold Dark Matter, scaled by a bias factor: omega_HAE=b^2 omega_DM where b=2.4^{+0.1}_{-0.2}. The corresponding 'characteristic' mass for the halos hosting HAEs is log(M_h/[M_sun/h])=11.7+/-0.1. Comparing to the latest semi-analytic GALFORM predictions for the evolution of HAEs in a LCDM cosmology, we find broad agreement with the observations, with GALFORM predicting a HAE correlation length of ~4 Mpc/h. Motivated by this agreement, we exploit the simulations to construct a parametric model of the halo occupation distribution (HOD) of HAEs, and use this to fit the observed clustering. Our best-fitting HOD can adequately reproduce the observed angular clustering of HAEs, yielding an effective halo mass and bias in agreement with that derived from the scaled omega_DM fit, but with the relatively small sample size the current data provide a poor constraint on the HOD. Our results support the broad picture that 'typical' (~L*) star-forming galaxies have been hosted by dark matter haloes with M_h<10^12 M_sun/h since z~2, but with a broad occupation distribution and clustering that is likely to be a strong function of luminosity.

We have discovered a 2.5 Mpc (projected) long filament of infrared-bright galaxies connecting two of the three ~5x10^14 Msun clusters making up the RCS 2319+00 supercluster at z=0.9. The filament is revealed in a deep Herschel Spectral and Photometric Imaging REceiver (SPIRE) map that shows 250-500um emission associated with a spectroscopically identified filament of galaxies spanning two X-ray bright cluster cores. We estimate that the total (8-1000um) infrared luminosity of the filament is Lir~5x10^12 Lsun, which, if due to star formation alone, corresponds to a total SFR 900 Msun/yr. We are witnessing the scene of the build-up of a >10^15 Msun cluster of galaxies, seen prior to the merging of three massive components, each of which already contains a population of red, passive galaxies that formed at z>2. The infrared filament demonstrates that significant stellar mass assembly is taking place in the moderate density, dynamically active circumcluster environments of the most massive clusters at high-redshift, and this activity is concomitant with the hierarchical build-up of large scale structure.

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

We present an application of unsupervised machine learning - the self-organised map (SOM) - as a tool for visualising, exploring and mining the catalogues of large astronomical surveys. Self-organisation culminates in a low-resolution representation of the 'topology' of a parameter volume, and this can be exploited in various ways pertinent to astronomy. Using data from the Cosmological Evolution Survey (COSMOS), we demonstrate two key astronomical applications of the SOM: (i) object classification and selection, using the example of galaxies with active galactic nuclei as a demonstration, and (ii) photometric redshift estimation, illustrating how SOMs can be used as totally empirical predictive tools. With a training set of ~3800 galaxies with z_spec<1, we achieve photometric redshift accuracies competitive with other (mainly template fitting) techniques that use a similar number of photometric bands (sigma(Dz)=0.03 with a ~2% outlier rate when using u*-band to 8um photometry). We also test the SOM as a photo-z tool using the PHoto-z Accuracy Testing (PHAT) synthetic catalogue of Hildebrandt et al. (2010), which compares several different photo-z codes using a common input/training set. We find that the SOM can deliver accuracies that are competitive with many of the established template-fitting and empirical methods. This technique is not without clear limitations, which are discussed, but we suggest it could be a powerful tool in the era of extremely large - 'petabyte' - databases where efficient data-mining is a paramount concern.

We present a new cluster detection algorithm designed for the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) survey but with generic application to any multiband data. The method makes no prior assumptions about the properties of clusters other than (a) the similarity in colour of cluster galaxies (the "red sequence") and (b) an enhanced projected surface density. The detector has three main steps: (i) it identifies cluster members by photometrically filtering the input catalogue to isolate galaxies in colour-magnitude space, (ii) a Voronoi diagram identifies regions of high surface density, (iii) galaxies are grouped into clusters with a Friends-of-Friends technique. Where multiple colours are available, we require systems to exhibit sequences in two colours. In this paper we present the algorithm and demonstrate it on two datasets. The first is a 7 square degree sample of the deep Sloan Digital Sky Survey equatorial stripe (Stripe 82), from which we detect 97 clusters with z<=0.6. Benefiting from deeper data, we are 100% complete in the maxBCG optically-selected cluster catalogue (based on shallower single epoch SDSS data) and find an additional 78 previously unidentified clusters. The second dataset is a mock Medium Deep Survey (MDS) Pan-STARRS catalogue, based on the Lambda-CDM model and a semi-analytic galaxy formation recipe. Knowledge of galaxy-halo memberships in the mock allows a quantification of algorithm performance. We detect 305 mock clusters in haloes with mass >10^13 solar masses at z<=0.6 and determine a spurious detection rate of <1%, consistent with tests on the Stripe 82 catalogue. The detector performs well in the recovery of model Lambda-CDM clusters. (abridged)

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

We present the results of a narrow-band (H2S1, 2.121/0.021um) imaging search with WFCAM/UKIRT for H-alpha emitters around several potential signposts of rare (10^-7-10^-8 Mpc^-3) over-dense regions at z=2.23: an over-density of QSOs (2QZ cluster), a powerful, high-redshift radio galaxy (HzRG), and a concentration of submillimetre galaxies (SMGs) and optically faint radio galaxies (OFRGs). In total, we detect 137 narrow-band emitter candidates down to emission-line fluxes of 0.5-1 x 10^-16 erg s^-1 cm^-2, across a total area of 0.56 sq. degrees (2.1 x 10^5 comoving Mpc at z=2.23) in these fields. The BzK colours of the emitters suggest that at least 80% of our sample are likely to be H-alpha emitters (HAEs) at z=2.23. This is one of the largest HAE samples known at z>2. We find modest (~3-sigma) local over-densities of emitters associated with all the three targets. In the 2QZ cluster field, the emitters show a striking filamentary structure connecting four of the z=2.23 QSOs extending over 30Mpc (comoving). In the HzRG and SMG/OFRG fields, the structures appear to be smaller and seen only in the vicinities of the targets. The K-band magnitudes and the H-alpha equivalent widths of the emitters are weakly correlated with the over-density of the emitters: emitters in over-dense region are more evolved systems compared to those in under-dense regions at z=2.23. We find several examples of extended HAEs in our target fields, including a striking example with a spatial extent of 7.5 arcsec (60 kpc at z=2.23) in the 2QZ field, suggesting that these are relatively common in high-density regions. We conclude that narrow-band H-alpha surveys are efficient routes to map over-dense regions at high-z and thus to understand the relation between the growth of galaxies and their surrounding large-scale structures.

We present IRAM Plateau de Bure interferometric detections of CO(1-0) emission from a 24um-selected sample of star-forming galaxies at z=0.4. The galaxies have PAH 7.7um-derived star formation rates of SFR~30-60 M_Sun/yr and stellar masses M*~10^{11} M_Sun. The CO(1-0) luminosities of the galaxies imply that the disks still contain a large reservoir of molecular gas, contributing ~20% of the baryonic mass, but have star-formation 'efficiencies' similar to local quiescent disks and gas-dominated disks at z~1.5-2. We reveal evidence that the average molecular gas fraction has undergone strong evolution since z~2, with f_gas ~ (1+z)^{2 +/- 0.5}. The evolution of f_gas encodes fundamental information about the relative depletion/replenishment of molecular fuel in galaxies, and is expected to be a strong function of halo mass. We show that the latest predictions for the evolution of the molecular gas fraction in semi-analytic models of galaxy formation within a LCDM Universe are supported by these new observations.

We present a catalogue of 4098 photometrically selected galaxy clusters with a median redshift <z> = 0.32 in the 270 square degree 'Stripe 82' region of the Sloan Digital Sky Survey (SDSS), covering the celestial equator in the Southern Galactic Cap (-50 < RA < 59 deg, |Dec| < 1.25 deg). Owing to the multi-epoch SDSS coverage of this region, the ugriz photometry is ~2 magnitudes deeper than single scans within the main SDSS footprint. We exploit this to detect clusters of galaxies using an algorithm that searches for statistically significant overdensities of galaxies in a Voronoi tessellation of the projected sky. 32% of the clusters have at least one member with a spectroscopic redshift from existing public data (SDSS Data Release 7, 2SLAQ & WiggleZ), and the remainder have a robust photometric redshift (accurate to ~5-9% at the median redshift of the sample). The weighted average of the member galaxies' redshifts provides a reasonably accurate estimate of the cluster redshift. The cluster catalogue is publicly available for exploitation by the community to pursue a range of science objectives. In addition to the cluster catalogue, we provide a linked catalogue of 18,295 V<21 mag quasar sight-lines with impact parameters within <3 Mpc of the cluster cores selected from the catalogue of Veron et al. (2010). The background quasars cover 0.25 < z < 2, where MgII absorption-line systems associated with the clusters are detectable in optical spectra.

Aimed at understanding the evolution of galaxies in clusters, the GLACE survey is mapping a set of optical lines ([OII]3727, [OIII]5007, Hbeta and Halpha/[NII] when possible) in several galaxy clusters at redshift around 0.40, 0.63 and 0.86, using the Tuneable Filters (TF) of the OSIRIS instrument (Cepa et al. 2005) at the 10.4m GTC telescope. This study will address key questions about the physical processes acting upon the infalling galaxies during the course of hierarchical growth of clusters. GLACE is already ongoing: we present some preliminary results on our observations of the galaxy cluster Cl0024+1654 at z = 0.395; on the other hand, GLACE@0.86 has been approved as ESO/GTC large project to be started in 2011.

We present the results of a panoramic (15 Mpc-scale) survey of the Cl0016+16 supercluster (z=0.55) using Spitzer Space Telescope MIPS 24um and Galaxy Evolution Explorer near-UV (2500A; NUV) imaging. The supercluster regions probed are characterised by several dense nodes connected by a pronounced intermediate-density filamentary structure. We have studied the mid-IR and NUV properties of potential cluster members within a Dz=0.1 photometric redshift slice, compared to an identical blank field selection. We have two main findings: (a) the star-formation rates of individual star-forming galaxies throughout the cluster are not significantly different to identically selected field galaxies, and (b) the cluster harbours pockets of 'accelerated' activity where galaxies have an enhanced probability of undergoing star formation. This observation could be explained in a simple model of 'pre-processing' of galaxies during cluster infall: galaxies in satellite groups have an increased chance of having star-formation triggered via gravitational tidal interactions compared to their counterparts in the field, but there is no environmental mechanism boosting the individual star-formation rates of galaxies. We estimate a lower-limit for the total star-formation rate of galaxies in the supercluster as ~850 Msun/yr (field corrected). If this rate is maintained over the typical infall time of a few Gyr, then the infall population could contribute ~1-2x10^12 Msun of stellar mass to the structure.

At z=0, clusters are primarily populated by red, elliptical and massive galaxies, while blue, spiral and lower-mass galaxies are common in low-density environments. Understanding how and when these differences were established is of absolute importance for our understanding of galaxy formation and evolution, but results at high-z remain contradictory. By taking advantage of the widest and deepest H-alpha narrow-band survey at z=0.84 over the COSMOS and UKIDSS UDS fields, probing a wide range of densities (from poor fields to rich groups and clusters, including a confirmed super-cluster with a striking filamentary structure), we show that the fraction of star-forming galaxies falls continuously from ~40% in fields to approaching 0% in rich groups/clusters. We also find that the median SFR increases with environmental density, at least up to group densities - but only for low and medium mass galaxies, and thus such enhancement is mass-dependent at z~1. The environment also plays a role in setting the faint-end slope (alpha) of the H-alpha luminosity function. Our findings provide a sharper view on galaxy formation and evolution and reconcile previously contradictory results at z~1: stellar mass is the primary predictor of star formation activity, but the environment also plays a major role.

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

We present an analysis of the molecular and atomic gas emission in the rest-frame far-infrared and sub-millimetre, from the lensed z=2.3 sub-millimetre galaxy SMM J2135-0102. We obtain very high signal-to-noise detections of 11 transitions from 3 species and limits on a further 20 transitions from 9 species. We use the 12CO, [CI] and HCN line strengths to investigate the gas mass, kinematic structure and interstellar medium (ISM) chemistry, and find strong evidence for a two-phase medium comprising a hot, dense, luminous component and an underlying extended cool, low-excitation massive component. Employing photo-dissociation region models we show that on average the molecular gas is exposed to a UV radiation field that is ~1000 x more intense than the Milky Way, with star-forming regions having a characteristic density of n~10^4 /cm^3. These conditions are similar to those found in local ULIRGs and in the central regions of typical starburst galaxies, even though the star formation rate is far higher in this system. The 12CO spectral line energy distribution and line profiles give strong evidence that the system comprises multiple kinematic components with different conditions, including temperature, and line ratios suggestive of high cosmic ray flux within clouds. We show that, when integrated over the galaxy, the gas and star-formation surface densities appear to follow the Kennicutt-Schmidt relation, although when compared to high-resolution sub-mm imaging, our data suggest that this relation breaks down on scales of <100pc. By virtue of the lens amplification, these observations uncover a wealth of information on the star formation and ISM at z~2.3 at a level of detail that has only recently become possible at z<0.1, and show the potential physical properties that will be studied in unlensed galaxies when ALMA is in full operation. (Abridged).

A space-based galaxy redshift survey would have enormous power in constraining dark energy and testing general relativity, provided that its parameters are suitably optimized. We study viable space-based galaxy redshift surveys, exploring the dependence of the Dark Energy Task Force (DETF) figure-of-merit (FoM) on redshift accuracy, redshift range, survey area, target selection, and forecast method. Fitting formulae are provided for convenience. We also consider the dependence on the information used: the full galaxy power spectrum P(k, P(k) marginalized over its shape, or just the Baryon Acoustic Oscillations (BAO). We find that the inclusion of growth rate information (extracted using redshift space distortion and galaxy clustering amplitude measurements) leads to a factor of ~ 3 improvement in the FoM, assuming general relativity is not modified. This inclusion partially compensates for the loss of information when only the BAO are used to give geometrical constraints, rather than using the full P(k) as a standard ruler. We find that a space-based galaxy redshift survey covering ~20,000 deg^2 over 0.5< z < 2 with \sigma_z/(1+z)<= 0.001 exploits a redshift range that is only easily accessible from space, extends to sufficiently low redshifts to allow both a vast 3-D map of the universe using a single tracer population, and overlaps with ground-based surveys to enable robust modeling of systematic effects. We argue that these parameters are close to their optimal values given current instrumental and practical constraints.

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

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

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

We report observations of Lyman Alpha Blob 1 (LAB1) in the SSA 22 protocluster region (z=3.09) with the integral-field spectrograph SAURON. We increased the signal-to-noise in the spectra by more than a factor three compared to our previous observations. This allows us to probe the structure of the LAB system in detail, examining its structure in the spatial and wavelength dimensions. We find that the emission from the system comes largely from five distinct blobs. Two of the emission regions are associated with Lyman Break Galaxies, while a third appears to be associated with a heavily obscured submillimeter galaxy. The fourth and fifth components do not appear to be associated with any galaxy despite the deep imaging that is available in this field. If we interpret wavelength shifts in the line centroid as velocity structure in the underlying gas, many of these emission systems show evidence of velocity shear. It remains difficult to distinguish between an underlying rotation of the gas and an outflow driven by the central object. We have examined all of the line profiles for evidence of strong absorption features. While several systems are better fitted by the inclusion of a weak absorption component, we do not see evidence for a large-scale coherent absorption feature such as that seen in LAB2.

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

Future galaxy redshift surveys aim to measure cosmological quantities from the galaxy power spectrum. A prime example is the detection of baryonic acoustic oscillations (BAOs), providing a standard ruler to measure the dark energy equation of state, w(z), to high precision. The strongest practical limitation for these experiments is how quickly accurate redshifts can be measured for sufficient galaxies to map the large-scale structure. A promising strategy is to target emission-line (i.e. star-forming) galaxies at high-redshift (z~0.5-2); not only is the space density of this population increasing out to z~2, but also emission-lines provide an efficient method of redshift determination. Motivated by the prospect of future dark energy surveys targeting H-alpha emitters at near-infrared wavelengths (i.e. z>0.5), we use the latest empirical data to model the evolution of the H-alpha luminosity function out to z~2, and thus provide predictions for the abundance of H-alpha emitters for practical limiting fluxes. We caution that the estimates presented in this work must be tempered by an efficiency factor, epsilon, giving the redshift success rate from these potential targets. For a range of practical efficiencies and limiting fluxes, we provide an estimate of nP_{0.2}, where n is the 3D galaxy number density and P_{0.2} is the galaxy power spectrum evaluated at k=0.2h/Mpc. Ideal surveys must provide nP_{0.2}>1 in order to balance shot-noise and cosmic variance errors. We show that a realistic emission-line survey (epsilon=0.5) could achieve nP_{0.2}=1 out to z~1.5 with a limiting flux of 10^{-16} erg/s/cm^{-2}. If the limiting flux is a factor 5 brighter, then this goal can only be achieved out to z~0.5, highlighting the importance of survey depth and efficiency in cosmological redshift surveys.

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

We have undertaken a deep (sigma~1.1 mJy) 1.1-mm survey of the z=0.54 cluster MS 0451.6-0305 using the AzTEC camera on the James Clerk Maxwell Telescope. We detect 36 sources with S/N>3.5 in the central 0.10 deg^2 and present the AzTEC map, catalogue and number counts. We identify counterparts to 18 sources (50%) using radio, mid-infrared, Spitzer IRAC and Submillimeter Array data. Optical, near- and mid-infrared spectral energy distributions are compiled for the 14 of these galaxies with detectable counterparts, which are expected to contain all likely cluster members. We then use photometric redshifts and colour selection to separate background galaxies from potential cluster members and test the reliability of this technique using archival observations of submillimetre galaxies. We find two potential MS 0451-03 members, which, if they are both cluster galaxies have a total star-formation rate (SFR) of ~100 solar masses per year -- a significant fraction of the combined SFR of all the other galaxies in MS 0451-03. We also examine the stacked rest-frame mid-infrared, millimetre and radio emission of cluster members below our AzTEC detection limit and find that the SFRs of mid-IR selected galaxies in the cluster and redshift-matched field populations are comparable. In contrast, the average SFR of the morphologically classified late-type cluster population is ~3 times less than the corresponding redshift-matched field galaxies. This suggests that these galaxies may be in the process of being transformed on the red-sequence by the cluster environment. Our survey demonstrates that although the environment of MS 0451-03 appears to suppress star-formation in late-type galaxies, it can support active, dust-obscured mid-IR galaxies and potentially millimetre-detected LIRGs.

We analyse a 100-ks Chandra observation of the powerful radio galaxy, 4C60.07 at z=3.79. We identify extended X-ray emission with Lx~10^45 erg/s across a ~90-kpc region around the radio galaxy. The energetics of this X-ray halo and its morphological similarity to the radio emission from the galaxy suggest that it arises from inverse Compton (IC) scattering, by relativistic electrons in the radio jets, of Cosmic Microwave Background photons and potentially far-infrared photons from the dusty starbursts around this galaxy. The X-ray emission has a similar extent and morphology to the Ly-alpha halo around the galaxy, suggesting that it may be ionising this halo. Indeed we find that the GHz-radio and X-ray and Ly-alpha luminosities of the halo around 4C60.07 are identical to those of 4C41.17 (also at z=3.8) implying that these three components are linked by a single physical process. This is only the second example of highly-extended IC emission known at z>3, but it underlines the potential importance of IC emission in the formation of the most massive galaxies at high redshifts. In addition, we detect two X-ray luminous active galactic nuclei (AGN) within ~30kpc of the radio galaxy. These two companion AGN imply that the radio and starburst activity in the radio galaxy is triggered through multiple mergers of massive progenitors on a short timescale, ~100Myrs. These discoveries demonstrate the wealth of information which sensitive X-ray observations can yield into the formation of massive galaxies at high redshifts.

We present X-ray point-source catalogs for a deep 400 ks Chandra ACIS-I exposure of the SSA22 field. The observations are centred on a z = 3.09 protocluster, which is populated by Lyman break galaxies (LBGs), Lyalpha emitters (LAEs), and extended Lyalpha-emitting blobs (LABs). The survey reaches ultimate (3 count) sensitivity limits of 5.7 X 10^-17 and 3.0 X 10^-16 ergs cm^-2 s^-1 for the 0.5-2 keV and 2-8 keV bands, respectively. These limits make SSA22 the fourth deepest extragalactic Chandra survey yet conducted, and the only one focused on a known high redshift structure. In total, we detect 297 X-ray point sources and identify one obvious bright extended X-ray source over a ~330 arcmin^2 region. In addition to our X-ray catalogs, we provide all available optical spectroscopic redshifts and near-infrared and mid-infrared photometry available for our sources. The basic X-ray and infrared properties of our Chandra sources indicate a variety of source types, although absorbed active galactic nuclei (AGNs) appear to dominate. In total, we have identified 12 X-ray sources (either via optical spectroscopic redshifts or LAE selection) at z =3.06-3.12 that are likely to be associated with the SSA22 protocluster. These sources have X-ray and multiwavelength properties that suggest they are powered by AGN with 0.5-8 keV luminosities in the range of 10^43-10^45 ergs s^-1 We have analysed the AGN fraction of sources in the protocluster as a function of local LAE source density and find suggestive evidence for a correlation between AGN fraction and local LAE source density (at the 96 per cent confidence level), implying that supermassive black hole growth at z = 3 is strongest in the highest density regions.

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

We present the results of a 400ks Chandra survey of 29 extended Ly-alpha emitting nebulae (Ly-alpha Blobs, LABs) in the z=3.09 proto-cluster in the SSA22 field. We detect luminous X-ray counterparts in five LABs, implying a large fraction of active galactic nuclei (AGN) in LABs, f_AGN = 17% down to L_2-32keV ~ 10^44 erg/s. All of the AGN appear to be heavily obscured, with spectral indices implying obscuring column densities of N_H > 10^23 cm^-2. The AGN fraction should be considered a lower limit, since several more LABs not detected with Chandra show AGN signatures in their mid-infrared emission. We show that the UV luminosities of the AGN are easily capable of powering the extended Ly-alpha emission via photo-ionization alone. When combined with the UV flux from a starburst component, and energy deposited by mechanical feedback, we demonstrate that 'heating' by a central source, rather than gravitational cooling is the most likely power source of LABs. We argue that all LABs could be powered in this manner, but that the luminous host galaxies are often just below the sensitivity limits of current instrumentation, or are heavily obscured. No individual LABs show evidence for extended X-ray emission, and a stack equivalent to a >9Ms exposure of an average LAB also yields no statistical detection of a diffuse X-ray component. The resulting diffuse X-ray/Ly-alpha luminosity limit implies there is no hot (T>10^7 K) gas component in these halos, and also rules out inverse Compton scattering of cosmic microwave background photons, or local far-infrared photons, as a viable power source for LABs.

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

We report on the identification of the highest redshift submm-selected source currently known: LESSJ033229.4-275619. This source was detected in the Large Apex BOlometer CAmera (LABOCA) Extended Chandra Deep Field South (ECDFS) Submillimetre Survey (LESS), a sensitive 870-um survey (~1.2-mJy rms) of the full 30'x30' ECDFS with the LABOCA camera on the Atacama Pathfinder EXperiment (APEX) telescope. The submm emission is identified with a radio counterpart for which optical spectroscopy provides a redshift of z=4.76. We show that the bolometric emission is dominated by a starburst with a star formation rate of ~1000 Msun/yr, although we also identify a moderate luminosity Active Galactic Nucleus (AGN) in this galaxy. Thus it has characteristics similar to those of z~2 submm galaxies (SMGs), with a mix of starburst and obscured AGN signatures. This demonstrates that ultraluminous starburst activity is not just restricted to the hosts of the most luminous (and hence rare) QSOs at z~5, but was also occurring in less extreme galaxies at a time when the Universe was less than 10% of its current age. Assuming that we are seeing the major phase of star formation in this galaxy, then we demonstrate that it would be identified as a luminous distant red galaxy at z~3 and that the current estimate of the space density of z>4 SMGs is only sufficient to produce ~10% of the luminous red galaxy population at these early times. However, this leaves open the possibility that some of these galaxies formed through less intense, but more extended star formation events. If the progenitors of all of the luminous red galaxies at z~3 go through an ultraluminous starburst at z>4 then the required volume density of z>4 SMGs will exceed that predicted by current galaxy formation models by more than an order of magnitude.

We present IRAM Plateau de Bure Interferometer 3mm observations of CO(1-0) emission in two 24um-selected starburst galaxies in the outskirts (2-3xR_virial) of the rich cluster Cl0024+16 (z=0.395). The galaxies' inferred far-infrared luminosities place them in the luminous infrared galaxy class (LIRGs, L_FIR>10^11 L_Sun), with star formation rates of ~60 M_Sun/yr. Strong CO(1-0) emission is detected in both galaxies, and we use the CO line luminosity to estimate the mass of cold molecular gas, M(H_2). Assuming M(H_2)/L'_CO = 0.8 M_Sun/(K km^-1 pc^2), we estimate M(H_2) = (5.4-9.1)x10^9 M_Sun for the two galaxies. We estimate the galaxies' dynamical masses from their CO line-widths, M_dyn~1-3x10^10 M_Sun, implying large cold gas fractions in the galaxies' central regions. At their current rates they will complete the assembly of M_Stars~10^10 M_Sun and double their stellar mass within as little as ~150Myr. If these galaxies are destined to evolve into S0s, then the short time-scale for stellar mass assembly implies that their major episode of bulge growth occurs while they are still in the cluster outskirts, long before they reach the core regions. Subsequent fading of the disc component relative to the stellar bulge after the gas reservoirs have been exhausted could complete the transformation of spiral-to-S0.

Swift-detected GRB 080307 showed an unusual smooth rise in its X-ray light-curve around 100 seconds after the burst, at the start of which the emission briefly softened. This `hump' has a longer duration than is normal for a flare at early times and does not demonstrate a typical flare profile. Using a two component power-law-to-exponential model, the rising emission can be modelled as the onset of the afterglow, something which is very rarely seen in Swift-X-ray light-curves. We cannot, however, rule out that the hump is a particularly slow early-time flare, or that it is caused by upscattered reverse shock electrons.

We present results from a new ultra-deep 400 ks Chandra observation of the SSA22 protocluster at z = 3.09. We have studied the X-ray properties of 234 z ~ 3 Lyman break galaxies (LBGs; protocluster and field) and 158 z = 3.09 Ly-alpha emitters (LAEs) in SSA22 to measure the influence of the high-density protocluster environment on the accretion activity of supermassive black holes (SMBHs) in these UV-selected star forming populations. We detect individually X-ray emission from active galactic nuclei (AGNs) in six LBGs and five LAEs; due to small overlap between the LBG and LAE source population, ten of these sources are unique. At least six and potentially eight of these sources are members of the protocluster. These sources have rest-frame 8-32 keV luminosities in the range of L_8-32 keV = (3-50) X 10^{43} ergs/s and an average observed-frame 2-8 keV to 0.5-2 keV band-ratio of ~0.8 (mean effective photon index of Gamma_eff = 1.1), suggesting significant absorption columns of N_H > 10^{22}-10^{24} cm^{-2}. We find that the fraction of LBGs and LAEs in the z = 3.09 protocluster harboring an AGN with L_8-32 keV > 3 X 10^{43} ergs/s is 9.5^{+12.7}_{-6.1}% and 5.1^{+6.8}_{-3.3}%, respectively. These AGN fractions are somewhat larger (by a mean factor of 6.1^{+10.3}_{-3.6}; significant at the 95% confidence level) than z ~ 3 sources found in lower-density "field" environments. Theoretical models imply that these results may be due to the presence of more actively growing and/or massive SMBHs in LBGs and LAEs within the protocluster compared to the field. Such a result is expected in a scenario where enhanced merger activity in the protocluster drives accelerated galaxy and SMBH growth at z > 2-3. (abridged)

We present the results of a Spitzer Infrared Spectrograph (IRS) survey of 24um-selected luminous infrared galaxies (LIRGs, L_IR > 10^11 L_sun) in the rich cluster Cl0024+16 at z=0.4. Optically, these LIRGs resemble unremarkable spiral galaxies with e(a)/e(c) spectral classifications and [Oii]-derived star formation rates (SFRs) of <2 M_sun/yr, generally indistinguishable from the 'quiescent' star forming population in the cluster. Our IRS spectra show that the majority of the 24um-detected galaxies exhibit polycyclic aromatic hydrocarbon (PAH) emission with implied SFRs ~30-60 M_sun/yr, with only one (<10%) of the sample displaying unambiguous evidence of an active galactic nucleus in the mid-infrared. This confirms the presence of a large population of obscured starburst galaxies in distant clusters, which comprise the bulk of the star formation occurring in these environments at z~0.5. We suggest that, although several mechanisms could be at play, these dusty starbursts could be the signature of an important evolutionary transition converting gas-rich spiral galaxies in distant clusters into the passive, bulge-dominated lenticular galaxies that become increasingly abundant in the cores of rich clusters in the ~4Gyr to the present day.

We present results from a near-infrared narrow-band survey of emission-line galaxies at z=2.23, using the Wide Field Camera on the United Kingdom Infrared Telescope. The H2S1 narrow-band filter (lambda_c=2.121um) we employ selects the H-alpha emission line redshifted to z=2.23, and is thus suitable for selecting 'typical' star forming galaxies and active galactic nuclei at this epoch. The pilot study was undertaken in the well studied Cosmological Evolution Survey field (COSMOS) and is already the largest near-infrared narrow-band survey at this depth, with a line flux limit of F(H-alpha)~10^-16 erg/s/cm^2 over 0.60 square degrees, probing ~220x10^3 Mpc^3 (co-moving) down to a limiting star formation rate of 30 M_sun/yr (3sigma). In this paper we present the results from our pilot survey and evaluate the H-alpha luminosity function and estimate the clustering properties of H-alpha emitters at z=2.23 from 55 detected galaxies. The integrated luminosity function is used to estimate the volume averaged star formation rate at z=2.23: SFRD = 0.17^{+0.16}_{-0.09} M_sun/yr/Mpc^3 for L(H-alpha) > 10^42 erg/s. For the first time, we use the H-alpha star-formation tracer to reliably constrain SFRD out to z=2.23 demonstrating the rapid increase in SRFD out to this redshift as well as confirming the flattening in SFRD between z~1-2. In addition to the luminosity distribution, we analyse the clustering properties of these galaxies. Using the 2-point angular correlation function, we estimate a real space correlation length of r_0= 4.2^{+0.4}_{-0.2} h^-1 Mpc. In comparison to models of clustering which take into account bias evolution, we estimate that these galaxies are hosted by dark matter halos of mass M_halo ~ 10^12 M_sun consistent with the progenitors of the Milky Way.

We present multi-object spectroscopy of galaxies in the immediate (Mpc-scale) environments of four low-power (L_1.4 GHz < 10^25 W/Hz) radio galaxies at z~0.5, selected from the Subaru/XMM-Newton Deep Field. We use the spectra to calculate velocity dispersions and central redshifts of the groups the radio galaxies inhabit, and combined with XMM-Newton (0.3-10 keV) X-ray observations investigate the L_X--sigma_v and T_X--sigma_v scaling relationships. All the radio galaxies reside in moderately rich groups -- intermediate environments between poor groups and rich clusters, with remarkably similar X-ray properties. We concentrate our discussion on our best statistical example that we interpret as a low-power (FRI) source triggered within a sub-group, which in turn is interacting with a nearby group of galaxies, containing the bulk of the X-ray emission for the system -- a basic scenario which can be compared to more powerful radio sources at both high (z>4) and low (z<0.1) redshifts. This suggests that galaxy-galaxy interactions triggered by group mergers may play an important role in the life-cycle of radio galaxies at all epochs and luminosities.

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

Using Spitzer Space Telescope IRAC (3.6-8um) and MIPS (24um) imaging, as well as Hubble Space Telescope optical observations, we identify the IRAC counterparts of the luminous power sources residing within the two largest and brightest Lyman-alpha emitting nebulae (LABs) in the SA22 protocluster at z=3.09 (LAB1 and LAB2). These sources are also both submillimeter galaxies (SMGs). From their rest-frame optical/near-infrared colors, we conclude that the SMG in LAB1 is likely to be starburst dominated and heavily obscured (Av~3). In contrast, LAB2 has excess rest-frame ~2um emission (over that expected from starlight) and hosts a hard X-ray source at the proposed location of the SMG, consistent with the presence of an active galactic nucleus (AGN). We conclude that LAB1 and LAB2 appear to have very different energy sources despite having similar Lyman-alpha spatial extents and luminosities, although it remains unclear whether on-going star-formation or periodic AGN heating is responsible for the extended Lyman-alpha emission. We find that the mid-infrared properties of the SMGs lying in LAB1 and LAB2 are similar to those of the wider SMG population, and so it is possible that extended Lyman-alpha haloes are a common feature of SMGs in general.

We report the results of a systematic near-infrared spectroscopic survey using the Subaru, VLT and Keck Telescopes of a sample of high redshift Ultra-luminous Infrared Galaxies (ULIRGs) mainly composed of submillimeter-selected galaxies. Our observations span the restframe optical range containing nebular emission lines such as H_beta, [OIII], and [OII], which are essential for making robust diagnostics of the physical properties of these ULIRGs. Using the H_alpha/H_beta emission line ratios, we derive internal extinction estimates for these galaxies similar to those of local ULIRGs:A_v~2.9+-0.5. Correcting the H_alpha estimates of the star formation rate for dust extinction using the Balmer decrement, results in rates which are consistent with those estimated from the far-infrared luminosity. The majority (>60%) of our sample show spectral features characteristic of AGN, with ~65% exhibiting broad Balmer emission lines. A proportion of these sources show relatively low [OIII]/H_beta line ratios, which are similar to those of Narrow Line Seyfert 1 galaxies suggesting small mass black holes which are rapidly growing. In the subsample of our survey with both [OIII] and hard X-ray coverage, at least ~60% show an excess of [OIII] emission, by a factor of 5-10 times, relative to the hard X-ray luminosity compared to the correlation between these two properties seen in Seyferts and QSOs locally. From our spectral diagnostics, we propose that the strong [OIII] emission in these galaxies arises from shocks in dense gaseous regions in this vigorously star-forming population.

(abridged) We present panoramic Spitzer MIPS 24um observations covering 9x9Mpc (25'x25') fields around two massive clusters, Cl0024+16 and MS0451-03, at z=0.39 and z=0.55. Our observations cover a very wide range of environments within these clusters, from high-density regions around the cores out to the turn-around radius. Cross-correlating the mid-infrared catalogs with deep optical and near-infrared imaging of these fields, we investigate the optical/near-infrared colors of the mid-infrared sources. We find excesses of mid-infrared sources with optical/near-infrared colors expected of cluster members in the two clusters and test this selection using spectroscopically confirmed 24um members. The much more significant excess is associated with Cl0024+16, whereas MS0451-03 has comparatively few mid-infrared sources. The mid-infrared galaxy population in Cl0024+16 appears to be associated with dusty star-forming galaxies (typically redder than the general cluster population by up to A_V~1-2 mags) rather than emission from dusty tori around active galactic nuclei (AGN) in early-type hosts. The inferred total-infrared star-formation rates in Cl0024+16 are typically >5x greater than those found from a similar Halpha survey, indicating significant obscured activity in the cluster population. We find evidence for strong evolution of the level of dust-obscured star-formation in dense environments out to z=0.5, analogous to the rise in fraction of optically-selected star-forming galaxies seen in clusters and the field out to similar redshifts. However, there are clearly significant cluster-to-cluster variations in the populations of mid-infrared sources, probably reflecting differences in the intracluster media and recent dynamical evolution of these systems.

We present the results from a submm survey of a sample of 23 giant Lya emitting nebulae in the overdensity at z=3.09 in the SA22 field. These objects, which have become known as Lya Blobs (LABs) have a diverse range of morphology and surface brightness, but the nature of their power source is unclear - with cooling flows and/or AGN/starburst ionised winds being possibilities. Using the SCUBA submm camera we measure the 850um flux of a sample of LABs, detecting four LABs at >3.5sigma individually, and a modest statistical detection of the full sample at about 3mJy. These fluxes correspond to bolometric luminosities in the ultraluminous regime, with star-formation rates of about 1e3 Msun/yr. We show there is a trend between Lya luminosity and bolometric output, which suggests that a galactic scale superwind generated from starbursts of age 10-100Myr may be responsible for the Lya emission. We estimate the star-formation rate density in SA22 to be >3 Msun/yr/Mpc^3 - greater than the field at this epoch, and note that there are now 7 submm galaxies in the SA22 structure, making this region the richest association of these intensely active galaxies. Finally we suggest that Lya haloes may be a common feature of the submm population in general, and have an important role in the heating and enrichment of the intergalactic medium.

We present deep wide-field K-band observations of six Hyperluminous Infrared Galaxies (HLIRGs) spanning a redshift range 0.44<z<1.55. The sample resides in a wide variety of environments, from the field to Abell 2 clusters, with a mean galaxy-HLIRG clustering amplitude of B_gh=190+/-45. The range in environments, and the mean clustering level, are both greater than those seen in local IR-luminous galaxies, from which we infer that the range of galaxy evolution processes driving IR-luminous galaxy evolution at z>0.5 is greater than locally, and includes mergers between gas-rich spiral galaxies in the field, but also includes encounters in clusters and hierarchical buildup. The similarity in the range of environments and mean clustering amplitude between our sample and QSOs over a similar redshift range is consistent with the interpretation where evolutionary connections between IR-luminous galaxies and QSOs are stronger at z>0.5 than locally, and that, at these redshifts, the processes that drive QSO evolution are similar to those that drive IR-luminous galaxy evolution. From comparison of the HLIRG and QSO host galaxies we further postulate that a larger fraction of IR-luminous galaxies pass through an optical QSO stage at z>0.5 than locally.

We present the K-band Hubble diagrams (K-z relations) of sub-mm-selected galaxies and hyperluminous galaxies (HLIRGs). We report the discovery of a remarkably tight K-z relation of HLIRGs, indistinguishable from that of the most luminous radiogalaxies. Like radiogalaxies, the HLIRG K-z relation at z<~3 is consistent with a passively evolving ~3L* instantaneous starburst starting from a redshift of z~10. In contrast, many sub-mm selected galaxies are >~2 magnitudes fainter, and the population has a much larger dispersion. We argue that dust obscuration and/or a larger mass range may be responsible for this scatter. The galaxies so far proved to be hyperluminous may have been biased towards higher AGN bolometric contributions than sub-mm-selected galaxies due to the 60 micron selection of some, so the location on the K-z relation may be related to the presence of the most massive AGN. Alternatively, a particular host galaxy mass range may be responsible for both extreme star formation and the most massive active nuclei.