Action-angle coordinates are an essential tool for understanding the properties of the six dimensional phase space involved in orbits of stars in galactic potentials. A new method, which does not require specific knowledge of a generating function, is described, implemented and tested that calculates the actions of an orbit in an arbitrary potential of an integrable Hamiltonian given a set of Cartesian phase space points. The method chooses between the simple harmonic oscillator and isochrone potentials to fit the data using a Levenberg-Marquardt routine. An average is taken over the angle coordinates by calculating volumes in phase space using the metric free FiEstAS algorithm. The perfect ellipsoidal potential, with actions chosen a priori, is used to test the output of the algorithm, giving some results that agree within 1%. Minimisation of a sampling error is discussed along with an identification of a source of noise in the data.

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.

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.

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.

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

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

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.

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

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

We present deep submillimetre photometry for 14 galaxies at z=0.5 that are hosts of type 1a supernovae, with the aim of examining the evolution of dust mass and extinction in normal galaxies. We combine these results with our previous observations of 17 z~0.5 SN1a hosts to look for any evolution in the dust content of normal galaxies between z=0 and z=0.5. The average observed frame 850micron flux of SN1a hosts in the full sample, excluding 2 bright individually detected objects, is 0.44 +/- 0.22 mJy. This flux level is consistent with there being little or no evolution in the dust content, or optical extinction, of normal galaxies from z=0 to z=0.5. One galaxy, the host of SN1996cf, is detected individually, and we also present a deep HST STIS image for this object. It appears to be an edge on disk system, similar to the submm bright host of SN1997ey. We thus examine the dust properties of these and one other individually detected object. 450-to-850 micron flux ratios and limits suggest that the dust in the two brightest submm sources, SN1996cf and SN1997ey, is cold, T ~20K, implying that they contain a substantial mass of dust ~10^9 Msun. The presence of two bright (F850 >7mJy) submm sources at z~0.5 in a sample of ostensibly normal galaxies is surprising, and has important implications. It supports the idea that a substantial part of the Cosmic Infrared Background (CIB) may be produced at z<1, while also suggesting that 'foreground' objects such as these may be a significant 'contaminant' in submm surveys. Finally, we examine the overall submm luminosity distribution at z=0.5 implied by our results, and conclude that either there is substantial evolution in the submm luminosity function from z=0 to 0.5, or our submm detected sources are somehow not representative of the bulk of galaxies at this redshift.

We discuss optical associations, spectral energy distributions and photometric redshifts for SWIRE sources in the ELAIS-N1 area and the Lockman Validation Field. The band-merged IRAC (3.6, 4.5, 5.8 and 8.0 mu) and MIPS (24, 70, 160 mu) data have been associated with optical UgriZ data from the INT Wide Field Survey in ELAIS-N1, and with our own optical Ugri data in Lockman-VF. The spectral energy distributions of selected ELAIS sources in N1 detected by SWIRE, most with spectroscopic redshifts, are modelled in terms of a simple set of galaxy and quasar templates in the optical and near infrared, and with a set of dust emission templates (cirrus, M82 starburst, Arp 220 starburst, and AGN dust torus) in the mid infrared. The optical data, together with the IRAC 3.6 and 4.5 mu data, have been used to determine photometric redshifts. For galaxies with known spectroscopic redshifts there is a notable improvement in the photometric redshift when the IRAC data are used, with a reduction in the rms scatter from 10% in (1+z) to 7%. The photometric redshifts are used to derive the 3.6 and 24 mu redshift distribution and to compare this with the predictions of models. For those sources with a clear mid infrared excess, relative to the galaxy starlight model used for the optical and near infrared, the mid and far infrared data are modelled in terms of the same dust emission templates. The proportions found of each template type are: cirrus 31%, M82 29%, Arp 220 10%, AGN dust tori 29%. The distribution of the different infrared sed types in the L_{ir}/L_{opt} versus L_{ir} plane, where L_{ir} and L_{opt} are the infrared and optical bolometric luminosities, is discussed.

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

We present deep sub-millimetre observations of seventeen galaxies at z=0.5, selected through being hosts of a type 1a supernova. Two galaxies are detected directly, and the sample is detected statistically with a mean 850 micron flux of 1.01mJy +/- 0.33mJy, which is 25% - 135% higher than locally. We infer that the mean value of A_v in normal galaxies at z=0.5 is comparable to or greater than the mean A_v in local normal galaxies, in agreement with galaxy chemical evolution models and indirect observational evidence. Scaling from the local value given by Rowan-Robinson (2003) gives a mean extinction at z=0.5 of A_v = 0.56 +/- 0.17. The dust in the brightest sub-mm object in our sample is best interpreted as normal `cirrus' dust similar to that seen locally. The detection rate of our sample suggests that some sources found in blank-field sub-mm surveys may not be high redshift starbursts, but rather cirrus galaxies at moderate redshifts and with lower star formation rates. Finally, an increase in host dust extinction with redshift may impact the cosmological results from distant supernova searches. This emphasizes the need to carefully monitor dust extinction when using type Ia supernovae to measure the cosmological parameters.

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

We present an extended analysis of the SCUBA observations of the Hubble Deep Field (HDF), expanding the areal coverage of the Hughes et al. 1998 study by a factor of ~1.8 and containing at least three further sources in addition to the five in that study. We also announce the public release of the reduced data products. The map is the deepest ever made in the sub-millimetre, obtained in excellent conditions (median 850um optical depth of 0.16). Two independent reductions were made, one with SURF and the other with a wholly algorithmic IDL analysis which we present in detail here. Of the three new sources, all appear to be at z>~0.9 and one is provisionally associated with an Extremely Red Object (I-K>5). There appears to be no significant cross-correlation signal between the 850um fluctuations and sources detected by ISOCAM, the VLA or Chandra, nor with Very Red Objects (I-K>4), nor quasars and quasar candidates in the HDF (notwithstanding a small number of individual weak candidate detections). This is consistent with interpretations where the 850um-selected galaxies are at higher redshifts than those currently probed by ISOCAM/VLA, and predominantly not Compton-thin AGN. There are only one or two compelling cases for the radio source being the sub-mm source. Nevertheless, most SCUBA-HDF point sources have a nearby radio source apparently well-separated from the sub-mm centroid.

(Abridged) We examine the power source of 41 local Ultraluminous Infrared Galaxies using archival infrared and optical photometry. We fit the observed Spectral Energy Distributions (SEDs) with starburst and AGN components; each component being drawn from a family of templates. We find all of the sample require a starburst, whereas only half require an AGN. In 90% of the sample the starburst provides over half the IR emission, with a mean fractional luminosity of 82%. When combined with other galaxy samples we find that starburst and AGN luminosities correlate over 6 decades in IR luminosity suggesting that a common factor governs both luminosities, plausibly the gas masses in the nuclear regions. We find that the mid-IR 7.7 micron line-continuum ratio is no indication of the starburst luminosity, or the fractional AGN luminosity, and therefore that this ratio is not a reliable diagnostic of the power source in ULIRGs. We propose that the scatter in the radio-IR correlation in ULIRGs is due to a skewed starburst IMF and/or relic relativistic electrons from a previous starburst, rather than contamination from an obscured AGN. We show that most ULIRGs undergo multiple starbursts during their lifetime, and by inference that mergers between more than two galaxies may be common amongst ULIRGs. Our results support the evolutionary model for ULIRGs proposed by Farrah et al 2001, where they can follow many different evolutionary paths of starburst and AGN activity in transforming merging spiral galaxies into elliptical galaxies, but that most do not go through an optical QSO phase. The lower level of AGN activity in our local sample than in z~1 HLIRGs implies that the two samples are distinct populations. We postulate that different galaxy formation processes at high-z are responsible for this difference.

We present maps, source lists, and number counts from the largest, unbiassed, extragalactic sub-mm survey so far undertaken with the SCUBA camera on the JCMT. Our maps cover 260 sq. arcmin, to a noise level S(850)=2.5 mJy/beam. We have reduced the data using both SURF, and our own pipeline which produces zero-footprint maps and noise images. The uncorrelated noise maps produced by the latter approach have allowed application of a maximum-likelihood method to measure the statistical significance of each peak, leading to properly quantified flux-density errors for all potential sources. We detect 19 sources with S/N > 4, 38 with S/N > 3.5, and 72 with S/N > 3. To assess completeness and the impact of source confusion we have applied our source extraction algorithm to a series of simulated images. The result is a new estimate of the sub-mm source counts in the flux-density range S(850)=5-15mJy, which we compare with other estimates, and with model predictions. Our estimate of the cumulative source count at S(850) > 8 mJy is 320 (+80,-100) per square degree. Assuming that the majority of sources have z > 1.5, the co-moving number density of high-z galaxies forming stars at a rate >1000 solar masses per year is 10^-5 per Mpc^3, with only a weak dependence on the precise redshift distribution. This number density corresponds to that of massive ellipticals with L > 3-4 L* at low redshift and is also the same as the co-moving number density of comparably massive, passively-evolving objects in the redshift band 1<z<2 inferred from recent surveys of extremely red objects. Thus the bright sub-mm sources can plausibly account for the formation of all present-day massive ellipticals. Improved z constraints, and a proper measurement of sub-mm clustering can refine or refute this picture.

We explore the relationship between the hard X-ray and sub-mm populations using deep Chandra observations of a large, contiguous SCUBA survey. In agreement with other recent findings, we confirm that the direct overlap is small. Of the 17 sub-mm sources detected in this field at 850 microns, only one is coincident with a Chandra source. The resulting limits imply that the majority of SCUBA sources are not powered by AGN, unless the central engine is obscured by Compton-thick material with a low (<1 per cent) scattered component. Furthermore, since Chandra detects only ~5 per cent of SCUBA sources, the typical obscuration would need to be almost isotropic. The X-ray upper limits are so strong that in most cases we can also rule out a starburst SED at low redshift, suggesting that the majority of SCUBA sources lie at z>1 even if they are purely starburst galaxies. Despite the low detection rate, we find evidence for strong angular clustering between the X-ray and sub-mm populations. The implication is that AGN and SCUBA sources trace the same large-scale structure but do not generally coincide. If bright sub-mm sources represent massive elliptical galaxies in formation, we suggest that (for a given galaxy) the major episode of star-formation must be distinct from the period of observable quasar activity.

The SCUBA 8-mJy survey is the largest submm extragalactic mapping survey undertaken to date, centred on the Lockman Hole and ELAIS N2 regions. Here, we present new 1.4-GHz imaging of the depth and resolution necessary to reliably identify radio counterparts for 18 of 30 submm sources. Armed with this greatly improved positional information, we present and analyse new optical, near-IR and XMM-Newton X-ray imaging to identify optical/IR host galaxies to half of the submm-selected sources in those fields. As many as 15% of the submm sources detected at 1.4 GHz are resolved by the 1.4'' beam and a further 25% have more than one radio counterpart, suggesting that radio and submm emission arise from extended starbursts and that interactions are common. We note that less than a quarter of the submm-selected sample would have been recovered by targeting optically faint radio sources, underlining the selective nature of such surveys. At least 60% of the radio-confirmed optical/IR host galaxies appear to be morphologically distorted; many are composite systems - red galaxies with relatively blue companions; just over half are found to be very red or extremely red; contrary to popular belief, most are sufficiently bright to be tackled with spectrographs on 8-m telescopes. We estimate the median redshift of the S(850um) ~ 8 mJy submm galaxy population: if the radio/far-IR correlation holds at high redshift, and our sample is unbiased, we derive a conservative median z >= 2.0, or >= 2.4 using spectral templates more representative of known submm galaxies [abridged].

We present the results of a multi-wavelength study of the 19 most significant sub-mm sources detected in the SCUBA 8-mJy survey. As described in Scott et al. (2001), this survey covers ~260 arcmin^2 using the sub-millimetre camera SCUBA, to a limiting source detection limit S(850um) ~ 8 mJy. One advantage of this relatively bright flux-density limit is that accurate astrometric positions are potentially achievable for every source using existing radio and/or mm-wave interferometers. However, an associated advantage is that SED-based redshift constraints should be more powerful than in fainter sub-mm surveys. Here we therefore exploit the parallel SCUBA 450um data, in combination with existing radio and ISO data at longer and shorter wavelengths to set constraints on the redshift of each source. We also analyse new and existing optical and near-infrared imaging of our SCUBA survey fields to select potential identifications consistent with these constraints. Our derived SED-based redshift constraints, and the lack of statistically significant associations with even moderately bright galaxies allow us to conclude that all 19 sources lie at z > 1, and at least half of them apparently lie at z > 2.

The IRAM Plateau de Bure mm interferometer and deep K-band imaging have been used to identify the brightest submm source detected in the Lockman field of the UK 8mJy SCUBA survey. The near infrared counterpart is an extended (20-30kpc), clumpy, and extremely red object. The spectral energy distribution suggests it to be a dusty star forming object at a redshift of about 3 (2-4). Its star formation rate and near-infrared properties are consistent with Lockman850.1 being a massive elliptical in formation.