The Lya line in the UV and the [CII] line in the FIR are widely used tools to identify galaxies and to obtain insights into ISM properties in the early Universe. By combining data obtained with ALMA in band 7 at ~ 320 GHz as part of the ALMA Large Program to INvestigate [CII] at Early Times (ALPINE) with spectroscopic data from DEIMOS at Keck, VIMOS and FORS2 at the VLT, we assembled a unique sample of 53 main-sequence star-forming galaxies at 4.4 < z < 6 in which we detect both the Lya line and the [CII]. We used [CII], observed with ALMA, as a tracer of the systemic velocity of the galaxies, and we find that 90% of the selected objects have Lya-[CII] velocity offsets in the range 0 < Dv_Lya-[CII] < 400 km/s, in line with the few measurements available so far in the early Universe, and significantly smaller than those observed at lower z. We observe ISM-[CII] offsets in the range -500 < Dv_ISM-[CII] < 0 km/s, in line with values at all redshifts. We find significant anticorrelations between Dv_Lya-[CII] and the Lya rest-frame equivalent width EW0(Lya) (or equivalently, the Lya escape fraction f_esc(Lya)). According to available models for the radiative transfer of Lya photons, the escape of Lya photons would be favored in galaxies with high outflow velocities, in agreement with our observations. The uniform shell model would also predict that the Lya escape in galaxies with slow outflows (0 < v_out < 300 km/s) is mainly determined by the neutral hydrogen column density (NHI), while the alternative model by Steidel+10 would favor a combination of NHI and covering fraction as driver of the Lya escape. We suggest that the observed increase in Lya escape that is observed between z~2 and z~6 is not due to a higher incidence of fast outflows at high redshift, but rather to a decrease in average NHI along the line of sight, or alternatively, a decrease in HI covering fraction. [abridged]

A significant fraction of high redshift star-forming disc galaxies are known to host giant clumps, whose formation, nature and role in galaxy evolution are yet to be understood. In this work we first present a new automated method based on deep neural networks to detect clumps in galaxy images and show that it is more sensitive and faster than previous proposed methods. We then use this method to systematically detect clumps in the rest-frame optical and UV images of a complete sample of $\sim1500$ star forming galaxies at $1<z<3$ in the CANDELS survey as well as in images from the VELA forward modeled zoom-in cosmological simulations. We show that observational effects have a dramatic impact on the derived clump properties leading to an overestimation of the clump mass up to a factor of 10, which highlights the importance of fair and direct comparisons between observations and simulations. After correcting for these effects with a simple mixture density network, we estimate that the clump stellar mass function follows a power-law down to the completeness limit ($10^{7}$ solar masses) with the majority of the clumps being less massive than $10^9$ solar masses. We find that between $20-40\%$ of $\sim10^{10}$ solar mass star-forming galaxies contain at least one massive ($M_c>10^{7}M_\odot$) off-centered clump but only $\sim2-5\%$ of the total galaxy stellar mass is in those clumps. We also show indications that the contribution of clumps to the stellar mass is more important in extended and low mass galaxies. The simulations explored in this work overall reproduce the shape of the observed clump stellar mass function when confronted under the same conditions although tend to lie in the lower limit of the confidence intervals of the observations. This agreement suggests that most of the observed clumps are formed in-situ through violent disk instabilities.

Using analytic modeling and simulations, we address the origin of an abundance of star-forming, clumpy, extended gas rings about massive central bodies in massive galaxies at $z \!<\! 4$. Rings form by high-angular-momentum streams and survive in galaxies of $M_{\rm star} \!>\! 10^{9.5-10} M_\odot$ where merger-driven spin flips and supernova feedback are ineffective. The rings survive after events of compaction to central nuggets. Ring longevity was unexpected based on inward mass transport driven by torques from violent disc instability. However, evaluating the torques from a tightly wound spiral structure, we find that the timescale for transport per orbital time is long and $\propto\! \delta_{\rm d}^{-3}$, with $\delta_{\rm d}$ the cold-to-total mass ratio interior to the ring. A long-lived ring forms when the ring transport is slower than its replenishment by accretion and the interior depletion by SFR, both valid for $\delta_{\rm d} \!<\! 0.3$. The central mass that lowers $\delta_{\rm d}$ is a compaction-driven bulge and/or dark matter, aided by the lower gas fraction at $z \!<\! 4$, provided that it is not too low. The ring is Toomre unstable for clump and star formation. The high-$z$ dynamic rings are not likely to arise form secular resonances or collisions. AGN feedback is not expected to affect the rings. Mock images of simulated rings through dust indicate qualitative consistency with observed rings about bulges in massive $z\!\sim\!0.5\!-\!3$ galaxies, in $H_{\alpha}$ and deep HST imaging. ALMA mock images indicate that $z\!\sim\!0.5\!-\!1$ rings should be detectable. We quote expected observable properties of rings and their central nuggets.

Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

Using the CANDELS photometric catalogs for the HST/ACS and WFC3, we identified massive evolved galaxies at $3 < z < 4.5$, employing three different selection methods. We find the comoving number density of these objects to be $\sim 2 \times 10^{-5}$ and $8 \times 10^{-6}Mpc^{-3}$ after correction for completeness for two redshift bins centered at $z=3.4, 4.7$. We quantify a measure of how much confidence we should have for each candidate galaxy from different selections and what are the conservative error estimates propagated into our selection. Then we compare the evolution of the corresponding number densities and their stellar mass density with numerical simulations, semi-analytical models, and previous observational estimates, which shows slight tension at higher redshifts as the models tend to underestimate the number and mass densities. By estimating the average halo masses of the candidates ($M_h \approx 4.2, 1.9, 1.3 \times 10^{12} M_\odot$ for redshift bins centered at $z=3.4, 4.1, 4.7$), we find them to be consistent with halos that were efficient in turning baryons to stars and were relatively immune to the feedback effects and on the verge of transition into hot-mode accretion. This can suggest the relative cosmological starvation of the cold gas followed by an overconsumption phase in which the galaxy consumes the available cold gas rapidly as one of the possible drivers for the quenching of the massive evolved population at high redshift.

The origin of the correlations between mass, morphology, quenched fraction, and formation history in galaxies is difficult to define, primarily due to the uncertainties in galaxy star-formation histories. Star-formation histories are better constrained for higher redshift galaxies, observed closer to their formation and quenching epochs. Here we use "non-parametric" star-formation histories and a nested sampling method to derive constraints on the formation and quenching timescales of quiescent galaxies at $0.7<z<2.5$. We model deep HST grism spectroscopy and photometry from the CLEAR (CANDELS Lyman$-\alpha$ Emission at Reionization) survey. The galaxy formation redshifts, $z_{50}$ (defined as the point where they had formed 50\% of their stellar mass) range from $z_{50}\sim 2$ (shortly prior to the observed epoch) up to $z_{50} \simeq 5-8$. \editone{We find that early formation redshifts are correlated with high stellar-mass surface densities, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) >$10.25, where $\Sigma_1$ is the stellar mass within 1~pkpc (proper kpc). Quiescent galaxies with the highest stellar-mass surface density, $\log\Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) > 10.25$, } show a \textit{minimum} formation redshift: all such objects in our sample have $z_{50} > 2.9$. Quiescent galaxies with lower surface density, $\log \Sigma_1 / (M_\odot\ \mathrm{kpc}^{-2}) = 9.5 - 10.25$, show a range of formation epochs ($z_{50} \simeq 1.5 - 8$), implying these galaxies experienced a range of formation and assembly histories. We argue that the surface density threshold $\log\Sigma_1/(M_\odot\ \mathrm{kpc}^{-2})>10.25$ uniquely identifies galaxies that formed in the first few Gyr after the Big Bang, and we discuss the implications this has for galaxy formation models.

The [CII] 158 micron line is one of the strongest IR emission lines, which has been shown to trace the star-formation rate (SFR) of galaxies in the nearby Universe and up to $z \sim 2$. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at $4.4 < z< 5.9$, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature we examine the relation between the [CII] luminosity and the SFR over the entire redshift range from $z \sim 4-8$. ALPINE galaxies, which are both detected in [CII] and dust continuum, show a good agreement with the local L([CII])-SFR relation. Galaxies undetected in the continuum with ALMA are found to be over-luminous in [CII], when the UV SFR is used. After accounting for dust-obscured star formation, by an amount SFR(IR)$\approx$SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L([CII])-SFR relation comparable to the local one. When [CII] non-detections are taken into account, the slope may be marginally steeper at high-z, although this is still somewhat uncertain. When compared in a homogeneous manner, the $z>6 $ [CII] measurements (detections and upper limits) do not behave very differently from the $z \sim 4-6$ data. We find a weak dependence of L([CII])/SFR on the Lyman-alpha equivalent width. Finally, we find that the ratio L([CII])/LIR $\sim (1-3) \times 10^{-3}$ for the ALPINE sources, comparable to that of "normal" galaxies at lower redshift. Our analysis, which includes the largest sample ($\sim 150$ galaxies) of [CII] measurements at $z>4$ available so far, suggests no or little evolution of the L([CII])-SFR relation over the last 13 Gyr of cosmic time.

In this paper, we extend the source detection in the GOODS-ALMA field (69 arcmin$^2$, rms sensitivity $\sigma$ $\simeq$ 0.18 mJy.beam$^{-1}$), to deeper levels than presented in Franco et al. (2018). Using positional information at 3.6 and 4.5 $\mu$m (from Spitzer-IRAC), we explore the presence of galaxies detected at 1.1 mm with ALMA below our original blind detection limit of 4.8-$\sigma$ at which the number of spurious sources starts to dominate over that of real sources. In this Supplementary Catalog, we find a total of 16 galaxies, including 2 galaxies with no counterpart in HST images (also known as optically-dark galaxies) down to a 5$\sigma$ limiting depth of H = 28.2 AB (HST/WFC3 F160W). This brings the total sample of GOODS-ALMA 1.1 mm sources to 35 galaxies. Galaxies in the new sample cover a wider dynamic range in redshift ($z$ = 0.65 - 4.73), are on average twice as large (1.3 vs 0.65 kpc) and and have lower stellar mass (M$_{\star}^{\rm SC}$ = 7.6$\times$10$^{10}$M$_\odot$ vs M$_{\star}^{\rm MC}$ = 1.2$\times$10$^{11}$M$_\odot$). Although exhibiting larger physical sizes, these galaxies have still far-infrared sizes significantly more compact than inferred from their optical emission. We show that the astrometry of the HST image does not only suffer from a global astrometric shift, as already discussed in previous papers, but also from local shifts. These distortions were artificially introduced in the process of building the mosaic of the GOODS-South HST image. By comparing the positions of almost 400 galaxies detected by HST, Pan-STARRS and ALMA, we create a distortion map which can be used to correct for these astrometric issues.

We investigate the properties of a sample of 35 galaxies, detected with ALMA at 1.1 mm in the GOODS-ALMA field (area of 69 arcmin$^2$, resolution = 0.60", RMS $\simeq$ 0.18 mJy beam$^{-1}$). Using the UV-to-radio deep multiwavelength coverage of the GOODS-South field, we fit the spectral energy distributions of these galaxies to derive their key physical properties. The galaxies detected by ALMA are among the most massive at $z$ = 2-4 (M$_{\star,med}$ = 8.5$ \times$ 10$^{10}$ M$_\odot$) and are either starburst or located in the upper part of the galaxy star-forming main sequence. A significant portion of our galaxy population ($\sim$ 40%), located at $z\sim$ 2.5-3, exhibits abnormally low gas fractions. The sizes of these galaxies, measured with ALMA, are compatible with the trend between $H$-band size and stellar mass observed for $z\sim2$ elliptical galaxies suggesting that they are building compact bulges. We show that there is a strong link between star formation surface density (at 1.1 mm) and gas depletion time: the more compact a galaxy's star-forming region is, the shorter its lifetime will be (without gas replenishment). The identified compact sources associated with relatively short depletion timescales ($\sim$100 Myr), are the ideal candidates to be the progenitors of compact elliptical galaxies at $z$ $\sim$ 2.

We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_{\lambda}$(Ly$\alpha$) across the range $-58 \unicode{xC5} \lesssim$ $W_{\lambda}$(Ly$\alpha$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=-18\unicode{xC5}$) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_{\lambda}$(Ly$\alpha$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

We present a robust method, weighted von Mises kernel density estimation, along with boundary correction to reconstruct the underlying number density field of galaxies. We apply this method to galaxies brighter than $\rm HST/F160w\le 26$ AB mag at the redshift range of $0.4\leq z \leq 5$ in the five CANDELS fields (GOODS-N, GOODS-S, EGS, UDS, and COSMOS). We then use these measurements to explore the environmental dependence of the star formation activity of galaxies. We find strong evidence of environmental quenching for massive galaxies ($\rm M \gtrsim 10^{11} \rm {M}_\odot$) out to $z\sim 3.5$ such that an over-dense environment hosts $\gtrsim 20\%$ more massive quiescent galaxies compared to an under-dense region. We also find that environmental quenching efficiency grows with stellar mass and reaches $\sim 60\%$ for massive galaxies at $z\sim 0.5$. The environmental quenching is also more efficient in comparison to the stellar mass quenching for low mass galaxies ($\rm M \lesssim 10^{10} \rm {M}_\odot$) at low and intermediate redshifts ($z\lesssim 1.2$). Our findings concur thoroughly with the "over-consumption" quenching model where the termination of cool gas accretion (cosmological starvation) happens in an over-dense environment and the galaxy starts to consume its remaining gas reservoir in depletion time. The depletion time depends on the stellar mass and could explain the evolution of environmental quenching efficiency with the stellar mass.

We report on the HST detection of the Lyman-continuum (LyC) radiation emitted by a galaxy at redshift z=3.794, dubbed Ion1 (Vanzella et al. 2012). The LyC from Ion1 is detected at rest-frame wavelength 820$\sim$890 \AA with HST WFC3/UVIS in the F410M band ($m_{410}=27.60\pm0.36$ magnitude (AB), peak SNR = 4.17 in a circular aperture with radius r = 0.12'') and at 700$\sim$830 \AA with the VLT/VIMOS in the U-band ($m_U = 27.84\pm0.19$ magnitude (AB), peak SNR = 6.7 with a r = 0.6'' aperture). A 20-hr VLT/VIMOS spectrum shows low- and high-ionization interstellar metal absorption lines, the P-Cygni profile of CIV and Ly$\alpha$ in absorption. The latter spectral feature differs from what observed in known LyC emitters, which show strong Ly$\alpha$ emission. An HST far-UV color map reveals that the LyC emission escapes from a region of the galaxy that is bluer than the rest, presumably because of lower dust obscuration. The F410M image shows that the centroid of the LyC emission is offset from the centroid of the non-ionizing UV emission by 0.12''$\pm$0.03'', corresponding to 0.85$\pm$0.21 kpc (physical), and that its morphology is likely moderately resolved. These morphological characteristics favor a scenario where the LyC photons produced by massive stars escape from low HI column-density "cavities" in the ISM, possibly carved by stellar winds and/or supernova. We also collect the VIMOS U-band images of a sample of 107 Lyman-break galaxies with spectroscopic redshifts at $3.40<z<3.95$, i.e. sampling the LyC, and stack them with inverse-variance weights. No LyC emission is detected in the stacked image, resulting in a 32.5 magnitude (AB) flux limit (1$\sigma$) and an upper limit of absolute LyC escape fraction $f_{esc}^{abs} < 0.63\%$. LyC emitters like Ion1 are very likely at the bright-end of the LyC luminosity function.

We study the efficiency of galactic feedback in the early Universe by stacking the [C II] 158 um emission in a large sample of normal star-forming galaxies at 4 < z < 6 from the ALMA Large Program to INvestigate [C II] at Early times (ALPINE) survey. Searching for typical signatures of outflows in the high-velocity tails of the stacked [C II] profile, we observe (i) deviations from a single-component Gaussian model in the combined residuals and (ii) broad emission in the stacked [C II] spectrum, with velocities of |v|<~ 500 km/s. The significance of these features increases when stacking the subset of galaxies with star formation rates (SFRs) higher than the median (SFRmed = 25 Msun/yr), thus confirming their star-formation-driven nature. The estimated mass outflow rates are comparable to the SFRs, yielding mass-loading factors of the order of unity (similarly to local star-forming galaxies), suggesting that star-formation-driven feedback may play a lesser role in quenching galaxies at z > 4. From the stacking analysis of the datacubes, we find that the combined [C II] core emission (|v|< 200 km/s) of the higher-SFR galaxies is extended on physical sizes of ~ 30 kpc (diameter scale), well beyond the analogous [C II] core emission of lower-SFR galaxies and the stacked far-infrared continuum. The detection of such extended metal-enriched gas, likely tracing circumgalactic gas enriched by past outflows, corroborates previous similar studies, confirming that baryon cycle and gas exchanges with the circumgalactic medium are at work in normal star-forming galaxies already at early epochs.

Supermassive black holes (SMBHs) are tightly correlated with their hosts but the origin of such connection remains elusive. To explore the cosmic build-up of this scaling relation, we present an empirically-motivated model that tracks galaxy and SMBH growth down to z=0. Starting from a random mass seed distribution at z=10, we assume that each galaxy evolves on the star-forming "main sequence" (MS) and each BH follows the recently-derived stellar mass (M$_{\star}$) dependent ratio between BH accretion rate and star formation rate, going as BHAR/SFR$\propto$M$_{\star}^{0.73[+0.22,-0.29]}$. Our simple recipe naturally describes the BH-galaxy build-up in two stages. At first, the SMBH lags behind the host that evolves along the MS. Later, as the galaxy grows in M$_{\star}$, our M$_{\star}$-dependent BHAR/SFR induces a super-linear BH growth, as M$_{\rm BH}$$\propto$M$_{\star}^{1.7}$. According to this formalism, smaller BH seeds increase their relative mass faster and earlier than bigger BH seeds, at fixed M$_{\star}$, thus setting along a gradually tighter M$_{\rm BH}$-M$_{\star}$ locus towards higher M$_{\star}$. Assuming reasonable values of the radiative efficiency $\epsilon \sim$0.1, our empirical trend agrees with both high-redshift model predictions and intrinsic M$_{\rm BH}$-M$_{\star}$ relations of local BHs. We speculate that the observed non-linear BH-galaxy build-up is reflected in a twofold behavior with dark matter halo mass (M$_{\rm DM}$), displaying a clear turnover at M$_{\rm DM}\sim$2$\times$10$^{12}$M$_{\odot}$. While Supernovae-driven feedback suppresses BH growth in smaller halos (BHAR/SFR$\propto$M$_{\rm DM}^{1.6}$), above the M$_{\rm DM}$ threshold cold gas inflows possibly fuel both BH accretion and star formation in a similar fashion (BHAR/SFR$\propto$M$_{\rm DM}^{0.3}$).

We present a VLT/X-Shooter spectroscopy of the Lyman continuum (LyC) emitting galaxy 'Ion2' at z=3.2121 and compare it to that of the recently discovered strongly lensed LyC-emitter at z=2.37, known as the 'Sunburst' arc. Three main results emerge from the X-Shooter spectrum: (a) the Lya has three distinct peaks with the central one at the systemic redshift, indicating a ionised tunnel through which both Lya and LyC radiation escape; (b) the large O32 oxygen index ([OIII]4959-5007 / [OII]3727-3729) of 9.18(-1.32/+1.82) is compatible to those measured in local (z~0.4) LyC leakers; (c) there are narrow nebular high-ionisation metal lines with \sigma_v < 20 km/s, which confirms the presence of young hot, massive stars. The HeII1640 appears broad, consistent with a young stellar component including Wolf-Rayet stars. Similarly, the Sunburst LyC-emitter shows a triple-peaked Lya profile and from VLT/MUSE spectroscopy the presence of spectral features arising from young hot and massive stars. The strong lensing magnification, (\mu > 20), suggests that this exceptional object is a gravitationally-bound star cluster observed at a cosmological distance, with a stellar mass M <~ 10^7 Msun and an effective radius smaller than 20 pc. Intriguingly, sources like Sunburst but without lensing magnification might appear as Ion2-like galaxies, in which unresolved massive star clusters dominate the ultraviolet emission. This work supports the idea that dense young star clusters can contribute to the ionisation of the IGM through holes created by stellar feedback.

We report the detection of [CII]158um emission from a system of three closely-separated sources in the COSMOS field at z~4.56, as part of the ALMA Large Program to INvestigate CII at Early times (ALPINE). The two dominant sources are closely associated, both spatially (1.6"~11kpc) and in velocity (~100km/s), while the third source is slightly more distant (2.8"~18kpc, ~300km/s). The second strongest source features a slight velocity gradient, while no significant velocity gradient is seen in the other two sources. Using the observed [CII] luminosities, we derive a total log(SFR_[CII]/[Msol/year])=2.8+/-0.2, which may be split into contributions of 59%, 31%, and 10% from the central, east, and west sources, respectively. Comparison of these [CII] detections to recent zoom-in cosmological simulations suggests an ongoing major merger. We are thus witnessing a system in a major phase of mass build-up by merging, including an on-going major merger and an upcoming minor merger, which is expected to end up in a single massive galaxy by z~2.5.

We present a WFC3 F160W ($H$-band) selected catalog in the CANDELS/GOODS-N field containing photometry from the ultraviolet (UV) to the far-infrared (IR), photometric redshifts and stellar parameters derived from the analysis of the multi-wavelength data. The catalog contains 35,445 sources over the 171 arcmin$^{2}$ of the CANDELS F160W mosaic. The 5$\sigma$ detection limits (within an aperture of radius 0\farcs17) of the mosaic range between $H=27.8$, 28.2 and 28.7 in the wide, intermediate and deep regions, that span approximately 50\%, 15\% and 35\% of the total area. The multi-wavelength photometry includes broad-band data from UV (U band from KPNO and LBC), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), near-to-mid IR (HST/WFC3 F105W, F125W, F140W and F160W, Subaru/MOIRCS Ks, CFHT/Megacam K, and \spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) and far IR (\spitzer/MIPS 24$\mu$m, HERSCHEL/PACS 100 and 160$\mu$m, SPIRE 250, 350 and 500$\mu$m) observations. In addition, the catalog also includes optical medium-band data (R$\sim50$) in 25 consecutive bands, $\lambda=500$ to 950~nm, from the SHARDS survey and WFC3 IR spectroscopic observations with the G102 and G141 grisms (R$\sim210$ and 130). The use of higher spectral resolution data to estimate photometric redshifts provides very high, and nearly uniform, precision from $z=0-2.5$. The comparison to 1,485 good quality spectroscopic redshifts up to $z\sim3$ yields $\Delta z$/(1+$z_{\rm spec}$)$=$0.0032 and an outlier fraction of $\eta=$4.3\%. In addition to the multi-band photometry, we release added-value catalogs with emission line fluxes, stellar masses, dust attenuations, UV- and IR- based star formation rates and rest-frame colors.

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work.

The Hubble Space Telescope has produced astonishing science over the past thirty years. Hubble's productivity can continue to soar for years to come provided some worn out components get upgraded. While powerful new ground-based and space telescopes are expected to come online over the next decade, none of them will have the UV capabilities that make Hubble a unique observatory. Without Hubble, progress in UV and blue optical astrophysics will be halted. Observations at these wavelengths are key for a range of unresolved astrophysics questions, ranging from the characterization of solar system planets to understanding interaction of galaxies with the intergalactic medium and the formation history of the universe. Hubble will remain our only source of high-angular resolution UV imaging and high-sensitivity UV spectroscopy for the next two decades, offering the ability for continued unique science and maximizing the science return from complementary observatories. Therefore, we recommend that NASA, ESA, and the private sector study the scientific merit, technical feasibility, and risk of a new servicing mission to Hubble to boost its orbit, fix aging components, and expand its instrumentation. Doing so would: 1) keep Hubble on its path to reach its unmet full potential, 2) extend the mission's lifetime past the next decade, which will maximize the synergy of Hubble with other upcoming facilities, and 3) enable and enhance the continuation of scientific discoveries in UV and optical astrophysics.

We present the first results from the deep and wide 5 GHz radio observations of the Great Observatories Origins Deep Survey (GOODS)-North ($\sigma=3.5 \; \mu Jy \; beam^{-1}$, synthesized beam size $\theta =$ 1.47 arcsec $\times$ 1.42 arcsec, and 52 sources over 109 arcmin$^{2}$) and GOODS-South ($\sigma=3.0 \; \mu Jy \; beam^{-1}$, $\theta=$0.98 arcsec $\times$ 0.45 arcsec, and 88 sources over 190 arcmin$^{2}$) fields using the Karl G. Jansky Very Large Array. We derive radio spectral indices {\alpha} between 1.4 and 5 GHz using the beam-matched images and show that the overall spectral index distribution is broad even when the measured noise and flux bias are considered. We also find a clustering of faint radio sources around $\alpha=0.8$, but only within $S_{5GHz} < 150 \; \mu Jy$. We demonstrate that the correct radio spectral index is important for deriving accurate rest frame radio power and analyzing the radio-FIR correlation, and adopting a single value of $\alpha=0.8$ leads to a significant scatter and a strong bias in the analysis of the radio-FIR correlation, resulting from the broad and asymmetric spectral index distribution. When characterized by specific star formation rates, the starburst population (58%) dominates the 5 GHz radio source population, and the quiescent galaxy population (30%) follows a distinct trend in spectral index distribution and the radio-FIR correlation. Lastly, we offer suggestions on sensitivity and angular resolution for future ultra-deep surveys designed to trace the cosmic history of star formation and AGN activity using radio continuum as a probe.

We present an ALMA survey of dust continuum emission in a sample of 70 galaxies in the redshift range z=2-5 selected from the CANDELS GOODS-S field. Multi-Epoch Abundance Matching (MEAM) is used to define potential progenitors of a z = 0 galaxy of stellar mass 1.5 10^11 M_sun. Gas masses are derived from the 850um luminosity. Ancillary data from the CANDELS GOODS-S survey are used to derive the gas mass fractions. The results at z<=3 are mostly in accord with expectations: The detection rates are 75% for the z=2 redshift bin, 50% for the z=3 bin and 0% for z>=4. The average gas mass fraction for the detected z=2 galaxies is f_gas = 0.55+/-0.12 and f_gas = 0.62+/-0.15 for the z=3 sample. This agrees with expectations for galaxies on the star-forming main sequence, and shows that gas fractions have decreased at a roughly constant rate from z=3 to z=0. Stacked images of the galaxies not detected with ALMA give upper limits to f_gas of <0.08 and <0.15, for the z=2 and z=3 redshift bins. None of our galaxies in the z=4 and z=5 sample are detected and the upper limit from stacked images, corrected for low metallicity, is f_gas<0.66. We do not think that lower gas-phase metallicities can entirely explain the lower dust luminosities. We briefly consider the possibility of accretion of very low-metallicity gas to explain the absence of detectable dust emission in our galaxies at z>4.

We establish a robust statistical description of the star-forming galaxy population at the end of cosmic HI reionization ($5.0\le{}z\le6.6$) from a large sample of 52 galaxies with spectroscopically confirmed redshifts from the VIMOS UltraDeep Survey. We identify galaxies with Ly$\alpha$ either in absorption or in emission, at variance with most spectroscopic samples in the literature where Ly$\alpha$ emitters dominate. We find that star-forming galaxies at these redshifts are distributed along a main sequence in the stellar mass vs. SFR plane. We report a flat evolution of the sSFR(z) in 3<z<6 compared to lower redshift measurements. UV-continuum slopes vary with luminosity, with a large dispersion. We determing UV and Ly$\alpha$ luminosity functions using V$_{max}$ method and use them to derive star formation rate densities (SFRD). We find that both UV-derived and Ly$\alpha$-derived SFRDs are in excellent agreement after correcting Ly$\alpha$ luminosity density for IGM absorption. Our new SFRD measurements at a mean redshift z=5.6 confirm the steep decline of the SFRD at z>2. The bright end of the Ly$\alpha$ luminosity function has a high number density, indicating a significant star formation activity concentrated in the brightest Ly$\alpha$ emitters (LAE) at these redshifts. LAE with EW>25\AA ~contribute to about 75\% of the total UV-derived SFRD. While our analysis favors a low dust content in 5.0<z<6.6, uncertainties on the dust extinction correction and associated degeneracies in spectral fitting will remain an issue to estimate the total SFRD until future survey extending spectroscopy to the NIR rest-frame spectral domain, e.g. with JWST.

We use deep \textit{Hubble Space Telescope} spectroscopy to constrain the metallicities and (\editone{light-weighted}) ages of massive ($\log M_\ast/M_\odot\gtrsim10$) galaxies selected to have quiescent stellar populations at $1.0<z<1.8$. The data include 12--orbit depth coverage with the WFC3/G102 grism covering $\sim$ $8,000<\lambda<11,500$~\AA\, at a spectral resolution of $R\sim 210$ taken as part of the CANDELS Lyman-$\alpha$ Emission at Reionization (CLEAR) survey. At $1.0<z<1.8$, the spectra cover important stellar population features in the rest-frame optical. We simulate a suite of stellar population models at the grism resolution, fit these to the data for each galaxy, and derive posterior likelihood distributions for metallicity and age. We stack the posteriors for subgroups of galaxies in different redshift ranges that include different combinations of stellar absorption features. Our results give \editone{light-weighted ages of $t_{z \sim 1.1}= 3.2\pm 0.7$~Gyr, $t_{z \sim 1.2}= 2.2\pm 0.6$~Gyr, $t_{z\sim1.3}= 3.1\pm 0.6$~Gyr, and $t_{z\sim1.6}= 2.0 \pm 0.6$~Gyr, \editone{for galaxies at $z\sim 1.1$, 1.2, 1.3, and 1.6. This} implies that most of the massive quiescent galaxies at $1<z<1.8$ had formed $>68$\% of their stellar mass by a redshift of $z>2$}. The posteriors give metallicities of \editone{$Z_{z\sim1.1}=1.16 \pm 0.29$~$Z_\odot$, $Z_{z\sim1.2}=1.05 \pm 0.34$~$Z_\odot$, $Z_{z\sim1.3}=1.00 \pm 0.31$~$Z_\odot$, and $Z_{z\sim1.6}=0.95 \pm 0.39$~$Z_\odot$}. This is evidence that massive galaxies had enriched rapidly to approximately Solar metallicities as early as $z\sim3$.

We present a 69 arcmin$^2$ ALMA survey at 1.1mm, GOODS-ALMA, matching the deepest HST-WFC3 H-band part of the GOODS-South field. We taper the 0"24 original image with a homogeneous and circular synthesized beam of 0"60 to reduce the number of independent beams - thus reducing the number of purely statistical spurious detections - and optimize the sensitivity to point sources. We extract a catalogue of galaxies purely selected by ALMA and identify sources with and without HST counterparts down to a 5$\sigma$ limiting depth of H=28.2 AB (HST/WFC3 F160W). ALMA detects 20 sources brighter than 0.7 mJy in the 0"60 tapered mosaic (rms sensitivity =0.18 mJy/beam) with a purity greater than 80%. Among these detections, we identify three sources with no HST nor Spitzer-IRAC counterpart, consistent with the expected number of spurious galaxies from the analysis of the inverted image; their definitive status will require additional investigation. An additional three sources with HST counterparts are detected either at high significance in the higher resolution map, or with different detection-algorithm parameters ensuring a purity greater than 80%. Hence we identify in total 20 robust detections. Our wide contiguous survey allows us to push further in redshift the blind detection of massive galaxies with ALMA with a median redshift of $z$=2.92 and a median stellar mass of M$_{\star}$ = 1.1 $\times 10^{11}$M$_\odot$. Our sample includes 20% HST-dark galaxies (4 out of 20), all detected in the mid-infrared with IRAC. The near-infrared based photometric redshifts of two of them $z\sim$4.3 and 4.8) suggest that these sources have redshifts $z$>4. At least 40% of the ALMA sources host an X-ray AGN, compared to 14% for other galaxies of similar mass and redshift. The wide area of our ALMA survey provides lower values at the bright end of number counts than single-dish telescopes

We present deep spectroscopic follow-up observations of the Bremer Deep Field (BDF) where the two $z\sim$7 bright Ly$\alpha$ emitters (LAE) BDF521 and BDF3299 were previously discovered by Vanzella et al. (2011) and where a factor of $\sim$3-4 overdensity of faint LBGs has been found by Castellano et al. (2016). We confirm a new bright Ly$\alpha$ emitter, BDF2195, at the same redshift of BDF521, $z=7.008$, and at only $\sim$90 kpc physical distance from it, confirming that the BDF area is likely an overdense, reionized region. A quantitative assessment of the Ly$\alpha$ fraction shows that the number of detected bright emitters is much higher than the average found at z$\sim$7, suggesting a high Ly$\alpha$ transmission through the inter-galactic medium (IGM). However, the line visibility from fainter galaxies is at odds with this finding, since no Ly$\alpha$ emission is found in any of the observed candidates with $M_{UV}>$-20.25. This discrepancy can be understood either if some mechanism prevents Ly$\alpha$ emission from fainter galaxies within the ionized bubbles from reaching the observer, or if faint galaxies are located outside the reionized area and bright LAEs are solely responsible for the creation of their own HII regions. A thorough assessment of the nature of the BDF region and of its sources of re-ionizing radiation will be made possible by JWST spectroscopic capabilities.

We report evidence of environmental quenching among galaxies at redshift ~ 2, namely the probability that a galaxy quenches its star formation activity is enhanced in the regions of space in proximity of other quenched, more massive galaxies. The effect is observed as strong clustering of quiescent galaxies around quiescent galaxies on angular scales \theta < 20 arcsec, corresponding to a proper(comoving) scale of 168 (502) kpc at z = 2. The effect is observed only for quiescent galaxies around other quiescent galaxies; the probability to find star-forming galaxies around quiescent or around star-forming ones is consistent with the clustering strength of galaxies of the same mass and at the same redshift, as observed in dedicated studies of galaxy clustering. The effect is mass dependent in the sense that the quenching probability is stronger for galaxies of smaller mass ($\rm{M_*<10^{10} Msun}$) than for more massive ones, i.e. it follows the opposite trend with mass relative to gravitational galaxy clustering. The spatial scale where the effect is observed suggests these environments are massive halos, in which case the observed effect would likely be satellite quenching. The effect is also redshift dependent in that the clustering strength of quiescent galaxies around other quiescent galaxies at z = 1.6 is ~ 1.7 times larger than that of the galaxies with the same stellar mass at z = 2.6. This redshift dependence allows for a crude estimate of the time scale of environmental quenching of low-mass galaxies, which is in the range 1.5 - 4 Gyr, in broad agreement with other estimates and with our ideas on satellite quenching.

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

We report on the serendipitous discovery of a z=4.0, M1500=-22.20 star-forming galaxy (Ion3) showing copious Lyman continuum (LyC) leakage (~60% escaping), a remarkable multiple peaked Lya emission, and significant Lya radiation directly emerging at the resonance frequency. This is the highest redshift confirmed LyC emitter in which the ionising and Lya radiation possibly share a common ionised cavity (with N_HI<10^17.2 cm^-2). Ion3 is spatially resolved, it shows clear stellar winds signatures like the P-Cygni NV1240 profile, and has blue ultraviolet continuum (\beta = -2.5 +/- 0.25, F_\lambda~ \lambda^\beta) with weak low-ionisation interstellar metal lines. Deep VLT/HAWKI Ks and Spitzer/IRAC 3.6um and 4.5um imaging show a clear photometric signature of the Halpha line with equivalent width of 1000A rest-frame emerging over a flat continuum (Ks-4.5um ~ 0). From the SED fitting we derive a stellar mass of 1.5x10^9 Msun, SFR of 140 Msun/yr and age of ~10 Myr, with a low dust extinction, E(B-V)< 0.1, placing the source in the starburst region of the SFR-M^* plane. Ion3 shows similar properties of another LyC emitter previously discovered (z=3.21, Ion2, Vanzella et al. 2016). Ion3 (and Ion2) represents ideal high-redshift reference cases to guide the search for reionising sources at z>6.5 with JWST.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

We use the deep CANDELS observations in the GOODS North and South fields to revisit the correlations between stellar mass ($M_*$), star--formation rate (SFR) and morphology, and to introduce a fourth dimension, the mass-weighted stellar age, in galaxies at $1.2<z<4$. We do this by making new measures of $M_*$, $SFR$, and stellar age thanks to an improved SED fitting procedure that allows various star formation history for each galaxy. Like others, we find that the slope of the Main Sequence (MS) of star formation in the $(M_*;SFR)$ plane bends at high mass. We observe clear morphological differences among galaxies across the MS, which also correlate with stellar age. At all redshifts, galaxies that are quenching or quenched, and thus old, have high $\Sigma_1$ (the projected density within the central 1 kpc), while younger, star-forming galaxies span a much broader range of $\Sigma_1$, which includes the high values observed for quenched galaxies, but also extends to much lower values. As galaxies age and quench, the stellar age and the dispersion of $\Sigma_1$ for fixed values of $M_{*}$ shows two different regimes, one, at the low--mass end, where quenching might be driven by causes external to the galaxies; the other, at the high--mass end, where quenching is driven by internal causes, very likely the mass given the low scatter of $\Sigma_1$ (mass quenching). We suggest that the monotonic increase of central density as galaxies grow is one manifestation of a more general phenomenon of structural transformation that galaxies undergo as they evolve.

Studying giant star-forming clumps in distant galaxies is important to understand galaxy formation and evolution. At present, however, observers and theorists have not reached a consensus on whether the observed "clumps" in distant galaxies are the same phenomenon that is seen in simulations. In this paper, as a step to establish a benchmark of direct comparisons between observations and theories, we publish a sample of clumps constructed to represent the commonly observed "clumps" in the literature. This sample contains 3193 clumps detected from 1270 galaxies at $0.5 \leq z < 3.0$. The clumps are detected from rest-frame UV images, as described in our previous paper. Their physical properties, e.g., rest-frame color, stellar mass (M*), star formation rate (SFR), age, and dust extinction, are measured by fitting the spectral energy distribution (SED) to synthetic stellar population models. We carefully test the procedures of measuring clump properties, especially the method of subtracting background fluxes from the diffuse component of galaxies. With our fiducial background subtraction, we find a radial clump U-V color variation, where clumps close to galactic centers are redder than those in outskirts. The slope of the color gradient (clump color as a function of their galactocentric distance scaled by the semi-major axis of galaxies) changes with redshift and M* of the host galaxies: at a fixed M*, the slope becomes steeper toward low redshift, and at a fixed redshift, it becomes slightly steeper with M*. Based on our SED-fitting, this observed color gradient can be explained by a combination of a negative age gradient, a negative E(B-V) gradient, and a positive specific star formation rate gradient of the clumps. We also find that the color gradients of clumps are steeper than those of intra-clump regions. [Abridged]

We analyze the CIII]-1908\AA emission properties in a sample of 3899 star-forming galaxies (SFGs) at 2<z<3.8 from the VIMOS Ultra-Deep Survey (VUDS). We find a median EW(CIII])=2.0$\pm$0.2 to 2.2$\pm$0.2\AA ~for the SFG population at 2<z<3 and 3<z<4, resp. About 24% of SFGs are showing EW(CIII])>3\AA, including 20% with 3<EW(CIII])<10\AA ~and 4% with strong emission EW(CIII])>10\AA. A significant fraction of 1.2% of SFGs presents strong CIII] emission 20<EW(CIII])<40\AA. This makes CIII] the second most-frequent emission line in the UV-rest spectra of SFGs after Lyman-$\alpha$. We find a large dispersion in the weak EW(CIII]) - EW(Ly\alpha) correlation, with galaxies showing strong CIII] and no Ly$\alpha$, and vice-versa. SFGs with 10<EW(CIII])<0\AA ~present strong emission lines including CIV-1549, HeII-1640, OIII-1664, but also weaker emission of highly ionized elements SiIV-1403, NIV-1485, NIII-1750, or SiIII-1888, indicating a hard radiation field. We present a broad range of observational evidence supporting the presence of AGN in the strong CIII] emitting population. As EW(CIII]) is rising, we identify powerful outflows with velocities up to 1014 km/s. The strongest CIII] emitters are preferentially located below the SFGs main sequence, with the SFR reduced by x2. In addition, the median stellar age of the strongest emitters is ~0.8 Gyr, three times that of galaxies with EW(CIII])<10\AA. Spectral line analysis presented in a joint paper by Nakajima et al. (2017) confirms that the strongest emitters require an AGN. We conclude that these properties are indicative of AGN feedback in SFGs at 2<z<3.8, contributing to star-formation quenching. We find that quenching timescales of 0.25-0.5x10^9 years are necessary for such AGN feedback to turn part of the star-forming galaxy population with Mstar>10^10 MSun at z~3 into the population of quiescent galaxies observed at redshift z~1-2.

We examine the fraction of massive ($M_{*}>10^{10} M_{\odot}$), compact star-forming galaxies (cSFGs) that host an active galactic nucleus (AGN) at $z\sim2$. These cSFGs are likely the direct progenitors of the compact quiescent galaxies observed at this epoch, which are the first population of passive galaxies to appear in large numbers in the early Universe. We identify cSFGs that host an AGN using a combination of Hubble WFC3 imaging and Chandra X-ray observations in four fields: the Chandra Deep Fields, the Extended Groth Strip, and the UKIDSS Ultra Deep Survey field. We find that $39.2^{+3.9}_{-3.6}$\% (65/166) of cSFGs at $1.4<z<3.0$ host an X-ray detected AGN. This fraction is 3.2 times higher than the incidence of AGN in extended star-forming galaxies with similar masses at these redshifts. This difference is significant at the $6.2\sigma$ level. Our results are consistent with models in which cSFGs are formed through a dissipative contraction that triggers a compact starburst and concurrent growth of the central black hole. We also discuss our findings in the context of cosmological galaxy evolution simulations that require feedback energy to rapidly quench cSFGs. We show that the AGN fraction peaks precisely where energy injection is needed to reproduce the decline in the number density of cSFGs with redshift. Our results suggest that the first abundant population of massive, quenched galaxies emerged directly following a phase of elevated supermassive black hole growth and further hints at a possible connection between AGN and the rapid quenching of star formation in these galaxies.

We present a detailed spectral analysis of the brightest Active Galactic Nuclei (AGN) identified in the 7Ms Chandra Deep Field South (CDF-S) survey over a time span of 16 years. Using a model of an intrinsically absorbed power-law plus reflection, with possible soft excess and narrow Fe K$\alpha$ line, we perform a systematic X-ray spectral analysis, both on the total 7Ms exposure and in four different periods with lengths of 2-21 months. With this approach, we not only present the power-law slopes, column densities $N_H$, observed fluxes, and absorption-corrected 2-10~keV luminosities $L_X$ for our sample of AGNs, but also identify significant spectral variabilities among them on time scales of years. We find that the $N_H$ variabilities can be ascribed to two different types of mechanisms, either flux-driven or flux-independent. We also find that the correlation between the narrow Fe line EW and $N_H$ can be well explained by the continuum suppression with increasing $N_H$. Accounting for the sample incompleteness and bias, we measure the intrinsic distribution of $N_H$ for the CDF-S AGN population and present re-selected subsamples which are complete with respect to $N_H$. The $N_H$-complete subsamples enable us to decouple the dependences of $N_H$ on $L_X$ and on redshift. Combining our data with that from C-COSMOS, we confirm the anti-correlation between the average $N_H$ and $L_X$ of AGN, and find a significant increase of the AGN obscured fraction with redshift at any luminosity. The obscured fraction can be described as $f_{obscured}\thickapprox 0.42\ (1+z)^{0.60}$.

We investigate the environmental quenching of galaxies, especially those with stellar masses (M*)$<10^{9.5} M_\odot$, beyond the local universe. Essentially all local low-mass quenched galaxies (QGs) are believed to live close to massive central galaxies, which is a demonstration of environmental quenching. We use CANDELS data to test {\it whether or not} such a dwarf QG--massive central galaxy connection exists beyond the local universe. To this purpose, we only need a statistically representative, rather than a complete, sample of low-mass galaxies, which enables our study to $z\gtrsim1.5$. For each low-mass galaxy, we measure the projected distance ($d_{proj}$) to its nearest massive neighbor (M*$>10^{10.5} M_\odot$) within a redshift range. At a given redshift and M*, the environmental quenching effect is considered to be observed if the $d_{proj}$ distribution of QGs ($d_{proj}^Q$) is significantly skewed toward lower values than that of star-forming galaxies ($d_{proj}^{SF}$). For galaxies with $10^{8} M_\odot < M* < 10^{10} M_\odot$, such a difference between $d_{proj}^Q$ and $d_{proj}^{SF}$ is detected up to $z\sim1$. Also, about 10\% of the quenched galaxies in our sample are located between two and four virial radii ($R_{Vir}$) of the massive halos. The median projected distance from low-mass QGs to their massive neighbors, $d_{proj}^Q / R_{Vir}$, decreases with satellite M* at $M* \lesssim 10^{9.5} M_\odot$, but increases with satellite M* at $M* \gtrsim 10^{9.5} M_\odot$. This trend suggests a smooth, if any, transition of the quenching timescale around $M* \sim 10^{9.5} M_\odot$ at $0.5<z<1.0$.

We present a 0.4-8$\mu$m multi-wavelength photometric catalog in the Extended Groth Strip (EGS) field. This catalog is built on the Hubble Space Telescope (HST) WFC3 and ACS data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), and it incorporates the existing HST data from the All-wavelength Extended Groth strip International Survey (AEGIS) and the 3D-HST program. The catalog is based on detections in the F160W band reaching a depth of F160W=26.62 AB (90% completeness, point-sources). It includes the photometry for 41457 objects over an area of $\approx 206$ arcmin$^2$ in the following bands: HST ACS F606W and F814W; HST WFC3 F125W, F140W and F160W; CFHT/Megacam $u^*$, $g'$, $r'$, $i'$ and $z'$; CFHT/WIRCAM $J$, $H$ and $K_\mathrm{S}$; Mayall/NEWFIRM $J1$, $J2$, $J3$, $H1$, $H2$, $K$; Spitzer IRAC $3.6\mu$m, $4.5\mu$m, $5.8\mu$m and $8.0\mu$m. We are also releasing value-added catalogs that provide robust photometric redshifts and stellar mass measurements. The catalogs are publicly available through the CANDELS repository.

We directly detect dust emission in an optically-detected, multiply-imaged galaxy lensed by the Frontier Fields cluster MACSJ0717.5+3745. We detect two images of the same galaxy at 1.1mm with the AzTEC camera on the Large Millimeter Telescope leaving no ambiguity in the counterpart identification. This galaxy, MACS071_Az9, is at z>4 and the strong lensing model (mu=7.5) allows us to calculate an intrinsic IR luminosity of 9.7e10 Lsun and an obscured star formation rate of 14.6 +/- 4.5 Msun/yr. The unobscured star formation rate from the UV is only 4.1 +/- 0.3 Msun/yr which means the total star formation rate (18.7 +/- 4.5 Msun/yr) is dominated (75-80%) by the obscured component. With an intrinsic stellar mass of only 6.9e9Msun, MACS0717_Az9 is one of only a handful of z>4 galaxies at these lower masses that is detected in dust emission. This galaxy lies close to the estimated star formation sequence at this epoch. However, it does not lie on the dust obscuration relation (IRX-beta) for local starburst galaxies and is instead consistent with the Small Magellanic Cloud (SMC) attenuation law. This remarkable lower mass galaxy showing signs of both low metallicity and high dust content may challenge our picture of dust production in the early Universe.

We report the detection of morphology dependent stellar age in massive quenched galaxies (QGs) at z~1.2. The sense of the dependence is that compact QGs are 0.5-2 Gyr older than normal-sized ones. The evidence comes from three different age indicators, Dn4000, Hdelta and fits to spectral synthesis models, applied to their stacked optical spectra. All age indicators consistently show that the stellar populations of compact QGs are older than their normally-sized counterparts. We detect weak [OII] emission in a fraction of QGs, and the strength of the line, when present, is similar between the two samples; however, compact galaxies exhibit significantly lower frequency of [OII] emission than normal ones. A fraction of both samples are individually detected in 7 Ms Chandra X-ray images (luminosities$\sim10^{40}$-$10^{41}$ erg/sec). 7 Ms stacks of non-detected galaxies show similarly low luminosities in the soft band only, consistent with a hot gas origin for the X-ray emission. While both [OII] emitters and non-emitters are also X-ray sources among normal galaxies, no compact galaxy with [OII] emission is an X-ray source, arguing against an AGN powering the line in compact galaxies. We interpret the [OII] properties as further evidence that compact galaxies are older and further along into the process of quenching star-formation and suppressing gas accretion. Finally, we argue that the older age of compact QGs is evidence of progenitor bias: compact QGs simply reflect the smaller sizes of galaxies at their earlier quenching epoch, with stellar density most likely having nothing directly to do with cessation of star-formation.

The reionization of the Universe is one of the most important topics of present day astrophysical research. The most plausible candidates for the reionization process are star-forming galaxies, which according to the predictions of the majority of the theoretical and semi-analytical models should dominate the HI ionizing background at z~3. We aim at measuring the Lyman continuum escape fraction, which is one of the key parameters to compute the contribution of star-forming galaxies to the UV background. We have used ultra-deep U-band imaging (U=30.2mag at 1sigma) by LBC/LBT in the CANDELS/GOODS-North field, as well as deep imaging in COSMOS and EGS fields, in order to estimate the Lyman continuum escape fraction of 69 star-forming galaxies with secure spectroscopic redshifts at 3.27<z<3.40 to faint magnitude limits (L=0.2L*, or equivalently M1500~-19). We have measured through stacks a stringent upper limit (<1.7% at 1sigma) for the relative escape fraction of HI ionizing photons from bright galaxies (L>L*), while for the faint population (L=0.2L*) the limit to the escape fraction is ~10%. We have computed the contribution of star-forming galaxies to the observed UV background at z~3 and we have found that it is not enough to keep the Universe ionized at these redshifts, unless their escape fraction increases significantly (>10%) at low luminosities (M1500>-19). We compare our results on the Lyman continuum escape fraction of high-z galaxies with recent estimates in the literature and discuss future prospects to shed light on the end of the Dark Ages. In the future, strong gravitational lensing will be fundamental to measure the Lyman continuum escape fraction down to faint magnitudes (M1500~-16) which are inaccessible with the present instrumentation on blank fields.

In this paper, we investigate the relationship between star formation and structure, using a mass-complete sample of 27,893 galaxies at $0.5<z<2.5$ selected from 3D-HST. We confirm that star-forming galaxies are larger than quiescent galaxies at fixed stellar mass (M$_{\star}$). However, in contrast with some simulations, there is only a weak relation between star formation rate (SFR) and size within the star-forming population: when dividing into quartiles based on residual offsets in SFR, we find that the sizes of star-forming galaxies in the lowest quartile are 0.27$\pm$0.06 dex smaller than the highest quartile. We show that 50% of star formation in galaxies at fixed M$_{\star}$ takes place within a narrow range of sizes (0.26 dex). Taken together, these results suggest that there is an abrupt cessation of star formation after galaxies attain particular structural properties. Confirming earlier results, we find that central stellar density within a 1 kpc fixed physical radius is the key parameter connecting galaxy morphology and star formation histories: galaxies with high central densities are red and have increasingly lower SFR/M$_{\star}$, whereas galaxies with low central densities are blue and have a roughly constant (higher) SFR/M$_{\star}$ at a given redshift. We find remarkably little scatter in the average trends and a strong evolution of $>$0.5 dex in the central density threshold correlated with quiescence from $z\sim0.7-2.0$. Neither a compact size nor high-$n$ are sufficient to assess the likelihood of quiescence for the average galaxy; rather, the combination of these two parameters together with M$_{\star}$ results in a unique quenching threshold in central density/velocity.

Deep observations are revealing a growing number of young galaxies in the first billion year of cosmic time. Compared to typical galaxies at later times, they show more extreme emission-line properties, higher star formation rates, lower masses, and smaller sizes. However, their faintness precludes studies of their chemical abundances and ionization conditions, strongly limiting our understanding of the physics driving early galaxy build-up and metal enrichment. Here we study a rare population of UV-selected, sub$-L^{*}$(z=3) galaxies at redshift 2.4$<z<$3.5 that exhibit all the rest-frame properties expected from primeval galaxies. These low-mass, highly-compact systems are rapidly-forming galaxies able to double their stellar mass in only few tens million years. They are characterized by very blue UV spectra with weak absorption features and bright nebular emission lines, which imply hard radiation fields from young hot massive stars. Their highly-ionized gas phase has strongly sub-solar carbon and oxygen abundances, with metallicities more than a factor of two lower than that found in typical galaxies of similar mass and star formation rate at $z\lesssim$2.5. These young galaxies reveal an early and short stage in the assembly of their galactic structures and their chemical evolution, a vigorous phase which is likely to be dominated by the effects of gas-rich mergers, accretion of metal-poor gas and strong outflows.

Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the CANDELS dataset. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected because the Universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$ Gyr). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms that operate on different timescales.

We present new results on [CII]158$\mu$ m emission from four galaxies in the reionization epoch. These galaxies were previously confirmed to be at redshifts between 6.6 and 7.15 from the presence of the Ly$\alpha$ emission line in their spectra. The Ly$\alpha$ emission line is redshifted by 100-200 km/s compared to the systemic redshift given by the [CII] line. These velocity offsets are smaller than what is observed in z~3 Lyman break galaxies with similar UV luminosities and emission line properties. Smaller velocity shifts reduce the visibility of Ly$\alpha$ and hence somewhat alleviate the need for a very neutral IGM at z~7 to explain the drop in the fraction of Ly$\alpha$ emitters observed at this epoch. The galaxies show [CII] emission with L[CII]=0.6-1.6 x10$^8 L_\odot$: these luminosities place them consistently below the SFR-L[CII] relation observed for low redshift star forming and metal poor galaxies and also below z =5.5 Lyman break galaxies with similar star formation rates. We argue that previous undetections of [CII] in z~7 galaxies with similar or smaller star formation rates are due to selection effects: previous targets were mostly strong Ly$\alpha$ emitters and therefore probably metal poor systems, while our galaxies are more representative of the general high redshift star forming population .

In this paper we investigate the impact of different star formation histories (SFHs) on the relation between stellar mass M$_{*}$ and star formation rate (SFR) using a sample of galaxies with reliable spectroscopic redshift zspec>2 drawn from the VIMOS Ultra-Deep Survey (VUDS). We produce an extensive database of dusty model galaxies, calculated starting from the new library of single stellar population (SSPs) models presented in Cassara' et al. 2013 and weighted by a set of 28 different SFHs based on the Schmidt function, and characterized by different ratios of the gas infall time scale $\tau_{infall}$ to the star formation efficiency $\nu$. The treatment of dust extinction and re-emission has been carried out by means of the radiative transfer calculation. The spectral energy distribution (SED) fitting technique is performed by using GOSSIP+, a tool able to combine both photometric and spectroscopic information to extract the best value of the physical quantities of interest, and to consider the Intergalactic Medium (IGM) attenuation as a free parameter. We find that the main contribution to the scatter observed in the $SFR-M_{*}$ plane is the possibility of choosing between different families of SFHs in the SED fitting procedure, while the redshift range plays a minor role. The majority of the galaxies, at all cosmic times, are best-fit by models with SFHs characterized by a high $\tau_{\rm infall}/\nu$ ratio. We discuss the reliability of the presence of a small percentage of dusty and highly star forming galaxies, in the light of their detection in the FIR.

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

We measure galaxy sizes on a sample of $\sim1200$ galaxies with confirmed spectroscopic redshifts $2 \leq z_{spec} \leq 4.5$ in the VIMOS Ultra Deep Survey (VUDS), representative of star-forming galaxies with $i_\mathrm{AB} \leq 25$. We first derive galaxy sizes applying a classical parametric profile fitting method using GALFIT. We then measure the total pixel area covered by a galaxy above a given surface brightness threshold, which overcomes the difficulty of measuring sizes of galaxies with irregular shapes. We then compare the results obtained for the equivalent circularized radius enclosing 100\% of the measured galaxy light $r_T^{100}$ to those obtained with the effective radius $r_{e,\mathrm{circ}}$ measured with GALFIT. We find that the sizes of galaxies computed with our non-parametric approach span a large range but remain roughly constant on average with a median value $r_T^{100}\sim2.2$ kpc for galaxies with $2<z<4.5$. This is in stark contrast with the strong downward evolution of $r_e$ with increasing redshift, down to sizes of $<1$ kpc at $z\sim4.5$. We analyze the difference and find that parametric fitting of complex, asymmetric, multi-component galaxies is severely underestimating their sizes. By comparing $r_T^{100}$ with physical parameters obtained through SED fitting we find that the star-forming galaxies that are the largest at any redshift are, on average, more massive and more star-forming. We discover that galaxies present more concentrated light profiles as we move towards higher redshifts. We interpret these results as the signature of several, possibly different, evolutionary paths of galaxies in their early stages of assembly, including major and minor merging or star-formation in multiple bright regions. (abridged)

This paper describes the first data release (DR1) of the VIMOS Ultra Deep Survey (VUDS). The DR1 includes all low-resolution spectroscopic data obtained in 276.9 arcmin2 of the CANDELS-COSMOS and CANDELS-ECFDS survey areas, including accurate spectroscopic redshifts z_spec and individual spectra obtained with VIMOS on the ESO-VLT. A total of 698 objects have a measured redshift, with 677 galaxies, two type-I AGN and a small number of 19 contaminating stars. The targets of the spectroscopic survey are selected primarily on the basis of their photometric redshifts to ensure a broad population coverage. About 500 galaxies have z_spec>2, 48 with z_spec>4, and the highest reliable redshifts reach beyond z_spec=6. This dataset approximately doubles the number of galaxies with spectroscopic redshifts at z>3 in these fields. We discuss the general properties of the sample in terms of the spectroscopic redshift distribution, the distribution of Lyman-alpha equivalent widths, and physical properties including stellar masses M_star and star formation rates (SFR) derived from spectral energy distribution fitting with the knowledge of z_spec. We highlight the properties of the most massive star-forming galaxies, noting the large range in spectral properties, with Lyman-alpha in emission or in absorption, and in imaging properties with compact, multi-component or pair morphologies. We present the catalogue database and data products. All data are publicly available and can be retrieved from a dedicated query-based database available at http://cesam.lam.fr/vuds.

Star-forming galaxies are considered to be the leading candidate sources that dominate the cosmic reionization at z>7, and the search for analogs at moderate redshift showing Lyman continuum (LyC) leakage is currently a active line of research. We have observed a star-forming galaxy at z=3.2 with Hubble/WFC3 in the F336W filter, corresponding to the 730-890A rest-frame, and detect LyC emission. This galaxy is very compact and also has large Oxygen ratio [OIII]5007/[OII]3727 (>=10). No nuclear activity is revealed from optical/near-infrared spectroscopy and deep multi-band photometry (including the 6Ms X-ray, Chandra). The measured escape fraction of ionizing radiation spans the range 50-100\%, depending on the IGM attenuation. The LyC emission is detected at S/N=10 with m(F336W)=27.57+/-0.11 and it is spatially unresolved, with effective radius R_e<200pc. Predictions from photoionization and radiative transfer models are in line with the properties reported here, indicating that stellar winds and supernova explosions in a nucleated star-forming region can blow cavities generating density-bounded conditions compatible with optically thin media. Irrespective to the nature of the ionizing radiation, spectral signatures of these sources over the entire electromagnetic spectrum are of central importance for their identification during the epoch of reionization, when the LyC is unobservable. Intriguingly, the Spitzer/IRAC photometric signature of intense rest-frame optical emissions ([OIII]+Hbeta) observed recently at z~7.5-8.5 is similar to what is observed in this galaxy. Only the James Webb Space Telescope will measure optical line ratios at z>7 allowing a direct comparison with lower redshift LyC emitters, as reported here.

We present the analysis of deep HST multi-band imaging of the BDF field specifically designed to identify faint companions around two of the few Ly-alpha emitting galaxies spectroscopically confirmed at z~7 (Vanzella et al. 2011). Although separated by only 4.4 proper Mpc these galaxies cannot generate HII regions large enough to explain visibility of their Ly-alpha line, thus requiring a population of fainter ionizing sources in their vicinity. We use deep HST and VLT-Hawk-I data to select z~7 Lyman break galaxies around the emitters. We select 6 new robust z~7 LBGs at Y~26.5-27.5 whose average spectral energy distribution is consistent with the objects being at the redshift of the close-by Ly-alpha emitters. The resulting number density of z~7 LBGs in the BDF field is a factor ~3-4 higher than expected in random pointings of the same size. We compare these findings with cosmological hydrodynamic plus radiative transfer simulations of a universe with a half neutral IGM: we find that indeed Ly-alpha emitter pairs are only found in completely ionized regions characterized by significant LBG overdensities. Our findings match the theoretical prediction that the first ionization fronts are generated within significant galaxy overdensities and support a scenario where faint, "normal" star-forming galaxies are responsible for reionization.

Aim: We aim to measure the LyC signal from a sample of sources in the Chandra deep field south. We collect star-forming galaxies (SFGs) and active galactic nuclei (AGN) with accurate spectroscopic redshifts, for which Hubble Space Telescope (HST) coverage and multi-wavelength photometry are available. Method: We selected a sample of about 200 sources at z~3. Taking advantage of HST resolution, we applied a careful cleaning procedure and rejected sources showing nearby clumps with different colours, which could be lower-z interlopers. Our clean sample consisted of 86 SFGs (including 19 narrow-band selected Lya emitters) and 8 AGN (including 6 detected in X-rays). We measured the LyC flux from aperture photometry in four narrow-band filters covering wavelengths below a 912 A rest frame (3.11<z<3.53). We estimated the ratio between ionizing (LyC flux) and 1400 A non-ionizing emissions for AGN and galaxies. Results: By running population synthesis models, we assume an average intrinsic L(1400 A)/L(900 A) ratio of 5 as the representative value for our sample. With this value and an average treatment of the lines of sight of the inter-galactic medium, we estimate the LyC escape fraction relative to the intrinsic value (fesc_rel(LyC)). We do not directly detect ionizing radiation from any individual SFG, but we are able to set a 1(2)sigma upper limit of fesc_rel(LyC)<12(24)%. This result is consistent with other non-detections published in the literature. No meaningful limits can be calculated for the sub-sample of Lya emitters. We obtain one significant direct detection for an AGN at z=3.46, with fesc_rel(LyC) = (72+/-18)%. Conclusions: Our upper limit on fescrel(LyC) implies that the SFGs studied here do not present either the physical properties or the geometric conditions suitable for efficient LyC-photon escape.

The aim of this paper is to investigate spectral and photometric properties of 854 faint ($i_{AB}$<~25 mag) star-forming galaxies (SFGs) at 2<z<2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes ($\beta$) as well as Ly$\alpha$ equivalent widths (EW). On average, the spectroscopically measured $\beta$ (-1.36$\pm$0.02), is comparable to the photometrically measured $\beta$ (-1.32$\pm$0.02), and has smaller measurement uncertainties. The positive correlation of $\beta$ with the Spectral Energy Distribution (SED)-based measurement of dust extinction, E$_{\rm s}$(B-V), emphasizes the importance of $\beta$ as an alternative dust indicator at high redshifts. To make a proper comparison, we divide these SFGs into three subgroups based on their rest-frame Ly$\alpha$ EW: SFGs with no Ly$\alpha$ emission (SFG$_{\rm N}$; EW$\le$0\AA), SFGs with Ly$\alpha$ emission (SFG$_{\rm L}$; EW$>$0\AA), and Ly$\alpha$ emitters (LAEs; EW$\ge$20\AA). The fraction of LAEs at these redshifts is $\sim$10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SFG$_{\rm N}$, SFG$_{\rm L}$ and LAE samples. For the luminosities probed here ($\sim$L$^*$), we find that galaxies with and without Ly$\alpha$ in emission have small but significant differences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant compared to the SFR and E$_{\rm s}$(B-V). [abridged]

We introduce a new color-selection technique to identify high-redshift, massive galaxies that are systematically missed by Lyman-break selection. The new selection is based on the H_{160} and IRAC 4.5um bands, specifically H - [4.5] > 2.25 mag. These galaxies, dubbed "HIEROs", include two major populations that can be separated with an additional J - H color. The populations are massive and dusty star-forming galaxies at z > 3 (JH-blue) and extremely dusty galaxies at z < 3 (JH-red). The 350 arcmin^2 of the GOODS-N and GOODS-S fields with the deepest HST/WFC3 and IRAC data contain 285 HIEROs down to [4.5] < 24 mag. We focus here primarily on JH-blue (z > 3) HIEROs, which have a median photometric redshift z ~4.4 and stellar massM_{*}~10^{10.6} Msun, and are much fainter in the rest-frame UV than similarly massive Lyman-break galaxies (LBGs). Their star formation rates (SFRs) reaches ~240 Msun yr^{-1} leading to a specific SFR, sSFR ~4.2 Gyr^{-1}, suggesting that the sSFRs for massive galaxies continue to grow at z > 2 but at a lower growth rate than from z=0 to z=2. With a median half-light radius of 2 kpc, including ~20% as compact as quiescent galaxies at similar redshifts, JH-blue HIEROs represent perfect star-forming progenitors of the most massive (M_{*} > 10^{11.2} Msun) compact quiescent galaxies at z ~ 3 and have the right number density. HIEROs make up ~60% of all galaxies with M_{*} > 10^{10.5} Msun identified at z > 3 from their photometric redshifts. This is five times more than LBGs with nearly no overlap between the two populations. While HIEROs make up 15-25% of the total SFR density at z ~ 4-5, they completely dominate the SFR density taking place in M_{*} >10^{10.5} Msun galaxies, and are therefore crucial to understanding the very early phase of massive galaxy formation.

[Abridged] We investigate the physical properties of a Lyman continuum emitter candidate at $z=3.212$ with photometric coverage from $U$ to MIPS 24$\mu$m band and VIMOS/VLT and MOSFIRE/Keck spectroscopy. Investigation of the UV spectrum confirms a direct spectroscopic detection of the Lyman continuum emission with $S/N>5$. Non-zero Ly$\alpha$ flux at the systemic redshift and high Lyman-$\alpha$ escape fraction suggest a low HI column density. The weak C and Si low-ionization absorption lines are also consistent with a low covering fraction along the line of sight. The [OIII]$\lambda\lambda4959,5007+\mathrm{H}\beta$ equivalent width is one of the largest reported for a galaxy at $z>3$ ($\mathrm{EW}([\mathrm{OIII}]\lambda\lambda4959,5007+\mathrm{H}\beta) \simeq 1600\AA$, rest-frame) and the NIR spectrum shows that this is mainly due to an extremely strong [OIII] emission. The large observed [OIII]/[OII] ratio ($>10$) and high ionization parameter are consistent with prediction from photoionization models in case of a density-bounded nebula scenario. Furthermore, the $\mathrm{EW}([\mathrm{OIII}]\lambda\lambda4959,5007+\mathrm{H}\beta)$ is comparable to recent measurements reported at $z\sim7-9$, in the reionization epoch. We also investigate the possibility of an AGN contribution to explain the ionizing emission but most of the AGN identification diagnostics suggest that stellar emission dominates instead. This source is currently the first high-$z$ example of a Lyman continuum emitter exhibiting indirect and direct evidences of a Lyman continuum leakage and having physical properties consistent with theoretical expectation from Lyman continuum emission from a density-bounded nebula.

The Lyman continuum (LyC) flux escaping from high-z galaxies into the IGM is a fundamental quantity to understand the physical processes involved in the reionization epoch. We have investigated a sample of star-forming galaxies at z~3.3 in order to search for possible detections of LyC photons escaping from galaxy halos. UV deep imaging in the COSMOS field obtained with the prime focus camera LBC at the LBT telescope was used together with a catalog of spectroscopic redshifts obtained by the VIMOS Ultra Deep Survey (VUDS) to build a sample of 45 galaxies at z~3.3 with L>0.5L*. We obtained deep LBC images of galaxies with spectroscopic redshifts in the interval 3.27<z<3.40 both in the R and deep U bands. A sub-sample of 10 galaxies apparently shows escape fractions>28% but a detailed analysis of their properties reveals that, with the exception of two marginal detections (S/N~2) in the U band, all the other 8 galaxies are most likely contaminated by the UV flux of low-z interlopers located close to the high-z targets. The average escape fraction derived from the stacking of the cleaned sample was constrained to fesc_rel<2%. The implied HI photo-ionization rate is a factor two lower than that needed to keep the IGM ionized at z~3, as observed in the Lyman forest of high-z QSO spectra or by the proximity effect. These results support a scenario where high redshift, relatively bright (L>0.5L*) star-forming galaxies alone are unable to sustain the level of ionization observed in the cosmic IGM at z~3. Star-forming galaxies at higher redshift and at fainter luminosities (L<<L*) can be the major contributors to the reionization of the Universe only if their physical properties are subject to rapid changes from z~3 to z~6-10. Alternatively, ionizing sources could be discovered looking for fainter sources among the AGN population at high-z.

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $\Delta\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet (NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer break of galaxies at z~0.8-3.4, and one of the breaks over the rest of the redshift range. Our comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, and improves the scatter at z<0.5 and at z>2. The panchromatic coverage of the UDF from the NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction available.

We present new results from our search for z~7 galaxies from deep spectroscopic observations of candidate z-dropouts in the CANDELS fields. Despite the extremely low flux limits achieved by our sensitive observations, only 2 galaxies have robust redshift identifications, one from its Lyalpha emission line at z=6.65, the other from its Lyman-break, i.e. the continuum discontinuity at the Lyalpha wavelength consistent with a redshift 6.42, but with no emission line. In addition, for 23 galaxies we present deep limits in the Lyalpha EW derived from the non detections in ultra-deep observations. Using this new data as well as previous samples, we assemble a total of 68 candidate z~7 galaxies with deep spectroscopic observations, of which 12 have a line detection. With this much enlarged sample we can place solid constraints on the declining fraction of Ly$\alpha$ emission in z~7 Lyman break galaxies compared to z~6, both for bright and faint galaxies. Applying a simple analytical model, we show that the present data favor a patchy reionization process rather than a smooth one.

We report the discovery of 31 low-luminosity (-14.5 > M_{AB}(B) > -18.8), extreme emission line galaxies (EELGs) at 0.2 < z < 0.9 identified by their unusually high rest-frame equivalent widths (100 < EW[OIII] < 1700 A) as part of the VIMOS Ultra Deep Survey (VUDS). VIMOS optical spectra of unprecedented sensitivity ($I_{AB}$ ~ 25 mag) along with multiwavelength photometry and HST imaging are used to investigate spectrophotometric properties of this unique sample and explore, for the first time, the very low stellar mass end (M* < 10^8 M$_{\odot}$) of the luminosity-metallicity (LZR) and mass-metallicity (MZR) relations at z < 1. Characterized by their extreme compactness (R50 < 1 kpc), low stellar mass and enhanced specific star formation rates (SFR/M* ~ 10^{-9} - 10^{-7} yr^{-1}), the VUDS EELGs are blue dwarf galaxies likely experiencing the first stages of a vigorous galaxy-wide starburst. Using T_e-sensitive direct and strong-line methods, we find that VUDS EELGs are low-metallicity (7.5 < 12+log(O/H) < 8.3) galaxies with high ionization conditions, including at least three EELGs showing HeII 4686A emission and four EELGs of extremely metal-poor (<10% solar) galaxies. The LZR and MZR followed by EELGs show relatively large scatter, being broadly consistent with the extrapolation toward low luminosity and mass from previous studies at similar redshift. However, we find evidences that galaxies with younger and more vigorous star formation -- as characterized by their larger EWs, ionization and sSFR -- tend to be more metal-poor at a given stellar mass.

We present the Swift X-ray Cluster Survey (SWXCS) catalog obtained using archival data from the X-ray telescope (XRT) on board the Swift satellite acquired from 2005 to 2012, extending the first release of the SWXCS. The catalog provides positions, soft fluxes, and, when possible, optical counterparts for a flux-limited sample of X-ray group and cluster candidates. We consider the fields with Galactic latitude |b| > 20 degree to avoid high HI column densities. We discard all of the observations targeted at groups or clusters of galaxies, as well as particular extragalactic fields not suitable to search for faint extended sources. We finally select ~3000 useful fields covering a total solid angle of ~400 degree^2. We identify extended source candidates in the soft-band (0.5-2keV) images of these fields using the software EXSdetect, which is specifically calibrated for the XRT data. Extensive simulations are used to evaluate contamination and completeness as a function of the source signal, allowing us to minimize the number of spurious detections and to robustly assess the selection function. Our catalog includes 263 candidate galaxy clusters and groups down to a flux limit of 7E-15 erg/cm^2/s in the soft band, and the logN-logS is in very good agreement with previous deep X-ray surveys. The final list of sources is cross-correlated with published optical, X-ray, and SZ catalogs of clusters. We find that 137 sources have been previously identified as clusters, while 126 are new detections. Currently, we have collected redshift information for 158 sources (60% of the entire sample). Once the optical follow-up and the X-ray spectral analysis of the sources are complete, the SWXCS will provide a large and well-defined catalog of groups and clusters of galaxies to perform statistical studies of cluster properties and tests of cosmological models.

We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$).

We aim to investigate the effect of the escaping ionizing radiation on the color selection of high redshift galaxies and identify candidate Lyman continuum (LyC) emitters. The intergalactic medium prescription of Inoue et al.(2014) and galaxy synthesis models of Bruzual&Charlot (2003) have been used to properly treat the ultraviolet stellar emission, the stochasticity of the intergalactic transmission and mean free path in the ionizing regime. Color tracks are computed by turning on/off the escape fraction of ionizing radiation. At variance with recent studies, a careful treatment of IGM transmission leads to no significant effects on the high-redshift broad-band color selection. The decreasing mean free path of ionizing photons with increasing redshift further diminishes the contribution of the LyC to broad-band colors. We also demonstrate that prominent LyC sources can be selected under suitable conditions by calculating the probability of a null escaping ionizing radiation. The method is applied to a sample of galaxies extracted from the GOODS-S field. A known LyC source at z=3.795 is successfully recovered as a LyC emitter candidate and another convincing candidate at z=3.212 is reported. A detailed analysis of the two sources (including their variability and morphology) suggests a possible mixture of stellar and non-stellar (AGN) contribution in the ultraviolet. Conclusions: Classical broad-band color selection of 2.5<z<4.5 galaxies does not prevent the inclusion of LyC emitters in the selected samples. Large fesc in relatively bright galaxies (L>0.1L*) could be favored by the presence of a faint AGN not easily detected at any wavelength. A hybrid stellar and non-stellar (AGN) ionizing emission could coexist in these systems and explain the tensions found among the UV excess and the stellar population synthesis models reported in literature.

Narrow-band imaging of the rest-frame Lyman continuum (LyC) of galaxies at z~3.1 has produced a large number of candidate LyC-emitting galaxies. These samples are contaminated by galaxies at lower redshift. To better understand LyC escape, we need an uncontaminated sample of galaxies that emit strongly in the LyC. Here we present deep Hubble imaging of five bright galaxies at z~3.1 that had previously been identified as candidate LyC-emitters with ground-based images. The WFC3 F336W images probe the LyC of galaxies at z>3.06 and provide an order-of-magnitude increase in spatial resolution over ground-based imaging. The non-ionizing UV images often show multiple galaxies (or components) within ~1'' of the candidate LyC emission seen from the ground. In each case, only one of the components is emitting light in the F336W filter, which would indicate LyC escape if that component is at z>3.06. We use Keck/NIRSPEC near-IR spectroscopy to measure redshifts of these components to distinguish LyC-emitters from foreground contamination. We find that two candidates are low redshift contaminants, one candidate had a previously misidentified redshift, and the other two cannot be confirmed as LyC-emitters. The level of contamination is consistent with previous estimates. For the galaxies with z>3.06, we derive strong 1 sigma limits on the relative escape fraction between 0.07 and 0.09. We still do not have a sample of definitive LyC-emitters, and a much larger study of low luminosity galaxies is required. The combination of high resolution imaging and deep spectroscopy is critical for distinguishing LyC-emitters from foreground contaminants.

We present the public release of the stellar mass catalogs for the GOODS-S and UDS fields obtained using some of the deepest near-IR images available, achieved as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) project. We combine the effort from ten different teams, who computed the stellar masses using the same photometry and the same redshifts. Each team adopted their preferred fitting code, assumptions, priors, and parameter grid. The combination of results using the same underlying stellar isochrones reduces the systematics associated with the fitting code and other choices. Thanks to the availability of different estimates, we can test the effect of some specific parameters and assumptions on the stellar mass estimate. The choice of the stellar isochrone library turns out to have the largest effect on the galaxy stellar mass estimates, resulting in the largest distributions around the median value (with a semi interquartile range larger than 0.1 dex). On the other hand, for most galaxies, the stellar mass estimates are relatively insensitive to the different parameterizations of the star formation history. The inclusion of nebular emission in the model spectra does not have a significant impact for the majority of galaxies (less than a factor of 2 for ~80% of the sample). Nevertheless, the stellar mass for the subsample of young galaxies (age < 100 Myr), especially in particular redshift ranges (e.g., 2.2 < z < 2.4, 3.2 < z < 3.6, and 5.5 < z < 6.5), can be seriously overestimated (by up to a factor of 10 for < 20 Myr sources) if nebular contribution is ignored.

Although giant clumps of stars are crucial to galaxy formation and evolution, the most basic demographics of clumps are still uncertain, mainly because the definition of clumps has not been thoroughly discussed. In this paper, we study the basic demographics of clumps in star-forming galaxies (SFGs) at 0.5<z<3, using our proposed physical definition that UV-bright clumps are discrete star-forming regions that individually contribute more than 8% of the rest-frame UV light of their galaxies. Clumps defined this way are significantly brighter than the HII regions of nearby large spiral galaxies, either individually or blended, when physical spatial resolution and cosmological dimming are considered. Under this definition, we measure the fraction of SFGs that contain at least one off-center clump (Fclumpy) and the contributions of clumps to the rest-frame UV light and star formation rate of SFGs in the CANDELS/GOODS-S and UDS fields, where our mass-complete sample consists of 3239 galaxies with axial ratio q>0.5. The redshift evolution of Fclumpy changes with the stellar mass (M*) of the galaxies. Low-mass (log(M*/Msun)<9.8) galaxies keep an almost constant Fclumpy of about 60% from z~3.0 to z~0.5. Intermediate-mass and massive galaxies drop their Fclumpy from 55% at z~3.0 to 40% and 15%, respectively, at z~0.5. We find that (1) the trend of disk stabilization predicted by violent disk instability matches the Fclumpy trend of massive galaxies; (2) minor mergers are a viable explanation of the Fclumpy trend of intermediate-mass galaxies at z<1.5, given a realistic observability timescale; and (3) major mergers are unlikely responsible for the Fclumpy trend in all masses at z<1.5. The clump contribution to the rest-frame UV light of SFGs shows a broad peak around galaxies with log(M*/Msun)~10.5 at all redshifts, possibly linked to the molecular gas fraction of the galaxies. (Abridged)

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $\sigma(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

We present an analysis of the deepest Herschel images in four major extragalactic fields GOODS-North, GOODS-South, UDS and COSMOS obtained within the GOODS-Herschel and CANDELS-Herschel key programs. The picture provided by 10497 individual far-infrared detections is supplemented by the stacking analysis of a mass-complete sample of 62361 star-forming galaxies from the CANDELS-HST H band-selected catalogs and from two deep ground-based Ks band-selected catalogs in the GOODS-North and the COSMOS-wide fields, in order to obtain one of the most accurate and unbiased understanding to date of the stellar mass growth over the cosmic history. We show, for the first time, that stacking also provides a powerful tool to determine the dispersion of a physical correlation and describe our method called "scatter stacking" that may be easily generalized to other experiments. We demonstrate that galaxies of all masses from z=4 to 0 follow a universal scaling law, the so-called main sequence of star-forming galaxies. We find a universal close-to-linear slope of the logSFR-logM* relation with evidence for a flattening of the main sequence at high masses (log(M*/Msun) > 10.5) that becomes less prominent with increasing redshift and almost vanishes by z~2. This flattening may be due to the parallel stellar growth of quiescent bulges in star-forming galaxies. Within the main sequence, we measure a non varying SFR dispersion of 0.3 dex. The specific SFR (sSFR=SFR/M*) of star-forming galaxies is found to continuously increase from z=0 to 4. Finally we discuss the implications of our findings on the cosmic SFR history and show that more than 2/3 of present-day stars must have formed in a regime dominated by the main sequence mode. As a consequence we conclude that, although omnipresent in the distant Universe, galaxy mergers had little impact in shaping the global star formation history over the last 12.5 Gyr.

The relation between the galaxy stellar mass M_star and the dark matter halo mass M_h gives important information on the efficiency in forming stars and assembling stellar mass in galaxies. We present the stellar mass to halo mass ratio (SMHR) measurements at redshifts 2<z<5, obtained from the VIMOS Ultra Deep Survey. We use halo occupation distribution (HOD) modelling of clustering measurements on ~3000 galaxies with spectroscopic redshifts to derive the dark matter halo mass M_h, and SED fitting over a large set of multi-wavelength data to derive the stellar mass M_star and compute the SMHR=M_star/M_h. We find that the SMHR ranges from 1% to 2.5% for galaxies with M_star=1.3x10^9 M_sun to M_star=7.4x10^9 M_sun in DM halos with M_h=1.3x10^{11} M_sun} to M_h=3x10^{11} M_sun. We derive the integrated star formation efficiency (ISFE) of these galaxies and find that the star formation efficiency is a moderate 6-9% for lower mass galaxies while it is relatively high at 16% for galaxies with the median stellar mass of the sample ~7x10^9 M_sun. The lower ISFE at lower masses may indicate that some efficient means of suppressing star formation is at work (like SNe feedback), while the high ISFE for the average galaxy at z~3 is indicating that these galaxies are efficiently building-up their stellar mass at a key epoch in the mass assembly process. We further infer that the average mass galaxy at z~3 will start experiencing star formation quenching within a few hundred millions years.

Quenched galaxies at z>2 are nearly all very compact relative to z~0, suggesting a physical connection between high stellar density and efficient, rapid cessation of star-formation. We present restframe UV spectra of Lyman-break galaxies (LBGs) at z~3 selected to be candidate progenitors of quenched galaxies at z~2 based on their compact restframe optical sizes and high surface density of star-formation. We compare their UV properties to those of more extended LBGs of similar mass and star formation rate (non-candidates). We find that candidate progenitors have faster ISM gas velocities and higher equivalent widths of interstellar absorption lines, implying larger velocity spread among absorbing clouds. Candidates deviate from the relationship between equivalent widths of Lyman-alpha and interstellar absorption lines in that their Lyman-alpha emission remains strong despite high interstellar absorption, possibly indicating that the neutral HI fraction is patchy such that Lyman-alpha photons can escape. We detect stronger CIV P-Cygni features (emission and absorption) and HeII emission in candidates, indicative of larger populations of metal rich Wolf-Rayet stars compared to non-candidates. The faster bulk motions, broader spread of gas velocity, and Lyman-alpha properties of candidates are consistent with their ISM being subject to more energetic feedback than non-candidates. Together with their larger metallicity (implying more evolved star-formation activity) this leads us to propose, if speculatively, that they are likely to quench sooner than non-candidates, supporting the validity of selection criteria used to identify them as progenitors of z~2 passive galaxies. We propose that massive, compact galaxies undergo more rapid growth of stellar mass content, perhaps because the gas accretion mechanisms are different, and quench sooner than normally-sized LBGs at these early epochs.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

Using new spectroscopic observations obtained as part of the VIMOS Ultra-Deep Survey (VUDS), we perform a systematic search for overdense environments in the early universe ($z>2$) and report here on the discovery of Cl J0227-0421, a massive protocluster at $z=3.29$. This protocluster is characterized by both the large overdensity of spectroscopically confirmed members, $\delta_{gal}=10.5\pm2.8$, and a significant overdensity in photometric redshift members. The halo mass of this protocluster is estimated, by a variety of methods, to be roughly $3\times10^{14}$ $\mathcal{M}_{\odot}$ at $z\sim3.3$, which, evolved to $z=0$ results in a halo mass rivaling or exceeding that of the Coma cluster. The properties of 19 spectroscopically confirmed member galaxies are compared with a large sample of VUDS/VVDS galaxies in lower density field environments at similar redshifts. We find tentative evidence for an excess of redder, brighter, and more massive galaxies within the confines of the protocluster relative to the field population, which suggests that we may be observing the beginning of environmentally-induced quenching. The properties of these galaxies are investigated, including a discussion of the brightest protocluster galaxy which appears to be undergoing vigorous coeval nuclear and starburst activity. The remaining member galaxies appear to have characteristics which are largely similar to the field population. Though we find weaker evidence of the suppression of the median star formation rates amongst and differences in stacked spectra of member galaxies with respect to the field, we defer any conclusions of these trends to future work with the ensemble of protostructures that are found in the full VUDS sample.

We investigate the evolution of galaxy clustering for galaxies in the redshift range 2.0<$z$<5.0 using the VIMOS Ultra Deep Survey (VUDS). We present the projected (real-space) two-point correlation function $w_p(r_p)$ measured by using 3022 galaxies with robust spectroscopic redshifts in two independent fields (COSMOS and VVDS-02h) covering in total 0.8 deg$^2$. We quantify how the scale dependent clustering amplitude $r_0$ changes with redshift making use of mock samples to evaluate and correct the survey selection function. Using a power-law model $\xi(r) = (r/r_0)^{-\gamma}$ we find that the correlation function for the general population is best fit by a model with a clustering length $r_0$=3.95$^{+0.48}_{-0.54}$ h$^{-1}$Mpc and slope $\gamma$=1.8$^{+0.02}_{-0.06}$ at $z$~2.5, $r_0$=4.35$\pm$0.60 h$^{-1}$Mpc and $\gamma$=1.6$^{+0.12}_{-0.13}$ at $z$~3.5. We use these clustering parameters to derive the large-scale linear galaxy bias $b_L^{PL}$, between galaxies and dark matter. We find $b_L^{PL}$ = 2.68$\pm$0.22 at redshift $z$~3 (assuming $\sigma_8$ = 0.8), significantly higher than found at intermediate and low redshifts. We fit an HOD model to the data and we obtain that the average halo mass at redshift $z$~3 is $M_h$=10$^{11.75\pm0.23}$ h$^{-1}$M$_{\odot}$. From this fit we confirm that the large-scale linear galaxy bias is relatively high at $b_L^{HOD}$ = 2.82$\pm$0.27. Comparing these measurements with similar measurements at lower redshifts we infer that the star-forming population of galaxies at $z$~3 should evolve into the massive and bright ($M_r$<-21.5) galaxy population which typically occupy haloes of mass $\langle M_h\rangle$ = 10$^{13.9}$ h$^{-1}$ $M_{\odot}$ at redshift $z$=0.

(arXiv abridged abstract) The observed UV rest-frame spectra of distant galaxies are the result of their intrinsic emission combined with absorption along the line of sight produced by the inter-galactic medium (IGM). Here we analyse the evolution of the mean IGM transmission Tr(Ly_alpha) and its dispersion along the line of sight for 2127 galaxies with 2.5<z<5.5 in the VIMOS Ultra Deep Survey (VUDS). We fit model spectra combined with a range of IGM transmission to the galaxy spectra using the spectral fitting algorithm GOSSIP+. We use these fits to derive the mean IGM transmission towards each galaxy for several redshift slices from z=2.5 to z=5.5. We find that the mean IGM transmission defined as Tr(Ly_alpha)=e^{-tau} (with tau the HI optical depth) is 79%, 69%, 59%, 55% and 46% at redshifts 2.75, 3,22, 3.70, 4.23, 4.77, respectively. We compare these results to measurements obtained from quasars lines of sight and find that the IGM transmission towards galaxies is in excellent agreement with quasar values up to redshift z~4. We find tentative evidence for a higher IGM transmission at z>= 4 compared to results from QSOs, but a degeneracy between dust extinction and IGM prevents to draw firm conclusions if the internal dust extinction for star-forming galaxies at z>4 takes a mean value significantly in excess of E(B-V)>0.15. Most importantly, we find a large dispersion of IGM transmission along the lines of sight towards distant galaxies with 68% of the distribution within 10 to 17% of the median value in delta z=0.5 bins, similar to what is found on the LOS towards QSOs. We demonstrate the importance of taking into account this large range of IGM transmission when selecting high redshift galaxies based on their colour properties (e.g. LBG or photometric redshift selection) or otherwise face a significant incompleteness in selecting high redshift galaxy populations.

We study the evolution of the star formation rate (SFR) - stellar mass (M_star) relation and specific star formation rate (sSFR) of star forming galaxies (SFGs) since a redshift z~5.5 using 2435 (4531) galaxies with highly reliable (reliable) spectroscopic redshifts in the VIMOS Ultra-Deep Survey (VUDS). It is the first time that these relations can be followed over such a large redshift range from a single homogeneously selected sample of galaxies with spectroscopic redshifts. The log(SFR) - log(M_star) relation for SFGs remains roughly linear all the way up to z=5 but the SFR steadily increases at fixed mass with increasing redshift. We find that for stellar masses M_star>3.2 x 10^9 M_sun the SFR increases by a factor ~13 between z=0.4 and z=2.3. We extend this relation up to z=5, finding an additional increase in SFR by a factor 1.7 from z=2.3 to z=4.8 for masses M_star > 10^10 M_sun. We observe a turn-off in the SFR-M_star relation at the highest mass end up to a redshift z~3.5. We interpret this turn-off as the signature of a strong on-going quenching mechanism and rapid mass growth. The sSFR increases strongly up to z~2 but it grows much less rapidly in 2<z<5. We find that the shape of the sSFR evolution is not well reproduced by cold gas accretion-driven models or the latest hydrodynamical models. Below z~2 these models have a flatter evolution (1+z)^{Phi} with Phi=2-2.25 compared to the data which evolves more rapidly with Phi=2.8+-0.2. Above z~2, the reverse is happening with the data evolving more slowly with Phi=1.2+-0.1. The observed sSFR evolution over a large redshift range 0<z<5 and our finding of a non linear main sequence at high mass both indicate that the evolution of SFR and M_star is not solely driven by gas accretion. The results presented in this paper emphasize the need to invoke a more complex mix of physical processes {abridge}

Aims. The aim of this work is to constrain the evolution of the fraction of Lya emitters among UV selected star forming galaxies at 2<z<6, and to measure the stellar escape fraction of Lya photons over the same redshift range. Methods. We exploit the ultradeep spectroscopic observations collected by the VIMOS Ultra Deep Survey (VUDS) to build an unique, complete and unbiased sample of 4000 spectroscopically confirmed star forming galaxies at 2<z<6. Our galaxy sample UV luminosities brighter than M* at 2<z<6, and luminosities down to one magnitude fainter than M* at 2<z<3.5. Results. We find that 80% of the star forming galaxies in our sample have EW0(Lya)<10A, and correspondingly fesc(Lya)<1%. By comparing these results with literature, we conclude that the bulk of the Lya luminosity at 2<z<6 comes from galaxies that are fainter in the UV than those we sample in this work. The strong Lya emitters constitute, at each redshift, the tail of the distribution of the galaxies with extreme EW0(Lya) and fesc(Lya) . This tail of large EW0 and fesc(Lya) becomes more important as the redshift increases, and causes the fraction of Lya with EW0> 25A to increase from 5% at z=2 to 30% at z=6, with the increase being relatively stronger beyond z=4. We observe no difference, for the narrow range of UV luminosities explored in this work, between the fraction of strong Lya emitters among galaxies fainter or brighter than M*, although the fraction for the FUV faint galaxies evolves faster, at 2<z<3.5, than for the bright ones. We do observe an anticorrelation between E(B-V) and fesc(Lya): generally galaxies with high fesc(Lya) have also small amounts of dust (and viceversa). However, when the dust content is low (E(B-V)<0.05) we observe a very broad range of fesc(Lya), ranging from 10^-3 to 1. This implies that the dust alone is not the only regulator of the amount of escaping Lya photons.

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

The observed deficit of strongly Lyman-alpha emitting galaxies at z>6.5 is attributed to either increasing neutral hydrogen in the intergalactic medium (IGM) and/or to the evolving galaxy properties. To investigate this, we have performed very deep near-IR spectroscopy of z>7 galaxies using MOSFIRE on the Keck-I Telescope. We measure the Lyman-alpha fraction at z~8 (combined photometric redshift peak at z=7.7) using two methods. First, we derived NLy{\alpha}/Ntot directly using extensive simulations to correct for incompleteness. Second, we used a Bayesian formalism (introduced by Treu et al. 2012) that compares the z>7 galaxy spectra to models of the Lyman-alpha equivalent width (WLy{\alpha}) distribution at z~6. We explored two simple evolutionary scenarios: smooth evolution where Lyman-alpha is attenuated in all galaxies by a constant factor (perhaps owing to processes from galaxy evolution or a slowly increasing IGM opacity), and patchy evolution where Lyman-alpha is blocked in some fraction of galaxies (perhaps due to the IGM being opaque along only some fraction of sightlines). The Bayesian formalism places stronger constraints compared with the direct method. Combining our data with that in the literature we find that at z~8 the Lyman-alpha fraction has dropped by a factor >3(84% confidence interval) using both the smooth and patchy scenarios compared to the z~6 values. Furthermore, we find a tentative evidence that the data favor the patchy scenario over smooth (with "positive" Bayesian evidence), extending trends observed at z~7 to higher redshift. If this decrease is a result of reionization as predicted by theory, then our data imply the volume averaged neutral hydrogen fraction in the IGM to be >0.3 suggesting that the reionization of the universe is in progress at z~8.

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

We use data taken as part of HST/WFC3 observations of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) to identify massive and evolved galaxies at 3<z<4.5. This is performed using the strength of the Balmer break feature at rest-frame 3648A, which is a diagnostic of the age of the stellar population in galaxies. Using WFC3 H-band selected catalog for the CANDELS GOODS-S field and deep multi-waveband photometry from optical (HST) to mid-infrared (Spitzer) wavelengths, we identify a population of old and evolved post-starburst galaxies based on the strength of their Balmer breaks (Balmer Break Galaxies- BBGs). The galaxies are also selected to be bright in rest-frame near-IR wavelengths and hence, massive. We identify a total of 16 BBGs. Fitting the spectral energy distribution (SED) of the BBGs show that the candidate galaxies have average estimated ages of ~800 Myr and average stellar masses of ~5x10^10 M_sun, consistent with being old and massive systems. Two of our BBG candidates are also identified by the criteria that is sensitive to star forming galaxies (LBG selection). We find a number density of ~3.2x10^-5 Mpc^-3 for the BBGs corresponding to a mass density of ~2.0x10^6 M_sun/Mpc^3 in the redshift range covering the survey. Given the old age and the passive evolution, it is argued that some of these objects formed the bulk of their mass only a few hundred million years after the Big Bang.

We determine the intrinsic, 3-dimensional shape distribution of star-forming galaxies at 0<z<2.5, as inferred from their observed projected axis ratios. In the present-day universe star-forming galaxies of all masses 1e9 - 1e11 Msol are predominantly thin, nearly oblate disks, in line with previous studies. We now extend this to higher redshifts, and find that among massive galaxies (M* > 1e10 Msol) disks are the most common geometric shape at all z < 2. Lower-mass galaxies at z>1 possess a broad range of geometric shapes: the fraction of elongated (prolate) galaxies increases toward higher redshifts and lower masses. Galaxies with stellar mass 1e9 Msol (1e10 Msol) are a mix of roughly equal numbers of elongated and disk galaxies at z~1 (z~2). This suggests that galaxies in this mass range do not yet have disks that are sustained over many orbital periods, implying that galaxies with present-day stellar mass comparable to that of the Milky Way typically first formed such sustained stellar disks at redshift z~1.5-2. Combined with constraints on the evolution of the star formation rate density and the distribution of star formation over galaxies with different masses, our findings imply that, averaged over cosmic time, the majority of stars formed in disks.

We measure new estimates for the galaxy stellar mass function and star formation rates for samples of galaxies at $z \sim 4,~5,~6~\&~7$ using data in the CANDELS GOODS South field. The deep near-infrared observations allow us to construct the stellar mass function at $z \geq 6$ directly for the first time. We estimate stellar masses for our sample by fitting the observed spectral energy distributions with synthetic stellar populations, including nebular line and continuum emission. The observed UV luminosity functions for the samples are consistent with previous observations, however we find that the observed $M_{UV}$ - M$_{*}$ relation has a shallow slope more consistent with a constant mass to light ratio and a normalisation which evolves with redshift. Our stellar mass functions have steep low-mass slopes ($\alpha \approx -1.9$), steeper than previously observed at these redshifts and closer to that of the UV luminosity function. Integrating our new mass functions, we find the observed stellar mass density evolves from $\log_{10} \rho_{*} = 6.64^{+0.58}_{-0.89}$ at $z \sim 7$ to $7.36\pm0.06$ $\text{M}_{\odot} \text{Mpc}^{-3}$ at $z \sim 4$. Finally, combining the measured UV continuum slopes ($\beta$) with their rest-frame UV luminosities, we calculate dust corrected star-formation rates (SFR) for our sample. We find the specific star-formation rate for a fixed stellar mass increases with redshift whilst the global SFR density falls rapidly over this period. Our new SFR density estimates are higher than previously observed at this redshift.

[Abridged] We characterise a massive proto-cluster at z=2.895 that we found in the COSMOS field using the spectroscopic sample of the VIMOS Ultra-Deep Survey (VUDS). This is one of the rare structures at z~3 not identified around AGNs or radio galaxies, so it is an ideal laboratory to study galaxy formation in dense environments. The structure comprises 12 galaxies with secure spectroscopic redshift in an area of 7'x8', in a z bin of Dz=0.016. The measured galaxy number overdensity is delta_g=12+/-2. This overdensity has total mass of M~8.1x10^(14)M_sun in a volume of 13x15x17 Mpc^3. Simulations indicate that such an overdensity at z~2.9 is a proto-cluster that will collapse in a cluster of total mass M~2.5x10^(15)M_sun at z=0. We compare the properties of the galaxies within the overdensity with a control sample at the same z but outside the overdensity. We did not find any statistically significant difference between the properties (stellar mass, SFR, sSFR, NUV-r, r-K) of the galaxies inside and outside the overdensity. The stacked spectrum of galaxies in the overdensity background shows a significant absorption feature at the wavelength of Lya redshifted at z=2.895 (lambda=4736 A), with a rest frame EW = 4+/- 1.4 A. Stacking only background galaxies without intervening sources at z~2.9 along their line of sight, we find that this absorption feature has a rest frame EW of 10.8+/-3.7 A, with a detection S/N of ~4. These EW values imply a high column density (N(HI)~3-20x10^(19)cm^(-2)), consistent with a scenario where such absorption is due to intervening cold gas streams, falling into the halo potential wells of the proto-cluster galaxies. However, we cannot exclude the hypothesis that this absorption is due to the diffuse gas within the overdensity.

We aim to determine the redshift of GDS-1408, the most solid z~7 galaxy candidate lying in the Hubble Ultra Deep Field. We have used all the VLT spectra of GDS-1408 collected by us and two other groups with FORS2 at VLT in the last five years, for a total integration time of 52hr. The combined spectrum is the deepest ever obtained of a galaxy in the Reionization epoch. We do not detect any emission line or continuum over the whole wavelength range, up to 10100A. Based on an accurate set of simulations, we are able to put a stringent upper limit of f(Lya)<3x10^(-18) erg/s/cm2 at 3-9 sigma in the explored wavelength range, corresponding to a rest-frame equivalent width EW<9A. Combining this limit with the SED modelling we refine the redshift to be z=6.82+/- 0.1 (1-sigma). The same SED fitting indicates that GDS-1408 is relatively extinct (A1600~1) with a dust corrected star formation rate of ~ 20 Msol/yr. The comparison between the un-attenuated equivalent width predicted by the case-B recombination theory and the observed upper limit, provides a limit on the effective Lya escape fraction of f_(esc)^(eff)(Lya)<8%. Even though we cannot rule out a major contribution of the inter/circum galactic medium in damping the line, a plausible interpretation is that GDS-1408 is moderately evolved and contains sufficient gas and dust to attenuate the Lya emission, before it reaches the intergalactic medium. The redshift confirmation of even the best z~7 candidates is very hard to achieve (unless the Lya or unusually strong rest-UV nebular emission lines are present) with the current generation of 8-10m class telescopes. We show that both JWST and E-ELT will be necessary to make decisive progresses. Currently, the increased redshift accuracy obtained with this kind of analysis makes ALMA an interesting option for the redshift confirmation.

We report on a study of the 2300-2600\AA FeII/FeII* multiplets in the rest-UV spectra of star-forming galaxies at 1.0<z<2.6 as probes of galactic-scale outflows. We extracted a mass-limited sample of 97 galaxies at z~1.0-2.6 from ultra-deep spectra obtained during the GMASS spetroscopic survey in the GOODS South field with the VLT and FORS2. We obtain robust measures of the rest equivalent width of the FeII absorption lines down to a limit of W_r>1.5 \AA and of the FeII* emission lines to W_r>0.5 \AA. Whenever we can measure the systemic redshift of the galaxies from the [OII] emission line, we find that both the FeII and MgII absorption lines are blueshifted, indicative that both species trace gaseous outflows. We also find, however, that the FeII gas has generally lower outflow velocity relative to that of MgII. We investigate the variation of FeII line profiles as a function of the radiative transfer properties of the lines, and find that transitions with higher oscillator strengths are more blueshifted in terms of both line centroids and line wings. We discuss the possibility that FeII lines are suppressed by stellar absorptions. The lower velocities of the FeII lines relative to the MgII doublet, as well as the absence of spatially extended FeII* emission in 2D stacked spectra, suggest that most clouds responsible for the FeII absorption lie close (3~4 kpc) to the disks of galaxies. We show that the FeII/FeII* multiplets offer unique probes of the kinematic structure of galactic outflows.

We investigate the star formation histories (SFHs) of high redshift (3 <~ z <~ 5) star-forming galaxies selected based on their rest-frame ultraviolet (UV) colors in the CANDELS/GOODS-S field. By comparing the results from the spectral-energy-distribution-fitting analysis with two different assumptions about the SFHs --- i.e., exponentially declining SFHs as well as increasing ones, we conclude that the SFHs of high-redshift star-forming galaxies increase with time rather than exponentially decline. We also examine the correlations between the star formation rates (SFRs) and the stellar masses. When the galaxies are fit with rising SFRs, we find that the trend seen in the data qualitatively matches the expectations from a semi-analytic model of galaxy formation. The mean specific SFR is shown to increase with redshift, also in agreement with the theoretical prediction. From the derived tight correlation between stellar masses and SFRs, we derive the mean SFH of star-forming galaxies in the redshift range of 3 <~ z <~ 5, which shows a steep power-law (with power alpha = 5.85) increase with time. We also investigate the formation timescales and the mean stellar population ages of these star-forming galaxies. Our analysis reveals that UV-selected star-forming galaxies have a broad range of the formation redshift. The derived stellar masses and the stellar population ages show positive correlation in a sense that more massive galaxies are on average older, but with significant scatter. This large scatter implies that the galaxies' mass is not the only factor which affects the growth or star formation of high-redshift galaxies.

We present the VIMOS Ultra Deep Survey (VUDS), a spectroscopic redshift survey of ~10.000 very faint galaxies to study the major phase of galaxy assembly 2<z<~6. The survey covers 1 deg^2 in 3 separate fields: COSMOS, ECDFS and VVDS-02h, with targets selection based on an inclusive combination of photometric redshifts and color properties. Spectra covering 3650<lambda<9350 A are obtained with VIMOS on the ESO-VLT with integration times of 14h. Here we present the survey strategy, the target selection, the data processing, as well as the redshift measurement process, emphasizing the specific methods adapted to this high redshift range. The spectra quality and redshift reliability are discussed, and we derive a completeness in redshift measurement of 91%, or 74% for the most reliable measurements, down to i_AB=25, and measurements are performed all the way down to i_AB=27. The redshift distribution of the main sample peaks at z=3-4 and extends over a large redshift range mainly in 2 < z < 6. At 3<z<5, the galaxies cover a large range of luminosities -23< M_U < -20.5, stellar mass 10^9 M_sun< M_star < 10^{11} M_sun, and star formation rates 1 M_sun/yr< SFR < 10^3 M_sun/yr. We discuss the spectral properties of galaxies using individual as well as stacked spectra. The comparison between spectroscopic and photometric redshifts as well as color selection demonstrate the effectiveness of our selection scheme. With ~6000 galaxies with reliable spectroscopic redshifts in 2<z<6 expected when complete, this survey is the largest at these redshifts and offers the opportunity for unprecedented studies of the star-forming galaxy population and its distribution in large scale structures during the major phase of galaxy assembly.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

We report the LBT/MODS1 spectroscopic confirmation of two images of faint Lyman alpha emitters at $z=6.4$ behind the Frontier Fields galaxy cluster MACSJ0717.5+3745. A wide range of lens models suggests that the two images are highly magnified, with a strong lower limit of mu>5. These are the faintest z>6 candidates spectroscopically confirmed to date. These may be also multiple images of the same z=6.4 source as supported by their similar intrinsic properties, but the lens models are inconclusive regarding this interpretation. To be cautious, we derive the physical properties of each image individually. Thanks to the high magnification, the observed near-infrared (restframe ultraviolet) part of the spectral energy distributions and Ly-alpha lines are well detected with S/N(m_1500)>~10 and S/N(Ly-alpha)~10-15. Adopting mu>5, the absolute magnitudes, M_1500, and Ly-alpha fluxes, are fainter than -18.7 and 2.8x10^(-18)erg/s/cm2, respectively. We find a very steep ultraviolet spectral slope beta=-3.0+/-0.5 (F_lambda=lambda^(beta)), implying that these are very young, dust-free and low metallicity objects, made of standard stellar populations or even extremely metal poor stars (age<~30Myr, E(B-V)=0 and metallicity 0.0-0.2 Z/Zsolar). The objects are compact (< 1 kpc^(2)), and with a stellar mass M* < 10^(8) M_solar. The very steep beta, the presence of the Ly-alpha line and the intrinsic FWHM (<300 km/s) of these newborn objects do not exclude a possible leakage of ionizing radiation. We discuss the possibility that such faint galaxies may resemble those responsible for cosmic reionization.

We analyze the star-forming and structural properties of 45 massive (log(M/Msun)>10) compact star-forming galaxies (SFGs) at 2<z<3 to explore whether they are progenitors of compact quiescent galaxies at z~2. The optical/NIR and far-IR Spitzer/Herschel colors indicate that most compact SFGs are heavily obscured. Nearly half (47%) host an X-ray bright AGN. In contrast, only about 10% of other massive galaxies at that time host AGNs. Compact SFGs have centrally-concentrated light profiles and spheroidal morphologies similar to quiescent galaxies, and are thus strikingly different from other SFGs. Most compact SFGs lie either within the SFR-M main sequence (65%) or below (30%), on the expected evolutionary path towards quiescent galaxies. These results show conclusively that galaxies become more compact before they lose their gas and dust, quenching star formation. Using extensive HST photometry from CANDELS and grism spectroscopy from the 3D-HST survey, we model their stellar populations with either exponentially declining (tau) star formation histories (SFHs) or physically-motivated SFHs drawn from semi-analytic models (SAMs). SAMs predict longer formation timescales and older ages ~2 Gyr, which are nearly twice as old as the estimates of the tau models. While both models yield good SED fits, SAM SFHs better match the observed slope and zero point of the SFR-M main sequence. Some low-mass compact SFGs (log(M/Msun)=10-10.6) have younger ages but lower sSFRs than that of more massive galaxies, suggesting that the low-mass galaxies reach the red sequence faster. If the progenitors of compact SFGs are extended SFGs, state-of-the-art SAMs show that mergers and disk instabilities are both able to shrink galaxies, but disk instabilities are more frequent (60% versus 40%) and form more concentrated galaxies. We confirm this result via high-resolution hydrodynamic simulations.

Out of several dozen z > 7 candidate galaxies observed spectroscopically, only five have been confirmed via Lyman-alpha emission, at z=7.008, 7.045, 7.109, 7.213 and 7.215. The small fraction of confirmed galaxies may indicate that the neutral fraction in the intergalactic medium (IGM) rises quickly at z > 6.5, as Lyman-alpha is resonantly scattered by neutral gas. However, the small samples and limited depth of previous observations makes these conclusions tentative. Here we report the results of a deep near-infrared spectroscopic survey of 43 z > 6.5 galaxies. We detect only a single galaxy, confirming that some process is making Lyman-alpha difficult to detect. The detected emission line at 1.0343 um is likely to be Lyman-alpha emission, placing this galaxy at a redshift z = 7.51, an epoch 700 million years after the Big Bang. This galaxy's colors are consistent with significant metal content, implying that galaxies become enriched rapidly. We measure a surprisingly high star formation rate of 330 Msol/yr, more than a factor of 100 greater than seen in the Milky Way. Such a galaxy is unexpected in a survey of our size, suggesting that the early universe may harbor more intense sites of star-formation than expected.

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

We use GOODS and CANDELS images to identify progenitors of massive (log M > 10 Msun) compact "early-type" galaxies (ETGs) at z~1.6. Since merging and accretion increase the size of the stellar component of galaxies, if the progenitors are among known star-forming galaxies, these must be compact themselves. We select candidate progenitors among compact Lyman-break galaxies at z~3 based on their mass, SFR and central stellar density and find that these account for a large fraction of, and possibly all, compact ETGs at z~1.6. We find that the average far-UV SED of the candidates is redder than that of the non-candidates, but the optical and mid-IR SED are the same, implying that the redder UV of the candidates is inconsistent with larger dust obscuration, and consistent with more evolved (aging) star-formation. This is in line with other evidence that compactness is a sensitive predictor of passivity among high-redshift massive galaxies. We also find that the light distribution of both the compact ETGs and their candidate progenitors does not show any extended "halos" surrounding the compact "core", both in individual images and in stacks. We argue that this is generally inconsistent with the morphology of merger remnants, even if gas-rich, as predicted by N-body simulations. This suggests that the compact ETGs formed via highly dissipative, mostly gaseous accretion of units whose stellar components are very small and undetected in the HST images, with their stellar mass assembling in-situ, and that they have not experienced any major merging until the epoch of observations at z~1.6.

We present results from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) photometric redshift methods investigation. In this investigation, the results from eleven participants, each using a different combination of photometric redshift code, template spectral energy distributions (SEDs) and priors, are used to examine the properties of photometric redshifts applied to deep fields with broad-band multi-wavelength coverage. The photometry used includes U-band through mid-infrared filters and was derived using the TFIT method. Comparing the results, we find that there is no particular code or set of template SEDs that results in significantly better photometric redshifts compared to others. However, we find codes producing the lowest scatter and outlier fraction utilize a training sample to optimize photometric redshifts by adding zero-point offsets, template adjusting or adding extra smoothing errors. These results therefore stress the importance of the training procedure. We find a strong dependence of the photometric redshift accuracy on the signal-to-noise ratio of the photometry. On the other hand, we find a weak dependence of the photometric redshift scatter with redshift and galaxy color. We find that most photometric redshift codes quote redshift errors (e.g., 68% confidence intervals) that are too small compared to that expected from the spectroscopic control sample. We find that all codes show a statistically significant bias in the photometric redshifts. However, the bias is in all cases smaller than the scatter, the latter therefore dominates the errors. Finally, we find that combining results from multiple codes significantly decreases the photometric redshift scatter and outlier fraction. We discuss different ways of combining data to produce accurate photometric redshifts and error estimates.

We present a UV-to-mid infrared multi-wavelength catalog in the CANDELS/GOODS-S field, combining the newly obtained CANDELS HST/WFC3 F105W, F125W, and F160W data with existing public data. The catalog is based on source detection in the WFC3 F160W band. The F160W mosaic includes the data from CANDELS deep and wide observations as well as previous ERS and HUDF09 programs. The mosaic reaches a 5$\sigma$ limiting depth (within an aperture of radius 0.17 arcsec) of 27.4, 28.2, and 29.7 AB for CANDELS wide, deep, and HUDF regions, respectively. The catalog contains 34930 sources with the representative 50% completeness reaching 25.9, 26.6, and 28.1 AB in the F160W band for the three regions. In addition to WFC3 bands, the catalog also includes data from UV (U-band from both CTIO/MOSAIC and VLT/VIMOS), optical (HST/ACS F435W, F606W, F775W, F814W, and F850LP), and infrared (HST/WFC3 F098M, VLT/ISAAC Ks, VLT/HAWK-I Ks, and Spitzer/IRAC 3.6, 4.5, 5.8, 8.0 $\mu$m) observations. The catalog is validated via stellar colors, comparison with other published catalogs, zeropoint offsets determined from the best-fit templates of the spectral energy distribution of spectroscopically observed objects, and the accuracy of photometric redshifts. The catalog is able to detect unreddened star-forming (passive) galaxies with stellar mass of 10^{10}M_\odot at a 50% completeness level to z$\sim$3.4 (2.8), 4.6 (3.2), and 7.0 (4.2) in the three regions. As an example of application, the catalog is used to select both star-forming and passive galaxies at z$\sim$2--4 via the Balmer break. It is also used to study the color--magnitude diagram of galaxies at 0<z<4.

We discuss the state of the assembly of the Hubble Sequence in the mix of bright galaxies at redshift 1.4< z \le 2.5 with a large sample of 1,671 galaxies down to H_{AB}~26, selected from the HST/ACS and WFC3 images of the GOODS--South field obtained as part of the GOODS and CANDELS observations. We investigate the relationship between the star formation properties and morphology using various parametric diagnostics, such as the Sersic light profile, Gini (G), M_{20}, Concentration (C), Asymmetry (A) and multiplicity parameters. Our sample clearly separates into massive, red and passive galaxies versus less massive, blue and star forming ones, and this dichotomy correlates very well with the galaxies' morphological properties. Star--forming galaxies show a broad variety of morphological features, including clumpy structures and bulges mixed with faint low surface brightness features, generally characterized by disky-type light profiles. Passively evolving galaxies, on the other hand, very often have compact light distribution and morphology typical of today's spheroidal systems. We also find that artificially redshifted local galaxies have a similar distribution with z~2galaxies in a G-M_{20} plane. Visual inspection between the rest-frame optical and UV images show that there is a generally weak morphological k-correction for galaxies at z~2, but the comparison with non-parametric measures show that galaxies in the rest-frame UV are somewhat clumpier than rest-frame optical. Similar general trends are observed in the local universe among massive galaxies, suggesting that the backbone of the Hubble sequence was already in place at z~2.

We present the multiwavelength - ultraviolet to mid-infrared - catalog of the UKIDSS Ultra-Deep Survey (UDS) field observed as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Based on publicly available data, the catalog includes: the CANDELS data from the Hubble Space Telescope (near-infrared WFC3 F125W and F160W data and visible ACS F606W and F814W data), u-band data from CFHT/Megacam, B, V, Rc, i' and z' band data from Subaru/Suprime-Cam, Y and Ks band data from VLT/HAWK-I, J, H and K bands data from UKIDSS (Data Release 8), and Spitzer/IRAC data (3.6, 4.5 from SEDS, 5.8 and 8.0um from SpUDS). The present catalog is F160W-selected and contains 35932 sources over an area of 201.7 square arcmin and includes radio and X-ray detected sources and spectroscopic redshifts available for 210 sources.

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

The mass assembly of galaxies can proceed through different physical processes. Here we report on the spectroscopic identification of close physical pairs of galaxies at redshifts 2 <z < 4 and discuss the impact of major mergers in shaping galaxies at these early cosmological times. We aim to identify and characterize close physical pairs of galaxies destined to merge and use their properties to infer the contribution of merging processes to the early mass assembly of galaxies. We search for galaxy pairs with a transverse separation r_p < 25 h^{-1} kpc and a velocity difference Delta v < 500 km s^{-1} using the VIMOS VLT Deep Survey (VVDS) and early data from the VIMOS Ultra Deep Survey (VUDS) that comprise a sample of 1111 galaxies with spectroscopic redshifts measurements at redshifts 1.8 < z < 4 in the COSMOS, ECDFS, and VVDS-02h fields. We analyse their spectra and associated visible and near-infrared photometry to assess the main properties of merging galaxies that have an average stellar mass M_{star}=2.3 x 10^{10} M_{sun} at these redshifts. Using the 12 physical pairs found in our sample we obtain a first estimate of the merger fraction at these rdshifts, fmerg ~(15-20)%. The pair separations indicate that these pairs will merge within 1 Gyr, on average, with each producing a more massive galaxy by the time the cosmic star formation peaks at z~1-2. From the average mass ratio between galaxies in the pairs, the stellar mass of the resulting galaxy after merging will be ~60% larger than the most massive galaxy in the pair before merging. We conclude that major merging of galaxy pairs is on-going at 2<z<4 and significantly contributing to the major assembly phase of galaxies at this early epoch.

We study the evolution of the number density, as a function of the size, of passive early-type galaxies with a wide range of stellar masses 10^10<M*/Msun<10^11.5) from z~3 to z~1, exploiting the unique dataset available in the GOODS-South field, including the recently obtained WFC3 images as a part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). In particular, we select a sample of 107 massive (M*>10^10 M_sun), passive (SSFR<10^-2 Gyr^-1) and morphologically spheroidal galaxies at 1.2<z<3, taking advantage of the panchromatic dataset available for GOODS, including VLT, CFHT, Spitzer, Chandra and HST ACS+WFC3 data. We find that at 1<z<3 the passively evolving early-type galaxies are the reddest and most massive objects in the Universe, and we prove that a correlation between mass, morphology, color and star-formation activity is already in place at that epoch. We measure a significant evolution in the mass-size relation of passive early-type galaxies (ETGs) from z~3 to z~1, with galaxies growing on average by a factor of 2 in size in a 3 Gyr timescale only. We witness also an increase in the number density of passive ETGs of 50 times over the same time interval. We find that the first ETGs to form at z>2 are all compact or ultra-compact, while normal sized ETGs (meaning ETGs with sizes comparable to those of local counterparts of the same mass) are the most common ETGs only at z<1. The increase of the average size of ETGs at 0<z<1 is primarily driven by the appearance of new large ETGs rather than by the size increase of individual galaxies.

We use multi-band spatially resolved photometry from the Cosmic Assembly Near-IR Deep Legacy Survey (CANDELS) in the 4 MSec Chandra Deep Field-South (CDF-S) to explore the nuclear and extended colors, color gradients and stellar populations of X-ray selected AGN host galaxies out to z=3. Based on a study of their central light, we develop X-ray based criteria to exclude objects with strong AGN contamination. We use stellar masses from the FIREWORKS database to understand and account for stellar mass selection effects, and carefully study, for the first time, the resolved host galaxy properties of AGNs at z~2 in their rest-frame optical light without substantial nuclear contamination. AGN hosts span a sizable range of stellar masses, colors and color gradients at these redshifts. Their colors, color gradients and stellar population properties are very similar to inactive galaxies of the same stellar mass. At z~1, we find a slightly narrower range in host colors compared to inactive galaxies, as well as hints of more recent star-formation. These differences are weaker or non-existent among AGN hosts at z~2. We discuss the importance of AGN driven feedback in the quenching of galaxies at z>~1 and speculate on possible evolution in the relationship between black hole accretion and the host galaxy towards high redshifts.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

The exploration of the relation between galaxy sizes and other physical parameters has provided important clues for understanding galaxy formation. We use the CANDELS Deep+Wide surveys in the GOODS-South, UDS and EGS fields, complemented by data from the HUDF09 program, to address the relation between size and luminosity at z\sim7. We select 153 z-band drop-out galaxies in six different fields characterized by a wide combination of depth and areal coverage, ideally suited to sample without biases the observed size-magnitude plane. Detailed simulations allow us to derive the completeness as a function of size and magnitude and to quantify measurements errors/biases. We find that the half light radius distribution function of z\sim7 galaxies fainter than J=26.6 is peaked at <0.1 arcsec (or equivalently 0.5 kpc proper), while at brighter magnitudes high-z galaxies are typically larger than ~0.15 arcsec. We also find a well defined size-luminosity relation, Rh\simL^1/2. We compute the Luminosity Function in the HUDF and P12HUDF fields, finding large spatial variation on the number density of faint galaxies. Adopting the size distribution and the size-luminosity relation found for faint galaxies at z=7, we derive a mean slope of -1.7\pm0.1 for the luminosity function of LBGs at this redshift. Using this LF, we find that the amount of ionizing photons cannot keep the Universe re-ionized if the IGM is clumpy (C_HII>3) and the Lyman continuum escape fraction of high-z LBGs is relatively low (f_esc<0.3). With future CANDELS data, we can put severe limits to the role of galaxies in the reionization of the Universe.

This paper studies the properties of kiloparsec-scale clumps in star-forming galaxies at z~2 through multi-wavelength broad band photometry. A sample of 40 clumps is identified through auto-detection and visual inspection from 10 galaxies with 1.5<z<2.5 in the Hubble Ultra Deep Field, where deep and high-resolution HST/WFC3 and ACS images enable us to resolve structures of z~2 galaxies down to kpc scale in the rest-frame UV and optical bands as well as to detect clumps toward the faint end. The physical properties of clumps are measured through fitting spatially resolved seven-band (BVizYJH) spectral energy distribution to models. On average, the clumps are blue and have similar median rest-frame UV--optical color as the diffuse components of their host galaxies, but the clumps have large scatter in their colors. Although the star formation rate (SFR)--stellar mass relation of galaxies is dominated by the diffuse components, clumps emerge as regions with enhanced specific SFRs, contributing individually ~10% and together ~50% of the SFR of the host galaxies. However, the contributions of clumps to the rest-frame UV/optical luminosity and stellar mass are smaller, typically a few percent individually and ~20% together. On average, clumps are younger by 0.2 dex and denser by a factor of 8 than diffuse components. Clump properties have obvious radial variations in the sense that central clumps are redder, older, more extincted, denser, and less active on forming stars than outskirts clumps. Our results are broadly consistent with a widely held view that clumps are formed through gravitational instability in gas-rich turbulent disks and would eventually migrate toward galactic centers and coalesce into bulges. Roughly 40% of the galaxies in our sample contain a massive clump that could be identified as a proto-bulge, which seems qualitatively consistent with such a bulge-formation scenario.

We study the evolution of galaxy rest-frame ultraviolet (UV) colors in the epoch 4 < z < 8. We use new wide-field near-infrared data in GOODS-S from the CANDELS, HUDF09 and ERS programs to select galaxies via photometric redshift measurements. Our sample consists of 2812 candidate galaxies at z > 3.5, including 113 at z = 7 to 8. We fit the observed spectral energy distribution to a suite of synthetic stellar population models, and measure the value of the UV spectral slope (beta) from the best-fit model spectrum. The median value of beta evolves significantly from -1.82 (+0.00,-0.04) at z = 4, to -2.37 (+0.26,-0.06) at z = 7. Additionally, we find that faint galaxies at z = 7 have beta = -2.68 (+0.39,-0.24) (~ -2.4 after correcting for observational bias); this is redder than previous claims in the literature, and does not require "exotic" stellar populations to explain their colors. This evolution can be explained by an increase in dust extinction, with the timescale consistent with low-mass AGB stars forming the bulk of the dust. We find no significant (< 2-sigma) correlation between beta and M_UV when measuring M_UV at a consistent rest-frame wavelength of 1500 A. This is particularly true at bright magnitudes, though our results do show evidence for a weak correlation at faint magnitudes when galaxies in the HUDF are considered separately, hinting that dynamic range in sample luminosities may play a role. We do find a strong correlation between beta and the stellar mass at all redshifts, in that more massive galaxies exhibit redder colors. The most massive galaxies in our sample have red colors at each redshift, implying that dust can build up quickly in massive galaxies, and that feedback is likely removing dust from low-mass galaxies at z > 7. Thus the stellar-mass - metallicity relation, previously observed up to z ~ 3, may extend out to z = 7 - 8.

We combine high-resolution HST/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (log(M*) > 10) galaxies at redshifts z = 1.4 - 3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5 - 3. At z > 2 most cSFGs have specific star-formation rates (sSFR = 10^-9 yr^-1) half that of typical, massive SFGs at the same epoch, and host X-ray luminous AGN 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2 - 3 and fade to cQGs by z = 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary scenarios of QG formation: an early (z > 2), fast-formation path of rapidly-quenched cSFGs that evolve into cQGs that later enlarge within the quiescent phase, and a slow, late-arrival (z < 2) path for SFGs to form QGs without passing through a compact state.

[Abridged] We present the results of new near-IR spectroscopic observations of passive galaxies at z>1.4 in a concentration of BzK-selected galaxies in the COSMOS field. The observations have been conducted with Subaru/MOIRCS, and have resulted in absorption lines and/or continuum detection for 18 out of 34 objects. This allows us to measure spectroscopic redshifts for a sample almost complete to K(AB)=21. COSMOS photometric redshifts are found in fair agreement overall with the spectroscopic redshifts, with a standard deviation of ~0.05; however, ~30% of objects have photometric redshifts systematically underestimated by up to ~25%. We show that these systematic offsets in photometric redshifts can be removed by using these objects as a training set. All galaxies fall in four distinct redshift spikes at z=1.43, 1.53, 1.67 and 1.82, with this latter one including 7 galaxies. SED fits to broad-band fluxes indicate stellar masses in the range of ~4-40x10^10Msun and that star formation was quenched ~1 Gyr before the cosmic epoch at which they are observed. The spectra of several individual galaxies have allowed us to measure their Hdelta_F and Dn4000 indices, which confirms their identification as passive galaxies, as does a composite spectrum resulting from the coaddition of 17 individual spectra. The effective radii of the galaxies have been measured on the HST/ACS F814W image, confirming the coexistence at these redshifts of passive galaxies which are substantially more compact than their local counterparts with others that follow the local size-stellar mass relation. For the galaxy with best S/N spectrum we were able to measure a velocity dispersion of 270+/-105 km/s, indicating that this galaxy lies closely on the virial relation given its stellar mass and effective radius.

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 6<z<8. These galaxies were selected from new deep near-infrared HST imaging from the CANDELS, HUDF09 and ERS programs. In contrast to the majority of previous analyses, which assume that the distribution of galaxy ultraviolet (UV) luminosities follows a Schechter distribution, and that the distribution continues to luminosities far below our observable limit, we investigate the contribution to reionization from galaxies which we can observe, free from these assumptions. We find that the observable population of galaxies can sustain a fully reionized IGM at z=6, if the average ionizing photon escape fraction (f_esc) is ~30%. A number of previous studies have measured UV luminosity densities at these redshifts that vary by 5X, with many concluding that galaxies could not complete reionization by z=6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high values of f_esc were present. The observed UV luminosity density from our observed galaxy samples at z=7-8 is not sufficient to maintain a fully reionized IGM unless f_esc>50%. Combining our observations with constraints on the emission rate of ionizing photons from Ly-alpha forest observations at z=6, we can constrain f_esc<34% (2-sigma) if the observed galaxies are the only contributors to reionization, or <13% (2-sigma) if the luminosity function extends to M_UV = -13. These escape fractions are sufficient to complete reionization by z=6. These constraints imply that the volume ionized fraction of the IGM becomes less than unity at z>7, consistent with a number of complementary reionization probes. If faint galaxies dominate reionization, future JWST observations will probe deep enough to see them, providing an indirect constraint on the ionizing photon escape fraction [abridged].

We use HST/WFC3 imaging from the CANDELS Multicycle Treasury Survey, in conjunction with the Sloan Digital Sky Survey, to explore the evolution of galactic structure for galaxies with stellar masses >3e10M_sun from z=2.2 to the present epoch, a time span of 10Gyr. We explore the relationship between rest-frame optical color, stellar mass, star formation activity and galaxy structure. We confirm the dramatic increase from z=2.2 to the present day in the number density of non-star-forming galaxies above 3e10M_sun reported by others. We further find that the vast majority of these quiescent systems have concentrated light profiles, as parametrized by the Sersic index, and the population of concentrated galaxies grows similarly rapidly. We examine the joint distribution of star formation activity, Sersic index, stellar mass, inferred velocity dispersion, and stellar surface density. Quiescence correlates poorly with stellar mass at all z<2.2. Quiescence correlates well with Sersic index at all redshifts. Quiescence correlates well with `velocity dispersion' and stellar surface density at z>1.3, and somewhat less well at lower redshifts. Yet, there is significant scatter between quiescence and galaxy structure: while the vast majority of quiescent galaxies have prominent bulges, many of them have significant disks, and a number of bulge-dominated galaxies have significant star formation. Noting the rarity of quiescent galaxies without prominent bulges, we argue that a prominent bulge (and perhaps, by association, a supermassive black hole) is an important condition for quenching star formation on galactic scales over the last 10Gyr, in qualitative agreement with the AGN feedback paradigm.

We perform a detailed analysis of the resolved colors and stellar populations of a complete sample of 323 star-forming galaxies at 0.5 < z < 1.5, and 326 star-forming galaxies at 1.5 < z < 2.5 in the ERS and CANDELS-Deep region of GOODS-South. Galaxies were selected to be more massive than 10^10 Msun and have specific star formation rates above 1/t_H. We model the 7-band optical ACS + near-IR WFC3 spectral energy distributions of individual bins of pixels, accounting simultaneously for the galaxy-integrated photometric constraints available over a longer wavelength range. We analyze variations in rest-frame color, stellar surface mass density, age, and extinction as a function of galactocentric radius and local surface brightness/density, and measure structural parameters on luminosity and stellar mass maps. We find evidence for redder colors, older stellar ages, and increased dust extinction in the nuclei of galaxies. Big star-forming clumps seen in star formation tracers are less prominent or even invisible on the inferred stellar mass distributions. Off-center clumps contribute up to ~20% to the integrated SFR, but only 7% or less to the integrated mass of all massive star-forming galaxies at z ~ 1 and z ~ 2, with the fractional contributions being a decreasing function of wavelength used to select the clumps. The stellar mass profiles tend to have smaller sizes and M20 coefficients, and higher concentration and Gini coefficients than the light distribution. Our results are consistent with an inside-out disk growth scenario with brief (100 - 200 Myr) episodic local enhancements in star formation superposed on the underlying disk. Alternatively, the young ages of off-center clumps may signal inward clump migration, provided this happens efficiently on the order of an orbital timescale.

We investigate how star-forming galaxies typically assemble their masses at high redshift. Using the deep multi-wavelength coverage of the GOODS dataset, we measure stellar mass of a large sample of star-forming galaxies at z~4 and 5, and make a robust determination of stellar mass function (SMF). We report a broad correlation between stellar mass and UV luminosity, such that more UV-luminous galaxies are, on average, more massive. However, the correlation has a substantial intrinsic scatter evidenced by a non-negligible number of UV-faint but massive galaxies. Furthermore, the low-mass end of the SMF does not rise as steeply as the UV luminosity function (alpha_UVLF} -(1.7-1.8), alpha_SMF -(1.3-1.4)) of the same galaxies. In a smooth formation scenario where star formation (SF) is sustained at the observed rates for a long time, these galaxies would have accumulated more mass (by a factor of ~3) than observed and therefore the SMF would mirror more closely that of the UVLF. The relatively shallow slope of the SMF is due to the fact that many of the UV-selected galaxies are not massive enough, and therefore are too faint in their rest-frame optical bands, to be detected in the current observations. Our results favor a episodic formation history in which SF in low-mass galaxies vary significantly over time, a scenario favored by galaxy clustering. Our findings for the UV-faint galaxies are in contrast with those found for more UV-luminous galaxies, which exhibit tighter SFR-M_star correlations. The discrepancy may suggest that galaxies at different luminosities may have different evolutionary paths.

A new set of color selection criteria (VJL) analogous with the BzK method is designed to select both star-forming galaxies (SFGs) and passively-evolving galaxies (PEGs) at 2.3<z<3.5 by using rest-frame UV--optical (V-J vs. J-L) colors. The criteria are thoroughly tested with theoretical stellar population synthesis models and real galaxies with spectroscopic redshifts to evaluate their efficiency and contamination. We apply the well-tested VJL criteria to the HST/WFC3 Early Release Science field and study the physical properties of selected galaxies. The redshift distribution of selected SFGs peaks at z~2.7, slightly lower than that of Lyman Break Galaxies at z~3. Comparing the observed mid-infrared fluxes of selected galaxies with the prediction of pure stellar emission, we find that our VJL method is effective at selecting massive dusty SFGs that are missed by the Lyman Break Technique. About half of the star formation in massive (M_{star}>10^{10}M_{Sun}) galaxies at 2.3<z<3.5 is contributed by dusty (extinction E(B-V)>0.4) SFGs, which however, only account for ~20% of the number density of massive SFGs. We also use the mid-infrared fluxes to clean our PEG sample, and find that galaxy size can be used as a secondary criterion to effectively eliminate the contamination of dusty SFGs. The redshift distribution of the cleaned PEG sample peaks at z~2.5. We find 6 PEG candidates at z>3 and discuss possible methods to distinguish them from dusty contamination. We conclude that at least part of our candidates are real PEGs at z~3, implying that this type of galaxies began to form their stars at z>5. We measure the integrated stellar mass density of PEGs at z~2.5 and set constraints on it at z>3. We find that the integrated stellar mass density grows by at least about factor of 10 in 1 Gyr at 3<z<5 and by another factor of 10 in next 3.5 Gyr (1<z<3).

We have analysed a sample of 25 extremely red H-[4.5]>4 galaxies, selected using 4.5 micron data from the Spitzer SEDS survey and deep H-band data from the Hubble Space Telescope (HST) CANDELS survey, over ~180 square arcmin of the UKIDSS Ultra Deep Survey (UDS) field. Our aim is to investigate the nature of this rare population of mid-infrared (mid-IR) sources that display such extreme near-to-mid-IR colours. Using up to 17-band photometry (U through 8.0 microns), we have studied in detail their spectral energy distributions, including possible degeneracies in the photometric redshift/internal extinction (zphot-Av) plane. Our sample appears to include sources of very different nature. Between 45% and 75% of them are dust-obscured, massive galaxies at 3<zphot<5. All of the 24 micron-detected sources in our sample are in this category. Two of these have S(24 micron)>300 microJy, which at 3<zphot<5 suggests that they probably host a dust-obscured active galactic nucleus (AGN). Our sample also contains four highly obscured (Av>5) sources at zphot<1. Finally, we analyse in detail two zphot~6 galaxy candidates, and discuss their plausibility and implications. Overall, our red galaxy sample contains the tip of the iceberg of a larger population of z>3 galaxies to be discovered with the future James Webb Space Telescope.

We use the spatially-resolved, multi-band photometry in the GOODS South field acquired by the CANDELS project to constrain the nature of candidate Lyman continuum (LyC) emitters at redshift z~3.7 identified using ultra-deep imaging below the Lyman limit (1-sigma limit of ~30 AB in a 2" diameter aperture). In 18 candidates, out of a sample of 19 with flux detected at >3-sigma level, the light centroid of the candidate LyC emission is offset from that of the LBG by up to 1.5". We fit the SED of the LyC candidates to spectral population synthesis models to measure photometric redshifts and the stellar population parameters. We also discuss the differences in the UV colors between the LBG and the LyC candidates, and how to estimate the escape fraction of ionizing radiation (f_esc) in cases, like in most of our galaxies, where the LyC emission is spatially offset from the host galaxy. In all but one case we conclude that the candidate LyC emission is most likely due to lower redshift interlopers. Based on these findings, we argue that the majority of similar measurements reported in the literature need further investigation before it can be firmly concluded that LyC emission is detected. Our only surviving LyC candidate is a LBG at z=3.795, which shows the bluest (B-V) color among LBGs at similar redshift, a stellar mass of M~2 x 10^9 Msun, weak interstellar absorption lines and a flat UV spectral slope with no Lya in emission. We estimate its f_esc to be in the range 25%-100%, depending on the dust and intergalactic attenuation.

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\sim1.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 \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 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.

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_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5\sigma\ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5\sigma\ point-source limit of H \gtrsim 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 report the results of the counterpart identification and a detailed analysis of the physical properties of the 48 sources discovered in our deep 1.1mm wavelength imaging survey of the GOODS-South field using the AzTEC instrument on the Atacama Submillimeter Telescope Experiment (ASTE). One or more robust or tentative counterpart candidate is found for 27 and 14 AzTEC sources, respectively, by employing deep radio continuum, Spitzer MIPS & IRAC, and LABOCA 870 micron data. Five of the sources (10%) have two robust counterparts each, supporting the idea that these galaxies are strongly clustered and/or heavily confused. Photometric redshifts and star formation rates (SFRs) are derived by analyzing UV-to-optical and IR-to-radio SEDs. The median redshift of z~2.6 is similar to other earlier estimates, but we show that 80% of the AzTEC-GOODS sources are at z>2, with a significant high redshift tail (20% at z>3.3). Rest-frame UV and optical properties of AzTEC sources are extremely diverse, spanning 10 magnitude in the i- and K-band photometry with median values of i=25.3 and K=22.6 and a broad range of red colour (i-K=0-6). These AzTEC sources are some of the most luminous galaxies in the rest-frame optical bands at z>2, with inferred stellar masses of (1-30) x 10^{10} solar masses and UV-derived star formation rates of SFR(UV) > 10-1000 solar masses per year. The IR-derived SFR, 200-2000 solar masses per year, is independent of redshift or stellar mass. The resulting specific star formation rates, SSFR = 1-100 per Gyr, are 10-100 times higher than similar mass galaxies at z=0, and they extend the previously observed rapid rise in the SSFR with redshift to z=2-5. These galaxies have a SFR high enough to have built up their entire stellar mass within their Hubble time. We find only marginal evidence for an AGN contribution to the near-IR and mid-IR SEDs. (abridged)

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

Using HST/WFC3 imaging taken as part of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS), we examine the role that major galaxy mergers play in triggering active galactic nuclei (AGN) activity at z~2. Our sample consists of 72 moderate-luminosity (Lx ~ 1E42-1E44 erg/s) AGN at 1.5<z<2.5 that are selected using the 4 Msec Chandra observations in the Chandra Deep Field South, the deepest X-ray observations to date. Employing visual classifications, we have analyzed the rest-frame optical morphologies of the AGN host galaxies and compared them to a mass-matched control sample of 216 non-active galaxies at the same redshift. We find that most of the AGN reside in disk galaxies (51.4%), while a smaller percentage are found in spheroids (27.8%). Roughly 16.7% of the AGN hosts have highly disturbed morphologies and appear to be involved in a major merger or interaction, while most of the hosts (55.6%) appear relatively relaxed and undisturbed. These fractions are statistically consistent with the fraction of control galaxies that show similar morphological disturbances. These results suggest that the hosts of moderate-luminosity AGN are no more likely to be involved in an ongoing merger or interaction relative to non-active galaxies of similar mass at z~2. The high disk fraction observed among the AGN hosts also appears to be at odds with predictions that merger-driven accretion should be the dominant AGN fueling mode at z~2, even at moderate X-ray luminosities. Although we cannot rule out that minor mergers are responsible for triggering these systems, the presence of a large population of relatively undisturbed disk-like hosts suggests that secular processes play a greater role in fueling AGN activity at z~2 than previously thought.

We report on the evolution of the number density and size of early-type galaxies from z~2 to z~0. We select a sample of 563 massive (M>10^{10} Msun), passively evolving (SSFR<10^{-2} Gyr^{-1}) and morphologically spheroidal galaxies at 0<z<2.5, using the panchromatic photometry and spectroscopic redshifts available in the GOODS fields. We combine ACS and WFC3 HST images to study the morphology of our galaxies in their optical rest-frame in the whole 0<z<2.5 range. We find that throughout the explored redshift range the passive galaxies selected with our criteria have weak morphological K-correction, with size being slightly smaller in the optical than in the UV rest-frame (by ~20 and ~10 at z>1.2 and z<1.2, respectively). We measure a significant evolution of the mass-size relation of early-type galaxies, with the fractional increment that is almost independent on the stellar mass. Early-type galaxies (ETGs) formed at z>1 appear to be preferentially small, and the evolution of the mass-size relation at z<1 is driven by both the continuous size growth of the compact galaxies and the appearance of new ETGs with large sizes. We also find that the number density of all passive early-type galaxies increases rapidly, by a factor of 5, from z~2 to z~1, and then more mildly by another factor of 1.5 from z~1 to z~0. We interpret these results as the evidence that the bulk of the ETGs are formed at 1<z<3 through a mechanism that leaves very compact remnants. At z<1 the compact ETGs grow gradually in size, becoming normal size galaxies, and at the same time new ETGs with normal-large sizes are formed.

We report the discovery of large amounts of cold (T ~ 10^4 K), chemically young gas in an overdensity of galaxies at redshift z ~ 1.6 in the Great Observatories Origins Deep Survey southern field (GOODS-S). The gas is identified thanks to the ultra-strong Mg II absorption features it imprints in the rest-frame UV spectra of galaxies in the background of the overdensity. There is no evidence that the optically-thick gas is part of any massive galaxy (i.e. M_star > 4x10^9 M_sun), but rather is associated with the overdensity; less massive and fainter galaxies (25.5 < z_850 < 27.5 mag) have too large an impact parameter to be causing ultra-strong absorption systems, based on our knowledge of such systems. The lack of corresponding Fe II absorption features, not detected even in co-added spectra, suggests that the gas is chemically more pristine than the ISM and outflows of star-forming galaxies at similar redshift, including those in the overdensity itself, and comparable to the most metal-poor stars in the Milky Way halo. A crude estimate of the projected covering factor of the high-column density gas (N_H >~ 10^20 cm-2) based on the observed fraction of galaxies with ultra-strong absorbers is C_F ~ 0.04. A broad, continuum absorption profile extending to the red of the interstellar Mg II absorption line by <~ 2000 km/s is possibly detected in two independent co-added spectra of galaxies of the overdensity, consistent with a large-scale infall motion of the gas onto the overdensity and its galaxies. Overall, these findings provides the first tentative evidence of accretion of cold, chemically young gas onto galaxies at high redshift, possibly feeding their star formation activity. The fact that the galaxies are members of a large structure, as opposed to field galaxies, might play a significant role in our ability to detect the accreting gas.

We take advantage of the sensitivity and resolution of Herschel at 100 and 160 micron to directly image the thermal dust emission and investigate the infrared luminosities, L(IR), and dust obscuration of typical star-forming (L*) galaxies at high redshift. Our sample consists of 146 UV-selected galaxies with spectroscopic redshifts 1.5<z<2.6 in the GOODS-North field. Supplemented with deep Very Large Array (VLA) and Spitzer imaging, we construct median stacks at the positions of these galaxies at 24, 100, and 160 micron, and 1.4 GHz. The comparison between these stacked fluxes and a variety of dust templates and calibrations implies that typical star-forming galaxies with UV luminosities L(UV)>1e10 Lsun at z~2 are luminous infrared galaxies (LIRGs) with a median L(IR)=(2.2+/-0.3)e11 Lsun. Typical galaxies at 1.5<z<2.6 have a median dust obscuration L(IR)/L(UV) = 7.1+/-1.1, which corresponds to a dust correction factor, required to recover the bolometric star formation rate (SFR) from the unobscured UV SFR, of 5.2+/-0.6. This result is similar to that inferred from previous investigations of the UV, H-alpha, 24 micron, radio, and X-ray properties of the same galaxies studied here. Stacking in bins of UV slope implies that L* galaxies with redder spectral slopes are also dustier, and that the correlation between UV slope and dustiness is similar to that found for local starburst galaxies. Hence, the rest-frame 30 and 50 micron fluxes validate on average the use of the local UV attenuation curve to recover the dust attenuation of typical star-forming galaxies at high redshift. In the simplest interpretation, the agreement between the local and high redshift UV attenuation curves suggests a similarity in the dust production and stellar and dust geometries of starburst galaxies over the last 10 billion years.

We present the final results from our ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 "z--band dropout" candidates selected from our VLT/Hawk-I imaging survey over three independent fields. In particular we report on two newly discovered galaxies at redshift ~6.7 in the NTT deep field: both galaxies show a Ly-alpha emission line with rest-frame EWs of the order 15-20 A and luminosities of 2-4 X 10^{42} erg/s. We also present the results of ultra-deep observations of a sample of i-dropout galaxies, from which we set a solid upper limit on the fraction of interlopers. Out of the 20 z-dropouts observed we confirm 5 galaxies at 6.7 < z < 7.1. This is systematically below the expectations drawn on the basis of lower redshift observations: in particular there is a significant lack of objects with intermediate Ly-alpha EWs (between 20 and 55 A). We conclude that the trend for the fraction of Ly-alpha emission in LBGs that is constantly increasing from z~3 to z~6 is most probably reversed from z~6 to z~7. Explaining the observed rapid change in the LAE fraction among the drop-out population with reionization requires a fast evolution of the neutral fraction of hydrogen in the Universe. Assuming that the Universe is completely ionized at z=6 and adopting the semi-analytical models of Dijkstra et al. (2011), we find that our data require a change of the neutral hydrogen fraction of the order Delta chi_{HI} ~ 0.6 in a time Delta z ~ 1, provided that the escape fraction does not increase dramatically over the same redshift interval.

In this paper, we present a derivation of the rest-frame 1400A luminosity function (LF) at redshift six from a new application of the maximum likelihood method by exploring the five deepest HST/ACS fields, i.e., the HUDF, two UDF05 fields, and two GOODS fields. We work on the latest improved data products, which makes our results more robust than those of previous studies. We use un-binned data and thereby make optimal use of the information contained in the dataset. We focus on the analysis to a magnitude limit where the completeness is larger than 50% to avoid possibly large errors in the faint end slope that are difficult to quantify. We also take into account scattering in and out of the dropout sample due to photometric errors by defining for each object a probability that it belongs to the dropout sample. We find the best fit Schechter parameters to the z~6 LF are: alpha = 1.87 +/- 0.14, M* = -20.25 +/- 0.23, and phi*=1.77^{+0.62}_{-0.49} * 10^{-3} Mpc^{-3}. Such a steep slope suggests that galaxies, especially the faint ones, are possibly the main sources of ionizing photons in the universe at redshift six. We also combine results from all studies at z~6 to reach an agreement in 95% confidence level that -20.45<M*<-20.05 and -1.90<alpha<-1.55. The luminosity density has been found not to evolve significantly between z~6 and z~5, but considerable evolution is detected from z~6 to z~3.

We report the detection of color gradients in six massive (stellar mass > 10^{10} M_{sun}) and passively evolving (specific SFR < 10^{-11}/yr) galaxies at redshift 1.3<z<2.5 identified in the HUDF using HST ACS and WFC3/IR images. After matching different PSFs, we obtain color maps and multi-band optical/near-IR photometry (BVizYJH) in concentric annuli, from the smallest resolved radial (~1.7 kpc) up to several times the H-band effective radius. We find that the inner regions of these galaxies have redder rest-frame UV--optical colors than the outer parts. The slopes of the color gradients mildly depend on the overall dust obscuration and rest-frame (U-V) color, with more obscured or redder galaxies having steeper color gradients. The z~2 color gradients are also steeper than those of local early-types. The gradient of a single parameter (age, extinction or metallicity) cannot fully explain the observed color gradients. Fitting spatially resolved HST seven-band photometry to stellar population synthesis models, we find that, regardless of assumptions for metallicity gradient, the redder inner regions of the galaxies have slightly higher dust obscuration than the bluer outer regions, although the magnitude depends on the assumed extinction law. The derived age gradient depends on the assumptions for metallicity gradient. We discuss the implications of a number of assumptions for metallicity gradient on the formation and evolution of these galaxies. We find that the evolution of the mass--size relationship from z~2 to z~0 cannot be driven by in--situ extended star formation, implying that accretion or merger is mostly responsible for the evolution. The lack of a correlation between color gradient and stellar mass argues against the metallicity gradient predicted by the monolithic collapse, which would require significant major mergers to evolve into the one observed at z~0. (Abridged)

The study of the Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at z=7 is important for ascertaining their role in the reionization of the Universe. We perform a detailed and critical analysis of the statistical and systematic errors in the z~7 LF determination: we have assembled a large sample of candidate LBGs at z~7 from different surveys, spanning a large variety of areas and depths. In particular, we have combined data from the deep (J<27.4) and ultradeep (J<29.2) surveys recently acquired with the new WFC3 NIR camera on HST, over the GOODS-ERS and the HUDF fields, with ground based surveys in wide and shallow areas from VLT and Subaru. We have used public ACS images in the z-band to select z-dropout galaxies, and other public data both in the blue (BVI) and in the red bands to reject possible low-redshift interlopers. We have compared our results with extensive simulations to quantify the observational effects of our selection criteria as well as the effects of photometric scatter, color selections or the morphology of the candidates. We have found that the number density of faint LBGs at z~7 is only marginally sensitive to the color selection adopted, but it is strongly dependent from the assumption made on the half light distributions of the simulated galaxies, used to correct the observed sample for incompleteness. The slope of the faint end of the LBGs LF has thus a rather large uncertainty, due to the unknown distribution of physical sizes of the z~7 LBGs. We conclude that galaxies at z~7 are unable to reionize the Universe unless there is a significant evolution in the clumpiness of the IGM or in the escape fraction of ionising photons or, alternatively, there is a large population of z~7 LBGs with large physical dimensions but still not detected by the present observations.

We measure the angular two-point correlation function of sub-millimeter galaxies (SMGs) from 1.1-millimeter imaging of the COSMOS field with the AzTEC camera and ASTE 10-meter telescope. These data yields one of the largest contiguous samples of SMGs to date, covering an area of 0.72 degrees^2 down to a 1.26 mJy/beam (1-sigma) limit, including 189 (328) sources with S/N greater than 3.5 (3). We can only set upper limits to the correlation length r_0, modeling the correlation function as a power-law with pre-assigned slope. Assuming existing redshift distributions, we derive 68.3% confidence level upper limits of r_0 < 6-8 h^-1 Mpc at 3.7 mJy, and r_0 < 11-12 h^-1 Mpc at 4.2 mJy. Although consistent with most previous estimates, these upper limits imply that the real r_0 is likely smaller. This casts doubts on the robustness of claims that SMGs are characterized by significantly stronger spatial clustering, (and thus larger mass), than differently selected galaxies at high-redshift. Using Monte Carlo simulations we show that even strongly clustered distributions of galaxies can appear unclustered when sampled with limited sensitivity and coarse angular resolution common to current sub-millimeter surveys. The simulations, however, also show that unclustered distributions can appear strongly clustered under these circumstances. From the simulations, we predict that at our survey depth, a mapped area of two degrees^2 is needed to reconstruct the correlation function, assuming smaller beam sizes of future surveys (e.g. the Large Millimeter Telescope's 6" beam size). At present, robust measures of the clustering strength of bright SMGs appear to be below the reach of most observations.

We report the spectroscopic confirmation of two Lyman break galaxies at redshift > 7. The galaxies were observed as part of an ultra-deep spectroscopic campaign with FORS2 at the ESO/VLT for the confirmation of z~7 ``z--band dropout'' candidates selected from our VLT/Hawk-I imaging survey. Both galaxies show a prominent emission line at 9735A and 9858A respectively: the lines have fluxes of ~ 1.6-1.2 x 10^(-17) erg/s/cm2 and exhibit a sharp decline on the blue side and a tail on the red side. The asymmetry is quantitatively comparable to the observed asymmetry in z ~ 6 Ly-alpha lines, where absorption by neutral hydrogen in the IGM truncates the blue side of the emission line profile. We carefully evaluate the possibility that the galaxies are instead at lower redshift and we are observing either [OII], [OIII] or H-alpha emission: however from the spectroscopic and the photometric data we conclude that there are no other plausible identifications, except for Ly-alpha at redshift >7, implying that these are two of the most robust redshift determination for galaxies in the reionization epoch. Based on their redshifts and broad--band photometry, we derive limits on the star formation rate and on the ultraviolet spectral slopes of the two galaxies. We argue that these two galaxies alone are unlikely to have ionized the IGM in their surroundings.

We present three bright z+ dropout candidates selected from deep Near-Infrared (NIR) imaging of the COSMOS 2 square degree field. All three objects match the 0.8-8um colors of other published z>7 candidates but are three magnitudes brighter, facilitating further study. Deep spectroscopy of two of the candidates covering 0.64-1.02um with Keck-DEIMOS and all three covering 0.94-1.10um and 1.52-1.80um with Keck-NIRSPEC detects weak spectral features tentatively identified as Ly-alpha at z=6.95 and z=7.69 in two of the objects. The third object is placed at z~1.6 based on a 24um and weak optical detection. A comparison with the spectral energy distributions of known z<7 galaxies, including objects with strong spectral lines, large extinction, and large systematic uncertainties in the photometry yields no objects with similar colors. However, the lambda>1um properties of all three objects can be matched to optically detected sources with photometric redshifts at z~1.8, so the non-detection in the i+ and z+ bands are the primary factors which favors a z>7 solution. If any of these objects are at z~7 the bright end of the luminosity function is significantly higher at z>7 than suggested by previous studies, but consistent within the statistical uncertainty and the dark matter halo distribution. If these objects are at low redshift, the Lyman-Break selection must be contaminated by a previously unknown population of low redshift objects with very strong breaks in their broad band spectral energy distributions and blue NIR colors. The implications of this result on luminosity function evolution at high redshift is discussed. We show that the primary limitation of z>7 galaxy searches with broad filters is the depth of the available optical data.

We explore methods to improve the estimates of star formation rates and mean stellar population ages from broadband photometry of high redshift star-forming galaxies. We use synthetic spectral templates with a variety of simple parametric star formation histories to fit broadband spectral energy distributions. These parametric models are used to infer ages, star formation rates and stellar masses for a mock data set drawn from a hierarchical semi-analytic model of galaxy evolution. Traditional parametric models generally assume an exponentially declining rate of star-formation after an initial instantaneous rise. Our results show that star formation histories with a much more gradual rise in the star formation rate are likely to be better templates, and are likely to give better overall estimates of the age distribution and star formation rate distribution of Lyman break galaxies. For B- and V-dropouts, we find the best simple parametric model to be one where the star formation rate increases linearly with time. The exponentially-declining model overpredicts the age by 100 % and 120 % for B- and V-dropouts, on average, while for a linearly-increasing model, the age is overpredicted by 9 % and 16 %, respectively. Similarly, the exponential model underpredicts star-formation rates by 56 % and 60 %, while the linearly-increasing model underpredicts by 15 % 22 %, respectively. For U-dropouts, the models where the star-formation rate has a peak (near z ~ 3) provide the best match for age -- overprediction is reduced from 110 % to 26 % -- and star-formation rate -- underprediction is reduced from 58 % to 22 %. We classify different types of star-formation histories in the semi-analytic models and show how the biases behave for the different classes. We also provide two-band calibration formulae for stellar mass and star formation rate estimations.

Distant galaxies show correlations between their current star-formation rates (SFRs) and stellar masses, implying that their star-formation histories (SFHs) are highly similar. Moreover, observations show that the UV luminosities and stellar masses grow from z=8 to 3, implying that the SFRs increase with time. We compare the cosmologically averaged evolution in galaxies at 3 < z < 8 at constant comoving number density, n = 2 x 10^-4 Mpc^-3. This allows us to study the evolution of stellar mass and star formation in the galaxy predecessors and descendants in ways not possible using galaxies selected at constant stellar mass or SFR, quantities that evolve strongly in time. We show that the average SFH of these galaxies increase smoothly from z=8 to 3 as SFR ~ t^alpha with alpha = 1.7 +/- 0.2. This conflicts with assumptions that the SFR is either constant or declines exponentially in time. We show that the stellar mass growth in these galaxies is consistent with this derived SFH. This provides evidence that the slope of the high-mass end of the IMF is approximately Salpeter unless the duty cycle of star formation is much less than unity. We argue that these relations follow from gas accretion (either through accretion or delivered by mergers) coupled with galaxy disk growth under the assumption that the SFR depends on the local gas surface density. This predicts that gas fractions decrease from z=8 to 3 on average as f_gas ~ (1+z)^0.9 for galaxies with this number density. The implied galaxy gas accretion rates at z > 4 are as fast and may even exceed the SFR: this is the "gas accretion epoch". At z < 4 the SFR overtakes the implied gas accretion rate, indicating a period where galaxies consume gas faster than it is acquired. At z < 3, galaxies with this number density depart from these relations implying that star formation and gas accretion are slowed at later times.

We use the deepest and the most comprehensive photometric data currently available for GOODS-South galaxies to measure their photometric redshifts. The photometry includes VLT/VIMOS (U-band), HST/ACS (F435W, F606W, F775W, and F850LP bands), VLT/ISAAC (J-, H-, and Ks-bands), and four Spitzer/IRAC channels (3.6, 4.5, 5.8, and 8.0 micron). The catalog is selected in the z-band (F850LP) and photometry in each band is carried out using the recently completed TFIT algorithm, which performs PSF matched photometry uniformly across different instruments and filters, despite large variations in PSFs and pixel scales. Photometric redshifts are derived using the GOODZ code, which is based on the template fitting method using priors. The code also implements "training" of the template SED set, using available spectroscopic redshifts in order to minimize systematic differences between the templates and the SEDs of the observed galaxies. Our final catalog covers an area of 153 sq. arcmin and includes photometric redshifts for a total of 32,505 objects. The scatter between our estimated photometric and spectroscopic redshifts is sigma=0.040 with 3.7% outliers to the full z-band depth of our catalog, decreasing to sigma=0.039 and 2.1% outliers at a magnitude limit m(z)<24.5. This is consistent with the best results previously published for GOODS-S galaxies, however, the present catalog is the deepest yet available and provides photometric redshifts for significantly more objects to deeper flux limits and higher redshifts than earlier works. Furthermore, we show that the photometric redshifts estimated here for galaxies selected as dropouts are consistent with those expected based on the Lyman break technique.

We present ultradeep optical spectroscopy obtained with FORS2 on VLT of seven Lyman-break galaxy (LBG) candidates at z>6.5 selected in the GOODS-S field from Hawk-I/VLT and WFC3/HST imaging. For one galaxy we detect a low significance emission line (S/N< 7), located at 9691.5 +/- 0.5A and with flux 3.4 x 10^(-18)erg/cm^2/s. If identified as Lyman alpha, it places the LBG at redshift z=6.972+/- 0.002, with a rest-frame equivalent width EW}=13A. Using Monte Carlo simulations and conservative EW distribution functions at 2<z<6, we estimate that the probability of observing no galaxies in our data with S/N>10 is ~ 2%, and that of observing only one galaxy out of seven with S/N=5 is ~4%, but these can be as small as ~1E-3, depending on the details of the EW distribution. We conclude that either a significant fraction of the candidates is not at high redshift or that some physical mechanism quenches the Lyman alpha emission emerging from the galaxies at z>6.5, abruptly reversing the trend of the increasing fraction of strong emitters with increasing redshift observed up to z~ 6.5. We discuss the possibility that an increasingly neutral intergalactic medium is responsible for such quenching.

We present the details and early results from a deep near-infrared survey utilising the NICMOS instrument on the Hubble Space Telescope centred around massive M_* > 10^11 M_0 galaxies at 1.7 < z < 2.9 found within the Great Observatories Origins Deep Survey (GOODS) fields. The GOODS NICMOS Survey (GNS) was designed to obtain deep F160W (H-band) imaging of 80 of these massive galaxies, as well as other colour selected objects such as Lyman-break drop-outs, BzK objects, Distant Red Galaxies, EROs, Spitzer Selected EROs, BX/BM galaxies, as well as sub-mm galaxies. We present in this paper details of the observations, our sample selection, as well as a description of features of the massive galaxies found within our survey fields. This includes: photometric redshifts, rest-frame colours, and stellar masses. We furthermore provide an analysis of the selection methods for finding massive galaxies at high redshifts, including colour selection, and how galaxy populations selected through different methods overlap. We find that a single colour selection method cannot locate all of the massive galaxies, with no one method finding more than 70 percent. We however find that the combination of these colour methods finds nearly all the massive galaxies, as selected by photometric redshifts with the exception of apparently rare blue massive galaxies. By investigating the rest-frame (U-B) vs. M_B diagram for these galaxies we furthermore show that there exists a bimodality in colour-magnitude space at z < 2, driven by stellar mass, such that the most massive galaxies are systematically red up to z~2.5, while lower mass galaxies tend to be blue. We also discuss the number densities for galaxies with stellar masses M_* > 10^11 M_0, whereby we find an increase of a factor of eight between z = 3 and z = 1.5, demonstrating that this is an epoch when massive galaxies establish most of their mass.

We present a study of the cosmic infrared background, which is a measure of the dust obscured activity in all galaxies in the Universe. We venture to isolate the galaxies responsible for the background at 1mm; with spectroscopic and photometric redshifts we constrain the redshift distribution of these galaxies. We create a deep 1.16mm map (sigma ~ 0.5mJy) by combining the AzTEC 1.1mm and MAMBO 1.2mm datasets in GOODS-N. This combined map contains 41 secure detections, 13 of which are new. By averaging the 1.16mm flux densities of individually undetected galaxies with 24um flux densities > 25uJy, we resolve 31--45 per cent of the 1.16mm background. Repeating our analysis on the SCUBA 850um map, we resolve a higher percentage (40--64 per cent) of the 850um background. A majority of the background resolved (attributed to individual galaxies) at both wavelengths comes from galaxies at z > 1.3. If the ratio of the resolved submillimeter to millimeter background is applied to a reasonable scenario for the origins of the unresolved submillimeter background, 60--88 per cent of the total 1.16mm background comes from galaxies at z > 1.3.

We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrow-band imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f_esc) of ionizing radiation for L >~ L*(z=3) Lyman Break Galaxies (LBGs) at redshift 3.4--4.5. Only one LBG, at redshift z=3.795, is detected in its Lyman continuum (LyC; S/N~5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R_eff=0.8, kpc physical). Three out of seven AGN are also detected in their LyC, including one at redshift z=3.951 and z850 = 26.1. From stacked data (LBGs) we set an upper limit to the average f_esc in the range 5%--20%, depending on the how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal and Gaussian) are explored. We find that the median f_esc is lower than ~6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f_esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f_esc, and estimates of f_esc from a sample of faint LBG from the HUDF (i775<28.5) are in broad quantitative agreement with such a scenario.

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a concept for an 8-meter to 16-meter UVOIR space observatory for launch in the 2025-2030 era. ATLAST will allow astronomers to answer fundamental questions at the forefront of modern astronphysics, including "Is there life elsewhere in the Galaxy?" We present a range of science drivers that define the main performance requirements for ATLAST (8 to 16 milliarcsec angular resolution, diffraction limited imaging at 0.5 {\mu}m wavelength, minimum collecting area of 45 square meters, high sensitivity to light wavelengths from 0.1 {\mu}m to 2.4 {\mu}m, high stability in wavefront sensing and control). We will also discuss the synergy between ATLAST and other anticipated future facilities (e.g., TMT, EELT, ALMA) and the priorities for technology development that will enable the construction for a cost that is comparable to current generation observatory-class space missions.

(Abridged) We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained - 30 observing hours with Hawk-I in the Y-band of two high galactic latitude fields. We combined the Y-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in -161 arcmin^2 of our survey to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a best-fit SED with photometric redshift z= 6.85 +0.20 -0.15 and an E(B-V)=0.05 +0.15 -0.05. We compute a binned estimate of the z ~ 7 LF and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than M_1500= -19.0, we derive a luminosity density rho_UV = 1.5^{+2.1}{-0.8} x 10^25 erg/s/Hz/Mpc^3, implying a decrease by a factor 3.5 from z=6 to z=6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I J and K band observations to derive an upper limit on the number density of M1500<~ -22.0 LBGs at z-8 (Y-dropouts).

We present the results of rest-frame, UV slitless spectroscopic observations of a sample of 32 z~0.7 Lyman break galaxy (LBG) analogs in the COSMOS field. The spectroscopic search was performed with the Solar Blind Channel (SBC) on Hubble Space Telescope. We report the detection of leaking Lyman continuum (LyC) radiation from an AGN-starburst composite. While we find no direct detections of LyC emission in the remainder of our sample, we achieve individual lower limits (3 sigma) of the observed non-ionizing UV to LyC flux density ratios, f_{nu}(1500A)/f_{nu}(830A) of 20 to 204 (median of 73.5) and 378.7 for the stack. Assuming an intrinsic Lyman break of 3.4 and an intergalactic medium (IGM) transmission of LyC photons along the line of sight to the galaxy of 85% we report an upper limit for the relative escape fraction in individual galaxies of 0.02 - 0.19 and a stacked 3 sigma upper limit of 0.01. We find no indication of a relative escape fraction near unity as seen in some LBGs at z~3. Our UV spectra achieve the deepest limits to date at any redshift for the escape fraction in individual sources. The contrast between these z~0.7 low escape fraction LBG analogs with z~3 LBGs suggests that either the processes conducive to high escape fractions are not being selected for in the z<1 samples or the average escape fraction is decreasing from z~3 to z~1. We discuss possible mechanisms which could affect the escape of LyC photons.

Aims: We investigate the nature of the source GDS J033218.92-275302.7at redshift ~ 5.56. Methods: The SED of the source is well sampled by 16 bands photometry, from UV-optical, near infrared and mid-infrared (MID-IR).The detection of signal in the MID-IR Spitzer/IRAC bands 5.8, 8.0 um -- where the nebular emission contribution is less effective -- suggests the presence of a Balmer break, signature of an underlying stellar population formed at earlier epochs. The optical spectrum shows a clear Lya emission line together with semi-forbidden NIV] 1483.3-1486.5 also in emission. Results: From the SED fitting and the Lya modelling it turns out that the source seems to have an evolved component with stellar mass of ~5 x10^(10) Msolar and age ~ 0.4 Gyrs, and a young component with an age of ~ 0.01 Gyrs and SFR in the range of 30-200 Msolar yr^(-1). The limits on the effective radius derived from the ACS/z850 and VLT/Ks bands indicate that this galaxy is denser than the local ones with similar mass. A relatively high nebular gas column density is favored from the Lya line modelling (NHI>=10^(21) cm^(-2)). A vigorous outflow (~ 450 km/s) has been measured from the optical spectrum,consistent with the Lya modelling. From ACS observations it turns out that the region emitting Lya photons is spatially compact and of the same order of the effective radius estimated at the ~1400A rest-frame wavelength, whose emission is dominated by the stellar continuum and/or AGN. The gas is blown out from the central region,but given the mass of the galaxy it is uncertain whether it will pollute the IGM to large distances. We argue that a burst of star formation in a dense gas environment is active (possibly containing hot and massive stars and/or a low luminosity AGN), superimposed to an already formed fraction of stellar mass (abridged).

We study the physical characteristics of galaxies at 6.3 < z < 8.6, selected from deep near-infrared imaging with the Wide Field Camera 3 (WFC3) on board the Hubble Space Telescope. Accounting for the photometric scatter using simulations, galaxies at z ~ 7 have bluer UV colors compared to typical local starburst galaxies at > 4 sigma confidence. Although these colors necessitate young ages (<100 Myr), low or zero dust attenuation, and low metallicities, these are explicable by normal (albeit unreddened) stellar populations, with no evidence for near-zero metallicities and/or top-heavy initial mass functions. The age of the Universe at these redshifts limits the amount of stellar mass in late-type populations, and the WFC3 photometry implies galaxy stellar masses ~ 10^8 - 10^9 Msol for Salpeter initial mass functions to a limiting magnitude of M_1500 ~ -18. The masses of ``characteristic'' (L*) z > 7 galaxies are smaller than those of L* Lyman break galaxies (LBGs) at lower redshifts, and are comparable to less evolved galaxies selected on the basis of their Lyman alpha emission at 3 < z < 6, implying that the 6.3 < z < 8.6 galaxies are the progenitors of more evolved galaxies at lower redshifts. We estimate that Lyman alpha emission is able to contribute to the observed WFC3 colors of galaxies at these redshifts, with an estimated typical line flux of ~ 10^-18 erg s^-1 cm^-2, roughly a factor of four below currently planned surveys. The integrated UV specific luminosity for the detected galaxies at z ~ 7 and z ~ 8 is within factors of a few of that required to reionize the IGM assuming low clumping factors, implying that in order to reionize the Universe galaxies at these redshifts have a high ( ~ 50%) escape fraction of Lyman continuum photons, possibly substantiated by the very blue colors of this population.

We present observations of a very massive galaxy at z=1.82 which show that its morphology, size, velocity dispersion and stellar population properties that are fully consistent with those expected for passively evolving progenitors of today's giant ellipticals. These findings are based on a deep optical rest-frame spectrum obtained with the Multi-Object InfraRed Camera and Spectrograph (MOIRCS) on the Subaru telescope of a high-z passive galaxy candidate (pBzK) from the COSMOS field, for which we accurately measure its redshift of z=1.8230 and obtain an upper limit on its velocity dispersion sigma_star<326 km/s. By detailed stellar population modeling of both the galaxy broad-band SED and the rest-frame optical spectrum we derive a star-formation-weighted age and formation redshift of t_sf~1-2 Gyr and z_form~2.5-4, and a stellar mass of M_star~(3-4)x10^{11} M_sun. This is in agreement with a virial mass limit of M_vir<7x10^{11}M_sun, derived from the measured sigma_star value and stellar half-light radius, as well as with the dynamical mass limit based on the Jeans equations. In contrast with previously reported super-dense passive galaxies at z~2, the present galaxy at z=1.82 appears to have both size and velocity dispersion similar to early-type galaxies in the local Universe with similar stellar mass. This suggests that z~2 massive and passive galaxies may exhibit a wide range of properties, then possibly following quite different evolutionary histories from z~2 to z=0.

We present near-IR images, obtained with the Hubble Space Telescope (HST) and the WFC3/IR camera, of six passive and massive galaxies at redshift 1.3<z<2.4 (SSFR<10^{-2} Gyr^{-1}; stellar mass M~10^{11} M_{sun}), selected from the Great Observatories Origins Deep Survey (GOODS). These images, which have a spatial resolution of ~1.5 kpc, provide the deepest view of the optical rest-frame morphology of such systems to date. We find that the light profile of these galaxies is regular and well described by a Sersic model with index typical of today's spheroids. Their size, however, is generally much smaller than today's early types of similar stellar mass, with four out of six galaxies having r_e ~ 1 kpc or less, in quantitative agreement with previous similar measures made at rest-frame UV wavelengths. The images reach limiting surface brightness mu~26.5 mag arcsec^{-2} in the F160W bandpass; yet, there is no evidence of a faint halo in the galaxies of our sample, even in their stacked image. We also find that these galaxies have very weak "morphological k-correction" between the rest-frame UV (from the ACS z-band), and the rest--frame optical (WFC3 H-band): the Sersic index, physical size and overall morphology are independent or only mildly dependent on the wavelength, within the errors.

[Abridged] We present the first results from a 1.1 mm confusion-limited map of the GOODS-S field taken with AzTEC on the ASTE telescope. We imaged a 270 sq. arcmin field to a 1\sigma depth of 0.48 - 0.73 mJy/beam, making this one of the deepest blank-field surveys at mm-wavelengths ever achieved. Although our GOODS-S map is extremely confused, we demonstrate that our source identification and number counts analyses are robust, and the techniques discussed in this paper are relevant for other deeply confused surveys. We find a total of 41 dusty starburst galaxies with S/N >= 3.5 within this uniformly covered region, where only two are expected to be false detections. We derive the 1.1mm number counts from this field using both a "P(d)" analysis and a semi-Bayesian technique, and find that both methods give consistent results. Our data are well-fit by a Schechter function model with (S', N(3mJy), \alpha) = (1.30+0.19 mJy, 160+27 (mJy/deg^2)^(-1), -2.0). Given the depth of this survey, we put the first tight constraints on the 1.1 mm number counts at S(1.1mm) = 0.5 mJy, and we find evidence that the faint-end of the number counts at S(850\mu m) < 2.0 mJy from various SCUBA surveys towards lensing clusters are biased high. In contrast to the 870 \mu m survey of this field with the LABOCA camera, we find no apparent under-density of sources compared to previous surveys at 1.1 mm. Additionally, we find a significant number of SMGs not identified in the LABOCA catalogue. We find that in contrast to observations at wavelengths < 500 \mu m, MIPS 24 \mu m sources do not resolve the total energy density in the cosmic infrared background at 1.1 mm, demonstrating that a population of z > 3 dust-obscured galaxies that are unaccounted for at these shorter wavelengths potentially contribute to a large fraction (~2/3) of the infrared background at 1.1 mm.

We have obtained deep Hubble Space Telescope far-UV images of 15 starburst galaxies at z~1.3 in the GOODS fields to search for escaping Lyman continuum photons. These are the deepest far-UV images m_{AB}=28.7, 3\sigma, 1" diameter) over this large an area (4.83 arcmin^2) and provide the best escape fraction constraints for any galaxy at any redshift. We do not detect any individual galaxies, with 3\sigma limits to the Lyman Continuum (~700 \AA) flux 50--149 times fainter (in f_nu) than the rest-frame UV (1500 \AA) continuum fluxes. Correcting for the mean IGM attenuation (factor ~2), as well as an intrinsic stellar Lyman Break (~3), these limits translate to relative escape fraction limits of f_{esc,rel}<[0.03,0.21]. The stacked limit is f_{esc,rel}(3\sigma)<0.02. We use a Monte Carlo simulation to properly account for the expected distribution of IGM opacities. When including constraints from previous surveys at z~1.3 we find that, at the 95% confidence level, no more than 8% of star--forming galaxies at z~1.3 can have relative escape fractions greater than 0.50. Alternatively, if the majority of galaxies have low, but non-zero, escaping Lyman Continuum, the escape fraction can not be more than 0.04. Both the stacked limits, and the limits from the Monte Carlo simulation suggest that the average ionizing emissivity (relative to non-ionizing UV emissivity) at z~1.3 is significantly lower than has been observed in Lyman Break Galaxies (LBGs) at z~3. If the ionizing emissivity of star-forming galaxies is in fact increasing with redshift, it would help to explain the high photoionization rates seen in the IGM at z>4 and reionization of the intergalactic medium at z>6. [Abridged]

We perform a deep search for galaxies in the redshift range 6.5<z<7.5, to measure the evolution of the number density of luminous galaxies in this redshift range and derive useful constraints on the evolution of their Luminosity Function. We present here the first results of an ESO Large Program, that exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained two Hawk-I pointings on the GOODS South field for a total of 32 observing hours, covering ~90 arcmin2. The images reach Y=26.7 mags for the two fields. We have used public ACS images in the z band to select z-dropout galaxies with the colour criteria Z-Y>1, Y-J<1.5 and Y-K<2. The other public data in the UBVRIJHK bands are used to reject possible low redshift interlopers. The output has been compared with extensive Monte Carlo simulations to quantify the observational effects of our selection criteria as well as the effects of photometric errors. We detect 7 high quality candidates in the magnitude range Y=25.5-26.7. This interval samples the critical range for M* at z>6 (M_1500 ~- 19.5 to -21.5). After accounting for the expected incompleteness, we rule out at a 99% confidence level a Luminosity Function constant from z=6 to z=7, even including the effects of cosmic variance. For galaxies brighter than M_1500=-19.0 we derive a luminosity density rho_UV=1.5^{+2.0}_{-0.9} 10^25 erg/s/Hz/Mpc3, implying a decrease by a factor 3.5 from z=6 to z~6.8. On the basis of our findings, we make predictions for the surface densities expected in future surveys surveys, based on ULTRA-VISTA, HST-WFC3 or JWST-NIRCam, evaluating the best observational strategy to maximise their impact.

We present results of our large-area survey for z'-band dropout galaxies at z=7 in a 1568 arcmin^2 sky area covering the SDF and GOODS-N fields. Combining our ultra-deep Subaru/Suprime-Cam z'- and y-band (lambda_eff=1um) images with legacy data of Subaru and HST, we have identified 22 bright z-dropout galaxies down to y=26, one of which has a spectroscopic redshift of z=6.96 determined from Lya emission. The z=7 luminosity function (LF) yields the best-fit Schechter parameters of phi*=0.69 +2.62/-0.55 x10^(-3) Mpc^(-3), Muv*=-20.10 +/-0.76 mag, and alpha=-1.72 +/-0.65, and indicates a decrease from z=6 at a >95% confidence level. This decrease is beyond the cosmic variance in our two fields, which is estimated to be a factor of <~2. We have found that the cosmic star formation rate density drops from the peak at z=2-3 to z=7 roughly by a factor of ~10 but not larger than ~100. A comparison with the reionization models suggests either that the Universe could not be totally ionized by only galaxies at z=7, or more likely that properties of galaxies at z=7 are different from those at low redshifts having, e.g., a larger escape fraction (>~0.2), and/or a flatter IMF. Our SDF z-dropout galaxies appear to form 60-Mpc long filamentary structures, and the z=6.96 galaxy with Lya emission is located at the center of an overdense region consisting of four UV bright dropout candidates, which might suggest an existence of a well-developed ionized bubble at z=7.

We provide a systematic measurement of the rest-frame UV continuum slope beta over a wide range in redshift (z~2-6) and rest-frame UV luminosity (0.1-2L*) to improve estimates of the SFR density at high redshift. We utilize the deep optical and infrared data (ACS/NICMOS) over the CDF-S and HDF-N GOODS fields, as well as the UDF for our primary UBVi "dropout" sample. We correct the observed distributions for selection biases and photometric scatter. We find that the UV-continuum slope of the most luminous galaxies is substantially redder at z~2-4 than it is at z~5-6. Lower luminosity galaxies are also found to be bluer than higher luminosity galaxies at z~2.5 and z~4. We do not find a large number of galaxies with beta's as red as -1 in our dropout selections at z~4, and particularly at z>~5, even though such sources could be readily selected from our data. This suggests that star-forming galaxies at z>~5 almost universally have very blue UV-continuum slopes, and that there are not likely to be a substantial number of dust-obscured galaxies at z>~5 that are missed in "dropout" searches. Using the same relation between UV-continuum slope and dust extinction as found to be appropriate at z~0 and z~2, we estimate the average dust extinction of galaxies as a function of redshift and UV luminosity in a consistent way. We find that the estimated dust extinction increases substantially with cosmic time for the most UV luminous galaxies, but remains small (<~2x) at all times for lower luminosity galaxies. Because these same lower luminosity galaxies dominate the luminosity density in the UV, the overall dust extinction correction remains modest at all redshifts. We include the contribution from ULIRGs in our SFR density estimates at z~2-6, but find that they contribute only ~20% of the total at z~2.5 and <~10% at z>~4.

We present deep imaging in the {\em U} band covering an area of 630 arcmin$^{2}$ centered on the southern field of the Great Observatories Origins Deep Survey (GOODS). The data were obtained with the VIMOS instrument at the ESO Very Large Telescope. The final images reach a magnitude limit $U_{lim} \approx 29.8$ (AB, 1$\sigma$, in a 1\arcsec radius aperture), and have good image quality, with full width at half maximum $\approx 0.8\arcsec$. They are significantly deeper than previous U--band images available for the GOODS fields, and better match the sensitivity of other multi--wavelength GOODS photometry. The deeper U--band data yield significantly improved photometric redshifts, especially in key redshift ranges such as $2<z<4$, and deeper color--selected galaxy samples, e.g., Lyman--break galaxies at $z\approx 3$. We also present the coaddition of archival ESO VIMOS R band data, with $R_{lim} \approx 29$ (AB, 1$\sigma$, 1\arcsec radius aperture), and image quality $\approx 0.75 \arcsec$. We discuss the strategies for the observations and data reduction, and present the first results from the analysis of the coadded images.

We present detailed morphological properties of Lyman alpha emitters (LAEs) at z~ 5.7 in the COSMOS field, based on {\it Hubble Space Telescope} Advanced Camera for Surveys (ACS) data. The ACS imaging in the F814W filter covered 85 LAEs of the 119 LAEs identified in the full two square degree field, and 47 LAEs of them are detected in the ACS images. Nearly half of them are spatially extended with a size larger than 0.15 arcsec (~0.88 kpc at z=5.7) up to 0.4 arcsec (~2.5 kpc at z=5.7). The others are nearly unresolved compact objects. Two LAEs show double-component structures, indicating interaction or merging of building components to form more massive galaxies. By stacking the ACS images of all the detected sources, we obtain a Sersic parameter of n~0.7 with a half-light radius of 0.13 arcsec (0.76 kpc), suggesting that the majority of ACS detected LAEs have not spheroidal-like but disk-like or irregular light profiles. Comparing ACS F814W magnitudes (I_814) with Subaru/Suprime-Cam magnitudes in the NB816, i', and z' bands, we find that the ACS imaging in the F814W band mainly probes UV continuum rather than Lyman alpha line emission. UV continuum sizes tend to be larger for LAEs with larger Ly\alpha emission regions as traced by the NB816 imaging. The non-detection of 38 LAEs in the ACS images is likely due to the fact that their surface brightness is even too low both in the UV continuum and Lyalpha emission. Estimating I_814 for the LAEs with ACS non-detection from the z' and NB816 magnitudes, we find that 16 of these are probably LAEs with a size larger than 0.15 arcsec in UV continuum. All these results suggest that our LAE sample contains systematically larger LAEs in UV continuum size than those previously studied at z~6.

The Advanced Technology Large-Aperture Space Telescope (ATLAST) is a set of mission concepts for the next generation of UVOIR space observatory with a primary aperture diameter in the 8-m to 16-m range that will allow us to perform some of the most challenging observations to answer some of our most compelling questions, including "Is there life elsewhere in the Galaxy?" We have identified two different telescope architectures, but with similar optical designs, that span the range in viable technologies. The architectures are a telescope with a monolithic primary mirror and two variations of a telescope with a large segmented primary mirror. This approach provides us with several pathways to realizing the mission, which will be narrowed to one as our technology development progresses. The concepts invoke heritage from HST and JWST design, but also take significant departures from these designs to minimize complexity, mass, or both. Our report provides details on the mission concepts, shows the extraordinary scientific progress they would enable, and describes the most important technology development items. These are the mirrors, the detectors, and the high-contrast imaging technologies, whether internal to the observatory, or using an external occulter. Experience with JWST has shown that determined competitors, motivated by the development contracts and flight opportunities of the new observatory, are capable of achieving huge advances in technical and operational performance while keeping construction costs on the same scale as prior great observatories.

We have carried out a search for galaxies at z ~ 7-10 in ~14.4 sq. arcmin of new NICMOS parallel imaging taken in the Great Observatories Origins Deep Survey (GOODS, 5.9 sq. arcmin), the Cosmic Origins Survey (COSMOS, 7.2 sq. arcmin), and SSA22 (1.3 sq. arcmin). These images reach 5 sigma sensitivities of J110 = 26.0-27.5 (AB), and combined they increase the amount of deep near-infrared data by more than 60% in fields where the investment in deep optical data has already been made. We find no z>7 candidates in our survey area, consistent with the Bouwens et al. (2008) measurements at z~7 and 9 (over 23 sq. arcmin), which predict 0.7 galaxies at z~7 and <0.03 galaxies at z~9. We estimate that 10-20% of z>7 galaxies are missed by this survey, due to incompleteness from foreground contamination by faint sources. For the case of luminosity evolution, assuming a Schecter parameterization with a typical phi* = 10^-3 Mpc^-3, we find M* > -20.0 for z~7 and M* > -20.7 for z~9 (68% confidence). This suggests that the downward luminosity evolution of LBGs continues to z~7, although our result is marginally consistent with the z~6 LF of Bouwens et al.(2006, 2007). In addition we present newly-acquired deep MMT/Megacam imaging of the z~9 candidate JD2325+1433, first presented in Henry et al. (2008). The resulting weak but significant detection at i' indicates that this galaxy is most likely an interloper at z~2.7.

We present a simple formalism to interpret two galaxy statistics, the UV luminosity function and two-point correlation functions for star-forming galaxies at z~4, 5, 6 in the context of LCDM cosmology. Both statistics are the result of how star formation takes place in DM halos, and thus are used to constrain how UV light depends on halo properties such as mass. The two measures were taken from the GOODS data, thus ideal for joint analysis. The two physical quantities we explore are the SF duty cycle, and the range of L_UV that a halo of mass M can have (mean and variance). The former addresses the typical duration of SF activity in halos while the latter addresses the averaged SF history and regularity of gas inflow into these systems. We explore various physical models consistent with data, and find the following: 1) the typical duration of SF observed in the data is <0.4 Gyr (1 sig), 2) the inferred scaling law between L_UV and halo mass M from the observed slope of the LFs is roughly linear at all redshifts, and 3) L_UV for a fixed halo mass decreases with time, implying that the SF efficiency (after dust extinction) is higher at earlier times. We explore several physical scenarios relating star formation to halo mass, but find that these scenarios are indistinguishable due to the limited range of halo mass probed by our data. In order to discriminate between different scenarios, we discuss constraining the bright-faint galaxy cross-correlation functions and luminosity-dependence of galaxy bias. (Abridged)

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

We present results of a survey for Ly alpha emitters at z=4.86 based on optical narrowband (lambda_c=7126 angstrom, Delta lambda=73 angstrom) and broadband (B, V, r', i', and z') observations of the Cosmic Evolution Survey (COSMOS) field using Suprime-Cam on the Subaru Telescope. We find 79 LAE candidates at z=4.86 over a contiguous survey area of 1.83 deg^2, down to the Ly alpha line flux of 1.47 x 10^-17 ergs s^-1 cm^-2. We obtain the Ly alpha luminosity function with a best-fit Schechter parameters of log L^*=42.9^+0.5_-0.3 ergs s^-1 and phi^* = 1.2^+8.0_-1.1 x 10^-4 Mpc^-3 for alpha=-1.5 (fixed). The two-point correlation function for our LAE sample is xi(r) = (r/4.4^+5.7_-2.9 Mpc)^-1.90+/-0.22. In order to investigate the field-to-field variations of the properties of Ly alpha emitters, we divide the survey area into nine tiles of 0.5^circ x 0.5^circ each. We find that the number density varies with a factor of ~ 2 from field to field with high statistical significance. However, we find no significant field-to-field variance when we divide the field into four tiles with 0.7^circ x 0.7^circ each. We conclude that at least 0.5 deg^2 survey area is required to derive averaged properties of LAEs at z~5, and our survey field is wide enough to overcome the cosmic variance.

We report on observations of Lyman-break galaxies (LBGs) selected from the Great Observatories Origins Deep Survey at mean redshift z~4, 5 and 6, obtained with FORS2 spectrograph at the ESO-VLT.This program has yielded spectroscopic identifications for 114 galaxies of which 51 are at z~4, 31 at z~5, and 32 at z~6. The adopted selection criteria are effective, identifying galaxies at the expected redshift with minimal foreground contamination. Once selection effects are properly accounted for, the rest-frame UV spectra of the higher-redshift LBGs appear to be similar to their counterparts at z~3. As at z~3, LBGs at z~4 and z~5 show Lya both in emission and in absorption; when in absorption, strong interstellar lines and red UV spectra are observed, a fact also measured at z~2 and 3. At z~6, sensitivity issues bias our sample towards galaxies with Lya in emission; nevertheless, these spectra appear to be similar to their lower-redshift counterparts. As in other studies at similar redshifts, we find clear evidence that brighter LBGs tend to have weaker Lya emission lines. At fixed rest-frame UV luminosity, the equivalent width of the Lya emission line is larger at higher redshifts. At all redshifts where the measurements can be reliably made, the redshift of the Lya emission line turns out to be larger than that of the interstellar absorption lines,with a median velocity difference DeltaV~400 km/s at z~4 and 5, consistent with results at lower redshifts. This shows that powerful, large-scale winds are common at high redshift.In general,there is no strong correlation between the morphology of the UV light and the spectroscopic properties. However, galaxies with deep interstellar absorption lines and strong Lya absorption appear to be more diffuse than galaxies with Lya in emission.(abridged)

We investigate the biases and uncertainties in estimates of physical parameters of high-redshift Lyman break galaxies (LBGs), such as stellar mass, mean stellar population age, and star formation rate (SFR), obtained from broad-band photometry. By combining LCDM hierarchical structure formation theory, semi-analytic treatments of baryonic physics, and stellar population synthesis models, we construct model galaxy catalogs from which we select LBGs at redshifts z ~ 3.4, 4.0, and 5.0. The broad-band spectral energy distributions (SEDs) of these model LBGs are then analysed by fitting galaxy template SEDs derived from stellar population synthesis models with smoothly declining SFRs. We compare the statistical properties of LBGs' physical parameters -- such as stellar mass, SFR, and stellar population age -- as derived from the best-fit galaxy templates with the intrinsic values from the semi-analytic model. We find some trends in these distributions: first, when the redshift is known, SED-fitting methods reproduce the input distributions of LBGs' stellar masses relatively well, with a minor tendency to underestimate the masses overall, but with substantial scatter. Second, there are large systematic biases in the distributions of best-fit SFRs and mean ages, in the sense that single-component SED-fitting methods underestimate SFRs and overestimate ages. We attribute these trends to the different star formation histories predicted by the semi-analytic models and assumed in the galaxy templates used in SED-fitting procedure, and to the fact that light from the current generation of star-formation can hide older generations of stars. These biases, which arise from the SED-fitting procedure, can significantly affect inferences about galaxy evolution from broadband photometry.

We present the serendipitous discovery of z=4.05 molecular gas CO emission lines with the IRAM Plateau de Bure interferometer coincident with GN20 and GN20.2, two luminous submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey North field (GOODS-N). These are among the most distant submillimeter-selected galaxies reliably identified through CO emission and also some of the most luminous known. In terms of CO to bolometric luminosities, stellar mass and star formation rates (SFRs), these newly discovered z>4 SMGs are similar to z~1.5-3 SMGs studied to date. These z~4 SMGs have much higher specific SFRs than typical B-band dropout Lyman break galaxies at the same redshift. The stellar mass-SFR correlation for normal galaxies does not seem to evolve much further, between z~2 and z~4. A significant z=4.05 spectroscopic redshift spike is observed in GOODS-N, and a strong spatial overdensity of B-band dropouts and IRAC selected z>3.5 galaxies appears to be centered on the GN20 and GN20.2 galaxies. This suggests a proto-cluster structure with total mass ~10^14 Msun. Using photometry at mid-IR, submm and radio wavelengths, we show that reliable photometric redshifts (Dz/(1+z)~0.1) can be derived for SMGs over 1<z<4. This new photometric redshift technique has been used to provide a first estimate of the space density of 3.5<z<6 hyper-luminous starburst galaxies, and to show that they contribute substantially to the SFR density at early epochs. Many of these high-redshift starbursts will be within reach of Herschel. We find that the radio to mid-IR flux density ratio can be used to select z>3.5 starbursts, regardless of their submm/mm emission [abridged].

We present an analysis of the radio properties of large samples of Lyman Break Galaxies (LBGs) at $z \sim 3$, 4, and 5 from the COSMOS field. The median stacking analysis yields a statistical detection of the $z \sim 3$ LBGs (U-band drop-outs), with a 1.4 GHz flux density of $0.90 \pm 0.21 \mu$Jy. The stacked emission is unresolved, with a size $< 1"$, or a physical size $< 8$kpc. The total star formation rate implied by this radio luminosity is $31\pm 7$ $M_\odot$ year$^{-1}$, based on the radio-FIR correlation in low redshift star forming galaxies. The star formation rate derived from a similar analysis of the UV luminosities is 17 $M_\odot$ year$^{-1}$, without any correction for UV dust attenuation. The simplest conclusion is that the dust attenuation factor is 1.8 at UV wavelengths. However, this factor is considerably smaller than the standard attenuation factor $\sim 5$, normally assumed for LBGs. We discuss potential reasons for this discrepancy, including the possibility that the dust attenuation factor at $z \ge 3$ is smaller than at lower redshifts. Conversely, the radio luminosity for a given star formation rate may be systematically lower at very high redshift. Two possible causes for a suppressed radio luminosity are: (i) increased inverse Compton cooling of the relativistic electron population due to scattering off the increasing CMB at high redshift, or (ii) cosmic ray diffusion from systematically smaller galaxies. The radio detections of individual sources are consistent with a radio-loud AGN fraction of 0.3%. One source is identified as a very dusty, extreme starburst galaxy (a 'submm galaxy').

High-redshift submillimetre-bright galaxies identified by blank field surveys at millimetre and submillimetre wavelengths appear in the region of the IRAC colour-colour diagrams previously identified as the domain of luminous active galactic nuclei (AGNs). Our analysis using a set of empirical and theoretical dusty starburst spectral energy distribution (SED) models shows that power-law continuum sources associated with hot dust heated by young (<100 Myr old), extreme starbursts at z>2 also occupy the same general area as AGNs in the IRAC colour-colour plots. A detailed comparison of the IRAC colours and SEDs demonstrates that the two populations are distinct from each other, with submillimetre-bright galaxies having a systematically flatter IRAC spectrum (>1 mag bluer in the observed [4.5]-[8.0] colour). Only about 20% of the objects overlap in the colour-colour plots, and this low fraction suggests that submillimetre galaxies powered by a dust-obscured AGN are not common. The red IR colours of the submillimetre galaxies are distinct from those of the ubiquitous foreground IRAC sources, and we propose a set of IR colour selection criteria for identifying SMG counterparts that can be used even in the absence of radio or Spitzer MIPS 24 micron data.

We report the spectroscopic confirmation of a sub-mm galaxy (SMG) at z=4.547 with an estimated L_IR=0.5-2.0x10^13 L_sun. The spectra, mid-IR, and X-ray properties indicate the bolometric luminosity is dominated by star formation at a rate of >1000M_sun per yr. Multiple, spatially separated components are visible in the Ly-Alpha line with an observed velocity difference of up to 380 km/sec and the object morphology indicates a merger. The best fit spectral energy distribution and spectral line indicators suggest the object is 2-8 Myr old and contains >10^10 M_sun of stellar mass. This object is a likely progenitor for the massive early type systems seen at z~2.

We present the first results of the VIsible Multiobject Spectrograph (VIMOS) ESO/GOODS program of spectroscopy of faint galaxies in the Chandra Deep Field South (CDF-S). The program complements the FORS2 ESO/GOODS campaign. 3312 spectra have been obtained in service mode with VIMOS at the ESO/VLT UT3. The VIMOS LR-Blue and MR grisms have been used to cover different redshift ranges. Galaxies at 1.8 < z < 3.5 have been observed in the GOODS VIMOS-LR-Blue campaign. Galaxies at z < 1 and Lyman Break Galaxies at z > 3.5 have been observed in the VIMOS MR survey. Here we report results for the first 6 masks (out of 10 total) that have been analyzed from each of the LR-Blue and MR grisms. Spectra of 2344 spectra have been extracted from these 6 LR-Blue masks and 968 from 6 MR masks. 33% of the LR-Blue and 18% of the MR spectra are serendipitous observations. We obtained 1481 redshifts in the LR-Blue campaign and 656 in the MR campaign for a total success rate of 63% and 68%, respectively, which increase to 70% and 75% when only the primary targets are considered. By complementing our VIMOS spectroscopic catalog with all existing spectroscopic redshifts publicly available in the CDF-S, we created a redshift master catalog. By comparing this redshift compilation with different photometric redshift catalogs we estimate the completeness level of the CDF-S spectroscopic coverage in several redshift bins. The completeness level is very high, > 60%, at z < 3.5, and it is very uncertain at higher redshift. The master catalog has been used also to estimate completeness and contamination levels of different galaxy photometric selection techniques, such as the BzK, the so called 'sub'-U-dropout and the drop-out methods and to identify large scale structures in the field.

Aims. We present the full data set of the spectroscopic campaign of the ESO/GOODS program in the GOODS-South field, obtained with the FORS2 spectrograph at the ESO/VLT. Method. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams, and comparing the results with public data. Results. Including the third part of the spectroscopic campaign (12 masks) to the previous work (26 masks, Vanzella et al. 2005, 2006), 1715 spectra of 1225 individual targets have been analyzed. The actual spectroscopic catalog provides 887 redshift determinations. The typical redshift uncertainty is estimated to be sigma(z) ~ 0.001. Galaxies have been selected adopting different color criteria and using photometric redshifts. The resulting redshift distribution typically spans two domains: from z=0.5 to 2 and z=3.5 to 6.3. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/

We report results from a search for massive and evolved galaxies at z>5 in the GOODS southern field. Combining HST ACS, VLT ISAAC and Spitzer IRAC photometric data, we develop a color selection technique to identify candidates for being evolved galaxies at high redshifts. The color selection is primarily based on locating the Balmer-break using the K- and 3.6micron bands. Stellar population synthesis models are fitted to the SEDs of these galaxies to identify the final sample. We find 11 candidates with photometric redshifts in the range 4.9 < z < 6.5, dominated by an old stellar population, with ages 0.2-1.0 Gyr, and stellar masses in the range (0.5 - 5) 10^{11} Msun. The majority of the stars in these galaxies were formed at z > 9. One candidate has a spectroscopically confirmed redshift, in good agreement with our photometric redshift. The galaxies are very compact, with half-light radii in the observed K-band smaller than ~2 kpc. Seven of the 11 candidates are also detected at 24micron with the MIPS instrument on Spitzer. The 24micron emission could be interpreted as PAH emission from a dusty starburst at z~2-3, however, it is also consistent with the presence of an obscured AGN at z>5. We estimate the completeness of the Balmer break galaxy sample to be ~40%. The comoving number density of galaxies with a stellar mass >10^{11} Msun, at an average redshift z=5.2, is 3.9 10^{-5} Mpc^{-3} (no-MIPS sample: 1.4 10^{-5} Mpc^{-3}). The corresponding stellar mass density is 8 10^{6} Msun/Mpc^3 (no-MIPS sample: 6.2 10^6 Msun/Mpc^3).

We present the first spatial clustering measurements of z~1, 24um-selected, star forming galaxies in the Great Observatories Origins Deep Survey (GOODS). The sample under investigation includes 495 objects in GOODS-South and 811 objects in GOODS-North selected down to flux densities of f_24>20 uJy and z_AB<23.5 mag, for which spectroscopic redshifts are available. The median redshift, IR luminosity and star formation rate (SFR) of the samples are z~0.8, L_IR~4.4 x 10^10 L_sun, and SFR~7.6 M_sun/yr, respectively. We measure the projected correlation function w(r_p) on scales of r_p=0.06-10 h^-1 Mpc, from which we derive a best fit comoving correlation length of r_0 = 4.0 +- 0.4 h^-1 Mpc and slope of gamma=1.5 +- 0.1 for the whole f_24>20uJy sample after combining the two fields. We find indications of a larger correlation length for objects of higher luminosity, with Luminous Infrared Galaxies (LIRGs, L_IR>10^11 L_sun) reaching r_0~5.1 h^-1 Mpc. This would imply that galaxies with larger SFRs are hosted in progressively more massive halos, reaching minimum halo masses of ~3 x 10^12 M_sun for LIRGs. We compare our measurements with the predictions from semi-analytic models based on the Millennium simulation. The variance in the models is used to estimate the errors in our GOODS clustering measurements, which are dominated by cosmic variance. The measurements from the two GOODS fields are found to be consistent within the errors. On scales of the GOODS fields, the real sources appear more strongly clustered than objects in the Millennium-simulation based catalogs, if the selection function is applied consistently. This suggests that star formation at z~0.5-1 is being hosted in more massive halos and denser environments than currently predicted by galaxy formation models.[truncated]

To derive a new H$\alpha$ luminosity function and to understand the clustering properties of star-forming galaxies at $z \approx 0.24$, we have made a narrow-band imaging survey for H$\alpha$ emitting galaxies in the HST COSMOS 2 square degree field. We used the narrow-band filter NB816 ($\lambda_c = 8150$ \AA, $\Delta \lambda = 120$ \AA) and sampled H$\alpha$ emitters with $EW_{\rm obs}(\rm H\alpha + [N\textsc{ii}]) > 12$ \AA in a redshift range between $z=0.233$ and $z=0.251$ corresponding to a depth of 70 Mpc. We obtained 980 H$\alpha$ emitting galaxies in a sky area of 5540 arcmin$^2$, corresponding to a survey volume of $3.1 \times 10^4 {\rm Mpc^3}$. We derive a H$\alpha$ luminosity function with a best-fit Schechter function parameter set of $\alpha = -1.35^{+0.11}_{-0.13}$, $\log\phi_* = -2.65^{+0.27}_{-0.38}$, and $\log L_* ({\rm erg s^{-1}}) = 41.94^{+0.38}_{-0.23}$. The H$\alpha$ luminosity density is $2.7^{+0.7}_{-0.6} \times 10^{39}$ ergs s$^{-1}$ Mpc$^{-3}$. After subtracting the AGN contribution (15 %) to the H$\alpha$ luminosity density, the star formation rate density is evaluated as $1.8^{+0.7}_{-0.4} \times 10^{-2}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. The angular two-point correlation function of H$\alpha$ emitting galaxies of $\log L({\rm H\alpha}) > 39.8$ is well fit by a power law form of $w(\theta) = 0.013^{+0.002}_{-0.001} \theta^{-0.88 \pm 0.03}$, corresponding to the correlation function of $\xi(r) = (r/1.9{\rm Mpc})^{-1.88}$. We also find that the H$\alpha$ emitters with higher H$\alpha$ luminosity are more strongly clustered than those with lower luminosity.

Examining a sample of massive galaxies at 1.4<z<2.5 with K_{Vega}<22 from the Great Observatories Origins Deep Survey, we compare photometry from Spitzer at mid- and far-IR, to submillimeter, radio and rest-frame ultraviolet wavelengths, to test the agreement between different tracers of star formation rates (SFRs) and to explore the implications for galaxy assembly. For z~2 galaxies with moderate luminosities(L_{8um}<10^{11}L_sun), we find that the SFR can be estimated consistently from the multiwavelength data based on local luminosity correlations. However,20--30% of massive galaxies, and nearly all those with L_{8um}>10^{11}L_sun, show a mid-IR excess which is likely due to the presence of obscured active nuclei, as shown in a companion paper. There is a tight and roughly linear correlation between stellar mass and SFR for 24um-detected galaxies. For a given mass, the SFR at z=2 was larger by a factor of ~4 and ~30 relative to that in star forming galaxies at z=1 and z=0, respectively. Typical ultraluminous infrared galaxies (ULIRGs) at z=2 are relatively 'transparent' to ultraviolet light, and their activity is long lived (~400 Myr), unlike that in local ULIRGs and high redshift submillimeter-selected galaxies. ULIRGs are the common mode of star formation in massive galaxies at z=2, and the high duty cycle suggests that major mergers are not the dominant trigger for this activity.Current galaxy formation models underpredict the normalization of the mass-SFR correlation by about a factor of 4, and the space density of ULIRGs by an orderof magnitude, but give better agreement for z>1.4 quiescent galaxies.

Approximately 20-30% of 1.4<z<2.5 galaxies with K<22 (Vega) detected with Spitzer MIPS at 24um show excess mid-IR emission relative to that expected based on the rates of star formation measured from other multiwavelength data.These galaxies also display some near-IR excess in Spitzer IRAC data, with a spectral energy distribution peaking longward of 1.6um in the rest frame, indicating the presence of warm-dust emission usually absent in star forming galaxies. Stacking Chandra data for the mid-IR excess galaxies yields a significant hard X-ray detection at rest-frame energies >6.2 keV. The stacked X-ray spectrum rises steeply at >10 keV, suggesting that these sources host Compton-thick Active Galactic Nuclei (AGNs) with column densities N_H~10^{24} cm^-2 and an average, unobscured X-ray luminosity L_{2-8 keV}~(1-4)x10^43 erg/s. Their sky density(~3200 deg^-2) and space density (~2.6x10^-4 Mpc^-3) are twice those of X-ray detected AGNs at z~2, and much larger than those of previously-known Compton thick sources at similar redshifts. The mid-IR excess galaxies are part of the long sought-after population of distant heavily obscured AGNs predicted by synthesis models of the X-ray background. The fraction of mid-IR excess objects increases with galaxy mass, reaching ~50-60% for M~10^11 M_sun, an effect likely connected with downsizing in galaxy formation. The ratio of theinferred black hole growth rate from these Compton-thick sources to the global star formation rate at z=2 is similar to the mass ratio of black holes to stars inlocal spheroids, implying concurrent growth of both within the precursors oftoday's massive galaxies.

We examine deep far-ultraviolet (1600 Angstrom) imaging of the Hubble Deep Field-North (HDFN) and the Hubble Ultra Deep Field (HUDF) to search for leaking Lyman continuum radiation from starburst galaxies at z~1.3. There are 21 (primarily sub-L*) galaxies with spectroscopic redshifts between 1.1<z<1.5 and none are detected in the far-UV. We fit stellar population templates to the galaxies' optical/near-infrared SEDs to determine the starburst age and level of dust attenuation, giving an accurate estimate of the intrinsic Lyman continuum ratio, f_1500/f_700, and allowing a conversion from f_700 limits to relative escape fractions. We show that previous high-redshift studies may have underestimated the amplitude of the Lyman Break, and thus the relative escape fraction, by a factor of ~2. Once the starburst age and intergalactic HI absorption are accounted for, 18 galaxies in our sample have limits to the relative escape fraction, f_esc,rel < 1.0 with some limits as low as f_esc,rel < 0.10 and a stacked limit of f_esc,rel < 0.08. This demonstrates, for the first time, that most sub-L* galaxies at high redshift do not have large escape fractions. When combined with a similar study of more luminous galaxies at the same redshift we show that, if all star-forming galaxies at z~1 have similar relative escape fractions, the value must be less than 0.14 (3 sigma). We also show that less than 20% (3 sigma) of star-forming galaxies at z~1 have relative escape fractions near unity. These limits contrast with the large escape fractions found at z~3 and suggest that the average escape fraction has decreased between z~3 and z~1. (Abridged)

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

We describe the details of the Hubble Space Telescope (HST) Advanced Camera for Surveys / Wide Field Channel (ACS/WFC) observations of the COSMOS field, including the data calibration and processing procedures. We obtained a total of 583 orbits of HST ACS/WFC imaging in the F814W filter, covering a field that is 1.64 square degrees in area, the largest contiguous field ever imaged with HST. The median exposure depth across the field is 2028 seconds (one HST orbit), achieving a limiting point-source depth AB(F814W) = 27.2 (5 sigma). We also present details about the astrometric image registration, distortion removal and image combination using MultiDrizzle, as well as motivating the choice of our final pixel scale (30 milliarcseconds per pixel), based on the requirements for weak lensing science. The final set of images are publicly available through the archive sites at IPAC and STScI, along with further documentation on how they were produced.

We have carried out a wide-field imaging survey for [OII]3727 emitting galaxies at z~1.2 in the HST COSMOS 2 square degree field using the Suprime-Cam on the Subaru Telescope. The survey covers a sky area of 6700 arcmin^2 in the COSMOS field, and a redshift range between 1.17 and 1.20 (Delta_z = 0.03), corresponding to a survey volume of 5.56*10^5 Mpc^3. We obtain a sample of 3176 [OII] emitting galaxies with observed emission-line equivalent widths greater than 26 AA. Since our survey tends to sample brighter [OII]3727 emitting galaxies, we also analyze a sample of fainter [OII]3727 emitting galaxies found in the Subaru Deep Field (SDF). We find an extinction-corrected [OII] luminosity density of 10^{40.35^+0.08_-0.06} ergs s^-1 Mpc-3, corresponding to star formation rate density of 0.32^+0.06_-0.04 M_sun yr-1 Mpc^-3 in the COSMOS field at z~1.2. This is the largest survey for [OII]3727 emitters beyond z=1 currently available.

We present the first identification of large-scale structures (LSS) at z $< 1.1$ in the Cosmic Evolution Survey (COSMOS). The structures are identified from adaptive smoothing of galaxy counts in the pseudo-3d space ($\alpha,\delta$,z) using the COSMOS photometric redshift catalog. The technique is tested on a simulation including galaxies distributed in model clusters and a field galaxy population -- recovering structures on all scales from 1 to 20\arcmin without {\it a priori} assumptions for the structure size or density profile. Our procedure makes {\bf no} {\it a priori} selection on galaxy spectral energy distribution (SED, for example the Red Sequence), enabling an unbiased investigation of environmental effects on galaxy evolution. The COSMOS photometric redshift catalog yields a sample of $1.5\times10^5$ galaxies with redshift accuracy, $\Delta z_{FWHM}/(1+z) \leq 0.1$ at z $< 1.1$ down to I$_{AB} \leq 25$ mag. Using this sample of galaxies, we identify 42 large-scale structures and clusters. abstract truncated for astroph 25 line limit -- see preprint

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

We present results from a narrow-band optical survey of a contiguous area of 1.95 deg^2, covered by the Cosmic Evolution Survey (COSMOS). Both optical narrow-band (lambda_c = 8150 AA and Delta_lambda = 120 AA) and broad-band (B, V, g', r', i', and z') imaging observations were performed with the Subaru prime-focus camera, Suprime-Cam on the Subaru Telescope. We provide the largest contiguous narrow-band survey, targetting Ly alpha emitters (LAEs) at z~5.7. We find a total of 119 LAE candidates at z~5.7. Over the wide-area covered by this survey, we find no strong evidence for large scale clustering of LAEs. We estimate a star formation rate (SFR) density of ~7*10^-4 M_sun yr^-1 Mpc^-3 for LAEs at z~5.7, and compare it with previous measurements.

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

We report on the discovery of a new potential galaxy threshing system in the COSMOS 2 square degree field using the prime-focus camera, Suprime-Cam, on the 8.2 m Subaru Telescope. This system consists of a giant elliptical galaxy with $M_V \approx -21.6$ and a tidally disrupted satellite galaxy with $M_V \approx -17.7$ at a photometric redshift of $z \approx 0.08$. This redshift is consistent with the spectroscopic redshift of 0.079 for the giant elliptical galaxy obtained from the Sloan Digital Sky Survey (SDSS) archive. The luminosity masses of the two galaxies are $3.7 \times 10^{12} \cal{M}_{\odot}$ and $3.1 \times 10^{9} \cal{M}_{\odot}$, respectively. The distance between the two galaxies is greater than 100 kpc. The two tidal tails emanating from the satellite galaxy extend over 150 kpc. This system would be the second well-defined galaxy threshing system found so far.

We measure photometric redshifts and spectral types for galaxies in the COSMOS survey. We use template fitting technique combined with luminosity function priors and with the option to simultaneously estimate dust extinction (i.e. E(B-V)) for each galaxy.Our estimated redshifts are accurate to i<25 and z~1.2. Using simulations with sampling and noise characteristics similar to those in COSMOS, the accuracy and reliability is estimated for the photometric redshifts as a function of the magnitude limits of the sample, S/N ratios and the number of bands used. From the simulations we find that the ratio of derived 95% confidence interval in the redshift probability distribution to the estimated photometric redshift (D95) can be used to identify and exclude the catastrophic failures in the photometric redshift estimates. We compare the derived redshifts with high-reliability spectroscopic redshifts for a sample of 868 normal galaxies with z < 1.2 from zCOSMOS. Considering different scenarios, depending on using prior, no prior and/or extinction, we compare the photometric and spectroscopic redshifts for this sample. This corresponds to an rms scatter of 0.031, with a small number of outliers (<2.5%). We also find good agreement (rms=0.10) between photometric and spectroscopic redshifts for Type II AGNs. We compare results from our photometric redshift procedure with three other independent codes and find them in excellent agreement. We show preliminary results, based on photometric redshifts for the entire COSMOS sample (to i < 25 mag.).

The Cosmic Evolution Survey (COSMOS) is designed to probe the correlated evolution of galaxies, star formation, active galactic nuclei (AGN) and dark matter (DM) with large-scale structure (LSS) over the redshift range z $> 0.5 $ to 6. The survey includes multi-wavelength imaging and spectroscopy from X-ray to radio wavelengths covering a 2 $\sq$\deg area, including HST imaging. Given the very high sensitivity and resolution of these datasets, COSMOS also provides unprecedented samples of objects at high redshift with greatly reduced cosmic variance, compared to earlier surveys. Here we provide a brief overview of the survey strategy, the characteristics of the major COSMOS datasets, and summarize the science goals.

We present deep optical imaging observations of 2 square degree area, covered by the Cosmic Evolution Survey (COSMOS), made by the prime-focus Camera (Supreme-Cam) on the 8.2m Subaru Telescope. Observations were done in six broad-band [B (4459.7 AA), g' (4723.1 AA), V (5483.8 AA), r' (6213.0 AA), i' (7640.8 AA), z' (8855.0 AA)], and one narrow-band (NB816) filters. A total of 10^6 galaxies were detected to i'~26.5 mag. These data, combined with observations at u* and K-band are used to construct the photometric catalogs for the COSMOS and to measure their photometric redshifts, multi-band spectral energy distributions, stellar masses and identification of high redshift candidates. This catalog provides multi-waveband data for scientific analysis of the COSMOS survey.

We present initial results from our ongoing campaign of spectroscopic identifications of Lyman--break galaxies (LBGs) at z ~ 4, 5 and 6 with FORS2 at the ESO VLT. The sample of LBGs is selected from the HST/ACS images in the GOODS southern field (CDF--S), and here we discuss the spectroscopic properties of 103 LBGs that we have cofirmed in the redshift interval 3.5 < z < 6.3, of which 34 at z>5 and 27 between 5.5 and 6.3. As for LBGs at lower redshifts, we find that LBGs with LyA in emission are characterized by weaker intestellar absorption lines and bluer UV continua than galaxies whose LyA is observed in absorption. No correlation between the LyA properties and the apparent elongation of the UV morphology of the galaxies is observed, suggesting that geometrical effects are not the cause of the spectral differences. A trend of increasing LyA equivalent width with decreasing luminosity is also observed. We also present new measures of the UV luinosity function of LBGs.

We describe the angular power spectrum of resolved sources at 3.6 microns (L-band) in Spitzer imaging data of the GOODS HDF-N, the GOODS CDF-S, and the NDWFS Bootes field in several source magnitude bins. We also measure angular power spectra of resolved sources in the Bootes field at K_S and J-bands using ground-based IR imaging data. In the three bands, J, K_S, and L, we detect the clustering of galaxies on top of the shot-noise power spectrum at multipoles between ell ~ 10^2 and 10^5. The angular power spectra range from the large, linear scales to small, non-linear scales of galaxy clustering, and in some magnitude ranges, show departure from a power-law clustering spectrum. We consider a halo model to describe clustering measurements and to establish the halo occup ation number parameters of IR bright galaxies at redshifts around one. We also extend our clustering results and completeness-corrected faint source number counts in GOODS fields to understand the underlying nature of unresolved sources responsible for IR background (IRB) anisotropies that were detected in deep Spitzer images. While these unresolved fluctuations were measured at sub-arcminute angular scales, if a high-redshift diffuse component associated with first galaxies exists in the IRB, then it's clustering properties are best studied with shallow, wide-field images that allow a measurement of the clustering spectrum from a few degrees to arcminute angular scales.

Using the multiwavelength photometric and spectroscopic data covering the CDF-S obtained within the GOODS, we investigate the rest-frame UV properties of galaxies to z~2.2, including the evolution of the luminosity function, the luminosity density, star formation rate (SFR) and galaxy morphology. We find a significant brightening (~ 1 mag) in the rest-frame 2800A characteristic magnitude (M*) over the redshift range 0.3<z<1.7 and no evolution at higher redshifts. The rest-frame 2800A luminosity density shows an increase by a factor ~4 to z~2.2. We estimate the SFR to z~2.2 from the 1500A and 2800A luminosities. The SFR derived from the 2800A luminosity density is almost factor two higher than that derived from the 1500A luminosities. Attributing this to differential dust extinction, we find that an E(B-V)=0.20 results in the same extinction corrected SFR from both 1500A and 2800A luminosities. We investigate the morphological composition of our sample by fitting Sersic profiles to the galaxy images at a fixed rest-frame wavelength of 2800A at 0.5<z<2.2. We find that the fraction of apparently bulge-dominated galaxies (Sersic index n>2.5) increases from ~10% at z~0.5 to ~30% at z~2.2. At the same time, we note that galaxies get bluer at increasing redshift. This suggests a scenario where an increased fraction of the star formation takes place in bulge-dominated systems at high redshift. This could be the evidence that the present day ellipticals are a result of assembly (i.e., mergers) of galaxies at z>1. Finally, we find that galaxy sizes for a luminosity selected sample evolves as (1+z)^{-1.1} between redshifts z=2.2 and z=1.1. This is consistent with previous measurements and suggests a similar evolution over the redshift range 0 < z < 6.

We have used the HST/ACS images to identify 4700 Lyman break galaxies (LBGs) in GOODS. We present the results from a parametric analysis of the 2-D surface brightness profiles, for 1333 LBGs at z > 2.5 with rest-frame UV(1600 Angstrom) AB magnitude < -20.5. Based on the Sersic index, n, which measures the profile shape, we find that about 40% of LBGs at z=3 have light profiles close to exponential, and only 30% have the high concentrations seen for spheroids. About 30% of LBGs appear to have multiple cores or disturbed morphologies suggestive of close pairs or on-going mergers. The fraction of spheroid-like (n > 2.5) LBGs decreases by about 15% from z = 5 to 3. A comparison of LBGs with the starburst galaxies at z = 1.2, shows that the fraction of spheroid-like profiles is about 20% higher among LBGs. The ellipticity distribution for LBGs exhibits a pronounced skew towards high ellipticities (> 0.5), which cannot be explained by morphologies similar to the local disks and spheroids viewed at random orientations. The peak of the distribution evolves toward lower ellipticities, from 0.7 at z = 4 to 0.5 at z = 3. At z = 1.2 the distribution is relatively flat as seen among the present-day galaxies. The dominance of elongated morphologies among LBGs suggests that in a significant fraction of them we may be witnessing star-formation in clumps along gas-rich filaments, or the earliest gas-rich bars that encompass essentially the entire visible galaxy. Similar features are found to be ubiquitous in hydrodynamical simulations in which galaxy formation at high redshifts occurs in filamentary inflows of dynamically cold gas within the dark matter halos, and involves gas- rich mergers.

Using the deep Spitzer IRAC observations of the Great Observatories Origins Deep Survey (GOODS), we study the stellar masses and star formation histories of galaxies at z ~ 6. Our study is based on the i-band dropout sample selected from the GOODS southern and northern fields (~ 330 arcmin^2 in total), several of which already have spectroscopic confirmations. In total, we derive stellar masses for 53 i-band dropouts that have robust IRAC detections. These galaxies have typical stellar masses of ~ 1e10 M_sun and typical ages of a couple of hundred million years, consistent with earlier results based on a smaller sample of z ~ 6 galaxies in the HUDF. We also study 79 i-band dropouts that are invisible in the IRAC data and find that they are typically less massive by a factor of ten. Based on our mass estimates, we derive a lower limit to the global stellar mass density at z ~ 6. Considering the range of systematic uncertainties in the derived stellar masses, this lower limit is 1.1--6.7e6 M_sun/Mpc^3 (co-moving), which is 0.2--1.1% of the present-day value. The prospect of detecting the progenitors of the most massive galaxies at yet higher redshifts is explored: a deep, wide-field near-IR survey using our current technology could possibly result in positive detections at z > 7. We also investigate the implication of our results for reionzation, and find that the progenitors of the galaxies comparable to those in our sample, even in the most optimized (probably unrealistic) scenario, cannot sustain the reionization for a period longer than ~ 2 million years. Thus most of the photons required for reionization must have been provided by other sources, such as the progenitors of the dwarf galaxies that are far below our current detection capability. (Abridged)

We present the second campaign of the ESO/GOODS program of spectroscopy of faint galaxies in the GOODS-South field. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams and comparing the results with public data. 807 spectra of 652 individua targets have been obtained in service mode with the FORS2 spectrograph at the ESO/VLT, providing 501 redshift determinations. The typical redshift uncertainty is estimated to be sigma_z ~ 0.001. Galaxies have been selected adopting three different color criteria and using the photometric redshifts.The resulting redshift distribution typically spans two redshift domains: from z=0.5 to 2 and z=3.5 to 6.2. In particular, 94 B435-,V606-,i775 "dropout" Lyman break galaxies have been observed, yielding redshifts for 65 objects in the interval 3.4<z<6.2. Three sources have been serendipitously discovered in the redshift interval 4.8<z<5.5. Together with the previous release, 930 sources have now been observed and 724 redshift determinations have been carried out. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/. Large scale structures are clearly detected at z=0.666, 0.734, 1.096, 1.221, 1.300, and 1.614. A sample of 34 sources with tilted [OII]3727 emission has been identified, 32 of them in the redshift range 0.9<z<1.5.

We present a compilation of multiwavelength data on different star formation indicators for a sample of nearby star forming galaxies. Here we discuss the observations, reductions and measurements of ultraviolet images obtained with STIS, on board the Hubble Space Telescope, ground-based Halpha, and VLA 8.46 GHz radio images. These observations are complemented with infrared fluxes, as well as large aperture optical radio and ultraviolet data from the literature. This database will be used in a forthcoming paper to compare star formation rates at different wavebands. We also present spectral energy distributions (SEDs) for those galaxies with at least one far-infrared measurements from ISO, longward of 100 um. These SEDs are divided in two groups, those which are dominated by the far-infrared emission, and those where the contribution from the far-infrared and optical emission is comparable. These SEDs are useful tools to study the properties of high redshift galaxies.

We present the results of a multiwavelength study of nearby galaxies, aimed at understanding the relation between the ultraviolet and far-infrared emission in star-forming galaxies. The dataset comprises new ultraviolet (from HST/STIS), ground-based Halpha, and radio continuum observations, together with archival infrared data (from IRAS and ISO). The local galaxies are used as benchmarks for comparison of the infrared-to-ultraviolet properties with two populations of high-redshift galaxies: the sub-millimeter star-forming galaxies detected by SCUBA and the ultraviolet-selected Lyman Break galaxies. In addition, the long wavelength baseline covered by the present data enables us to compare the star formation rates (SFRs) derived from the observed ultraviolet, Halpha, infrared, and radio luminosities, and to gauge the impact of dust opacity in the local galaxies. We also derive a new calibration for the non-thermal part of the radio SFR estimator, based on the comparison of 1.4 GHz measurements with a new estimator of the bolometric luminosity of the star forming regions. We find that more actively star forming galaxies show higher dust opacities, in line with previous results. We find that the local star-forming galaxies have lower F(205um)/F(UV) ratio, by two-three orders of magnitude than the submillimeter-selected galaxies, and may have similar or somewhat higher F(205um)/F(UV) than Lyman Break Galaxies. The F(205um)/F(UV) ratio of the local galaxy population may be influenced by the cool dust emission in the far-infrared heated by non-ionizing stellar populations, which may be reduced or absent in the LBGs.

We report on the angular correlation function of Lyman-break galaxies (LBGs) at z~4 and 5 from deep samples obtained from the Great Observatories Deep Origins Survey (GOODS). Similar to LBGs at z~3, the shape of w(theta) of the GOODS LBGs is well approximated by a power-law with slope beta~0.6 at angular separation theta > 10 arcsec. The clustering strength of z~4, 5 LBGs also depends on the rest-frame UV luminosity, with brighter galaxies more strongly clustered than fainter ones, implying a general correlation between halos' mass and LBGs' star-formation rate. At smaller separations, w(theta) of deep samples significantly exceeds the extrapolation of the large-scale power-law fit, implying enhanced spatial clustering at scales r < 1 Mpc. We also find that bright LBGs statistically have more faint companions on scales theta < 20 arcsec than fainter ones, showing that the enhanced small-scale clustering is very likely due to sub-structure, namely the fact that massive halos can host multiple galaxies. A simple model for the halo occupation distribution and the CDM halo mass function reproduce well the observed w(theta). The scaling relationship of the clustering strength with volume density and with redshift is quantitatively consistent with that of CDM halos. A comparison of the clustering strength of three samples of equal luminosity limit at z ~ 3, 4 and 5 shows that the LBGs at z~5 are hosted in halos about one order of magnitude less massive than those in the lower redshift bins, suggesting that star-formation was more efficient at higher-redshift.

We investigate the properties of massive galaxies at z=1-3.5 using HST observations, ground-based near-IR imaging, and Spitzer Space Telescope observations at 3-24 micron. We identify 153 distant red galaxies (DRGs) with J-K > 2.3 mag (Vega) in the southern GOODS field. This sample is approximately complete in stellar mass for passively evolving galaxies above 10^11 solar masses and z < 3. The galaxies identified by this selection are roughly split between objects whose optical and near-IR rest-frame light is dominated by evolved stars combined with ongoing star formation, and galaxies whose light is dominated by heavily reddened starbursts. Very few of the galaxies (< 10%) have no indication of current star formation. Using SFR estimates that include the reradiated IR emission, the DRGs at z=1.5-3 with stellar masses > 10^11 solar masses have specific SFRs (SFRs per unit stellar mass) ranging from 0.2 to 10 Gyr^-1, with a mean value of ~2.4 Gyr^-1. The DRGs with stellar masses > 10^11 solar masses and 1.5 < z < 3 have integrated specific SFRs greater the global value over all galaxies. In contrast, we find that galaxies at z = 0.3-0.75 with these stellar masses have integrated specific SFRs less than the global value, and more than an order of magnitude lower than that for massive DRGs at z = 1.5-3. At z < 1, lower-mass galaxies dominate the overall cosmic mass assembly. This suggests that the bulk of star formation in massive galaxies occurs at early cosmic epochs and is largely complete by z~1.5. [Abridged]

We present results from a search for high-redshift J--band ``dropout'' galaxies in the portion of the GOODS southern field that is covered by extremely deep imaging from the Hubble Ultradeep Field (HUDF).Using observations at optical, near-infrared and mid-infrared wavelengths from the Hubble and Spitzer Space Telescopes and the ESO-VLT, we search for very massive galaxies at high redshifts and find one particularly remarkable candidate. Its spectral energy distribution is consistent with a galaxy at z ~ 6.5 and a stellar mass of 6x10e11 M(sun) (for a Salpeter IMF). We interpret a prominent photometric break between the near-infrared and Spitzer bandpasses as the 3646A Balmer discontinuity. The best-fitting models have low reddening and ages of several hundred Myr, placing the formation of the bulk of the stars at z > 9. Alternative models of dusty galaxies at z ~ 2.5 are possible but provide significantly poorer fits. The object is detected with Spitzer at 24 micron. This emission originats from an obscured active nucleus or star formation. We present optical and near-infrared spectroscopy which has, thus far, failed to detect any spectral features. This helps limit the solution in which the galaxy is a starburst or active galaxy at z ~ 2.5. If the high-redshift interpretation is correct, this object would be an example of a galaxy that formed by a process strongly resembling traditional models of monolithic collapse, in a way which a very large mass of stars formed within a remarkably short period of time, at very high redshift.

We examine the possible reionization of the intergalactic medium (IGM) by the source UDF033238.7-274839.8 (hereafter HUDF-JD2), which was discovered in deep {\it HST}/VLT/{\it Spitzer} images obtained as part of the Great Observatory Origins Deep Survey and {\it Hubble} Ultra-Deep Field projects. Mobasher et al (2005) have identified HUDF-JD2 as a massive ($\sim6\times10^{11}M_\odot$) post-starburst galaxy at redshift z$\gtrsim6.5$. We find that HUDF-JD2 may be capable of reionizing its surrounding region of the Universe, starting the process at a redshift as high as z$\approx 15 \pm5$.

We apply a new approach to quantifying galaxy morphology and identifying galaxy mergers to the rest-frame far-ultraviolet images of 82 z ~ 4 Lyman break galaxies (LBGs) and 55 1.2 < z < 1.8 emission-line galaxies in the GOODS and Ultra Deep Fields. We compare the distributions of the Gini coefficient (G), second-order moment of the brightest 20% of galaxy light (M20), and concentration (C) for high-redshift and low-redshift galaxies with average signal to noise per pixel > 2.5 and Petrosian radii > 0.3 arcsec. Ten of the 82 LBGs have M20 >= -1.1 and possess bright double or multiple nuclei, implying a major-merger fraction of star-forming galaxies ~ 10-25% at M_{FUV} < -20, depending on our incompleteness corrections. Galaxies with bulge-like morphologies (G >= 0.55, M20 < -1.6) make up ~ 30% of the z ~ 4 LBG sample, while the remaining ~ 50% have G and M20 values higher than expected for smooth bulges and disks and may be star-forming disks, minor-mergers or post-mergers. The star-forming z ~ 1.5 galaxy sample has a morphological distribution which is similar to the UDF z ~ 4 LBGs, with an identical fraction of major-merger candidates but fewer spheroids. The observed morphological distributions are roughly consistent with current hierarchical model predictions for the major-merger rates and minor-merger induced starbursts at z ~ 1.5 and ~4. We also examine the rest-frame FUV-NUV and FUV-B colors as a function of morphology and find no strong correlations at either epoch.

We use data from the first epoch of observations with the IRAC/Spitzer for the GOODS to detect and study a collection of LBGs at z ~ 6 to 5 in the HUDF, six of which have spectroscopic confirmation. At these redshifts, IRAC samples rest-frame optical light in the range 0.5 to 0.8 um, where the effects of dust extinction are smaller and the sensitivity to light from evolved stars is greater than at shorter, rest-frame UV observable from the ground or with the HST. As such, it provides useful constraints on the ages and masses of these galaxies' stellar populations. We find that the SEDs for many of these galaxies are best fitted by models of stellar populations with masses of a few x 1e10 M_sun, and with ages of a few hundred Myrs, values quite similar to those derived for typical LBGs at z ~ 3. When the universe was only 1 Gyr old, some galaxies had already formed a mass of stars approaching that of the present-day Milky Way, and that they started forming those stars at z > 7, and in some cases much earlier. We find that the lower limits to the space density for galaxies in this mass range are consistent with predictions from recent hydrodynamic simulations of structure formation in a LCDM universe. All objects in our samples are consistent with having solar metallicity, suggesting that they might have already been significantly polluted by metals. The values for dust reddening derived from the model fitting are low or zero, and we find that some of the galaxies have rest-frame UV colors that are even bluer than those predicted by the stellar population models to which we compare them. These colors might be attributed to the presence of very massive stars (> 100 M_sun), or by weaker intergalactic HI absorption than what is commonly assumed.

We investigate the multi-wavelength emission of BzK selected star forming galaxies at z~2 in the Great Observatories Origins Deep Survey (GOODS) North region. Most (82%) of the sources are individually detected at 24um in the Spitzer MIPS imaging, and one fourth (26%) in the VLA radio data. Significant detections of the individually undetected objects are obtained through stacking in the radio, submm and X-ray domains. The typical star forming galaxy with stellar mass ~10^{11}Mo at z=2 is an Ultra-luminous Infrared Galaxy (ULIRG), with L_IR ~ 1-2x10^{12}Lo and star formation rate SFR 200-300Mo/yr, implying a comoving density of ULIRGs at z=2 at least 3 orders of magnitude above the local one. SFRs derived from the reddening corrected UV luminosities agree well, on average, with the longer wavelength estimates. The high 24um detection rate suggests a relatively large duty cycle for the BzK star forming phase, consistently with the available independent measurements of the space density of passively evolving galaxies at z>1.4. If the IMF at z=2 is similar to the local one, and in particular is not a top-heavy IMF, this suggests that a substantial fraction of the high mass tail (>10^{11}Mo) of the galaxy stellar mass function was completed by z~1.4.

Using Hubble Space Telescope and ground-based U through K- band photometry from the Great Observatories Origins Deep Survey (GOODS), we measure the evolution of the luminosity function and luminosity density in the rest-frame optical (UBR) to z ~ 2, bridging the poorly explored ``redshift desert'' between z~1 and z~2. We also use deep near-infrared observations to measure the evolution in the rest-frame J-band to z~1. Compared to local measurements from the SDSS, we find a brightening of the characteristic magnitude, (M*), by ~2.1, \~0.8 and ~0.7 mag between z=0.1 and z=1.9, in U, B, and R bands, respectively. The evolution of M* in the J-band is in the opposite sense, showing a dimming between redshifts z=0.4 and z=0.9. This is consistent with a scenario in which the mean star formation rate in galaxies was higher in the past, while the mean stellar mass was lower, in qualitative agreement with hierarchical galaxy formation models. We find that the shape of the luminosity function is strongly dependent on spectral type and that there is strong evolution with redshift in the relative contribution from the different spectral types to the luminosity density. We find good agreement in the luminosity function derived from an R-selected and a K-selected sample at z~1, suggesting that optically selected surveys of similar depth (R < 24) are not missing a significant fraction of objects at this redshift relative to a near-infrared-selected sample. We compare the rest-frame B-band luminosity functions from z~0--2 with the predictions of a semi-analytic hierarchical model of galaxy formation, and find qualitatively good agreement. In particular, the model predicts at least as many optically luminous galaxies at z~1--2 as are implied by our observations.

We have studied the evolution in the morphologies, sizes, stellar-masses, colors, and internal color dispersion (ICD) of galaxies at z=1 and 2.3, using a near-IR, flux-limited catalog for the HDF-N. At z=1 most luminous galaxies have morphologies of early-to-mid Hubble-types, and many show transformations between their rest-frame UV-optical morphologies. Galaxies at z=2.3 have compact and irregular morphologies with no clearly evident Hubble-sequence candidates. The mean galaxy size grows from z=2.3 to 1 by 40%, and the density of galaxies larger than 3 kpc increases by 7 times. At z=1, the size-luminosity distribution is broadly consistent with that of local galaxies, with passive evolution. However, galaxies at z=2.3 are smaller than the large present-day galaxies, and must continue to grow in size and stellar mass. We have measured the galaxies' UV-optical ICD, which quantifies differences in morphology and the relative amount of on-going star-formation. The mean and scatter in galaxies' total colors and ICD increase from z=2.3 to 1. At z=1 many galaxies with large ICD are spirals, with a few irregular systems. Few z=2.3 galaxies have high ICD, and those that do are actively merging. We interpret this as evidence for the presence of older and more diverse stellar populations at z=1 that are not generally present at z>2. We conclude that the star-formation histories of galaxies at z>2 are dominated by discrete, recurrent bursts, which quickly homogenize the galaxies' stellar content, and are possibly associated with mergers. The increase in the stellar-population diversification by z<1.4 implies that merger-induced starbursts occur less frequently than at higher redshifts, and more quiescent star-forming modes dominate. This transition coincides with the emergence of Hubble-sequence galaxies. [Abridged]

We constrain the X-ray emission properties of Lyman break galaxies (LBGs) at z ~ 3-6 using the 2 Ms Chandra Deep Field-North and 1 Ms Chandra Deep Field-South. Large samples of LBGs were discovered using HST as part of the Great Observatories Origins Deep Survey (GOODS). Deep optical and X-ray imaging over the GOODS fields have allowed us to place the most significant constraints on the X-ray properties of LBGs to date. Mean X-ray properties of 449, 1734, 629, and 247 LBGs with z ~ 3, 4, 5, and 6, respectively, were determined using stacking techniques. When stacked, we detect X-ray emission from LBGs at z ~ 3 (7 sigma) and from an optically bright subset (brightest 25%) of LBGs at z ~ 4 (3 sigma); the latter is the highest redshift detection yet for "normal" galaxies in the X-ray band. The effective exposure times for these stacked observations are 0.7 and 0.5 Gs, respectively. The derived average rest-frame 2.0-8.0 keV luminosities are 1.5 and 1.4 x 10^41 erg s^-1, respectively. X-ray emission from these LBGs is likely due to high mass X-ray binaries (HMXBs) and Type II supernovae; the corresponding star formation rates are 85-240 M_sol yr^-1. The X-ray to B-band mean luminosity ratio (Lx/LB) at z~3 is somewhat elevated with respect to that measured for starburst galaxies in the local Universe. When stacking full samples of LBGs at z ~ 4, 5, and 6 we do not obtain significant detections and derive rest-frame 2.0-8.0 keV luminosity upper limits of 0.9, 2.8, and 7.1 x 10^41 erg s^-1, respectively. These upper limits constrain any widespread AGN activity in these objects to be modest at best. Furthermore, we find that ~ 0.5% of our LBGs from z ~ 3-6 are detected individually in the X-ray band. Continued ...

We present the first results of the ESO/GOODS program of spectroscopy of faint galaxies in the Chandra Deep Field South (CDF-S). 399 spectra of 303 unique targets have been obtained in service mode with the FORS2 spectrograph at the ESO/VLT, providing 234 redshift determinations (the median of the redshift distribution is at 1.04). The typical redshift uncertainty is estimated to be sig(z) ~ 0.001. Galaxies have been color selected in a way that the resulting redshift distribution typically spans from z=0.5 to 2. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/ Large scale structure is clearly detected at z ~ 0.67, 0.73, 1.10 and 1.61. Three Lyman-break galaxies have also been included as targets and are confirmed to have redshifts z=4.800, 4.882 and 5.828. In a few cases, we observe clear [OII]3727 rotation curves, even at the relatively low resolution (R = 860) of the present observations. Assuming that the observed velocity structure is due to dynamically-relaxed rotation, this is an indication of large galactic masses (few times 10^(11) solar masses) at z ~ 1.

We present results from the Hubble Higher-z Supernova Search, the first space-based open field survey for supernovae (SNe). In cooperation with the Great Observatories Origins Deep Survey, we have used the Hubble Space Telescope with the Advanced Camera for Surveys to cover 300 square arcmin in the area of the Chandra Deep Field South and the Hubble Deep Field North on five separate search epochs (separated by 45 day intervals) to a limiting magnitude of z'=26. These deep observations have allowed us to discover 42 SNe in the redshift range 0.2 < z < 1.6. As these data span a large range in redshift, they are ideal for testing the validity of Type Ia supernova progenitor models with the distribution of expected ``delay times,'' from progenitor star formation to SN Ia explosion, and the SN rates these models predict. Through a Bayesian maximum likelihood test, we determine which delay-time models best reproduce the redshift distribution of SNe Ia discovered in this survey. We find that models that require a large fraction of ``prompt'' (less than 2 Gyr) SNe Ia poorly reproduce the observed redshift distribution and are rejected at 95% confidence. We find that Gaussian models best fit the observed data for mean delay times in the range of 3 to 4 Gyr.

We use a sample of 42 supernovae detected with the Advanced Camera for Surveys on-board the Hubble Space Telescope as part of the Great Observatories Origins Deep Survey to measure the rate of core collapse supernovae to z~0.7 and type Ia supernovae to z~1.6. This significantly increases the redshift range where supernova rates have been estimated from observations. The rate of core collapse supernovae can be used as an independent probe of the cosmic star formation rate. Based on the observations of 17 core collapse supernovae, we measure an increase in the core collapse supernova rate by a factor of 1.6 in the range 0.3<z<0.7, and an overall increase by a factor of 7 to z~0.7 in comparison to the local core collapse supernova rate. The increase in the rate in this redshift range in consistent with recent measurements of the star formation rate derived from UV-luminosity densities and IR datasets. Based on 25 type Ia supernovae, we find a SN Ia rate that is a factor 3-5 higher at z~1 compared to earlier estimates at lower redshifts (z<0.5), implying that the type Ia supernova rate traces a higher star formation rate at redshifts z>1 compared to low redshift. At higher redshift (z>1), we find a suggested decrease in the type Ia rate with redshift. This evolution of the Ia rate with redshift is consistent with a type Ia progenitor model where there is a substantial delay between the formation of the progenitor star and the explosion of the supernova. Assuming that the type Ia progenitor stars have initial main sequence masses 3-8 M_Sun, we find that 5-7% of the available progenitors explode as type Ia supernovae.

We have discovered 16 Type Ia supernovae (SNe Ia) with the Hubble Space Telescope (HST) and have used them to provide the first conclusive evidence for cosmic deceleration that preceded the current epoch of cosmic acceleration. These objects, discovered during the course of the GOODS ACS Treasury program, include 6 of the 7 highest-redshift SNe Ia known, all at z>1.25, and populate the Hubble diagram in unexplored territory. The luminosity distances to these and 170 previous SNe Ia are provided. A purely kinematic interpretation of the SN Ia sample provides evidence at the > 99% confidence level for a transition from deceleration to acceleration or similarly, strong evidence for a cosmic jerk. Using a simple model of the expansion history, the transition between the two epochs is constrained to be at z=0.46 +/- 0.13. The data are consistent with the cosmic concordance model of Omega_M ~ 0.3, Omega_Lambda~0.7 (chi^2_dof=1.06), and are inconsistent with a simple model of evolution or dust as an alternative to dark energy. For a flat Universe with a cosmological constant. When combined with external flat-Universe constraints we find w=-1.02 + 0.13 - 0.19 (and $<-0.76 at the 95% confidence level) for an assumed static equation of state of dark energy, P = w\rho c^2. Joint constraints on both the recent equation of state of dark energy, $w_0$, and its time evolution, dw/dz, are a factor of ~8 more precise than its first estimate and twice as precise as those without the SNe Ia discovered with HST. Our constraints are consistent with the static nature of and value of w expected for a cosmological constant (i.e., w_0 = -1.0, dw/dz = 0), and are inconsistent with very rapid evolution of dark energy. We address consequences of evolving dark energy for the fate of the Universe.

We examine the evolution of the sizes and number densities of disk galaxies using the high resolution images obtained by the Great Observatories Origins Deep Survey (GOODS) with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). The multiwavelength (B,V,i,z) images allow us to classify galaxies based on their rest-frame B-band morphologies out to redshift z = 1.25. In order to minimize the effect of selection biases, we confine our analysis to galaxies which occupy the region of magnitude-size plane where the survey is almost 90% complete at all redshifts. The observed size distribution is consistent with a log-normal distribution as seen for the disk galaxies in the local Universe and does not show any significant evolution over the redshift range 0.25 < z < 1.25. We find that the number densities of disk galaxies remains fairly constant over this redshift range, although a modest evolution by a factor of four may be possible within the 2-sigma uncertainties.

We have conducted a deep multi-color imaging survey of 0.2 degrees^2 centered on the Hubble Deep Field North (HDF-N). We shall refer to this region as the Hawaii-HDF-N. Deep data were collected in U, B, V, R, I, and z' bands over the central 0.2 degrees^2 and in HK' over a smaller region covering the Chandra Deep Field North (CDF-N). The data were reduced to have accurate relative photometry and astrometry across the entire field to facilitate photometric redshifts and spectroscopic followup. We have compiled a catalog of 48,858 objects in the central 0.2 degrees^2 detected at 5 sigma significance in a 3" aperture in either R or z' band. Number counts and color-magnitude diagrams are presented and shown to be consistent with previous observations. Using color selection we have measured the density of objects at 3<z<7. Our multi-color data indicates that samples selected at z>5.5 using the Lyman break technique suffer from more contamination by low redshift objects than suggested by previous studies.

The integrated colors of distant galaxies provide a means for interpreting the properties of their stellar content. Here, we use rest-frame UV-to-optical colors to constrain the spectral-energy distributions and stellar populations of color-selected, B-dropout galaxies at z ~ 4 in the Great Observatories Origins Deep Survey. We combine the ACS data with ground-based near-infrared images, which extend the coverage of galaxies at z ~ 4 to the rest-frame B-band. We observe a color-magnitude trend in the rest-frame m(UV) - B versus B diagram for the z ~ 4 galaxies that has a fairly well-defined "blue-envelope", and is strikingly similar to that of color-selected, U-dropout galaxies at z ~ 3. We also find that although the co-moving luminosity density at rest-frame UV wavelengths (1600 Angstroms) is roughly comparable at z ~ 3 and z ~ 4, the luminosity density at rest-frame optical wavelengths increases by about one-third from z ~ 4 to z ~ 3. Although the star-formation histories of individual galaxies may involve complex and stochastic events, the evolution in the global luminosity density of the UV-bright galaxy population corresponds to an average star-formation history with a star-formation rate that is constant or increasing over these redshifts. This suggests that the evolution in the luminosity density corresponds to an increase in the stellar-mass density of more than 33%.

We have measured the rest--frame lambda~1500 Ang comoving specific luminosity density of star--forming galaxies at redshift 3.5<z<6.5 from deep images taken with the Hubble Space Telescope (HST)and the Advanced Camera for Surveys (ACS), obtained as part of the Great Observatories Origins Deep Survey (GOODS). We used color selection criteria to construct samples of star--forming galaxies at redshifts z~4, 5 and 6, identified by the signature of the 912 Ang Lyman continuum discontinuity and Lyman-alpha forest blanketing in their rest--frame UV colors (Lyman--break galaxies). The ACS samples cover ~0.09 square degree, and are also relatively deep, reaching between 0.2 and 0.5 L_3^*, depending on the redshift, where $L_3^*$ is the characteristic UV luminosity of Lyman--break galaxies at z~3. The specific luminosity density of Lyman--break galaxies appears to be nearly constant with redshift from z~3 to z~6, although the measure at z~6 remains relatively uncertain, because it depends on the accurate estimate of the faint counts of the z~6 sample. If Lyman--break galaxies are fair tracers of the cosmic star formation activity, our results suggest that at z~6 the universe was already producing stars as vigorously as it did near its maximum several Gyr later, at 1<~z<~3. Thus, the onset of large--scale star formation in the universe is to be sought at around z~6 or higher, namely at less than ~7% of the current cosmic age.

Using U'- through Ks-band imaging data in the GOODS-South field, we construct a large, complete sample of 275 ``extremely red objects'' (EROs; K_s<22.0, R-K_s>3.35; AB), all with deep HST/ACS imaging in B_435, V_606, i_775, and z_850, and well-calibrated photometric redshifts. Quantitative concentration and asymmetry measurements fail to separate EROs into distinct morphological classes. We therefore visually classify the morphologies of all EROs into four broad types: ``Early'' (elliptical-like), ``Late'' (disk galaxies), ``Irregular'' and ``Other'' (chain galaxies and low surface brightness galaxies), and calculate their relative fractions and comoving space densities. For a broad range of limiting magnitudes and color thresholds, the relative number of early-type EROs is approximately constant at 33-44%, and the comoving space densities of Early- and Late-type EROs are comparable. Mean rest-frame spectral energy distributions (SEDs) at wavelengths between 0.1 and 1.2 um are constructed for all EROs. The SEDs are extremely similar in their range of shapes, independent of morphological type. The implication is that any differences between the broad-band SEDs of Early-type EROs and the other types are relatively subtle, and there is no robust way of photometrically distinguishing between different morphological types with usual optical/near-infrared photometry.

This Special Issue of the Astrophysical Journal Letters is dedicated to presenting initial results from the Great Observatories Origins Deep Survey (GOODS) that are primarily, but not exclusively, based on multi--band imaging data obtained with the {\it Hubble Space Telescope} ({\it HST}) and the Advanced Camera for Surveys (ACS). The survey covers roughly 320 square arcminutes in the ACS F435W, F606W, F814W, and F850LP bands, divided into two well-studied fields. Existing deep observations from the {\it Chandra} X-ray Observatory (\cxo) and ground-based facilities are supplemented with new, deep imaging in the optical and near-infrared from the European Southern Observatory (ESO) and from the Kitt Peak National Observatory (KPNO). Deep observations with the {\it Space Infrared Telescope Facility} (SIRTF) are scheduled. Reduced data from all facilities are being released worldwide within three to six months of acquisition. Together, this data set provides two deep reference fields for studies of distant normal and active galaxies, supernovae, and faint stars in our own galaxy. This paper serves to outline the survey strategy and describe the specific data that have been used in the accompanying letters, summarizing the reduction procedures and sensitivity limits.

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later.

We use extensive multi-wavelength photometric data from the Great Observatories Origins Deep Survey (GOODS) to estimate photometric redshifts for a sample of 434 galaxies with spectroscopic redshifts in the Chandra Deep Field South. Using the Bayesian method, which incorporates redshift/magnitude priors, we estimate photometric redshifts for galaxies in the range 18 < R (AB) < 25.5, giving an rms scatter of 0.11. The outlier fraction is < 10%, with the outlier-clipped rms being 0.047. We examine the accuracy of photometric redshifts for several, special sub--classes of objects. The results for extremely red objects are more accurate than those for the sample as a whole, with rms of 0.051 and very few outliers (3%). Photometric redshifts for active galaxies, identified from their X-ray emission, have a dispersion of 0.104, with 10% outlier fraction, similar to that for normal galaxies. Employing a redshift/magnitude prior in this process seems to be crucial in improving the agreement between photometric and spectroscopic redshifts.

We report early results on galaxies at z~6, selected from Hubble Space Telescope imaging for the Great Observatories Origins Deep Survey. Spectroscopy of one object with the Advanced Camera for Surveys grism and from the Keck and VLT observatories a shows a strong continuum break and asymmetric line emission, identified as Lya at z=5.83. We detect only five spatially extended, z~6 candidates with signal-to-noise ratios > 10, two of which have spectroscopic confirmation. This is many fewer than would be expected if galaxies at z=6 had the same luminosity function as those at z=3. There are many fainter candidates, but we expect substantial contamination from foreground interlopers and spurious detections. Our best estimates favor a z=6 galaxy population with fainter luminosities, higher space density, and similar co-moving ultraviolet emissivity to that at z=3, but this depends critically on counts at fluxes fainter than those reliably probed by the current data.

Cosmic variance is the uncertainty in observational estimates of the volume density of extragalactic objects such as galaxies or quasars arising from the underlying large-scale density fluctuations. This is often a significant source of uncertainty, especially in deep galaxy surveys, which tend to cover relatively small areas. We present estimates of the relative cosmic variance for one-point statistics (i.e. number densities) for typical scales and volumes sampled by the Great Observatories Origins Deep Survey (GOODS). We use two approaches: for objects with a known two-point correlation function that is well-approximated by a power law, one can use the standard analytic formalism to calculate the cosmic variance (in excess of shot noise). We use this approach to estimate the cosmic variance for several populations that are being studied in the GOODS program: Extremely Red Objects (ERO) at z~1, and Lyman Break Galaxies (LBG) at z~3 and z~4, using clustering information for similar populations in the literature. For populations with unknown clustering, one can use predictions from Cold Dark Matter theory to obtain a rough estimate of the variance as a function of number density. We present a convenient plot which allows one to use this approach to read off the cosmic variance for a population with a known mean redshift and estimated number density. We conclude that for the volumes sampled by GOODS, cosmic variance is a significant source of uncertainty for strongly clustered objects (~40-60% for EROs) and less serious for less clustered objects, ~10-20% for LBGs.

We investigate several fundamental properties of z ~ 4 Lyman-break galaxies by comparing observations with the predictions of a semi-analytic model based on the Cold Dark Matter theory of hierarchical structure formation. We use a sample of B_{435}-dropouts from the Great Observatories Origins Deep Survey, and complement the ACS optical B_{435}, V_{606}, i_{775}, and z_{850} data with the VLT ISAAC J, H, and K_{s} observations. We extract B_{435}-dropouts from our semi-analytic mock catalog using the same color criteria and magnitude limits that were applied to the observed sample. We find that the i_{775} - K_{s} colors of the model-derived and observed B_{435}-dropouts are in good agreement. However, we find that the i_{775}-z_{850} colors differ significantly, indicating perhaps that either too little dust or an incorrect extinction curve have been used. Motivated by the reasonably good agreement between the model and observed data we present predictions for the stellar masses, star formation rates, and ages for the z ~ 4 Lyman-break sample. We find that according to our model the color selection criteria used to select our z ~ 4 sample surveys 67% of all galaxies at this epoch down to z_{850} < 26.5. We find that our model predicts a roughly 40% mass build-up between the z \~ 4 and z ~ 3 epochs for the UV rest-frame L* galaxies. Furthermore, according to our model, at least 50% of the total stellar mass resides in relatively massive UV-faint objects that fall below our observational detection limit.

Multi-band near-infrared images of twelve BL Lacertae objects were obtained with the 2.5m telescope at the Las Campanas Observatory in order to determine the properties of their underlying host galaxies. Resolved emission was clearly detected in eight of the lowest redshift targets (up to z~0.3), and was modeled with a de Vaucouleurs r^{1/4} surface brightness law. We find that the morphologies match the elliptical galaxy profiles well, and that the BL Lac objects reside in large and luminous, but otherwise normal hosts -- consistent with previous studies done predominantly at optical wavelengths. The median absolute K-band magnitude of the galaxies in this study is -26.2, the average half-light radius is 4.2 +/- 2.3 kpc, and their average integrated R-K color is 2.7 +/- 0.3 mag. These are well within the range of values measured previously in the H-band by Kotilainen et al. and Scarpa et al. in a comparable number of targets. Taking their data together with our results, we find a best-fit K-band Kormendy relation of \mu_{e} = 4.3 log_{10} (r_{e}/kpc) + 14.2 mag arcsec^{-2}, virtually identical to that obtained for normal ellipticals. Finally, the near-infrared colors determined for five galaxies (average J-K = 0.8 +/- 0.3 mag) are the first such measurements for BL Lac hosts, and match those expected from old stellar populations at the BL Lac redshifts.