The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015–2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14 000 deg² of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

We exploit the deep and extended far-IR data sets (at 70, 100 and 160 μm) of the Herschel Guaranteed Time Observation (GTO) PACS Evolutionary Probe (PEP) Survey, in combination with the Herschel Multi-tiered Extragalactic Survey data at 250, 350 and 500 μm, to derive the evolution of the rest-frame 35-, 60-, 90- and total infrared (IR) luminosity functions (LFs) up to z ∼ 4. We detect very strong luminosity evolution for the total IR LF (L_IR ∝ (1 + z)^{3.55 ± 0.10} up to z ∼ 2, and ∝ (1 + z)^{1.62 ± 0.51} at 2 < z ≲ 4) combined with a density evolution (∝(1 + z)^{−0.57 ± 0.22} up to z ∼ 1 and ∝ (1 + z)^{−3.92 ± 0.34} at 1 < z ≲ 4). In agreement with previous findings, the IR luminosity density (ρ_IR) increases steeply to z ∼ 1, then flattens between z ∼ 1 and z ∼ 3 to decrease at z ≳ 3. Galaxies with different spectral energy distributions, masses and specific star formation rates (SFRs) evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to ρ_IR. Galaxies occupying the well-established SFR–stellar mass main sequence (MS) are found to dominate both the total IR LF and ρ_IR at all redshifts, with the contribution from off-MS sources (≥0.6 dex above MS) being nearly constant (∼20 per cent of the total ρ_IR) and showing no significant signs of increase with increasing z over the whole 0.8 < z < 2.2 range. Sources with mass in the range 10 ≤ log(M/M_⊙) ≤ 11 are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z (≳2) LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the Super Massive Black Holes (SMBH) grows and is obscured by dust (possibly triggered by a major merging event), is followed by an AGN-dominated phase, then evolving towards a local elliptical. On the other hand, moderately star-forming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role within the entire set of Herschel surveys, and the field selection that includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, and EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high-redshift galaxy populations, thus testing and superseding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

We present results from the deepest Herschel-Photodetector Array Camera and Spectrometer (PACS) far-infrared blank field extragalactic survey, obtained by combining observations of the Great Observatories Origins Deep Survey (GOODS) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and theconstruction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3σ depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 μm, respectively, and resolve ~75% of the cosmic infrared background at 100 μm and 160 μm into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths, and to determine the infrared luminosity function of galaxies down to L IR = 1011 L ⊙ at z ~ 1 and L IR = 1012 L ⊙ at z ~ 2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3σ) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

Aims. We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE). These sources, detected in COSMOS and ECDFS, have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z ≃ 6. Methods. We looked for counterparts of the ALMA sources in all the available multi-wavelength (from HST to VLA) and photometric redshift catalogues. We also made use of deeper UltraVISTA and Spitzer source lists and maps to identify optically dark sources with no matches in the public catalogues. We used the sources with estimated redshifts to derive the 250 μm rest-frame and total infrared (8–1000 μm) LFs from z ≃ 0.5 to 6. Results. Our ALMA blind survey (860 μm flux density range: ∼0.3–12.5 mJy) allows us to further push the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint end of the infrared LF, showing little evolution from z ≃ 2.5 to z ≃ 6, and a “flat” slope up to the highest redshifts (i.e. 4.5 < z < 6). The SFRD obtained by integrating the luminosity function remains almost constant between z ≃ 2 and z ≃ 6, and significantly higher than the optical or ultra-violet derivations, showing a significant contribution of dusty galaxies and obscured star formation at high-z. About 14% of all the ALPINE serendipitous continuum sources are found to be optically and near-infrared (near-IR) dark (to a depth Ks ∼ 24.9 mag). Six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [C II] emitters at z ≃ 5. The six HST+near-IR dark galaxies with mid-IR counterparts are found to contribute about 17% of the total SFRD at z ≃ 5 and to dominate the high-mass end of the stellar mass function at z > 3.