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.

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.

Extragalactic background light (EBL) anisotropy traces variations in the total production of photons over cosmic history, and may contain faint, extended components missed in galaxy point source surveys. Infrared EBL fluctuations have been attributed to primordial galaxies and black holes at the epoch of reionization (EOR), or alternately, intra-halo light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report new EBL anisotropy measurements from a specialized sounding rocket experiment at 1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely explained by IHL emission. The measured fluctuations are associated with an EBL intensity that is comparable to the background from known galaxies measured through number counts, and therefore a substantial contribution to the energy contained in photons in the cosmos.At near-infrared wavelengths, where the large zodiacal light foreground complicates absolute photometry measurements, the extragalactic background light (EBL) may be best accessed by anisotropy measurements. On large angular scales, fluctuations are produced by the clustering of galaxies, which is driven by the underlying distribution of dark matter. EBL anisotropy measurements can probe emission from epoch of reionization (EOR) galaxies (1-3) and directcollapse black holes (4) that formed during the EOR before the universe was fully ionized by exploiting the distinctive Lyman cutoff feature in the rest-frame ultraviolet (UV), thus probing the UV luminosity density at high redshifts (5). However, large-scale fluctuations may also arise from the intrahalo light (IHL) created by stars stripped from their parent galaxies during tidal interactions (6) at redshift z < 3. A multi-wavelength fluctuation analysis can distinguish among these scenarios and constrain the EOR star formation rate.A search for such background components must carefully account for fluctuations produced 2 by known galaxy populations. Linear galaxy clustering is an important contribution to fluctuations on scales much larger than galaxies themselves. On fine scales, the variation in the number of galaxies produces predominantly Poissonian fluctuations, with an amplitude that depends on the luminosity distribution. Anisotropy measurements suppress foreground galaxy fluctuations by masking known galaxies from an external catalog.The first detections of infrared fluctuations in excess of the contribution from known galaxies with the Spitzer Space Telescope (7-9) were interpreted as arising from a population of faint first-light galaxies at z > 7. The Hubble Space Telescope was used at shorter wavelengths (10) to carry out a fluctuation study in a small deep field but did not report fluctuations in excess of known galaxy populations. Measurements with the AKARIsatellite (11) show excess fluctuations with a blue spectrum rapidly rising from 4.1μm to 2.4μm. Fluctuation measurements in a large survey fi...

Unresolved anisotropies of the cosmic near-infrared background radiation are expected to have contributions from the earliest galaxies during the epoch of reionization and from faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to pinpoint conclusively the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from subarcminute to one-degree angular scales, and find the clustering amplitude to be larger than predicted by the models based on the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. However, a physical mechanism that increases the clustering amplitude is needed. Motivated by recent results related to the extended stellar light profile in dark-matter haloes, we consider the possibility that the fluctuations originate from intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 per cent relative to the total luminosity in dark-matter haloes of 10(9) to 10(12) solar masses at redshifts of about 1 to 4.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer (RT) self‐consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure RT tests (i.e. for fixed, non‐evolving density fields). This is the second paper of the Cosmological Radiative Transfer Comparison Project, in which we compare nine independent RT codes directly coupled to gas dynamics on three relatively simple astrophysical hydrodynamics problems: (i) the expansion of an H ii region in a uniform medium, (ii) an ionization front in a 1/r2 density profile with a flat core and (iii) the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The absolute intensity and the anisotropies of the near-IR background provide an observational probe of the first-light galaxies and their spatial distribution. We estimate the extragalactic background light intensity during reionization by accounting for the stellar and nebular emission from first-light galaxies. We require the UV photon density from these galaxies to generate a reionization history that is consistent with the optical depth to electron scattering from cosmic microwave background measurements. We also require the bright-end luminosity function of galaxies in our models to reproduce the measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 0.1 to 0.4 nW m −2 sr −1 at the peak of its spectrum at ∼ 1.1 µm. We also discuss the anisotropy power spectrum of the near-IR background using a halo model to describe the galaxy distribution. We compare our predictions for the anisotropy power spectrum to existing measurements from deep near-IR imaging data from Spitzer/IRAC, Hubble/NICMOS, and AKARI. The predicted rms fluctuations at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing measurements. While strong arguments have been made that the measured fluctuations do not have an origin involving faint low-redshift galaxies, we find that measurements in the literature are also incompatible with galaxies present during the era of reionization. The measured near-IR background anisotropies remain unexplained with an unknown origin. Subject headings: cosmology: theory -diffuse radiation -intergalactic medium -large scale structure of universe * −1.35

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multiwavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcmin-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 μm. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at redshifts greater than 8 to be . This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point-source detection level in current surveys.

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (fNL) and fourth-order (gNL and τNL) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of τNL and gNL with WMAP 5-year data, our study provides the first forecasts for future constraints on gNL. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio ANL = τNL/(6fNL/5) 2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to lmax = 600 we constrain the cubic order non-Gaussianity parameters τNL, and gNL and find −7.4 < gNL/10 5 < 8.2 and −0.6 < τNL/10 4 < 3.3 improving the previous COBE-based limit on τNL < 10 8 nearly four orders of magnitude with WMAP.