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.

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 present the optical spectra of a sample of 34 SWIRE-CDFS sources observed with EFOSC2 on the ESO 3.6m Telescope. We have used the spectra and spectroscopic redshifts to validate our photometric redshift codes and SED template fitting methods. 12 of our sources are Infrared Luminous Galaxies. Of these, five belong to the class of ULIRGs and one to the class of HLIRGs with evidence of both an AGN and starburst component contributing to their extreme infrared luminosity for 3, starburst contributing for 1 and AGN contributing for 2 of them.

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.

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

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.

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.

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 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 submillimetre observations of the z=3.653 quasar SDSS160705+533558 together with data in the optical and infrared. The object is unusually bright in the far-IR and submm with an IR luminosity of ~10^14 L_sun. We ascribe this luminosity to a combination of AGN and starburst emission, with the starburst forming stars at a rate of a few thousand solar masses per year. Submillimetre Array (SMA) imaging observations with a resolution ~1" show that the submm (850 micron) emission is extended on scales of 10--35kpc and is offset from the optical position by ~10 kpc. This morphology is dissimilar to that found in submm galaxies, which are generally un- or marginally resolved on arcsecond scales, or submm-luminous AGN where the AGN lies at the peak of the submm or molecular emission. The simplest explanation is that the object is in the early stages of a merger between a gas rich galaxy, which hosts the starburst, and a gas-poor AGN-host galaxy, which is responsible for the quasar emission. It is also possible that jet induced star formation might contribute to the unusual morphology.

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.

We present the SWIRE Photometric Redshift Catalogue, 1025119 redshifts of unprecedented reliability and accuracy. Our method is based on fixed galaxy and QSO templates applied to data at 0.36-4.5 mu, and on a set of 4 infrared emission templates fitted to infrared excess data at 3.6-170 mu. The code involves two passes through the data, to try to optimize recognition of AGN dust tori. A few carefully justified priors are used and are the key to supression of outliers. Extinction, A_V, is allowed as a free parameter. We use a set of 5982 spectroscopic redshifts, taken from the literature and from our own spectroscopic surveys, to analyze the performance of our method as a function of the number of photometric bands used in the solution and the reduced chi^2. For 7 photometric bands the rms value of (z_{phot}-z_{spec})/(1+z_{spec}) is 3.5%, and the percentage of catastrophic outliers is ~1%. We discuss the redshift distributions at 3.6 and 24 mu. In individual fields, structure in the redshift distribution corresponds to clusters which can be seen in the spectroscopic redshift distribution. 10% of sources in the SWIRE photometric redshift catalogue have z >2, and 4% have z>3, so this catalogue is a huge resource for high redshift galaxies. A key parameter for understanding the evolutionary status of infrared galaxies is L_{ir}/L_{opt}, which can be interpreted as the specific star-formation rate for starbursts. For dust tori around Type 1 AGN, L_{tor}/L_{opt} is a measure of the torus covering factor and we deduce a mean covering factor of 40%.

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

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

This paper presents galaxy source counts at 24 microns in the six Spitzer Wide-field InfraRed Extragalactic (SWIRE) fields. The source counts are compared to counts in other fields, and to model predictions that have been updated since the launch of Spitzer. This analysis confirms a very steep rise in the Euclidean-normalized differential number counts between 2 mJy and 0.3 mJy. Variations in the counts between fields show the effects of sample variance in the flux range 0.5-10 mJy, up to 100% larger than Poisson errors. Nonetheless, a "shoulder" in the normalized counts persists at around 3 mJy. The peak of the normalized counts at 0.3 mJy is higher and narrower than most models predict. In the ELAIS N1 field, the 24 micron data are combined with Spitzer-IRAC data and five-band optical imaging, and these bandmerged data are fit with photometric redshift templates. Above 1 mJy the counts are dominated by galaxies at z less than 0.3. By 300 microJy, about 25% are between z ~ 0.3-0.8, and a significant fraction are at z ~ 1.3-2. At low redshifts the counts are dominated by spirals, and starbursts rise in number density to outnumber the spirals' contribution to the counts below 1 mJy.

We investigate the clustering of galaxies selected in the 3.6 micron band of the Spitzer Wide-area Infrared Extragalactic (SWIRE) legacy survey. The angular two-point correlation function is calculated for eleven samples with flux limits of S_3.6 > 4-400 mujy, over an 8 square degree field. The angular clustering strength is measured at >5-sigma significance at all flux limits, with amplitudes of A=(0.49-29)\times10^{-3} at one degree, for a power-law model, A\theta^{-0.8}. We estimate the redshift distributions of the samples using phenomological models, simulations and photometric redshifts, and so derive the spatial correlation lengths. We compare our results with the GalICS (Galaxies In Cosmological Simulations) models of galaxy evolution and with parameterized models of clustering evolution. The GalICS simulations are consistent with our angular correlation functions, but fail to match the spatial clustering inferred from the phenomological models or the photometric redshifts. We find that the uncertainties in the redshift distributions of our samples dominate the statistical errors in our estimates of the spatial clustering. At low redshifts (median z<0.5) the comoving correlation length is approximately constant, r_0=6.1\pm0.5h^{-1} Mpc, and then decreases with increasing redshift to a value of 2.9\pm0.3h^{-1} Mpc for the faintest sample, for which the median redshift is z=1. We suggest that this trend can be attributed to a decrease in the average galaxy and halo mass in the fainter flux-limited samples, corresponding to changes in the relative numbers of early- and late-type galaxies. However, we cannot rule out strong evolution of the correlation length over 0.5<z<1.

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from common interlopers. The search yields 100 z~3 (U-dropout) QSO candidates with 19<r'<22 over 11.7 deg^2 in the ELAIS-N1 (EN1) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared Extragalactic (SWIRE) Legacy Survey. The z~3 selection is reliable, with spectroscopic follow-up of 10 candidates confirming they are all QSOs at 2.83<z<3.44. We find that our z~4$ (g'-dropout) sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs at 3.50<z<3.89. Detailed simulations show our z~3 completeness to be ~80-90% from 3.0<z<3.5, significantly better than the ~30-80% completeness of the SDSS at these redshifts. The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end slope of beta=-1.42 +- 0.15, consistent with values measured at lower redshift. Therefore, we see no evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS QSO sample, we have now directly measured the space density of QSOs responsible for ~70% of the QSO UV luminosity density at z~3. We derive a maximum rate of HI photoionization from QSOs at z~3.2, Gamma = 4.8x10^-13 s^-1, about half of the total rate inferred through studies of the Ly-alpha forest. Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI in the intergalactic medium at z~3.

We study the ultraviolet to far-infrared (hereafter UV-to-IR) SEDs of a sample of intermediate redshift (0.2 < z < 0.7) UV-selected galaxies from the ELAIS-N1 and ELAIS-N2 fields by fitting a multi-wavelength dataset to a library of GRASIL templates. Star formation related properties of the galaxies are derived from the library of models by using the Bayesian statistics. We find a decreasing presence of galaxies with low attenuation and low total luminosity as redshift decreases, which does not hold for high total luminosity galaxies. In addition the dust attenuation of low mass galaxies increases as redshift decreases, and this trend seems to disappear for galaxies with M* > 10^11 M_sun. This result is consistent with a mass dependent evolution of the dust to gas ratio, which could be driven by a mass dependent efficiency of star formation in star forming galaxies. The specific star formation rates (SSFR) decrease with increasing stellar mass at all redshifts, and for a given stellar mass the SSFR decreases with decreasing redshift. The differences in the slope of the M*--SSFR relation found between this work and others at similar redshift could be explained by the adopted selection criteria of the samples which, for a UV selected sample, favours blue, star forming galaxies.

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

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

We present measurements of the spatial clustering of ultraluminous infrared galaxies in two redshift intervals, 1.5<z<2.0 and 2<z<3. Both samples cluster strongly, with r_0=14.40+/-1.99 h^-1 Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1 Mpc for the 1.5<z<2.0 sample, making them among the most biased galaxies at these epochs. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13 solar masses. We infer that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0.

We investigate claims that recent deep X-ray surveys are detecting starbursts at cosmologically interesting redshifts (z=0-1). We combine X-ray data from the 2Ms CDF-North with multi-wavelength observations from GOODS to build the SEDs (UV, optical, IR) of X-ray sources in this field. These are fit with model templates providing an estimate of the total IR luminosity of each source. We then exploit the correlation between IR and X-ray luminosities for star-forming galaxies, established in the local Universe, to select sources that are dominated by star-formation. This approach is efficient in discriminating normal galaxies from AGN over a wide range of SFRs. The resulting sample consists of 45 X-ray selected star-forming systems at a median redshift z~0.5, the majority of which (60%) are LIRGs or ULIRGs. This sample is least affected by incompleteness and AGN contamination and is well suited for cosmological studies. We quantify the X-ray evolution of these sources by constructing their differential, dN/dS, and comparing them with evolving luminosity function models. The results are consistent with luminosity evolution of the form (1+z)^{p} with p~2.4. This is similar to the evolution rate of star-forming galaxies selected at other wavelengths, suggesting that the deep X-ray surveys, are indeed finding the starburst galaxy population that drives the rapid evolution of the global SFR density in the range z~0-1. Our analysis also reveals a separate population of IR-faint X-ray sources. These include old galaxies but also systems that are X-ray luminous for their stellar mass compared to local E/S0. We argue that these may be post-starbursts that will, over time, become fainter at X-rays and will eventually evolve into early-type systems (e.g. E/S0).

We report results from a Spitzer GO-1 program of IRS spectroscopy of a large sample of Luminous Infrared Galaxies and quasars selected from the European Large Area ISO Survey (ELAIS) which have a wide multiwavelength coverage, including ISOCAM, ISOPHOT, IRAC and MIPS (from SWIRE), and optical photometry. We present the sample selection and results from the IRS spectroscopy.

We have carried out a deep (21.5<R<23.9) systematic survey of galaxies selected from SWIRE using GMOS to target X-ray sources, AGN and galaxies with photometric redshifts greater than 1. We have obtained a total of 198 redshifts, of which 132 have secure redshifts in the range 0.12<z<1.21. This sample, combined with the 111 redshifts targeted with WIYN in the same field and the 1 square degree survey of ELAIS-N1 with ACIS-I, will allow us to calibrate our photometric redshift techniques; characterise the star formation history of the universe between 0.5<z<1.2 where galaxy formation activity peaks; and explore the connection between AGN and star-formation activity.

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 construct rest-frame luminosity functions at 3.6, 4.5, 5.8, 8 and 24 microns over the redshift range 0<z<2 for galaxies and 0<z<4 for optical QSOs, using optical and infrared data from the Spitzer Wide-area InfraRed Extragalactic survey. The 3.6 and 4.5 micron galaxy LFs show evidence for moderate positive luminosity evolution up to z~1.5, consistent with the passive ageing of evolved stellar populations. Their comoving luminosity density was found to evolve passively, gradually increasing out to z~0.5-1 but flattening, or even declining, at higher redshift. Conversely, the 24 micron galaxy LF, which is more sensitive to obscured star formation and/or AGN activity, undergoes strong positive evolution, with the derived IR energy density and SFR density proportional to (1+z)^gamma with gamma=4.5^{+0.7}_{-0.6} and the majority of this evolution occurring since z~1. Optical QSOs, however, show positive luminosity evolution in all bands, out to the highest redshifts (3<z<4). Modelling as L* proportional to (1+z)^gamma gave gamma=1.3^{+0.1}_{-0.1} at 3.6 micron, gamma=1.0^{+0.1}_{-0.1} at 4.5 micron and stronger evolution at the longer wavelengths (5.8, 8 and 24 micron), of gamma~3. Comparison of the galaxy LFs to predictions from a semi-analytic model based on CDM indicate that an IMF skewed towards higher mass star formation in bursts compared to locally is preferred. As a result the currently inferred massive star formation rates in distant sub-mm sources may require substantial downwards revision.

We fit a parametric model comprising a mixture of multi-dimensional Gaussian functions to the 3.6 to 8um colour and optical photo-z distribution of galaxy populations in the ELAIS-N1 and Lockman Fields of SWIRE. For 16,698 sources in ELAIS-N1 we find our data are best modelled (in the sense of the Bayesian Information Criterion) by the sum of four Gaussian distributions or modes (C_a, C_b, C_c and C_d). We compare the fit of our empirical model with predictions from existing semi-analytic and phenomological models. We infer that our empirical model provides a better description of the mid-infrared colour distribution of the SWIRE survey than these existing models. This colour distribution test is thus a powerful model discriminator and complementary to comparisons of number counts. We use our model to provide a galaxy classification scheme and explore the nature of the galaxies in the different modes of the model. C_a consists of dusty star-forming systems such as ULIRG's. Low redshift late-type spirals are found in C_b, where PAH emission dominates at 8um. C_c consists of dusty starburst systems at intermediate redshifts. Low redshift early-type spirals and ellipticals dominate C_d. We thus find a greater variety of galaxy types than one can with optical photometry alone. Finally we develop a new technique to identify unusual objects, and find a selection of outliers with very red IRAC colours. These objects are not detected in the optical, but have very strong detections in the mid-infrared. These sources are modelled as dust-enshrouded, strongly obscured AGN, where the high mid-infrared emission may either be attributed to dust heated by the AGN or substantial star-formation. These sources have z_ph ~ 2-4, making them incredibly infrared luminous, with a L_IR ~ 10^(12.6-14.1) L_sun.

We present measurements of the spatial clustering of galaxies with stellar masses >10^11Msun, infrared luminosities >10^12 Lsun, and star formation rates >200Msun per year in two redshift intervals; 1.5<z<2.0 and 2<z<3. Both samples cluster very strongly, with spatial correlation lengths of r_0=14.40+/-1.99 h^-1Mpc for the 2<z<3 sample, and r_0=9.40+/-2.24 h^-1Mpc for the 1.5<z<2.0 sample. These clustering amplitudes are consistent with both populations residing in dark matter haloes with masses of ~6x10^13Msun, making them among the most biased galaxies at these epochs. We infer, from this and previous results, that a minimum dark matter halo mass is an important factor for all forms of luminous, obscured activity in galaxies at z>1, both starbursts and AGN. Adopting plausible models for the growth of DM haloes with redshift, then the haloes hosting the 2<z<3 sample will likely host the richest clusters of galaxies at z=0, whereas the haloes hosting the 1.5<z<2.0 sample will likely host poor to rich clusters at z=0. We conclude that ULIRGs at z>1 signpost stellar buildup in galaxies that will reside in clusters at z=0, with ULIRGs at increasing redshifts signposting the buildup of stars in galaxies that will reside in increasingly rich clusters.

We use a simple optical/infrared photometric selection of high redshift QSOs which identifies a Lyman Break in the optical and requires a red IRAC color to distinguish QSOs from common interlopers. We find 100 U-dropout (z=3) QSO candidates with r'<22 within 11.2 deg^2 in the ELAIS-N1 & ELAIS-N2 fields in the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Legacy Survey. Spectroscopy of 10 candidates shows that they are all QSOs with 2.83<z<3.44. We use detailed simulations which incorporate variations in QSO SEDs, IGM transmission and imaging depth to derive a completeness of 85-90% between 3.0<z<3.4. The resulting luminosity function is two magnitudes fainter than SDSS and, when combined with those data, gives a faint end slop \beta = 1.62 \pm 0.18, consistent with measurements at z<2 and steeper than initial measurements at the same redshift.

We discuss spectral energy distributions, photometric redshifts, redshift distributions, luminosity functions, source-counts and the far infrared to optical luminosity ratio for sources in the SWIRE Legacy Survey. The spectral energy distributions of selected SWIRE sources 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 5%. While further spectroscopic data are needed to confirm this result, the prospect of determining good photometric redshifts for the 2 million extragalactic objects in SWIRE is excellent. 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. Source-counts at 24, 70 and 160 mu are discussed, and luminosity functions at 3.6 and 24 mu are presented.

Using the large multi-wavelength data set in the chandra/SWIRE Survey (0.6 square degrees in the Lockman Hole), we show evidence for the existence of highly obscured (Compton-thick) AGN, estimate a lower limit to their surface density and characterize their multi-wavelength properties. Two independent selection methods based on the X-ray and infrared spectral properties are presented. The two selected samples contain 1) 5 X-ray sources with hard X-ray spectra and column densities > 10^24 cm-2, and 2) 120 infrared sources with red and AGN-dominated infrared spectral energy distributions (SEDs). We estimate a surface density of at least 25 Compton-thick AGN per square degree detected in the infrared in the chandra/SWIRE field of which ~40% show distinct AGN signatures in their optical/near-infrared SEDs, the remainings being dominated by the host-galaxy emission. Only ~33% of all Compton-thick AGN are detected in the X-rays at our depth (F(0.3-8 keV)>10^-15 erg/cm2/s. We report the discovery of two sources in our sample of Compton-thick AGN, SWIRE_J104409.95+585224.8 (z=2.54) and SWIRE_J104406.30+583954.1 (z=2.43), which are the most luminous Compton-thick AGN at high-z currently known. The properties of these two sources are discussed in detail with an analysis of their spectra, SEDs, luminosities and black-hole masses.

We present ground-based optical and Spitzer infrared imaging observations of the interacting galaxy UGC 10214, the "Tadpole Galaxy" (z = 0.0310), focusing on the star formation activity in the nuclear, disk, spiral arms and tidal tail regions. The major findings of this study are that the Tadpole is actively forming stars in the main disk outside of the nucleus and in the tidal plume, with an estimated mean star formation rate of ~2 to 4 M_sun/yr. The most prominent sites of mid-infrared emission define a "ring" morphology that, combined with the overall morphology of the system, suggest the interaction may belong to the rare class of off-center collisional ring systems that form both shock-induced rings of star formation and tidal plumes. The nuclear emission is solely powered by older stars, with little evidence for ongoing star formation at the center of the Tadpole. Extra-nuclear star formation accounts for >50% of the total star formation in the disk and spiral arms, featuring infrared-bright 'hot spots' that exhibit strong PAH emission, whose band strength is comparable to that of late-type star-forming disk galaxies. The tidal tail, which extends 2 arcmin (~75 kpc) into the intergalactic medium, is populated by super massive star clusters likely triggered by the galaxy-galaxy interaction that has distorted UGC 10214 into its current "tadpole" shape.

We report the detection of a 3.6sigma 350micron-selected source in the Bootes Deep Field. The source, the first Short-wavelength Submillimeter-selected Galaxy (SSG 1), was discovered as part of a blank field extragalactic survey using the 350micron-optimised Submillimeter High Angular Resolution Camera (SHARC II) at the Caltech Submillimeter Observatory. With multiwavelength photometry from NOAO-NDWFS (R and I band), FLAMEX (J and K_s), Spitzer (IRAC and MIPS) and the Westerbork 1.4GHz Deep Survey (radio upper limit), we are able to constrain the photometric redshift using different methods, all of which suggest a redshift of close to 1. In the absence of long-wavelength submillimeter data we use SED templates to infer that this source is an ultraluminous infrared galaxy (ULIRG) with a dust temperature of 30+/-5 K, occupying a region of luminosity-temperature space shared by modarate redshift ISO-selected ULIRGs (rather than high redshift SCUBA-selected SMGs). SHARC II can thus select SMGs with moderately "warm'' dust that might be missed in submillimeter surveys at longer wavelengths.

We report on a search for the optical counterparts of 175 um - selected sources from the Far-Infrared Background (FIRBACK) survey in the European Large Area ISO Survey (ELAIS) N2 field. Applying a likelihood ratio technique to optical catalogues from the Isaac Newton Telescope - Wide Field Survey (INT--WFS), we found optical identifications for 33 out of 55 FIRBACK sources in this field. These were then reassessed in the light of associations with the ELAIS final catalogue for the the N2 field, to yield a final set of 31 associations. We have investigated the nature of this population through a comparison of their observed spectral energy distributions with predictions from radiative transfer models which simulate the emission from both cirrus and starburst components. We find the far-infrared sources to be 80 per cent star bursting galaxies with their starburst component at a high optical depth. The resulting SEDs were used to estimate far-infrared luminosities, star formation rates, dust temperatures and dust masses. The N2 FIRBACK population is found to consist of four suspected ULIRGs, a number of LIRGs and a population of low redshift quiescently star forming galaxies. We also discuss the implications of these results for current evolutionary models.

Photometric redshifts can be routinely obtained to accuracies of better than 0.1 in Delta(z)/(1+z). However, the issue of dust extinction is one that has still not been well quantified. In this paper the success of two template-fitting photometric redshift codes (IMPZ and HYPERZ) at reliably returning Av in addition to redshift is explored. New data on the CNOC2 spectroscopic sample of 0.2<z<0.7 galaxies are presented. These data allow us to estimate Av values from the observed Balmer decrements. We also investigate whether the empirical value of gamma=0.44, the ratio between gas- and star- derived extinction, as determined by Calzetti (2001), is necessarily the best value for this sample. When comparing the two codes to the Balmer-derived Av, a correlation between the photometrically derived Av, Phot-Av, and the Balmer-Av is found. The correlation is improved when the empirical value of gamma=0.44 is allowed to vary. From least-squares-fitting the minimum in the reduced chi^2 distribution is found for gamma~0.25 +/- 0.2. For the sample of galaxies here, the factor of two difference in covering factor implied by the Calzetti ratio is found to be plausible. The CNOC2 galaxies with detected Balmer lines have some preference for an increased covering factor difference, which would perhaps imply they are undergoing more rapid, `bursty' star formation than the galaxies Calzetti used in her derivation.

We present the Final Analysis of the European Large Area ISO Survey (ELAIS) 15 $\mu$m observations, carried out with the ISOCAM instrument on board the Infrared Space Observatory (ISO). The data reduction method, known as LARI method, is based on a mathematical model of the detector's behaviour and was specifically designed for the detection of faint sources in ISO-CAM/PHOT data. The method is fully interactive and leads to very reliable and complete source lists. The resulting catalogue includes 1923 sources detected with $S/N > 5$ in the {0.5 -- 100 mJy} flux range and over an area of 10.85 {deg$^2$} split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from ISO data. This paper presents the concepts underlying the data reduction method together with its latest enhancements. The data reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.

We exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at least 2 IRAC bands and 988 sources detected at 24micron with MIPS brighter than 0.1 mJy. Of these, 102 IRAC and 59 MIPS sources have Chandra counterparts, out of a total of 122 X-ray sources present in the area with S(0.5-8 kev)>10^(-15) erg/cm^2/s. We have constructed SEDs for each source using data from the 4 IRAC wavebands, Chandra fluxes, and optical follow-up data in the wavebands U, g', r', i', Z, and H. We fit a number of spectral templates to the SEDs at optical and infrared wavelengths to determine photometric redshifts and spectral categories, and also make use of diagnostics based on the X-ray luminosities, hardness ratios, X-ray to infrared spectral slopes and optical morphologies. Although we have spectroscopic redshifts for only a minority of the Chandra sources, the available SEDs constrain the redshifts for most of the sample sources, which turn out to be typically at 0.5<z<2. We find that 39% of the Chandra sources are dominated by type-1 AGN emission, 23% display optical/IR spectra typical of type-2 AGNs, while the remaining 38% fraction show starburst-like or even normal galaxy spectra. Since we prove that all these galaxies are dominated by AGN emission in X-rays this brings the fraction of type-1 AGNs to be 80% of the type-2: even assuming that all the Chandra sources undetected by Spitzer are type-2 AGNs, the type-1 fraction would exceed 1/3 of the total population (abridged).

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.

In a recent paper Pello et al. reported observations of a faint galaxy, gravitationally lensed by the galaxy cluster Abell 1835. Deep J-band spectroscopy revealed a weak emission line near 1.34 microns, detected in two spectra with different central wavelengths. The line was interpreted as Lyman-alpha at redshift z=10.0. This interpretation is supported by the broad-band photometric spectral energy distribution, and by the location of the galaxy close to the lens critical line for this redshift. We have reanalysed the two spectra, just released from the data archive. Our analysis includes allowance for wavelength shifts due to transverse drift of the object in the slit. We do not detect a significant emission line at the reported location, or nearby, at either grating setting, nor in the combined spectrum. We provide a possible explanation for the reported detection as due to spurious positive flux introduced in the sky-subtraction stage as a result of variable hot pixels. We provide our final reduced 2D frame, and corresponding error array.

We present a mid-infrared analysis of 35 quasars with spectroscopic redshifts selected from the Spitzer Wide-area InfraRed Extragalactic Survey (SWIRE). We discuss their optical and mid-infrared (MIR) colors, and show that these quasars occupy well defined regions in MIR color-color space. We examine the issue of type-I AGN candidate selection in detail and propose new selection methods based on mid-IR colors. The available multi-band data allows us to construct two new, well-sampled quasar templates, covering wavelengths from the ultraviolet to the MIR.

We present the 90 micron luminosity function of the Final Analysis of the European Large Area ISO Survey (ELAIS), extending the sample size of our previous analysis (paper IV) by about a factor of 4. Our sample extends to z=1.1, around 50 times the comoving volume of paper IV, and 10^{7.7} < h^{-2}L/Lsun < 10^{12.5}. From our optical spectroscopy campaigns of the northern ELAIS 90 mircon survey (7.4 deg^2 in total, to S(90um)>70mJy), we obtained redshifts for 61% of the sample (151 redshifts) to B<21 identified at 7 microns, 15 microns, 20cm or with bright (B<18.5) optical identifications. The selection function is well-defined, permitting the construction of the 90 micron luminosity function of the Final Analysis catalogue in the ELAIS northern fields, which is in excellent agreement with our Preliminary Analysis luminosity function in the ELAIS S1 field from paper IV. The luminosity function is also in good agreement with the IRAS-based prediction of Serjeant & Harrison (2004), which if correct requires luminosity evolution of (1+z)^{3.4 +/- 1.0} for consistency with the source counts. This implies an evolution in comoving volume averaged star formation rate at z<~1 consistent with that derived from rest-frame optical and ultraviolet surveys.

We present a combined galaxy-quasar approach to template-fitting photometric redshift techniques and show the method to be a powerful one. The code (ImpZ) is presented, developed and applied to two spectroscopic redshift catalogues, namely the Isaac Newton Telescope Wide Angle Survey ELAIS N1 and N2 fields and the Chandra Deep Field North. In particular, optical size information is used to improve the redshift determination. The success of the code is shown to be very good with Delta z/(1+z) constrained to within 0.1 for 92 per cent of the galaxies in our sample. The extension of template-fitting to quasars is found to be reasonable with Delta z/(1+z) constrained to within 0.25 for 68 per cent of the quasars in our sample. Various template extensions into the far-UV are also tested.

We characterize the SWIRE galaxy populations in the SWIRE validation field within the Lockman Hole, based on the 3.6-24$\mu$ Spitzer data and deep U,g',r',r' optical imaging within an area ~1/3 sq. deg for ~16,000 Spitzer-SWIRE sources. The entire SWIRE survey will discover over 2.3 million galaxies at 3.6$\mu$m and almost 350,000 at 24$\mu$m; ~70,000 of these will be 5-band 3.6-24$\mu$ detections. The colors cover a broad range, generally well represented by redshifted spectral energy distributions of known galaxy populations, however significant samples of unusually blue objects in the [3.6-4.5]$\mu$m color are found, as well as many objects very red in the 3.6-24$\mu$m mid-IR. Nine of these are investigated and are interpreted as star-forming systems, starbursts and AGN from z=0.37 to 2.8, with luminosities from L$_{IR}$=10$^{10.3}$ to 10$^{13.7}$ L$_{\odot}$

We present the final band-merged ELAIS catalogue at 6.7, 15, 90, and 175 $\mu$m, and the associated data at u,g,r,i,z,J,H,K, and 20cm. The origin of the survey, infrared and radio observations, data-reduction and optical identifications are briefly reviewed, and a summary of the area covered, and completeness limit for each infrared band is given. A detailed discussion of the band-merging and optical association strategy is given. The total catalogues consists of 2860 sources. For extragalactic sources observed in 3 or more infrared bands, colour-colour diagrams are presented and discussed in terms of the contributing infrared populations. Spectral energy distributions are shown for selected sources and compared with cirrus, M82 and Arp220 starburst, and AGN dust torus models. Spectroscopic redshifts are tabulated, where available. For the N1 and N2 areas, the INT ugriz Wide Field Survey, permits photometric redshifts to be estimated for galaxies and quasars. These agree well with the spectroscopic redshifts, within the uncertainty of the photometric method. The redshift distribution is given for selected ELAIS bands and colour-redshift diagrams are discussed. There is a high proportion of ultraluminous infrared galaxies in the ELAIS Catalogue ($> 10 %$ of 15 $\mu$m sources), many with Arp220-like colours. 10 hyperluminous infrared galaxies and 10 EROs are found in the survey. The large numbers of ultraluminous galaxies imply very strong evolution in the star-formation rate between z = 0 and 1.