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

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

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

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

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

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

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.

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

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

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

We use one of the deepest spectroscopic samples of broad line active galactic nuclei (AGN) currently available, extracted from the VIMOS VLT Deep Survey (VVDS), to compute MgII and CIV virial masses estimate of 120 super-massive black holes in the redshift range 1.0<z<1.9 and 2.6<z<4.3. We find that the mass-luminosity relation shows considerably enhanced dispersion towards low AGN luminosities (log L_bol ~ 45). At these luminosities, there is a substantial fraction of black holes accreting far below their Eddington limit (L_bol/L_Edd < 0.1), in marked contrast to what is generally found for AGN of higher luminosities. We speculate that these may be AGN on the decaying branch of their light-curves, well past their peak activity. This would agree with recent theoretical predictions of AGN evolution. In the electronic Appendix of this paper we publish an update of the VVDS type-1 AGN sample, including the first and most of the second epoch observations. This sample contains 298 objects of which 168 are new.

We investigate the redshift and luminosity evolution of the galaxy colour-density relation using the data from the First Epoch VIMOS-VLT Deep Survey (VVDS) on scales of R=5 h^(-1)Mpc up to redshift z ~ 1.5. While at lower redshift we confirm the existence of a steep colour-density relation, with the fraction of the reddest(/bluest) galaxies of the same luminosity increasing(/decreasing) as a function of density, this trend progressively disappears in the highest redshift bins investigated.Our results suggest the existence of an epoch (more remote for brighter galaxies) characterized by the absence of the colour-density relation on the R=5 h^(-1)Mpc scales investigated. The rest frame u*-g' colour-magnitude diagram shows a bimodal pattern in both low and high density environments up to redshift z ~ 1.5. We find that the bimodal distribution is not universal but strongly depends upon environment. Both the colour-density and colour-magnitude-density relations, on the R=5 h^(-1)Mpc scales, appear to be a transient, cumulative product of genetic and environmental factors that have been operating over at least a period of 9 Gyr. These findings support an evolutionary scenario in which star formation/gas depletion processes are accelerated in more luminous objects and in high density environments: star formation activity is progressively shifting with cosmic time towards lower luminosity galaxies (downsizing), and out of high density environments.

We investigate the redshift and luminosity evolution of the galaxy colour-density relation using the data from the First Epoch VIMOS-VLT Deep Survey (VVDS). The size (6582 galaxies), depth (I_AB<=24) and redshift sampling rate of the survey enable us to reconstruct the 3D galaxy environment on relatively local scales (R=5 Mpc) up to z~1.5. Particular attention has been devoted to calibrate a density reconstruction scheme, which factors out survey selection effects and reproduces in an unbiased way the underlying `real' galaxy environment. While at lower redshift we confirm the existence of a steep colour-density relation, with the fraction of the reddest(/bluest) galaxies of the same luminosity increasing(/decreasing) as a function of density, this trend progressively disappears in the highest redshift bins investigated. The rest frame u*-g' colour-magnitude diagram shows a bimodal pattern in both low and high density environments up to z~1.5. We find that the bimodal distribution is not universal but strongly depends upon environment: at lower redshifts the colour-magnitude diagrams in low and high density regions are significantly different while the progressive weakening of the colour-density relation causes the two bimodal distributions to nearly mirror each other in the highest redshift bin investigated. Both the colour-density and the colour-magnitude-density relations appear to be a transient, cumulative product of genetic and environmental factors operating over at least a period of 9 Gyr. These findings support an evolutionary scenario in which star formation/gas depletion processes are accelerated in more luminous objects and in high density environments: star formation activity is shifting with cosmic time towards lower luminosity (downsizing), and out of high density environments.

We present the type-1 active galactic nuclei (AGN) sample extracted from the VIMOS VLT Deep Survey first observations of 21000 spectra in 1.75 square degree. This sample, which is purely magnitude limited, free of morphological or color selection biases, contains 130 broad line AGN (BLAGN) spectra with redshift up to 5. Our data are divided into a wide (Iab < 22.5) and a deep (Iab < 24) subsample containing 56 and 74 objects respectively. Because of its depth and selection criteria, this sample is uniquely suited to study the population of faint type-1 AGN. Our measured surface density (~ 472 +- 48 BLAGN per square degree with Iab < 24) is significantly higher than that of any other optically selected sample of BLAGN with spectroscopic confirmation. By applying a morphological and color analysis to our AGN sample we find that: (1)~23% of the AGN brighter than Iab=22.5 are classified as extended; this percentage increases to ~42% for those with z < 1.6; (2) a non-negligible fraction of our BLAGN are lying close to the color space area occupied by stars in u*-g' versus g'-r' color-color diagram. This leads us to the conclusion that classical optical ultraviolet preselection technique, if employed at such deep magnitudes (Iab=22.5) in conjuction with a preselection of point-like sources, can miss miss up to ~35% of the AGN population. Finally, we present a composite spectrum of our sample of objects. While the continuum shape is very similar to that of the SDSS composite at short wavelengths, it is much redder than it at lambda > 3000 A. We interpret this as due to significant contamination from emission of the host galaxies, as expected from the faint absolute magnitudes sampled by our survey.

From the first epoch observations of the VVDS up to z=1.5 we have derived luminosity functions (LF) of different spectral type galaxies. The VVDS data, covering ~70% of the life of the Universe, allow for the first time to study from the same sample and with good statistical accuracy the evolution of the LFs by galaxy type in several rest frame bands from a purely magnitude selected sample. The magnitude limit of the VVDS allows the determination of the faint end slope of the LF with unprecedented accuracy. Galaxies have been classified in four spectral classes, using their colours and redshift, and LFs have been derived in the U, B, V, R and I rest frame bands from z=0.05 to z=1.5. We find a significant steepening of the LF going from early to late types. The M* parameter is significantly fainter for late type galaxies and this difference increases in the redder bands. Within each of the galaxy spectral types we find a brightening of M* with increasing redshift, ranging from =< 0.5 mag for early type galaxies to ~1 mag for the latest type galaxies, while the slope of the LF of each spectral type is consistent with being constant with redshift. The LF of early type galaxies is consistent with passive evolution up to z~1.1, while the number of bright early type galaxies has decreased by ~40% from z~0.3 to z~1.1. We also find a strong evolution in the normalization of the LF of latest type galaxies, with an increase of more than a factor 2 from z~0.3 to z~1.3: the density of bright late type galaxies in the same redshift range increases of a factor ~6.6. These results indicate a strong type-dependent evolution and identifies the latest spectral types as responsible for most of the evolution of the UV-optical luminosity function out to z=1.5.

We have computed the evolution of the rest-frame B-band luminosity function (LF) for bulge and disk-dominated galaxies since z=1.2. We use a sample of 605 spectroscopic redshifts with I_{AB}<24 in the Chandra Deep Field South from the VIMOS-VLT Deep Survey, 3555 galaxies with photometric redshifts from the COMBO-17 multi-color data, coupled with multi-color HST/ACS images from the Great Observatories Origin Deep Survey. We split the sample in bulge- and disk-dominated populations on the basis of asymmetry and concentration parameters measured in the rest-frame B-band. We find that at z=0.4-0.8, the LF slope is significantly steeper for the disk-dominated population (\alpha=-1.19 \pm 0.07) compared to the bulge-dominated population (\alpha=-0.53 \pm 0.13). The LF of the bulge-dominated population is composed of two distinct populations separated in rest-frame color: 68% of red (B-I)_{AB}>0.9 and bright galaxies showing a strongly decreasing LF slope \alpha=+0.55 \pm 0.21, and 32% of blue (B-I)_{AB}<0.9 and more compact galaxies which populate the LF faint-end. We observe that red bulge-dominated galaxies are already well in place at z~1, but the volume density of this population is increasing by a factor 2.7 between z~1 and z~0.6. It may be related to the building-up of massive elliptical galaxies in the hierarchical scenario. In addition, we observe that the blue bulge-dominated population is dimming by 0.7 magnitude between z~1 and z~0.6. Galaxies in this faint and more compact population could possibly be the progenitors of the local dwarf spheroidal galaxies.

We present photometric redshifts for an uniquely large and deep sample of 522286 objects with i'_{AB}<25 in the Canada-France Legacy Survey ``Deep Survey'' fields, which cover a total effective area of 3.2 deg^2. We use 3241 spectroscopic redshifts with 0<z<5 from the VIMOS VLT Deep Survey as a calibration to derive these photometric redshifts. We devise a robust calibration method which removes systematic trends in the photometric redshifts and significantly reduces the fraction of catastrophic errors. We use our unique spectroscopic sample to present a detailed assessment of the robustness of the photometric redshift sample. For a sample selected at i'_{AB}<24, we reach a redshift accuracy of \sigma_{\Delta z/(1+z)}=0.037 with \eta=3.7% of catastrophic error. The reliability of our photometric redshifts is lower for fainter objects: we find \sigma_{\Delta z/(1+z)}=0.029, 0.043 and \eta=1.7%, 5.4% for samples selected at i'_{AB}=17.5-22.5 and 22.5-24 respectively. We find that the photometric redshifts of starburst galaxies in our sample are less reliable: although these galaxies represent only 18% of the spectroscopic sample they are responsible for 54% of the catastrophic errors. We find an excellent agreement between the photometric and the VVDS spectroscopic redshift distributions at i'_{AB}<24. Finally, we compare the redshift distributions of i' selected galaxies on the four CFHTLS deep fields, showing that cosmic variance is already present on fields of 0.8 deg^2.

We measure the evolution of clustering for galaxies with different spectral types from 6495 galaxies with 17.5<=I_AB<=24 and measured spectroscopic redshift in the first epoch VIMOS-VLT Deep Survey. We classify our sample into 4 classes, based on the fit of well-defined galaxy spectral energy distributions on observed multi-color data. We measure the projected function wp(rp) and estimate the best-fit parameters for a power-law real-space correlation function. We find the clustering of early-spectral-type galaxies to be markedly stronger than that of late-type galaxies at all redshifts up to z<=1.2. At z~0.8, early-type galaxies display a correlation length r_0=4.8+/-0.9h^{-1}Mpc, while late types have r_0=2.5+/-0.4h^{-1}Mpc. The clustering of these objects increases up to r_0=3.42+/-0.7h^{-1}Mpc for z~1.4. The relative bias between early- and late-type galaxies within our magnitude-limited survey remains approximately constant with b~1.7-1.8 from z~=0.2 up to z~=1, with indications for a decrease at z>1.2, due to the growth in clustering of the star-forming population. We find similar results when splitting the sample into `red' and `blue' galaxies using the observed color bi-modality. When compared to the expected linear growth of mass fluctuations, a natural interpretation of these observations is that: (a) the assembly of massive early type galaxies is already mostly complete in the densest dark matter halos at z~=1; (b) luminous late-type galaxies are located in higher-density, more clustered regions of the Universe at z~=1.5 than at present, indicating that star formation activity is progressively increasing, going back in time, in the higher-density peaks that today are mostly dominated by old galaxies.

We measure the evolution of the galaxy Luminosity Function as a function of large-scale environment up to z=1.5 from the VIMOS-VLT Deep Survey (VVDS) first epoch data. The 3D galaxy density field is reconstructed using a sample of 6582 galaxies with 17.5 < I_{AB} < 24 and measured spectroscopic redshifts. We split the sample in four redshift bins up to z=1.5 and in under-dense and over-dense environments according to the average density contrast \delta=0. There is a strong dependence of the Luminosity Function (LF) with large-scale environment up to z=1.2: the LF shape is observed to have a steeper slope in under-dense environments. We find a continuous brightening of \Delta M* ~0.6 mag from z=0.25 to z=1.5 both in under-dense and over-dense environments. The rest-frame B-band luminosity density continuously increases in under-dense environments from z=0.25 to z=1.5 whereas its evolution in over-dense environments presents a peak at z~0.9. We interpret the peak by a complex interplay between the decrease of the star formation rate and the increasing fraction of galaxies at \delta>0 due to hierarchical growth of structures. As the environmental dependency of the LF shape is already present at least up to z=1.2, we therefore conclude that either the shape of the LF is imprinted very early on in the life of the Universe, a `nature' process, or that `nurture' physical processes shaping up environment relation have already been efficient earlier than a look-back time corresponding to 30% of the current age of the Universe.

With new generation spectrographs integral field spectroscopy is becoming a widely used observational technique. The Integral Field Unit of the VIsible Multi-Object Spectrograph on the ESO-VLT allows to sample a field as large as 54" x 54" covered by 6400 fibers coupled with micro-lenses. We are presenting here the methods of the data processing software developed to extract the astrophysical signal of faint sources from the VIMOS IFU observations. We focus on the treatment of the fiber-to-fiber relative transmission and the sky subtraction, and the dedicated tasks we have built to address the peculiarities and unprecedented complexity of the dataset. We review the automated process we have developed under the VIPGI data organization and reduction environment (Scodeggio et al. 2005), along with the quality control performed to validate the process. The VIPGI-IFU data processing environment is available to the scientific community to process VIMOS-IFU data since November 2003.

In this paper we present a new deep, wide-field near-infrared imaging survey. Our J- and K-band observations in four separate fields complement optical BVRI, ultraviolet and spectroscopic observations undertaken as part of the VIMOS-VLT deep survey (VVDS). In total, our survey spans ~400arcmis^2. Our catalogues are reliable in all fields to at least Kvega~20.75 and Jvega~21.50 (defined as the magnitude where object contamination is less than 10% and completeness greater than 90%). Taken together these four fields represents a unique combination of depth, wavelength coverage and area. We describe the complete data reduction process and outline a comprehensive series of tests carried out to characterise the reliability of the final catalogues. We compare the statistical properties of our catalogues with literature compilations. We find that our J- and K-selected galaxy counts are in good agreement with previously published works, as are our (J-K) versus K colour-magnitude diagrams. Stellar number counts extracted from our fields are consistent with a synthetic model of our galaxy. Using the location of the stellar locus in colour-magnitude space and the measured field-to-field variation in galaxy number counts we demonstrate that the absolute accuracy of our photometric calibration is at the 5% level or better. Finally, an investigation of the angular clustering of K- selected extended sources in our survey displays the expected scaling behaviour with limiting magnitude, with amplitudes in each magnitude bin in broad agreement with literature values.

We investigate the evolution of the galaxy luminosity function from the VIMOS-VLT Deep Survey (VVDS) from the present to z=2 in five (U, B, V, R and I) rest-frame band-passes. We use the first epoch VVDS deep sample of 11,034 spectra selected at 17.5 <= I_{AB} <= 24.0, on which we apply the Algorithm for Luminosity Function (ALF), described in this paper. We observe a substantial evolution with redshift of the global luminosity functions in all bands. From z=0.05 to z=2, we measure a brightening of the characteristic magnitude M* included in the magnitude range 1.8-2.5, 1.7-2.4, 1.2-1.9, 1.1-1.8 and 1.0-1.6 in the U, B, V, R and I rest-frame bands, respectively. We confirm this differential evolution of the luminosity function with rest-frame wavelength, from the measurement of the comoving density of bright galaxies (M < M*(z=0.1)). This density increases by a factor of around 2.6, 2.2, 1.8, 1.5, 1.5 between z=0.05 and z=1 in the U, B, V, R, I bands, respectively. We also measure a possible steepening of the faint-end slope of the luminosity functions, with \Delta\alpha ~ -0.3 between z=0.05 and z=1, similar in all bands.

We are presenting in this paper a detailed account of the methods used to compute the three-dimensional two-point galaxy correlation function in the VIMOS-VLT deep survey (VVDS). We investigate how instrumental selection effects and observational biases affect the measurements and identify the methods to correct them. We quantify the accuracy of our correction method using an ensemble of fifty mock galaxy surveys generated with the GalICS semi-analytic model of galaxy formation which incorporate the same selection biases and tiling strategy as the real data does. We demonstrate that we are able to recover the real-space two-point correlation function xi(s) to an accuracy better than 10% on scales larger than 1 h^{-1} Mpc, and of about 30% on scales below 1 h^{-1} Mpc, with the sampling strategy used for the first epoch VVDS data. The projected correlation function w_p(r_p) is recovered with an accuracy better than 10% on all scales 0.1 <= r <= 10 h^{-1} Mpc. There is a tendency for a small but systematic under-estimate of the correlation length derived from w_p(r_p) of 6% on average, remaining after our correction process. The large number of simulated surveys allows us to provide a reliable estimate of the cosmic error on the measurements of the correlation length r_0, of about 15-20% for the first epoch VVDS observation (Le Fevre et al.2004, astro-ph/0409133). The error estimation and measurement techniques outlined in this paper are being used in several studies which investigate in detail the clustering properties of galaxies in the VVDS data.

This paper presents the evolution of the clustering of the main population of galaxies from z=2.1 to z=0.2, from the first epoch VIMOS VLT Deep Survey (VVDS), a magnitude limited sample with 17.5<=I_{AB}<=24. We have computed the correlation functions \xi(r_p,\pi) and w_p(r_p), and the correlation length r_0(z), for the VVDS-02h and VVDS-CDFS fields, for a total of 7155 galaxies in a 0.61deg^2 area. We find that the correlation length in this sample stays roughly constant from z=0.5 to z=1.1, with r_0(z)=2.5-2.8 h^{-1} Mpc (comoving), for galaxies comparable in luminosity to the local 2dFGRS and SDSS samples, indicating that the amplitude of the correlation function was ~2.5x lower at z~1 than observed locally. The correlation length in our lowest redshift bin z=[0.2,0.5] is r_0=2.4 h^{-1} Mpc, lower than for any other population at the same redshift, indicating the low clustering of very low luminosity galaxies, 1.5 magnitudes fainter than in the 2dFGRS or SDSS. The correlation length is increasing to r_0~3.0 h^{-1} Mpc at higher redshifts z=[1.3,2.1], as we are observing increasingly brighter galaxies, comparable to galaxies with MB_AB=-20.5 locally. We compare our measurement to the DEEP2 measurements in the range z=[0.7,1.35] \citep{coil} on the population selected applying the same magnitude and color selection criteria as in their survey, and find comparable results. The slowly varying clustering of VVDS galaxies as redshift increases is markedly different from the predicted evolution of the clustering of dark matter, indicating that bright galaxies are already tracing the large scale structures emerging from the dark matter distribution 9-10 billion years ago, a supporting evidence for a strong evolution of the galaxy vs. dark matter bias.

This paper presents the ``First Epoch'' sample from the VIMOS VLT Deep Survey (VVDS). The VVDS goals, observations, data reduction with VIPGI, and redshift measurement with KBRED are discussed. Data have been obtained with the VIsible Multi Object Spectrograph (VIMOS) on the ESO-VLT UT3, allowing to observe ~600 slits simultaneously at R~230. A total of 11564 objects have been observed in the VVDS-02h and VVDS-CDFS Deep fields over a total area of 0.61deg^2, selected solely on the basis of apparent magnitude 17.5 <=I_{AB} <=24. The VVDS covers the redshift range 0 < z <= 5. It is successfully going through the ``redshift desert'' 1.5<z<2.2, while the range 2.2<z<2.7 remains of difficult access because of the VVDS wavelength coverage.A total of 9677 galaxies have a redshift measurement, 836 are stars, 90 are AGNs, and a redshift could not be measured for 961 objects. There are 1065 galaxies with a measured redshift z>1.4. The survey reaches a redshift measurement completeness of 78% overall (93% including less reliable objects), with a spatial sampling of the population of galaxies of 25% and ~30% in the VVDS-02h and VVDS-CDFS. The redshift accuracy measured from repeated observations with VIMOS and comparison to other surveys is ~276km/s. From this sample we present for the first time the redshift distribution of a magnitude limited spectroscopic sample down to IAB=24. The redshift distribution has a median of z=0.62, z=0.65, z=0.70, and z=0.76, for magnitude limited samples with IAB<=22.5, 23, 23.5, and 24. A high redshift tail above redshift 2 and up to redshift 5 becomes readily apparent for IAB>23.5, probing the bright star forming population of galaxies. This sample provides an unprecedented dataset to study galaxy evolution over 90% of the life of the universe

In this paper we present the optical and near-infrared identifications of the 1054 radio sources detected in the 20cm deep radio survey down to a 5sigma flux limit of about 80 microJy obtained with the VLA in the VIMOS VLT Deep Survey VVDS-02h deep field. Using U,B,V,R,I and K data, we identified 718 radio sources (~74% of the whole sample). The photometric redshift analysis shows that, in each magnitude bin, the radio sample has a higher median photometric redshift than the whole optical sample, while the median (V-I) color of the radio sources is redder than the median color of the whole optical sample. These results suggest that radio detection is preferentially selecting galaxies with higher intrinsic optical luminosity. From the analysis of the optical properties of the radio sources as function of the radio flux, we found that while about 35% of the radio sources are optically unidentified in the higher radio flux bin (S> 1.0 mJy), the percentage of unidentified sources decreases to about 25% in the faintest bins (S< 0.5 mJy). The median I magnitude for the total sample of radio sources,i.e. including also the unidentified ones, is brighter in the faintest radio bins than in the bin with higher radio flux. This suggests that most of the faintest radio sources are likely to be associated to relatively lower radio luminosity objects at relatively modest redshift, rather than radio-powerful, AGN type objects at high redshift.

The VIMOS VLT Deep Survey (VVDS), designed to measure 150,000 galaxy redshifts, requires a dedicated data reduction and analysis pipeline to process in a timely fashion the large amount of spectroscopic data being produced. This requirement has lead to the development of the VIMOS Interactive Pipeline and Graphical Interface (VIPGI), a new software package designed to simplify to a very high degree the task of reducing astronomical data obtained with VIMOS, the imaging spectrograph built by the VIRMOS Consortium for the European Southern Observatory, and mounted on Unit 3 (Melipal) of the Very Large Telescope (VLT) at Paranal Observatory (Chile). VIPGI provides the astronomer with specially designed VIMOS data reduction functions, a VIMOS-centric data organizer, and dedicated data browsing and plotting tools, that can be used to verify the quality and accuracy of the various stages of the data reduction process. The quality and accuracy of the data reduction pipeline are comparable to those obtained using well known IRAF tasks, but the speed of the data reduction process is significantly increased, thanks to the large set of dedicated features. In this paper we discuss the details of the MOS data reduction pipeline implemented in VIPGI, as applied to the reduction of some 20,000 VVDS spectra, assessing quantitatively the accuracy of the various reduction steps. We also provide a more general overview of VIPGI capabilities, a tool that can be used for the reduction of any kind of VIMOS data.

The VIMOS VLT Deep Survey (VVDS) is underway to study the evolution of galaxies, large scale structures and AGNs, from the measurement of more than 100000 spectra of faint objects. We present here the results from the first epoch observations of more than 20000 spectra. The main challenge of the program, the redshift measurements, is described, in particular entering the ``redshift desert'' in the range $1.5<z<3$ for which only very weak features are observed in the observed wavelength range. The redshift distribution of a magnitude limited sample brighter than $I_{AB}=24$ is presented for the first time, showing a peak at a low redshift $z\sim0.7$, and a tail extending all the way above $z=4$. The evolution of the luminosity function out to $z=1.5$ is presented, with the LF of blue star forming galaxies carrying most of the evolution, with L* changing by more than two magnitudes for this sub-sample.

We present a comparison between two optical cluster finding methods: a matched filter algorithm using galaxy angular coordinates and magnitudes, and a percolation algorithm using also redshift information. We test the algorithms on two mock catalogues. The first mock catalogue is built by adding clusters to a Poissonian background, while the other is derived from N-body simulations. Choosing the physically most sensible parameters for each method, we carry out a detailed comparison and investigate advantages and limits of each algorithm, showing the possible biases on final results. We show that, combining the two methods, we are able to detect a large part of the structures, thus pointing out the need to search for clusters in different ways in order to build complete and unbiased samples of clusters, to be used for statistical and cosmological studies. In addition, our results show the importance of testing cluster finding algorithms on different kinds of mock catalogues to have a complete assessment of their behaviour.

We have conducted a deep survey (r.m.s noise 17 microJy) with the Very Large Array (VLA) at 1.4 GHz, with a resolution of 6 arcsec, of a 1 square degree region included in the VIRMOS VLT Deep Survey. In the same field we already have multiband photometry down to I(AB)=25, and spectroscopic observations will be obtained during the VIRMOS VLT survey. The homogeneous sensitivity over the whole field has allowed to derive a complete sample of 1054 radio sources (5 sigma limit). We give a detailed description of the data reduction and of the analysis of the radio observations, with particular care to the effects of clean bias and bandwidth smearing, and of the methods used to obtain the catalogue of radio sources. To estimate the effect of the resolution bias on our observations we have modelled the effective angular-size distribution of the sources in our sample and we have used this distribution to simulate a sample of radio sources. Finally we present the radio count distribution down to 0.08 mJy derived from the catalogue. Our counts are in good agreement with the best fit derived from earlier surveys, and are about 50 % higher than the counts in the HDF. The radio count distribution clearly shows, with extremely good statistics, the change in the slope for the sub-mJy radio sources.

Research on the core of medium and high z clusters of galaxies can derive great benefits from integral field spectroscopy, and a key role in this respect will be played by the Integral Field Unit (IFU) being developed by the VIRMOS Consortium as part of the VIMOS spectrograph. After a brief technical outline of the instrument, some of the problems to be addressed by data reduction techniques are described, and the scientific issues to which VIMOS IFU is likely to give its major contribution are pointed out.