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 used wide area surveys over 39 deg$^2$ by the HerMES collaboration, performed with the Herschel Observatory SPIRE multi-wavelength camera, to estimate the low-redshift, $0.02<z<0.5$, monochromatic luminosity functions (LFs) of galaxies at 250, 350 and 500$\,\mu$m. SPIRE flux densities were also combined with Spitzer photometry and multi-wavelength archival data to perform a complete SED fitting analysis of SPIRE detected sources to calculate precise k-corrections, as well as the bolometric infrared (8-1000$\,\mu$m) luminosity functions and their low-$z$ evolution from a combination of statistical estimators. Integration of the latter prompted us to also compute the local luminosity density (LLD) and the comoving star formation rate density (SFRD) for our sources, and to compare them with theoretical predictions of galaxy formation models. The luminosity functions show significant and rapid luminosity evolution already at low redshifts, $0.02<z<0.2$, with L$_{IR}^* \propto (1+z)^{6.0\pm0.4}$ and $\Phi_{IR}^* \propto (1+z)^{-2.1\pm0.4}$, L$_{250}^* \propto (1+z)^{5.3\pm0.2}$ and $\Phi_{250}^* \propto (1+z)^{-0.6\pm0.4}$ estimated using the IR bolometric and the 250$\,\mu$m LFs respectively. Converting our IR LD estimate into an SFRD assuming a standard Salpeter IMF and including the unobscured contribution based on the UV dust-uncorrected emission from local galaxies, we estimate a SFRD scaling of SFRD$_0+0.08 z$, where SFRD$_0\simeq (1.9\pm 0.03)\times 10^{-2} [\mathrm{M}_\odot\,\mathrm{Mpc}^{-3}]$ is our total SFRD estimate at $z\sim0.02$.

We separate the extragalactic radio source population above ~50 uJy into active galactic nuclei (AGN) and star-forming sources. The primary method of our approach is to fit the infrared spectral energy distributions (SEDs), constructed using Spitzer/IRAC and MIPS and Herschel/SPIRE photometry, of 380 radio sources in the Extended Chandra Deep Field-South. From the fitted SEDs, we determine the relative AGN and star-forming contributions to their infrared emission. With the inclusion of other AGN diagnostics such as X-ray luminosity, Spitzer/IRAC colours, radio spectral index and the ratio of star-forming total infrared flux to k-corrected 1.4 GHz flux density, qIR, we determine whether the radio emission in these sources is powered by star formation or by an AGN. The majority of these radio sources (60 per cent) show the signature of an AGN at some wavelength. Of the sources with AGN signatures, 58 per cent are hybrid systems for which the radio emission is being powered by star formation. This implies that radio sources which have likely been selected on their star formation have a high AGN fraction. Below a 1.4 GHz flux density of 1 mJy, along with finding a strong contribution to the source counts from pure star-forming sources, we find that hybrid sources constitute 20-65 per cent of the sources. This result suggests that hybrid sources have a significant contribution, along with sources that do not host a detectable AGN, to the observed flattening of the source counts at ~1mJy for the extragalactic radio source population.

The Herschel Multi-tiered Extragalactic Survey (HerMES) has identified large numbers of dusty star-forming galaxies (DSFGs) over a wide range in redshift. A detailed understanding of these DSFGs is hampered by the limited spatial resolution of Herschel. We present 870um 0.45" resolution imaging from the Atacama Large Millimeter/submillimeter Array (ALMA) of 29 HerMES DSFGs with far-infrared (FIR) flux densities in between the brightest of sources found by Herschel and fainter DSFGs found in ground-based sub-millimeter (sub-mm) surveys. We identify 62 sources down to the 5-sigma point-source sensitivity limit in our ALMA sample (sigma~0.2mJy), of which 6 are strongly lensed (showing multiple images) and 36 experience significant amplification (mu>1.1). To characterize the properties of the ALMA sources, we introduce and make use of uvmcmcfit, a publicly available Markov chain Monte Carlo analysis tool for interferometric observations of lensed galaxies. Our lens models tentatively favor intrinsic number counts for DSFGs with a steep fall off above 8mJy at 880um. Nearly 70% of the Herschel sources comprise multiple ALMA counterparts, consistent with previous research indicating that the multiplicity rate is high in bright sub-mm sources. Our ALMA sources are located significantly closer to each other than expected based on results from theoretical models as well as fainter DSFGs identified in the LABOCA ECDFS Submillimeter Survey. The high multiplicity rate and low projected separations argue in favor of interactions and mergers driving the prodigious emission from the brightest DSFGs as well as the sharp downturn above S_880=8mJy.

Dust-Obscured galaxies (DOGs) are bright 24 um-selected sources with extreme obscuration at optical wavelengths. They are typically characterized by a rising power-law continuum of hot dust (T_D ~ 200-1000K) in the near-IR indicating that their mid-IR luminosity is dominated by an an active galactic nucleus (AGN). DOGs with a fainter 24 um flux display a stellar bump in the near-IR and their mid-IR luminosity appears to be mainly powered by dusty star formation. Alternatively, it may be that the mid-IR emission arising from AGN activity is dominant but the torus is sufficiently opaque to make the near-IR emission from the AGN negligible with respect to the emission from the host component. In an effort to characterize the astrophysical nature of the processes responsible for the IR emission in DOGs, this paper exploits Herschel data (PACS + SPIRE) on a sample of 95 DOGs within the COSMOS field. We derive a wealth of far-IR properties (e.g., total IR luminosities; mid-to-far IR colors; dust temperatures and masses) based on SED fitting. Of particular interest are the 24 um-bright DOGs (F_24um > 1mJy). They present bluer far-IR/mid-IR colors than the rest of the sample, unveiling the potential presence of an AGN. The AGN contribution to the total 8-1000 um flux increases as a function of the rest-frame 8 um-luminosity irrespective of the redshift. This confirms that faint DOGs (L_8 um< 10^12 L_sun) are dominated by star-formation while brighter DOGs show a larger contribution from an AGN.

[Abridged] We study the evolution of the radio spectral index and far-infrared/radio correlation (FRC) across the star-formation rate-stellar masse (i.e. SFR-M*) plane up to z 2. We start from a M*-selected sample of galaxies with reliable SFR and redshift estimates. We then grid the SFR-M* plane in several redshift ranges and measure the infrared luminosity, radio luminosity, radio spectral index, and ultimately the FRC index (i.e. qFIR) of each SFR-M*-z bin. The infrared luminosities of our SFR-M*-z bins are estimated using their stacked far-infrared flux densities inferred from observations obtained with Herschel. Their radio luminosities and radio spectral indices (i.e. alpha, where Snu nu^-alpha) are estimated using their stacked 1.4GHz and 610MHz flux densities from the VLA and GMRT, respectively. Our far-infrared and radio observations include the most widely studied blank extragalactic fields -GOODS-N/S, ECDFS, and COSMOS- covering a sky area of 2deg^2. Using this methodology, we constrain the radio spectral index and FRC index of star-forming galaxies with M*>10^10Msun and 0<z<2.3. We find that alpha^1.4GHz_610MHz does not evolve significantly with redshift or with the distance of a galaxy with respect to the main sequence (MS) of the SFR-M* plane (i.e. Delta_log(SSFR)_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). Instead, star-forming galaxies have a radio spectral index consistent with a canonical value of 0.8, which suggests that their radio spectra are dominated by non-thermal optically thin synchrotron emission. We find that qFIR displays a moderate but statistically significant redshift evolution as qFIR(z)=(2.35+/-0.08)*(1+z)^(-0.12+/-0.04), consistent with some previous literature. Finally, we find no significant correlation between qFIR and Delta_log(SSFR)_MS, though a weak positive trend, as observed in one of our redshift bins, cannot be firmly ruled out using our dataset.

We present the most complete study to date of the X-ray emission from star-formation in high redshift (median z=0.7; z<1.5), IR-luminous (L_IR=10^10-10^13 L_sun) galaxies detected by Herschel's PACS and SPIRE instruments. For our purpose we take advantage of the deepest X-ray data to date, the Chandra deep fields (North and South). Sources which host AGN are removed from our analysis by means of multiple AGN indicators. We find an AGN fraction of 18+/-2 per cent amongst our sample and note that AGN entirely dominate at values of log[L_X/L_IR]>-3 in both hard and soft X-ray bands. From the sources which are star-formation dominated, only a small fraction are individually X-ray detected and for the bulk of the sample we calculate average X-ray luminosities through stacking. We find an average soft X-ray to infrared ratio of log[L_SX/L_IR]=-4.3 and an average hard X-ray to infrared ratio of log[L_HX/L_IR]=-3.8. We report that the X-ray/IR correlation is approximately linear through the entire range of L_IR and z probed and, although broadly consistent with the local (z<0.1) one, it does display some discrepancies. We suggest that these discrepancies are unlikely to be physical, i.e. due to an intrinsic change in the X-ray properties of star-forming galaxies with cosmic time, as there is no significant evidence for evolution of the L_X/L_IR ratio with redshift. Instead they are possibly due to selection effects and remaining AGN contamination. We also examine whether dust obscuration in the galaxy plays a role in attenuating X-rays from star-formation, by investigating changes in the L_X/L_IR ratio as a function of the average dust temperature. We conclude that X-rays do not suffer any measurable attenuation in the host galaxy.

We present a study of the infrared properties for a sample of seven spectroscopically confirmed submillimeter galaxies at $z>$4.0. By combining ground-based near-infrared, Spitzer IRAC and MIPS, Herschel SPIRE, and ground-based submillimeter/millimeter photometry, we construct their Spectral Energy Distributions (SED) and a composite model to fit the SEDs. The model includes a stellar emission component at $\lambda_{\rm rest} <$ 3.5$ \mu$m; a hot dust component peaking at $\lambda_{rest} \sim$ 5$\,\mu$m; and cold dust component which becomes significant for $\lambda_{\rm rest} >$ 50$\,\mu$m. Six objects in the sample are detected at 250 and 350$ \mu$m. The dust temperatures for the sources in this sample are in the range of 40$-$80 K, and their $L_{\rm FIR}$ $\sim$ 10$^{13}$ L$_{\odot}$ qualifies them as Hyper$-$Luminous Infrared Galaxies (HyperLIRGs). The mean FIR-radio index for this sample is around $< q > = 2.2$ indicating no radio excess in their radio emission. Most sources in the sample have 24$ \mu$m detections corresponding to a rest-frame 4.5$ \mu$m luminosity of Log$_{10}$(L$_{4.5}$ / L$_{\odot}$) = 11 $\sim$ 11.5. Their L$_{\rm 4.5}$/$L_{\rm FIR}$ ratios are very similar to those of starburst dominated submillimeter galaxies at $z \sim$ 2. The $L_{\rm CO}-L_{\rm FIR}$ relation for this sample is consistent with that determined for local ULIRGs and SMGs at $z \sim$ 2. We conclude that submillimeter galaxies at $z >$ 4 are hotter and more luminous in the FIR, but otherwise very similar to those at $z \sim$ 2. None of these sources show any sign of the strong QSO phase being triggered.

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/

The promise of multi-wavelength astronomy has been tempered by the large disparity in sensitivity and resolution between different wavelength regimes. Here we present a statistical approach which attempts to overcome this by fitting parametric models directly to image data. Specifically, we fit a model for the radio luminosity function (LF) of star-forming galaxies to pixel intensity distributions at 1.4 GHz coincident with near-IR selected sources in COSMOS. Taking a mass-limited sample in redshift bins across the range $0<z<4$ we are able to fit the radio LF with ~0.2 dex precision in the key parameters (e.g. Phi*,L*). Good agreement is seen between our results and those using standard methods at radio and other wavelengths. Integrating our luminosity functions to get the star formation rate density we find that galaxies with a stellar mass greater than $10^{9.5}\,$M$_{\odot}$ contribute at least 50 per cent of cosmic star formation at since $z=4$. The scalability of our approach is empirically estimated, with the precision in LF parameter estimates found to scale with the number of sources in the stack as $\sqrt{N}$. This type of approach will be invaluable in the multi-wavelength analysis of upcoming surveys with the SKA pathfinder facilities; LOFAR, ASKAP and MeerKAT.

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.

The Herschel Multi-tiered Extragalactic Survey (HerMES) is the largest Guaranteed Time Key Programme on the Herschel Space Observatory. With a wedding cake survey strategy, it consists of nested fields with varying depth and area totalling ~380 deg^2. In this paper, we present deep point source catalogues extracted from Herschel-SPIRE observations of all HerMES fields, except for the later addition of the 270 deg^2 HeLMS field. These catalogues constitute the second Data Release (DR2) made in October 2013. A subset of these catalogues, which consists of bright sources extracted from Herschel-SPIRE observations completed by May 1, 2010 (covering ~ 74 deg^2) were released earlier in the first extensive Data Release (DR1) in March 2012. Two different methods are used to generate the point source catalogues, the SUSSEXtractor (SXT) point source extractor used in two earlier data releases (EDR and EDR2) and a new source detection and photometry method. The latter combines an iterative source detection algorithm, StarFinder (SF), and a De-blended SPIRE Photometry (DESPHOT) algorithm. We use end-to-end Herschel-SPIRE simulations with realistic number counts and clustering properties to characterise basic properties of the point source catalogues, such as the completeness, reliability, photometric and positional accuracy. Over 500, 000 catalogue entries in HerMES fields (except HeLMS) are released to the public through the HeDAM website (http://hedam.oamp.fr/herMES).

The potential for Planck to detect clusters of dusty, star-forming galaxies at z greater than 1 is tested by examining the Herschel-SPIRE images of Planck Early Release Compact Source Catalog (ERCSC) sources lying in fields observed by the HerMES survey. Of the 16 Planck sources that lie in the roughly 90 sq. deg. examined, we find that twelve are associated with single bright Herschel sources. The remaining four are associated with overdensities of Herschel sources, making them candidate clusters of dusty, starforming galaxies. We use complementary optical and NIR data for these clumps to test this idea, and find evidence for the presence of galaxy clusters in all four cases. We use photometric redshifts and red sequence galaxies to estimate the redshifts of these clusters, finding that they range from 0.8 to 2.3. These redshifts imply that the Herschel sources in these clusters, which contribute to the detected Planck flux, are forming stars very rapidly, with typical total cluster star formation rates greater than 1000Msun per yr. The high redshift clusters discovered in these observations are used to constrain the epoch of cluster galaxy formation, finding that the galaxies in our clusters are 1 to 1.5 Gy old at z about 1 to 2. Prospects for the discovery of further clusters of dusty galaxies are discussed, using not only all sky Planck surveys, but also deeper, smaller area, Herschel surveys.

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

We present a method for selecting $z>4$ dusty, star forming galaxies (DSFGs) using Herschel/SPIRE 250/350/500 $\mu m$ flux densities to search for red sources. We apply this method to 21 deg$^2$ of data from the HerMES survey to produce a catalog of 38 high-$z$ candidates. Follow-up of the first 5 of these sources confirms that this method is efficient at selecting high-$z$ DSFGs, with 4/5 at $z=4.3$ to $6.3$ (and the remaining source at $z=3.4$), and that they are some of the most luminous dusty sources known. Comparison with previous DSFG samples, mostly selected at longer wavelengths (e.g., 850 $\mu m$) and in single-band surveys, shows that our method is much more efficient at selecting high-$z$ DSFGs, in the sense that a much larger fraction are at $z>3$. Correcting for the selection completeness and purity, we find that the number of bright ($S_{500\,\mu m} \ge 30$ mJy), red Herschel sources is $3.3 \pm 0.8$ deg$^{-2}$. This is much higher than the number predicted by current models, suggesting that the DSFG population extends to higher redshifts than previously believed. If the shape of the luminosity function for high-$z$ DSFGs is similar to that at $z\sim2$, rest-frame UV based studies may be missing a significant component of the star formation density at $z=4$ to $6$, even after correction for extinction.

We quantify the fraction of the cosmic infrared background (CIB) that originates from galaxies identified in the UV/optical/near-infrared by stacking 81,250 (~35.7 arcmin^2) K-selected sources (K_AB < 24.0), split according to their rest-frame U - V vs. V - J colors into 72,216 star-forming and 9,034 quiescent galaxies, on maps from Spitzer/MIPS (24um), Herschel/SPIRE (250, 350, 500um), Herschel/PACS (100, 160um), and AzTEC (1100um). The fraction of the CIB resolved by our catalog is (69 $\pm$ 15)% at 24um, (78 $\pm$ 17)% at 70um, (58 $\pm$ 13)% at 100um, (78 $\pm$ 18)% at 160um, (80 $\pm$ 17)% at 250um, (69 $\pm$ 14)% at 350um, (65 $\pm$ 12)% at 500um, and (45 $\pm$ 8)% at 1100um. Of that total, about 95% originates from star-forming galaxies, while the remaining 5% is from apparently quiescent galaxies. The CIB at $\lambda$ < 200um is sourced predominantly from galaxies at z < 1, while at $\lambda$ > 200um the bulk originates from 1 < z < 2. Galaxies with stellar masses log(M/ M_sun)=9.5-11 are responsible for the majority of the CIB, with those in the log(M/ M_sun)=9.5-10 contributing mostly at $\lambda$ < 250um, and those in the log(M/ M_sun)=10.5-11 bin dominating at $\lambda$ > 350um. The contribution from galaxies in the log(M/ M_sun)=9.0-9.5 and log(M/ M_sun)=11.0-12.0 stellar mass bins contribute the least, both of order 5%, although the highest stellar-mass bin is a significant contributor to the luminosity density at z > 2. The luminosities of the galaxies responsible for the CIB shifts from a combination of "normal" and luminous infrared galaxies (LIRGs) at $\lambda$ < 160um, to LIRGs at 160um < $\lambda$ < 500um, to finally LIRGs and ultra-luminous infrared galaxies (ULIRGs) at $\lambda$ > 500um. Stacking analyses were performed with SIMSTACK (available at http://www.astro.caltech.edu/~viero/viero_homepage/toolbox.html) which accounts for possible biases due to clustering.

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

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

We combine Herschel PACS and SPIRE maps of the full 2 deg^2 COSMOS field with existing multi-wavelength data to obtain template and model-independent optical-to-far-infrared spectral energy distributions (SEDs) for 4,218 Herschel-selected sources with log(L_IR/L_sun) = 9.4-13.6 and z = 0.02-3.54. Median SEDs are created by binning the optical to far-infrared (FIR) bands available in COSMOS as a function of infrared luminosity. Herschel probes rest-frame wavelengths where the bulk of the infrared radiation is emitted, allowing us to more accurately determine fundamental dust properties of our sample of infrared luminous galaxies. We find that the SED peak wavelength (lambda_peak) decreases and the dust mass (M_dust) increases with increasing total infrared luminosity (L_IR). In the lowest infrared luminosity galaxies (log(L_IR/L_sun) = 10.0-11.5), we see evidence of Polycyclic Aromatic Hydrocarbons (PAH) features (lambda=7-9 um), while in the highest infrared luminosity galaxies (L_IR > 10^12 L_sun) we see an increasing contribution of hot dust and/or power-law emission, consistent with the presence of heating from an active galactic nucleus (AGN). We study the relationship between stellar mass and star formation rate of our sample of infrared luminous galaxies and find no evidence that Herschel-selected galaxies follow the SFR/M_* "main sequence" as previously determined from studies of optically selected, star-forming galaxies. Finally, we compare the mid-infrared (MIR) to FIR properties of our infrared luminous galaxies using the previously defined diagnostic, IR8 = L_IR / L_8, and find that galaxies with L_IR > 10^11.3 L_sun tend to systematically lie above (~3-5) the IR8 "infrared main sequence", suggesting either suppressed PAH emission or an increasing contribution from AGN heating.

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.

Dust-obscured galaxies (DOGs) are a UV-faint, IR-bright galaxy population that reside at z~2 and are believed to be in a phase of dusty star-forming and AGN activity. We present far-IR observations of a complete sample of DOGs in the 2 deg2 of COSMOS. The 3077 DOGs have <z>=1.9+/-0.3 and are selected from 24um and r+ observations using a color cut of r+ - [24]>=7.5 (AB mag) and S24>=100uJy. Based on the near-IR SEDs, 47% are star-formation dominated and 10% are AGN-dominated. We use SPIRE far-IR photometry from HerMES to calculate the IR luminosity and characteristic dust temperature for the 1572 (51%) DOGs that are detected at 250um (>=3{\sigma}). For the remaining 1505 (49%) that are undetected, we perform a median stacking analysis to probe fainter luminosities. Detected and undetected DOGs have average IR luminosities of (2.8+/-0.4) x 1012 LSun and (0.77+-0.08)x10^12LSun, and dust temperatures of 34+/-7 K and 37+/-3 K, respectively. The IR luminosity function of DOGs with S24>=100uJy is calculated, using far-IR observations and stacking. DOGs contribute 10-30% to the total star formation rate density of the Universe at z=1.5-2.5, dominated by 250um detected and bump DOGs. For comparison, DOGs contribute 30% to the star-formation rate density of the Universe for all galaxies at z=1.5-2.5 with S24>=100uJy. DOGs have a large scatter about the star-formation main sequence and their specific star-formation rates show that the observed phase of star-formation could be responsible for their observed stellar mass at z~2.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We have observed four massive galaxy clusters with the SPIRE instrument on the Herschel Space Observatory and measure a deficit of surface brightness within their central region after subtracting sources. We simulate the effects of instrumental sensitivity and resolution, the source population, and the lensing effect of the clusters to estimate the shape and amplitude of the deficit. The amplitude of the central deficit is a strong function of the surface density and flux distribution of the background sources. We find that for the current best fitting faint end number counts, and excellent lensing models, the most likely amplitude of the central deficit is the full intensity of the cosmic infrared background (CIB). Our measurement leads to a lower limit to the integrated total intensity of the CIB of I(250 microns) > 0.69_(-0.03)^(+0.03) (stat.)_(-0.06)^(+0.11) (sys.) MJy/sr, with more CIB possible from both low-redshift sources and from sources within the target clusters. It should be possible to observe this effect in existing high angular resolution data at other wavelengths where the CIB is bright, which would allow tests of models of the faint source component of the CIB.

In this work we explore the impact of the presence of an active galactic nucleus (AGN) on the mid- and far-infrared (IR) properties of galaxies as well as the effects of simultaneous AGN and starburst activity in these same galaxies. To do this we apply a multi-component, multi-band spectral synthesis technique to a sample of 250 micron selected galaxies of the Herschel Multi-tiered Extragalactic Survey, with IRS spectra available for all galaxies. Our results confirm that the inclusion of the IRS spectra plays a crucial role in the spectral analysis of galaxies with an AGN component improving the selection of the best-fit hot dust model (torus). We find a correlation between the obscured star formation rate (SFR) derived from the IR luminosity of the starburst component, SFR_IR and SFR_PAH, derived from the luminosity of the PAH features, L_PAH, with SFR_FIR taking higher values than SFR_PAH. The correlation is different for AGN- and starburst-dominated objects. The ratio of L_PAH to that of the starburst component, L_PAH/L_SB, is almost constant for AGN-dominated objects but decreases with increasing L_SB for starburst-dominated objects. SFR_FIR increases with the accretion luminosity, L_acc, with the increase less prominent for the very brightest, unobscured AGN-dominated sources. We find no correlation between the masses of the hot and cold dust components. We interpret this as a non-constant fraction of gas driven by the gravitational effects to the AGN while the starburst is ongoing. We also find no evidence of the AGN affecting the temperature of the cold dust component, though this conclusion is mostly based on objects with a non-dominant AGN component. We conclude that our findings do not provide evidence that the presence of AGN affects the star formation process in the host galaxy, but rather that the two phenomena occur simultaneously over a wide range of luminosities.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

We review our knowledge of the most basic properties of the AGN obscuring region - its location, scale, symmetry, and mean covering factor - and discuss new evidence on the distribution of covering factors in a sample of ~9000 quasars with WISE, UKIDSS, and SDSS photometry. The obscuring regions of AGN may be in some ways more complex than we thought - multi-scale, not symmetric, chaotic - and in some ways simpler - with no dependence on luminosity, and a covering factor distribution that may be determined by the simplest of considerations - e.g. random misalignments.

Using Herschel data from the deepest SPIRE and PACS surveys (HerMES and PEP) in COSMOS and GOODS (N+S), we examine the dust properties of IR-luminous (L_IR>10^10 L_sun) galaxies at 0.1<z<2 and determine how these evolve with cosmic time. The unique angle of this work is the rigorous analysis of survey selection effects, making this the first study of the star-formation-dominated, IR-luminous population within a framework almost entirely free of selection biases. We find that IR-luminous galaxies have SEDs with broad far-IR peaks characterised by cool/extended dust emission and average dust temperatures in the 25-45K range. Hot (T>45K) SEDs and cold (T<25K), cirrus-dominated SEDs are rare, with most sources being within the range occupied by warm starbursts such as M82 and cool spirals such as M51. We observe a luminosity-temperature (L-T) relation, where the average dust temperature of log[L_IR/L_sun]=12.5 galaxies is about 10K higher than that of their log[L_IR/L_sun]=10.5 counterparts. However, although the increased dust heating in more luminous systems is the driving factor behind the L-T relation, the increase in dust mass and/or starburst size with luminosity plays a dominant role in shaping it. Our results show that the dust conditions in IR-luminous sources evolve with cosmic time: at high redshift, dust temperatures are on average up to 10K lower than what is measured locally. This is manifested as a flattening of the L-T relation, suggesting that (U)LIRGs in the early Universe are typically characterised by a more extended dust distribution and/or higher dust masses than local equivalent sources. Interestingly, the evolution in dust temperature is luminosity dependent, with the fraction of LIRGs with T<35K showing a 2-fold increase from z~0 to z~2, whereas that of ULIRGs with T<35K shows a 6-fold increase.

We have constructed an extended halo model (EHM) which relates the total stellar mass and star-formation rate (SFR) to halo mass (M_h). An empirical relation between the distribution functions of total stellar mass of galaxies and host halo mass, tuned to match the spatial density of galaxies over 0<z<2 and the clustering properties at z~0, is extended to include two different scenarios describing the variation of SFR on M_h. We also present new measurements of the redshift evolution of the average SFR for star-forming galaxies of different stellar mass up to z=2, using data from the Herschel Multi-tiered Extragalactic Survey (HerMES) for infrared-bright galaxies. Combining the EHM with the halo accretion histories from numerical simulations, we trace the stellar mass growth and star-formation history in halos spanning a range of masses. We find that: (1) The intensity of the star-forming activity in halos in the probed mass range has steadily decreased from z~2 to 0; (2) At a given epoch, halos in the mass range between a few times 10^{11} M_Sun and a few times 10^{12} M_Sun are the most efficient at hosting star formation; (3) The peak of SFR density shifts to lower mass halos over time; (4) Galaxies that are forming stars most actively at z~2 evolve into quiescent galaxies in today's group environments, strongly supporting previous claims that the most powerful starbursts at z~2 are progenitors of today's elliptical galaxies.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

(abridged) Far-infrared Herschel photometry from the PEP and HerMES programs is combined with ancillary datasets in the GOODS-N, GOODS-S, and COSMOS fields. Based on this rich dataset, we reproduce the restframe UV to FIR ten-colors distribution of galaxies using a superposition of multi-variate Gaussian modes. The median SED of each mode is then fitted with a modified version of the MAGPHYS code that combines stellar light, emission from dust heated by stars and a possible warm dust contribution heated by an AGN. The defined Gaussian grouping is also used to identify rare sources. The zoology of outliers includes Herschel-detected ellipticals, very blue z~1 Ly-break galaxies, quiescent spirals, and torus-dominated AGN with star formation. Out of these groups and outliers, a new template library is assembled, consisting of 32 SEDs describing the intrinsic scatter in the restframe UV-to-submm colors of infrared galaxies. This library is tested against L(IR) estimates with and without Herschel data included, and compared to eight other popular methods often adopted in the literature. When implementing Herschel photometry, these approaches produce L(IR) values consistent with each other within a median absolute deviation of 10-20%, the scatter being dominated more by fine tuning of the codes, rather than by the choice of SED templates. Finally, the library is used to classify 24 micron detected sources in PEP GOODS fields. AGN appear to be distributed in the stellar mass (M*) vs. star formation rate (SFR) space along with all other galaxies, regardless of the amount of infrared luminosity they are powering, with the tendency to lie on the high SFR side of the "main sequence". The incidence of warmer star-forming sources grows for objects with higher specific star formation rates (sSFR), and they tend to populate the "off-sequence" region of the M*-SFR-z space.

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 combine far-infrared photometry from Herschel (PEP/HERMES) with deep mid-infrared spectroscopy from Spitzer to investigate the nature and the mass assembly history of a sample of 31 Luminous and Ultraluminous Infrared Galaxies at z~1 and 2 selected in GOODS-S with 24 $\mu$m fluxes between 0.2 and 0.5 mJy. We model the data with a self-consistent physical model (GRASIL) which includes a state-of-the-art treatment of dust extinction and reprocessing. We find that all of our galaxies appear to require massive populations of old (>1 Gyr) stars and, at the same time, to host a moderate ongoing activity of SF (SFR < 100 M$_{\odot}$/yr). The bulk of the stars appear to have been formed a few Gyr before the observation in essentially all cases. Only five galaxies of the sample require a recent starburst superimposed on a quiescent star formation history (SFH). We also find discrepancies between our results and those based on optical-only SED fitting for the same objects; by fitting their observed Spectral Energy Distributions with our physical model we find higher extinctions (by $\Delta$A_{V} ~ 0.81 and 1.14) and higher stellar masses (by $\Delta$Log(M*) ~ 0.16 and 0.36 dex) for z~1 and z~2 (U)LIRGs, respectively. The stellar mass difference is larger for the most dust obscured objects. We also find lower SFRs than those computed from L_{IR} using the Kennicutt relation due to the significant contribution to the dust heating by intermediate-age stellar populations through 'cirrus' emission (~73% and ~66% of total L_{IR} for z~1 and z~2 (U)LIRGs, respectively).

We investigate the range of covering factors (determined from the ratio of IR to UV/optical luminosity) seen in luminous type 1 quasars using a combination of data from the WISE, UKIDSS and SDSS surveys. Accretion disk (UV/optical) and obscuring dust (IR) luminosities are measured via the use of a simple three component SED model. We use these estimates to investigate the distribution of covering factors and its relationship to both accretion luminosity and IR SED shape. The distribution of covering factors (f_C) is observed to be log-normal, with a bias-corrected mean of <log10 f_C>=-0.41 and standard deviation of 0.2. The fraction of IR luminosity emitted in the near-IR (1--5 micron) is found to be high (~40 per cent), and strongly dependant on covering factor.

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

We use data from the Disc Emission via a Bias-free Reconnaissance in the Infrared/Submillimetre (DEBRIS) survey, taken at 100 um with the Photoconductor Array Camera and Spectrometer instrument on board the Herschel Space Observatory, to make a cosmic variance independent measurement of the extragalactic number counts. These data consist of 323 small-area mapping observations performed uniformly across the sky, and thus represent a sparse sampling of the astronomical sky with an effective coverage of ~2.5 deg^2. We find our cosmic variance independent analysis to be consistent with previous count measurements made using relatively small area surveys. Furthermore, we find no statistically significant cosmic variance on any scale within the errors of our data. Finally, we interpret these results to estimate the probability of galaxy source confusion in the study of debris discs.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

We present spectroscopic observations for a sample of 36 Herschel-SPIRE 250-500um selected galaxies (HSGs) at 2<z<5 from the Herschel Multi-tiered Extragalactic Survey (HerMES). Redshifts are confirmed as part of a large redshift survey of Herschel-SPIRE-selected sources covering ~0.93deg^2 in six extragalactic legacy fields. Observations were taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). Precise astrometry, needed for spectroscopic follow-up, is determined by identification of counterparts at 24um or 1.4GHz using a cross-identification likelihood matching method. Individual source luminosities range from log(L_IR/Lsun)=12.5-13.6 (corresponding to star formation rates 500-9000Msun/yr, assuming a Salpeter IMF), constituting some of the most intrinsically luminous, distant infrared galaxies yet discovered. We present both individual and composite rest-frame ultraviolet spectra and infrared spectral energy distributions (SEDs). The selection of these HSGs is reproducible and well characterized across large areas of sky in contrast to most z>2 HyLIRGs in the literature which are detected serendipitously or via tailored surveys searching only for high-z HyLIRGs; therefore, we can place lower limits on the contribution of HSGs to the cosmic star formation rate density at (7+-2)x10^(-3)Msun/yr h^3Mpc^(-3) at z~2.5, which is >10% of the estimated total star formation rate density (SFRD) of the Universe from optical surveys. The contribution at z~4 has a lower limit of 3x10^(-3)Msun/yr h^3 Mpc^(-3), ~>20% of the estimated total SFRD. This highlights the importance of extremely infrared-luminous galaxies with high star formation rates to the build-up of stellar mass, even at the earliest epochs.

We present Keck spectroscopic observations and redshifts for a sample of 767 Herschel-SPIRE selected galaxies (HSGs) at 250, 350, and 500um, taken with the Keck I Low Resolution Imaging Spectrometer (LRIS) and the Keck II DEep Imaging Multi-Object Spectrograph (DEIMOS). The redshift distribution of these SPIRE sources from the Herschel Multitiered Extragalactic Survey (HerMES) peaks at z=0.85, with 731 sources at z<2 and a tail of sources out to z~5. We measure more significant disagreement between photometric and spectroscopic redshifts (<delta_z>/(1+z)>=0.29) than is seen in non-infrared selected samples, likely due to enhanced star formation rates and dust obscuration in infrared-selected galaxies. We estimate that the vast majority (72-83%) of z<2 Herschel-selected galaxies would drop out of traditional submillimeter surveys at 0.85-1mm. We estimate the luminosity function and implied star-formation rate density contribution of HSGs at z<1.6 and find overall agreement with work based on 24um extrapolations of the LIRG, ULIRG and total infrared contributions. This work significantly increased the number of spectroscopically confirmed infrared-luminous galaxies at z>>0 and demonstrates the growing importance of dusty starbursts for galaxy evolution studies and the build-up of stellar mass throughout cosmic time. [abridged]

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

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 galaxies with log(M_*/M_sun) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deepest H-band Hubble Space Telescope survey of its type prior to the installation of WFC3. We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Survey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median log(M_*/M_sun) = 10.8 and z = 2.0, associated with 250 um sources detected with SNR > 3. We derive mean total IR luminosity log L_IR (L_sun) = 12.36 +/- 0.05 and corresponding star formation rate SFR_(IR+UV) = (280 +/- 40) M_sun/yr for these objects, and find them to have mean dust temperature T_dust ~ 35 K. We find that the SFR derived from the far-IR photometry combined with UV-based estimates of unobscured SFR for these galaxies is on average more than a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, although we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis to measure the mean SFR in bins of stellar mass. We obtain detections at the 2-4 sigma level at SPIRE wavelengths for samples with log(M_*/M_sun) > 10. In contrast to the Herschel detected GNS galaxies, we find that estimates of SFR_(IR+UV) for the stacked samples are comparable to those derived from extinction corrected UV emission, although the uncertainties are large. We find evidence for an increasing fraction of dust obscured star formation with stellar mass, finding SFR_IR/SFR_UV \propto M_*^{0.7 +/- 0.2}, which is likely a consequence of the mass--metallicity relation.

We report Herschel SPIRE (250, 350, and 500 micron) detections of 32 quasars with redshifts 0.5 < z < 3.6 from the Herschel Multi-tiered Extragalactic Survey. These sources are from a MIPS 24 micron flux-limited sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field permit construction of the rest-frame Spectral Energy Distributions (SEDs)from ultraviolet to the mid-infrared for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show dust temperatures in the range of 25K to 60K, with a mean of 34K. The FIR luminosities range from 10^{11.3} to 10^{13.5} Lsun, qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3--20 micron, rest-frame), and the bolometric luminosities derived using the 5100 A index may be underestimated for these FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities for quasars with higher near-infrared luminosities.

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

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

We use the Spitzer Wide-area InfraRed Extragalactic Legacy Survey (SWIRE) to explore the specific star-formation activity of galaxies and their evolution near the peak of the cosmic far-infrared (FIR) background at 70 and 160um. We use a stacking analysis to determine the mean FIR properties of well defined subsets of galaxies at flux levels well below the FIR catalogue detection limits of SWIRE and other Spitzer surveys. We tabulate the contribution of different subsets of galaxies to the FIR background at 70um and 160um. These long wavelengths provide a good constraint on the bolometric, obscured emission. The large area provides good constraints at low z and in finer redshift bins than previous work. At all redshifts we find that the specific FIR Luminosity (sLFIR) decreases with increasing mass, following a trend L_FIR/M* propto M_* ^beta with beta =-0.38\pm0.14. This is a more continuous change than expected from the {Delucia2007} semi-analytic model suggesting modifications to the feedback prescriptions. We see an increase in the sLFIR by about a factor of ~100 from 0<z<2 and find that the sLFIR evolves as (1+z)^alpha with alpha=4.4\pm0.3 for galaxies with 10.5 < log M*/Msun < 12. This is considerably steeper than the {Delucia2007} semi-analytic model (alpha \sim 2.5). When separating galaxies into early and late types on the basis of the optical/IR spectral energy distributions we find that the decrease in sLFIR with stellar mass is stronger in early type galaxies (beta ~ -0.46), while late type galaxies exhibit a flatter trend (beta \sim -0.15). The evolution is strong for both classes but stronger for the early type galaxies. The early types show a trend of decreasing strength of evolution as we move from lower to higher masses while the evolution of the late type galaxies has little dependence on stellar mass. We suggest that in late-type galaxies we are seeing a consistently declining sSFR..

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We describe the generation of single-band point source catalogues from submillimetre Herschel-SPIRE observations taken as part of the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey (HerMES). Flux densities are found by means of peak-finding and the fitting of a Gaussian point-response function. With highly-confused images, careful checks must be made on the completeness and flux density accuracy of the detected sources. This is done by injecting artificial sources into the images and analysing the resulting catalogues. Measured flux densities at which 50 per cent of injected sources result in good detections at (250, 350, 500) {\mu}m range from (11.6, 13.2, 13.1) mJy to (25.7, 27.1, 35.8) mJy, depending on the depth of the observation (where a `good' detection is taken to be one with positional offset less than one full-width half-maximum of the point-response function, and with the measured flux density within a factor of 2 of the flux density of the injected source). This paper acts as a reference for the 2010 July HerMES public data release.

For the first time, we investigate the X-ray/infrared (IR) correlation for star-forming galaxies at z~1, using SPIRE submm data from the recently-launched Herschel Space Observatory and deep X-ray data from the 2Ms Chandra deep field north (CDFN) survey. We examine the X-ray/IR correlation in the soft X-ray (SX, 0.5-2 keV) and hard X-ray (HX, 2-10 keV) bands by comparing our z~1 SPIRE-detected star-forming galaxies (SFGs) to equivalently IR-luminous (L_IR >10^10 L_sun) samples in the local/low redshift Universe. Our results suggest that the X-ray/IR properties of the SPIRE SFGs are on average similar to those of their local counterparts, as we find no evidence for evolution in the L_SX/L_IR and L_HX/L_IR ratios with redshift. We note however, that at all redshifts, both L_SX/L_IR and L_HX/L_IR are strongly dependent on IR luminosity, with luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs,L_IR >10^11 L_sun) having up to an order of magnitude lower values than normal infrared galaxies (L_IR <10^11 L_sun). We derive a L_SX-L_IR relation and confirm the applicability of an existing L_HX-L_IR relation for both local and distant LIRGs and ULIRGs, consistent with a scenario where X-ray luminosity is correlated with the star-formation rate (SFR).

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

As part of the Herschel Multi-tiered Extragalactic Survey we have investigated the rest-frame far-infrared (FIR) properties of a sample of more than 4800 Lyman Break Galaxies (LBGs) in the Great Observatories Origins Deep Survey North field. Most LBGs are not detected individually, but we do detect a sub-sample of 12 objects at 0.7 < z < 1.6 and one object at z ~ 2.0. The ones detected by Herschel SPIRE have redder observed NUV-U and U-R colors than the others, while the undetected ones have colors consistent with average LBGs at z > 2.5. The UV-to-FIR spectral energy distributions of the objects detected in the rest-frame FIR are investigated using the code CIGALE to estimate physical parameters. We find that LBGs detected by SPIRE are high mass, luminous infrared galaxies. It appears that LBGs are located in a triangle-shaped region in the A_FUV vs. Log L_FUV diagram limited by A_FUV=0 at the bottom and by a diagonal following the temporal evolution of the most massive galaxies from the bottom-right to the top-left of the diagram. This upper envelop can be used as upper limits for the UV dust attenuation as a function of L_FUV}. The limits of this region are well explained using a closed-box model, where the chemical evolution of galaxies produces metals, which in turn lead to higher dust attenuation when the galaxies age.

We report the detection of CO 5-4, 3-2, and 1-0 emission in the strongly lensed, Herschel/SPIRE-selected submillimeter galaxy (SMG) HLSW-01 at z=2.9574+/-0.0001, using the Plateau de Bure Interferometer, the Combined Array for Research in Millimeter-wave Astronomy, and the Green Bank Telescope. The observations spatially resolve the molecular gas into four lensed images with a maximum separation of ~9", and reveal the internal gas dynamics in this system. We derive lensing-corrected CO line luminosities of L'(CO 1-0) = (4.17+/-0.41), L'(CO 3-2) = (3.96+/-0.20) and L'(CO 5-4) = (3.45+/-0.20) x 10^10 (mu_L/10.9)^-1 Kkm/s pc^2, corresponding to luminosity ratios of r_31 = 0.95+/-0.10, r_53 = 0.87+/-0.06, and r_51 = 0.83+/-0.09. This suggests a total molecular gas mass of Mgas = 3.3 x 10^10 (alpha_CO/0.8) (mu_L/10.9)^-1 Msun. The gas mass, gas mass fraction, gas depletion timescale, star formation efficiency, and specific star formation rate are typical for an SMG. The velocity structure of the gas reservoir suggests that the brightest two lensed images are dynamically resolved projections of the same dust-obscured region in the galaxy that are kinematically offset from the unresolved fainter images. The resolved kinematics appear consistent with the complex velocity structure observed in major, `wet' (i.e., gas-rich) mergers. Major mergers are commonly observed in SMGs, and are likely to be responsible for fueling their intense starbursts at high gas consumption rates. This study demonstrates the level of detail to which galaxies in the early universe can be studied by utilizing the increase in effective spatial resolution and sensitivity provided by gravitational lensing.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

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

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

(abridged) We present the first study of the farIR properties of high redshift, radio-selected ULIRGs using deep observations obtained with SPIRE from the Herschel Multi-tiered Extragalactic Survey (HerMES). These galaxies span a large range of 850um fluxes from submillimetre-luminous ~10mJy "SCUBA galaxies -- SMGs" to ~1.5mJy from stacked SCUBA non-detections, thus likely representing a complete distribution of ULIRG spectral energy distributions. From Keck spectroscopic surveys in the Lockman-North field we identified a sample of 31 SMGs and 37 submillimetre-faint, optically-faint radio galaxies (OFRGs), all with radio-inferred IR luminosities >10^12 Lsun. These galaxies were cross-identified with SPIRE 250, 350 and 500um catalogs based on fluxes extracted at 24um positions in the SWIRE survey, yielding a sample of more than half of the galaxies well detected in at least two of the SPIRE bandpasses. By fitting greybody dust models to the SPIRE photometry together with SCUBA 850um measurements, we infer dust temperatures and far-infrared luminosities. The OFRGs detected by SPIRE have median <T_d>= 41+-5 K and the SMGs have <T_d>= 34+-5 K, both in reasonable agreement with previous (pre-Herschel) estimates, reaffirming that the local FIR/radio correlation holds (at least for this subset of high-z ULIRGs) at high redshift. Our observations firstly confirm that a substantial fraction of OFRGs exhibit large infrared luminosities corresponding to SFRs of ~400 Msun/yr. The SPIRE observations secondly confirm the higher dust temperatures for these OFRGs than similarly selected SMGs, consistent with early predictions of the submm-faint radio populations. Our observations also clearly confirm the large infrared luminosities of most SMGs selected with S850um>5 mJy and radio and strong 24um detections, corresponding to SFRs of ~700 Msun/yr.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

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

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

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

We present new radial velocities from AAOmega on the Anglo-Australian Telescope for 307 galaxies (b_J < 19.5) in the region of the rich cluster Abell 1386. Consistent with other studies of galaxy clusters that constitute sub-units of superstructures, we find that the velocity distribution of A1386 is very broad (21,000--42,000 kms^-1, or z=0.08--0.14) and complex. The mean redshift of the cluster that Abell designated as number 1386 is found to be ~0.104. However, we find that it consists of various superpositions of line-of-sight components. We investigate the reality of each component by testing for substructure and searching for giant elliptical galaxies in each and show that A1386 is made up of at least four significant clusters or groups along the line of sight whose global parameters we detail. Peculiar velocities of brightest galaxies for each of the groups are computed and found to be different from previous works, largely due to the complexity of the sky area and the depth of analysis performed in the present work. We also analyse A1386 in the context of its parent superclusters: Leo A, and especially the Sloan Great Wall. Although the new clusters may be moving toward mass concentrations in the Sloan Great Wall or beyond, many are most likely not yet physically bound to it.

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

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

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

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

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

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

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

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

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

One of the biggest problems faced by current and next-generation astronomical surveys is trying to produce large numbers of accurate cross identifications across a range of wavelength regimes with varying data quality and positional uncertainty. Until recently simple spatial "nearest neighbour" associations have been sufficient for most applications. However as advances in instrumentation allow more sensitive images to be made the rapid increase in the source density has meant that source confusion across multiple wavelengths is a serious problem. The field of far-IR and sub-mm astronomy has been particularly hampered by such problems. The poor angular resolution of current sub-mm and far-IR instruments is such that in a lot of cases there are multiple plausible counterparts for each source at other wavelengths. Here we present a new automated method of producing associations between sources at different wavelengths using a combination of spatial and SED information set in the Bayesian framework presented by Budavari & Szalay (2008). Testing of the technique is performed on both simulated catalogues of sources from GaLICS and real data from multi-wavelength observations of the SXDF. It is found that a single figure of merit, the Bayes factor, can be effectively used to describe the confidence in the match. Further applications of this technique to future Herschel datasets are discussed.

[abridged] We present a statistical exploration of the parameter space of the De Lucia and Blaizot version of the Munich semi-analytic model built upon the millennium dark matter simulation. This is achieved by applying a Monte Carlo Markov Chain method to constrain the 6 free parameters that define the stellar and black-hole mass functions at redshift zero. The model is tested against three different observational data sets, including the galaxy K-band luminosity function, B-V colours, and the black hole-bulge mass relation, separately and combined, to obtain mean values, confidence limits and likelihood contours for the best fit model. Using each observational data set independently, we discuss how the SA model parameters affect each galaxy property and to what extent the correlations between them can lead to improved understandings of the physics of galaxy formation. When all the observations are combined, we find reasonable agreement between the majority of the previously published parameter values and our confidence limits. However, the need to suppress dwarf galaxy formation requires the strength of the supernova feedback to be significantly higher in our best-fit solution than in previous work. To balance this, we require the feedback to become ineffective in halos of lower circular velocity than before, so as to permit the formation of sufficient high-luminosity galaxies: unfortunately, this leads to an excess of galaxies around L*. Although the best-fit is formally consistent with the data, there is no region of parameter space that reproduces the shape of galaxy luminosity function across the whole magnitude-range. We discuss modifications to the semi-analytic model that might simultaneously improve the fit to the observed luminosity function and reduce the reliance on excessive supernova feedback in small halos.

The advent of the Spitzer Space Telescope has revolutionized our understanding of the history of star formation and galaxy mass assembly in the Universe. The tremendous leap in sensitivity from previous mid-to-far IR missions has allowed Spitzer to perform deeper, and wider, surveys than previously possible at these wavelengths. In this brief review I highlight some of the key results to come out of these surveys, and the implications these have for current models of galaxy formation and evolution.

We have performed IRS low resolution 5-12 micron spectroscopy on a sample of galaxies selected to be at 3 distinct post-starburst evolutionary stages based on their optical spectral indices. The resulting IRS spectra show distinctive PAH emission line structures at 6.2, 7.7, 8.6 and 11.3 micron and little silicate absorption, indicative of ongoing star formation. However the PAH inter-line ratios, in particular the 11.3/6.2 micron and 7.7/6.2 micron ratio, show large variations. These variations are found to correlate with both time since the most recent starburst and AGN activity. We speculate that the evolution observed in these PAH ratios is related to an increase in AGN activity with time since star burst.

We present the final spectroscopic QSO catalogue from the 2dF-SDSS LRG and QSO (2SLAQ) Survey. This is a deep, 18<g<21.85 (extinction corrected), sample aimed at probing in detail the faint end of the broad line AGN luminosity distribution at z<2.6. The candidate QSOs were selected from SDSS photometry and observed spectroscopically with the 2dF spectrograph on the Anglo-Australian Telescope. This sample covers an area of 191.9 deg^2 and contains new spectra of 16326 objects, of which 8764 are QSOs, and 7623 are newly discovered (the remainder were previously identified by the 2QZ and SDSS surveys). The full QSO sample (including objects previously observed in the SDSS and 2QZ surveys) contains 12702 QSOs. The new 2SLAQ spectroscopic data set also contains 2343 Galactic stars, including 362 white dwarfs, and 2924 narrow emission line galaxies with a median redshift of z=0.22. We present detailed completeness estimates for the survey, based on modelling of QSO colours, including host galaxy contributions. This calculation shows that at g~21.85 QSO colours are significantly affected by the presence of a host galaxy up to redshift z~1 in the SDSS ugriz bands. In particular we see a significant reddening of the objects in g-i towards fainter g-band magnitudes. This reddening is consistent with the QSO host galaxies being dominated by a stellar population of age at least 2-3 Gyr. The full catalogue, including completeness estimates, is available on-line at http://www.2slaq.info/

We present an analysis of the small-to-intermediate scale clustering of samples of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey and the 2dF-SDSS LRG and QSO (2SLAQ) survey carefully matched to have the same rest-frame colours and luminosity. We study the spatial two-point auto-correlation function in both redshift-space and real-space of a combined sample of over 10,000 LRGs, which represent the most massive galaxies in the universe with stellar masses > 10^11 h^-1 M_sun and space densities 10^-4 h^-3 Mpc^-3. We find no significant evolution in the amplitude r_0 of the correlation function with redshift, but do see a slight decrease in the slope with increasing redshift over 0.19 < z < 0.55 and scales of 0.32 < r < 32 h^-1 Mpc. We compare our measurements with the predicted evolution of dark matter clustering and use the halo model to interpret our results. We find that our clustering measurements are inconsistent (>99.9% significance) with a passive model whereby the LRGs do not merge with one another; a model with a merger rate of 7.5 +/- 2.3% from z = 0.55 to z = 0.19 (i.e. an average rate of 2.4% Gyr^-1) provides a better fit to our observations. Our clustering and number density measurements are consistent with the hypothesis that the merged LRGs were originally central galaxies in different haloes which, following the merger of these haloes, merged to create a single Brightest Cluster Galaxy. In addition, we show that the small-scale clustering signal constrains the scatter in halo merger histories. When combined with measurements of the luminosity function, our results suggest that this scatter is sub-Poisson. While this is a generic prediction of hierarchical models, it has not been tested before.

We present a clustering analysis of Luminous Red Galaxies (LRGs) using nearly 9 000 objects from the final catalogue of the 2dF-SDSS LRG And QSO (2SLAQ) Survey. We measure the redshift-space two-point correlation function, xi(s), at the mean LRG redshift of z=0.55. A single power-law fits the deprojected correlation function, xi(r), with a correlation length of r_0=7.45+-0.35 Mpc and a power-law slope of gamma=1.72+-0.06 in the 0.4<r<50 Mpc range. But it is in the LRG angular correlation function that the strongest evidence for non-power-law features is found where a slope of gamma=-2.17+-0.07 is seen at 1<r<10 Mpc with a flatter gamma=-1.67+-0.03 slope apparent at r<~1 Mpc scales. We use the simple power-law fit to the galaxy xi(r) to model the redshift space distortions in the 2-D redshift-space correlation function, xi(sigma,pi). We fit for the LRG velocity dispersion, w_z, Omega_m and beta, where beta=Omega_m^0.6/b and b is the linear bias parameter. We find values of w_z=330kms^-1, Omega_m= 0.10+0.35-0.10 and beta=0.40+-0.05. These high redshift results, which incorporate the Alcock-Paczynski effect and the effects of dynamical infall, start to break the degeneracy between Omega_m and beta found in low-redshift galaxy surveys. This degeneracy is further broken by introducing an additional external constraint, the value of beta(z=0.1)=0.45 from 2dFGRS, and then considering the evolution of clustering from z~0 to z_LRG~0.55. With these combined methods we find Omega_m(z=0)=0.30+-0.15 and beta(z=0.55)=0.45+-0.05. Assuming these values, we find a value for b(z=0.55)=1.66+-0.35. We show that this is consistent with a simple ``high peaks'' bias prescription which assumes that LRGs have a constant co-moving density and their clustering evolves purely under gravity. [ABRIDGED]

We have combined optical data from the 2dF-SDSS Luminous Red Galaxy and QSO (2SLAQ) redshift survey with radio measurements from the 1.4 GHz VLA FIRST and NVSS surveys to identify a volume-limited sample of 391 radio galaxies at redshift 0.4<z<0.7. By determining an accurate radio luminosity function for early-type galaxies in this redshift range, we can investigate the cosmic evolution of the radio-galaxy population over a wide range in radio luminosity. The low-power radio galaxies in our LRG sample (those with 1.4 GHz radio luminosities in the range 10^{24} to 10^{25} W/Hz, corresponding to FR I radio galaxies in the local universe) undergo significant cosmic evolution over the redshift range 0<z<0.7, consistent with pure luminosity evolution of the form (1+z)^k where k=2.0+/-0.3. Our results appear to rule out (at the 6-7 sigma level) models in which low-power radio galaxies undergo no cosmic evolution. The most powerful radio galaxies in our sample (with radio luminosities above 10^{26} W/Hz) may undergo more rapid evolution over the same redshift range. The evolution seen in the low-power radio-galaxy population implies that the total energy input into massive early-type galaxies from AGN heating increases with redshift, and was roughly 50% higher at z~0.55 (the median redshift of the 2SLAQ LRG sample) than in the local universe.

The 2dF-SDSS LRG and QSO (2SLAQ) survey is a new survey of distant Luminous Red Galaxies (LRGs) and faint quasars selected from the Sloan Digital Sky Survey (SDSS) multi-color photometric data and spectroscopically observed using the 2dF instrument on the Anglo-Australian Telescope (AAT). In total, the 2SLAQ survey has measured over 11000 LRG redshifts, covering 180deg^2 of SDSS imaging data, from 87 allocated nights of AAT time. Over 90% of these galaxies are within the range 0.45<z<0.7 and have luminosities consistent with >=3L*. When combined with the lower redshift SDSS LRGs, the evolution in the luminosity function of these LRGs is fully consistent with that expected from a simple passive (luminosity) evolution model. This observation suggests that at least half of the LRGs seen at z~0.2 must already have more than half their stellar mass in place by z~0.6, i.e., our observations are inconsistent with a majority of LRGs experiencing a major merger in the last 6 Gyrs. However, some "frosting" (i.e., minor mergers) has taken place with ~5% of LRGs showing some evidence of recent and/or on-going star-formation, but it only contributes ~1% of their stellar mass.

We describe the construction of MegaZ-LRG, a photometric redshift catalogue of over one million luminous red galaxies (LRGs) in the redshift range 0.4 < z < 0.7 with limiting magnitude i < 20. The catalogue is selected from the imaging data of the Sloan Digital Sky Survey Data Release 4. The 2dF-SDSS LRG and Quasar (2SLAQ) spectroscopic redshift catalogue of 13,000 intermediate-redshift LRGs provides a photometric redshift training set, allowing use of ANNz, a neural network-based photometric-redshift estimator. The rms photometric redshift accuracy obtained for an evaluation set selected from the 2SLAQ sample is sigma_z = 0.049 averaged over all galaxies, and sigma_z = 0.040 for a brighter subsample (i < 19.0). The catalogue is expected to contain ~5 per cent stellar contamination. The ANNz code is used to compute a refined star/galaxy probability based on a range of photometric parameters; this allows the contamination fraction to be reduced to 2 per cent with negligible loss of genuine galaxies. The MegaZ-LRG catalogue is publicly available on the World Wide Web from http://www.2slaq.info .

We present a detailed investigation into the recent star formation histories of 5,697 Luminous Red Galaxies (LRGs) based on the Hdelta (4101A) and [OII] (3727A) lines. LRGs are luminous (L>3L*), galaxies which have been selected to have photometric properties consistent with an old, passively evolving stellar population. For this study we utilise LRGs from the recently completed 2dF-SDSS LRG and QSO survey (2SLAQ). Equivalent widths of the Hdelta and [OII] lines are measured and used to define three spectral types, those with only strong Hdelta absorption (k+a), those with strong [OII] in emission (em) and those with both (em+a). All other LRGs are considered to have passive star formation histories. The vast majority of LRGs are found to be passive (~80 per cent), however significant numbers of k+a (2.7 per cent), em+a (1.2 per cent) and em LRGs (8.6 per cent) are identified. An investigation into the redshift dependence of the fractions is also performed. A sample of SDSS MAIN galaxies with colours and luminosities consistent with the 2SLAQ LRGs is selected to provide a low redshift comparison. While the em and em+a fractions are consistent with the low redshift SDSS sample, the fraction of k+a LRGs is found to increase significantly with redshift. This result is interpreted as an indication of an increasing amount of recent star formation activity in LRGs with redshift. By considering the expected life time of the k+a phase, the number of LRGs which will undergo a k+a phase can be estimated. A crude comparison of this estimate with the predictions from semi-analytic models of galaxy formation shows that the predicted level of k+a and em+a activity is not sufficient to reconcile the predicted mass growth for massive early-types in a hierarchical merging scenario.

We present a spectroscopic survey of almost 15,000 candidate intermediate-redshift Luminous Red Galaxies (LRGs) brighter than i=19.8, observed with 2dF on the Anglo-Australian Telescope. The targets were selected photometrically from the Sloan Digital Sky Survey (SDSS) and lie along two narrow equatorial strips covering 180 sq deg. Reliable redshifts were obtained for 92% of the targets and the selection is very efficient: over 90% have redshifts between 0.45 and 0.8. More than 80% of the ~11,000 red galaxies have pure absorption-line spectra consistent with a passively-evolving old stellar population. The redshift, photometric and spatial distributions of the LRGs are described. The 2SLAQ data will be released publicly from mid-2006, providing a powerful resource for observational cosmology and the study of galaxy evolution.

We present new measurements of the luminosity function (LF) of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS) and the 2dF-SDSS LRG and Quasar (2SLAQ) survey. We have carefully quantified, and corrected for, uncertainties in the K and evolutionary corrections, differences in the colour selection methods, and the effects of photometric errors, thus ensuring we are studying the same galaxy population in both surveys. Using a limited subset of 6326 SDSS LRGs (with 0.17<z<0.24) and 1725 2SLAQ LRGs (with 0.5 <z<0.6), for which the matching colour selection is most reliable, we find no evidence for any additional evolution in the LRG LF, over this redshift range, beyond that expected from a simple passive evolution model. This lack of additional evolution is quantified using the comoving luminosity density of SDSS and 2SLAQ LRGs, brighter than M_r - 5logh = -22.5, which are 2.51+/-0.03 x 10^-7 L_sun Mpc^-3 and 2.44+/-0.15 x 10^-7 L_sun Mpc^-3 respectively (<10% uncertainty). We compare our LFs to the COMBO-17 data and find excellent agreement over the same redshift range. Together, these surveys show no evidence for additional evolution (beyond passive) in the LF of LRGs brighter than M_r - 5logh = -21 (or brighter than L*). We test our SDSS and 2SLAQ LFs against a simple ``dry merger'' model for the evolution of massive red galaxies and find that at least half of the LRGs at z=0.2 must already have been well-assembled (with more than half their stellar mass) by z=0.6. This limit is barely consistent with recent results from semi-analytical models of galaxy evolution.

We present a catalogue of galaxies in Abell 3653 from observations made with the 2dF spectrograph at the Anglo-Australian Telescope. Of the 391 objects observed, we find 111 are bone-fide members of Abell 3653. We show that the cluster has a velocity of cz = 32214 +/- 83 km/s (z=0.10738 +/- 0.00027), with a velocity dispersion typical of rich, massive clusters of sigma_{cz} = 880^{+66}_{-54}. We find that the cD galaxy has a peculiar velocity of 683 +/- 96 km/s in the cluster restframe - some 7sigma away from the mean cluster velocity, making it one of the largest and most significant peculiar velocities found for a cD galaxy to date. We investigate the cluster for signs of substructure, but do not find any significant groupings on any length scale. We consider the implications of our findings on cD formation theories.

We have used the 2dF instrument on the AAT to obtain redshifts of a sample of z<3, 18.0<g<21.85 quasars selected from SDSS imaging. These data are part of a larger joint programme: the 2dF-SDSS LRG and QSO Survey (2SLAQ). We describe the quasar selection algorithm and present the resulting luminosity function of 5645 quasars in 105.7 deg^2. The bright end number counts and luminosity function agree well with determinations from the 2dF QSO Redshift Survey (2QZ) data to g\sim20.2. However, at the faint end the 2SLAQ number counts and luminosity function are steeper than the final 2QZ results from Croom et al. (2004), but are consistent with the preliminary 2QZ results from Boyle et al. (2000). Using the functional form adopted for the 2QZ analysis, we find a faint end slope of beta=-1.78+/-0.03 if we allow all of the parameters to vary and beta=-1.45+/-0.03 if we allow only the faint end slope and normalization to vary. Our maximum likelihood fit to the data yields 32% more quasars than the final 2QZ parameterization, but is not inconsistent with other g>21 deep surveys. The 2SLAQ data exhibit no well defined ``break'' but do clearly flatten with increasing magnitude. The shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that derived from type I quasars found in hard X-ray surveys. [Abridged]