The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8-1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z > 1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10 9 M to z ≈ 2, reaching the knee of the ultraviolet luminosity function of galaxies to z ≈ 8. The survey covers approximately 800 arcmin 2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H = 27.7 mag) covers ∼125 arcmin 2 within Great Observatories Origins Deep Survey (GOODS)-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (Extended Groth Strip, COSMOS, and Ultra-deep Survey) and covers the full area to a 5σ pointsource limit of H 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding-cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

We present a physically motivated model for the early coevolution of massive spheroidal galaxies and active nuclei at their centers. Within dark matter halos, forming at the rate predicted by the canonical hierarchical clustering scenario, the gas evolution is controlled by gravity, radiative cooling, and heating by feedback from supernovae and from the growing active nucleus. Supernova heating is increasingly effective with decreasing binding energy in slowing down the star formation and in driving gas outflows. The more massive protogalaxies virializing at earlier times are thus the sites of the faster star formation. The correspondingly higher radiation drag fastens the angular momentum loss by the gas, resulting in a larger accretion rate onto the central black hole. In turn, the kinetic energy carried by outflows driven by active nuclei can unbind the residual gas, thus halting both the star formation and the black hole growth, in a time again shorter for larger halos. For the most massive galaxies the gas unbinding time is short enough for the bulk of the star formation to be completed before Type Ia supernovae can substantially increase the Fe abundance of the interstellar medium, thus accounting for the -enhancement seen in the largest galaxies. The feedback from supernovae and from the active nucleus also determines the relationship between the black hole mass and the mass, or the velocity dispersion, of the host galaxy, as well as the black hole mass function. In both cases the model predictions are in excellent agreement with the observational data. Coupling the model with GRASIL (Silva et al. 1998), the code computing in a selfconsistent way the chemical and spectrophotometric evolution of galaxies over a very wide wavelength interval, we have obtained predictions in excellent agreement with observations for a number of observables that proved to be extremely challenging for all the current semianalytic models, including the submillimeter counts and the corresponding redshift distributions, and the epoch-dependent K-band luminosity function of spheroidal galaxies.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at z ∼ 1.5 − 8, and to study Type Ia SNe beyond z > 1.5. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers ∼ 125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of ∼ 800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

This paper presents new deep and wide narrow-band surveys undertaken with UKIRT, Subaru and the VLT; a unique combined effort to select large, robust samples of Hα star-forming galaxies at z = 0.40, 0.84, 1.47 and 2.23 (corresponding to look-back times of 4.2, 7.0, 9.2 and 10.6 Gyrs) in a uniform manner over ∼ 2 deg 2 in the COSMOS and UDS fields. The deep multi-epoch Hα surveys reach a matched 3 σ flux limit of ≈ 3 M yr −1 out to z = 2.2 for the first time, while the wide area and the coverage over two independent fields allow to greatly overcome cosmic variance and assemble by far the largest samples of Hα emitters. Catalogues are presented for a total of 1742, 637, 515 and 807 Hα emitters, robustly selected at z = 0.40, 0.84, 1.47 and 2.23, respectively, and used to determine the Hα luminosity function and its evolution. The faint-end slope of the Hα luminosity function is found to be α = −1.60 ± 0.08 over z = 0 − 2.23, showing no significant evolution. The characteristic luminosity of SF galaxies, L * Hα , evolves significantly as log L * Hα (z) = 0.45z + log L * z=0 . This is the first time Hα has been used to trace SF activity with a single homogeneous survey at z = 0.4−2.23. Overall, the evolution seen with Hα is in good agreement with the evolution seen using inhomogeneous compilations of other tracers of star formation, such as FIR and UV, jointly pointing towards the bulk of the evolution in the last 11 Gyrs being driven by a statistically similar SF population across cosmic time, but with a strong luminosity increase from z ∼ 0 to z ∼ 2.2. Our uniform analysis allows to derive the Hα star formation history of the Universe (SFRH), showing a clear rise up to z ∼ 2.2, for which the simple parametrisation log 10 ρ SFR = −2.1(1 + z) −1 is valid over 80 per cent of the age of the Universe. The results reveal that both the shape and normalisation of the Hα SFRH are consistent with the measurements of the stellar mass density growth, confirming that our Hα SFRH is tracing the bulk of the formation of stars in the Universe for z < 2.23. The star formation activity over the last ∼ 11 Gyrs is responsible for producing ∼ 95 per cent of the total stellar mass density observed locally, with half of that being assembled in 2 Gyrs between z = 1.2-2.2, and the other half in 8 Gyrs (since z < 1.2). If the star-formation rate density continues to decline with time in the same way as seen in the past ∼ 11 Gyrs, then the stellar mass density of the Universe will reach a maximum which is only 5 per cent higher than the present-day value.