We examine the massive colliding cluster El Gordo, one of the most massive clusters at high redshift. We use a free-form lensing reconstruction method that avoids making assumptions about the mass distribution. We use data from the RELICS program and identify new multiply lensed system candidates. The new set of constraints and free-form method provides a new independent mass estimate of this intriguing colliding cluster. Our results are found to be consistent with earlier parametric models, indirectly confirming the assumptions made in earlier work. By fitting a double gNFW profile to the lens model, and extrapolating to the virial radius, we infer a total mass for the cluster of $M_{200c}=(1.08^{+0.65}_{-0.12})\times10^{15}$M$_{\odot}$. We estimate the uncertainty in the mass due to errors in the photometric redshifts, and discuss the uncertainty in the inferred virial mass due to the extrapolation from the lens model. We also find in our lens map a mass overdensity corresponding to the large cometary tail of hot gas, reinforcing its interpretation as a large tidal feature predicted by hydrodynamical simulations that mimic El Gordo. Finally, we discuss the observed relation between the plasma and the mass map, finding that the peak in the projected mass map may be associated with a large concentration of colder gas, exhibiting possible star formation. El Gordo is one of the first clusters that will be observed with JWST, which is expected to unveil new high redshift lensed galaxies around this interesting cluster, and provide a more accurate estimation of its mass.

We present a characterization of the physical properties of a sample of 35 securely-detected, dusty galaxies in the deep ALMA 1.2-mm image obtained as part of the ALMA Spectroscopic Survey in the {\it Hubble} Ultra Deep Field (ASPECS) Large Program. This sample is complemented by 26 additional sources identified via an optical/infrared source positional prior. Using their well-characterized spectral energy distributions, we derive median stellar masses and star formation rates (SFR) of $4.8\times10^{10}~M_\odot$ and 30 $M_\odot$ yr$^{-1}$, and interquartile ranges of $(2.4-11.7)\times10^{10}~M_\odot$ and $20-50~M_\odot$ yr$^{-1}$, respectively. We derive a median spectroscopic redshift of 1.8 with an interquartile range $1.1-2.6$, significantly lower than sub-millimeter galaxies detected in shallower, wide-field surveys. We find that 59\%$\pm$13\%, 6\%$\pm$4\%, and 34\%$\pm$9\% of our sources are within, above and below $\pm0.4$ dex from the SFR- stellar mass relation or main-sequence (MS), respectively. The ASPECS galaxies closely follow the SFR-molecular gas mass relation and other previously established scaling relations, confirming a factor of five increase of the gas-to-stellar mass ratio from $z=0.5$ to $z=2.5$, and a mild evolution of the gas depletion timescales with a typical value of 0.7 Gyr at $z=1-3$. ASPECS galaxies located significantly below the MS, a poorly exploited parameter space, have low gas-to-stellar mass ratios of $\sim0.1-0.2$ and long depletion timescales $>1$ Gyr. Galaxies along the MS dominate the cosmic density of molecular gas at all redshifts. Systems above the main sequence have an increasing contribution to the total gas reservoirs from $z<1$ to $z=2.5$, while the opposite is found for galaxies below the MS.

While low-luminosity galaxies dominate number counts at all redshifts, their contribution to cosmic Reionization is poorly understood due to a lack of knowledge of their physical properties. We isolate a sample of 35 z~4-5 continuum-faint Lyman-alpha emitters from deep VLT/MUSE spectroscopy and directly measure their Halpha emission using stacked Spitzer/IRAC Ch. 1 photometry. Based on Hubble Space Telescope imaging, we determine that the average UV continuum magnitude is fainter than -16 (~0.01 L_star), implying a median Lyman-alpha equivalent width of 249 Angstroms. By combining the Halpha measurement with the UV magnitude we determine the ionizing photon production efficiency, xi_ion, a first for such faint galaxies. The measurement of log (xi_ion [Hz/erg]) = 26.28 (+0.28; -0.40) is in excess of literature measurements of both continuum- and emission line-selected samples, implying a more efficient production of ionizing photons in these lower-luminosity, Lyman-alpha-selected systems. We conclude that this elevated efficiency can be explained by stellar populations with metallicities between 4e-4 and 0.008, with light-weighted ages less than 3 Myr.

We use Hubble Space Telescope (HST) imaging and near-infrared spectroscopy from Keck/MOSFIRE to study the sub-structure around the progenitor of a Milky Way-mass galaxy in the Hubble Frontier Fields (HFF). Specifically, we study an $r_e = 40^{+70}_{-30}$pc, $M_{\star} \sim 10^{8.2} M_{\odot}$ rest-frame ultra-violet luminous "clump" at a projected distance of $\sim$100~pc from a $M_{\star} \sim 10^{9.8}$M$_{\odot}$ galaxy at $z = 2.36$ with a magnification $\mu = 5.21$. We measure the star formation history of the clump and galaxy by jointly modeling the broadband spectral energy distribution from HST photometry and H$\alpha$ from MOSFIRE spectroscopy. Given our inferred properties (e.g., mass, metallicity, dust) of the clump and galaxy, we explore scenarios in which the clump formed \emph{in-situ} (e.g., a star forming complex) or \emph{ex-situ} (e.g., a dwarf galaxy being accreted). If it formed \emph{in-situ}, we conclude that the clump is likely a single entity as opposed to a aggregation of smaller star clusters, making it one of the most dense star clusters cataloged. If it formed \emph{ex-situ}, then we are witnessing an accretion event with a 1:40 stellar mass ratio. However, our data alone are not informative enough to distinguish between \emph{in-situ} and \emph{ex-situ} scenarios to a high level of significance. We posit that the addition of high-fidelity metallicity information, such as [OIII]4363\AA, which can be detected at modest S/N with only a few hours of JWST/NIRSpec time, may be a powerful discriminant. We suggest that studying larger samples of moderately lensed sub-structures across cosmic time can provide unique insight into the hierarchical formation of galaxies like the Milky Way.

We present the results from the 1.2 mm continuum image obtained as part of the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS). The 1.2 mm continuum image has a size of 2.9 (4.2) arcmin$^2$ within a primary beam response of 50% (10%) and a rms value of $9.3\thinspace{\rm\mu Jy\thinspace beam^{-1}}$. We detect 35 sources at high significance (Fidelity $\geq0.5$), 32 of these have well characterized near-infrared HST counterparts. We estimate the 1.2 mm number counts to flux levels of $<30\thinspace{\rm\mu Jy}$ in two different ways: we first use the detected sources to constrain the number counts and find a significant flattening of the counts below $S_\nu \sim 0.1$ mJy. In a second approach, we constrain the number counts by using a probability of deflection statistics (P(D)) analysis. For this latter approach, we describe new methods to accurately measure the noise in interferometric imaging (employing jack-knifing in the cube and in the visibility plane). This independent measurement confirms the flattening of the number counts. Our analysis of the differential number counts shows that we are detecting $\sim$93% ($\sim$100% if we include the lower fidelity detections) of the total continuum dust emission associated to galaxies in the HUDF. The ancillary data allows us to study the dependence of the 1.2 mm number counts on redshift ($z=0-4$), galaxy dust mass (${\rm M}_{\rm dust}=10^{7}-10^{9}{\rm M}_{\odot}$), stellar mass (${\rm M}_{*}=10^{9}-10^{12}{\rm M}_{\odot}$), and star-formation rate (${\rm SFR}=1-1000\thinspace{\rm M}_{\odot}\thinspace{\rm yr^{-1}}$). In an accompanying paper we show that the number counts are crucial to constrain galaxy evolution models and the understanding of star-forming galaxies at high redshift.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

The ionizing photon production efficiency, $\xi_{ion}$, is a critical parameter that provides a number of physical constraints to the nature of the early Universe, including the contribution of galaxies to the timely completion of the reionization of the Universe. Here we use KECK/MOSFIRE and ZFOURGE multi-band photometric data to explore the $\xi_{ion}$ of a population of galaxies at $z\sim2$ with $log_{10}(M_*/M_\odot)\sim9.0-11.5$. Our 130 \Halpha\ detections show a median $log_{10}(\xi_{ion}[Hz/erg])$ of $24.8\pm0.5$ when dust corrected using a Calzetti et al. (2000) dust prescription. Our values are typical of mass/magnitude selected $\xi_{ion}$ values observed in the $z\sim2$ Universe. Using BPASSv2.2.1 and Starburst99 stellar population models with simple parametric star-formation-histories (SFH), we find that even with models that account for effects of stellar evolution with binaries/stellar rotation, model galaxies at $log_{10}(\xi_{ion}[Hz/erg])\lesssim25.0$ have low H$\alpha$ equivalent widths (EW) and redder colors compared to our $z\sim2$ observed sample. We find that introducing star-bursts to the SFHs resolve the tension with the models, however, due to the rapid time evolution of $\xi_{ion}$, H$\alpha$ EWs, and rest-frame optical colors, our Monte Carlo simulations of star-bursts show that random distribution of star-bursts in evolutionary time of galaxies are unlikely to explain the observed distribution. Thus, either our observed sample is specially selected based on their past SFH or stellar models require additional mechanisms to reproduce the observed high UV luminosity of galaxies for a given production rate of hydrogen ionizing photons.

The Hubble Frontier Fields offer an exceptionally deep window into the high-redshift universe, covering a substantially larger area than the Hubble Ultra-Deep field at low magnification and probing 1--2 mags deeper in exceptional high-magnification regions. We aim to leverage Atacama Large Millimetre Array (ALMA) band 6 ($\approx$263\,GHz) mosaics in the central portions of five Frontier Fields to characterize the infrared (IR) properties of $1582$ ultraviolet (UV)-selected Lyman-Break Galaxies (LBGs) at redshifts of $z {\sim}$2--8. We investigated individual and stacked fluxes and IR excess (IRX) values of the LBG sample as functions of stellar mass ($\mathrm{M}_{\bigstar}$), redshift, UV luminosity and slope $\beta$, and lensing magnification. Two (2) LBG candidates were individually detected above a significance of $4.1{-}\sigma$, while stacked samples of the remaining LBG candidates yielded no significant detections. We investigated our detections and upper limits in the context of the IRX-$\mathrm{M}_{\bigstar}$ and IRX-$\beta$ relations, probing at least one dex lower in stellar mass than past studies have done. Our upper limits exclude substantial portions of parameter space and they are sufficiently deep in a handful of cases to create mild tension with the typically assumed attenuation and consensus relations. We observe a clear and smooth trend between $\mathrm{M}_{\bigstar}$ and $\beta$, which extends to low masses and blue (low) $\beta$ values, consistent with expectations from previous works.

We present 16 new ultrabright $H_{AB}\lesssim25$ galaxy candidates at z~8 identified over the COSMOS/UltraVISTA field. The new search takes advantage of the deepest-available ground-based optical and near-infrared observations, including the DR3 release of UltraVISTA and full-depth Spitzer/IRAC observations from the SMUVS and SPLASH programs. Candidates are selected using Lyman-break criteria, combined with strict optical non-detection and SED-fitting criteria, minimizing contamination by low-redshift galaxies and low-mass stars. HST/WFC3 coverage from the DASH program reveals that one source evident in our ground-based near-IR data has significant substructure and may actually correspond to 3 separate z~8 objects, resulting in a sample of 18 galaxies, 10 of which seem to be fairly robust (with a >97% probability of being at z>7). The UV-continuum slope $\beta$ for the bright z~8 sample is $\beta=-2.2\pm0.6$, bluer but still consistent with that of similarly bright galaxies at z~6 ($\beta=-1.55\pm0.17$) and z~7 ($\beta=-1.75\pm0.18$). Their typical stellar masses are 10$^{9.1^{+0.5}_{-0.4}}M_{\odot}$, with the SFRs of $32^{+44}_{-32}M_{\odot}$/year, specific SFR of $4^{+8}_{-4}$ Gyr$^{-1}$, stellar ages of $\sim22^{+69}_{-22}$\,Myr, and low dust content A$_V=0.15^{+0.30}_{-0.15}$ mag. Using this sample we constrain the bright end of the z~8 UV luminosity function (LF). When combined with recent empty field LF estimates at z~8-9, the resulting z~8 LF can be equally well represented by either a Schechter or a double power-law (DPL) form. Assuming a Schechter parameterization, the best-fit characteristic magnitude is $M^*= -20.95^{+0.30}_{-0.35}$ mag with a very steep faint end slope $\alpha=-2.15^{+0.20}_{-0.19}$. These new candidates include amongst the brightest yet found at these redshifts, 0.5-1.0 mag brighter than found over CANDELS, providing excellent targets for follow-up studies.

This manuscript describes the public release of the Hubble Legacy Fields (HLF) project photometric catalog for the extended GOODS-South region from the Hubble Space Telescope (HST) archival program AR-13252. The analysis is based on the version 2.0 HLF data release that now includes all ultraviolet (UV) imaging, combining three major UV surveys. The HLF data combines over a decade worth of 7475 exposures taken in 2635 orbits totaling 6.3 Msec with the HST Advanced Camera for Surveys Wide Field Channel (ACS/WFC) and the Wide Field Camera 3 (WFC3) UVIS/IR Channels in the greater GOODS-S extragalactic field, covering all major observational efforts (e.g., GOODS, GEMS, CANDELS, ERS, UVUDF and many other programs; see Illingworth et al 2019, in prep). The HLF GOODS-S catalogs include photometry in 13 bandpasses from the UV (WFC3/UVIS F225W, F275W and F336W filters), optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters), to near-infrared (WFC3/IR F098M, F105W, F125W, F140W and F160W filters). Such a data set makes it possible to construct the spectral energy distributions (SEDs) of objects over a wide wavelength range from high resolution mosaics that are largely contiguous. Here, we describe a photometric analysis of 186,474 objects in the HST imaging at wavelengths 0.2--1.6$\mu$m. We detect objects from an ultra-deep image combining the PSF-homogenized and noise-equalized F850LP, F125W, F140W and F160W images, including Gaia astrometric corrections. SEDs were determined by carefully taking the effects of the point-spread function in each observation into account. All of the data presented herein are available through the HLF website (https://archive.stsci.edu/prepds/hlf/).

We present the first measurements of the Lyman-continuum photon production efficiency $\xi_{\textrm{ion,0}}$ at $z\sim4$-5 for galaxies fainter than 0.2 $L^*$ ($-$19 mag). $\xi_{\textrm{ion,0}}$ quantifies the production rate of ionizing photons with respect to the UV luminosity density assuming a fiducial escape fraction of zero. Extending previous measurements of $\xi_{\textrm{ion,0}}$ to the faint population is important, as ultra-faint galaxies are expected to contribute the bulk of the ionizing emissivity. We probe $\xi_{\textrm{ion,0}}$ to such faint magnitudes by taking advantage of 200-hour depth Spitzer/IRAC observations from the GREATS program and $\approx$300 3<$z$<6 galaxies with spectroscopic redshifts from the MUSE GTO Deep + Wide programs. Stacked IRAC [3.6]$-$[4.5] colors are derived and used to infer the H$\alpha$ rest-frame equivalent widths, which range from 403\r{A} to 2818\r{A}. The derived $\xi_{\textrm{ion,0}}$ is $\log_{10}(\xi_{\textrm{ion,0}} / \textrm{Hz erg}^{-1}) = 25.36 \pm 0.08$ over $-$20.5 < M$_{\textrm{UV}}$ < $-$17.5, similar to those derived for brighter galaxy samples at the same redshift and therefore suggesting that $\xi_{\textrm{ion}}$ shows no strong dependence on $M_{UV}$. The $\xi_{\textrm{ion,0}}$ values found in our sample imply that the Lyman-continuum escape fraction for $M_{\textrm{UV}} \approx -19$ star-forming galaxies cannot exceed $\approx$8-20\% in the reionization era.

We present the results obtained with VLT/MUSE on the faint-end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The precise aim of the present study is to further constrain the abundance of Lyman-alpha emitters (LAEs) by taking advantage of the magnification provided by lensing clusters. We blindly selected a sample of 156 LAEs, with redshifts between $2.9 \le z \le 6.7$ and magnification-corrected luminosities in the range $ 39 \lesssim \log L_{Ly_{\alpha}}$ [erg s$^{-1}$] $\lesssim 43$. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure. To properly take into account the individual differences in detection conditions (including lensing configurations, spatial and spectral morphologies) when computing the LF, a new method based on the 1/Vmax approach was implemented. The LAE LF has been obtained in four different redshift bins with constraints down to $\log L_{Ly_{\alpha}} = 40.5$. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis including data from the literature to complete the present sample a steep faint-end slope was measured varying from $\alpha = -1.69^{+0.08}_{-0.08}$ to $\alpha = -1.87^{+0.12}_{-0.12}$ between the lowest and the highest redshift bins. The contribution of the LAE population to the star formation rate density at $z \sim 6$ is $\lesssim 50$% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. (abridged)

We analyze the interstellar medium properties of a sample of sixteen bright CO line emitting galaxies identified in the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) Large Program. This CO$-$selected galaxy sample is complemented by a couple of additional CO line emitters in the UDF that are identified based on their MUSE optical spectroscopic redshifts. The ASPECS CO$-$selected galaxies cover a larger range of star-formation rates and stellar masses compared to literature CO emitting galaxies at $z>1$ for which scaling relations have been established previously. Most of ASPECS CO-selected galaxies follow these established relations in terms of gas depletion timescales and gas fractions as a function of redshift, as well as the star-formation rate-stellar mass relation (`galaxy main sequence'). However, we find that $\sim30\%$ of the galaxies (5 out of 16) are offset from the galaxy main sequence at their respective redshift, with $\sim12\%$ (2 out of 16) falling below this relationship. Some CO-rich galaxies exhibit low star-formation rates, and yet show substantial molecular gas reservoirs, yielding long gas depletion timescales. Capitalizing on the well-defined cosmic volume probed by our observations, we measure the contribution of galaxies above, below, and on the galaxy main sequence to the total cosmic molecular gas density at different lookback times. We conclude that main sequence galaxies are the largest contributor to the molecular gas density at any redshift probed by our observations (z$\sim$1$-$3). The respective contribution by starburst galaxies above the main sequence decreases from z$\sim$2.5 to z$\sim$1, whereas we find tentative evidence for an increased contribution to the cosmic molecular gas density from the passive galaxies below the main sequence.

We present the Super Eight galaxies - a set of very luminous, high-redshift ($7.1<z<8.0$) galaxy candidates found in Brightest of Reionizing Galaxies (BoRG) Survey fields. The original sample includes eight galaxies that are $Y$-band dropout objects with $H$-band magnitudes of $m_H<25.5$. Four of these objects were originally reported in Calvi et al. 2016. Combining new Hubble Space Telescope (HST) WFC3/F814W imaging and $Spitzer$ IRAC data with archival imaging from BoRG and other surveys, we explore the properties of these galaxies. Photometric redshift fitting places six of these galaxies in the redshift range of $7.1<z<8.0$, resulting in three new high-redshift galaxies and confirming three of the four high-redshift galaxy candidates from Calvi et al. 2016. We calculate the half-light radii of the Super Eight galaxies using the HST F160W filter and find that the Super Eight sizes are in line with typical evolution of size with redshift. The Super Eights have a mean mass of log(M$_*$/M$_\odot$) $\sim10$, which is typical for sources in this luminosity range. Finally, we place our sample on the UV $z\sim8$ luminosity function and find that the Super Eight number density is consistent with other surveys in this magnitude and redshift range.

We report the results of an expanded search for z~9-10 candidates over the ~883 arcmin^2 CANDELS+ERS fields. This study adds 147 arcmin^2 to the search area we consider over the CANDELS COSMOS, UDS, and EGS fields, while expanding our selection to include sources with bluer J_{125}-H_{160} colors than our previous J_{125}-H_{160}>0.5 mag selection. In searching for new z~9-10 candidates, we make full use of all available HST, Spitzer/IRAC, and ground-based imaging data. As a result of our expanded search and use of broader color criteria, 3 new candidate z~9-10 galaxies are identified. We also find again the z=8.683 source previously confirmed by Zitrin+2015. This brings our sample of probable z~9-11 galaxy candidates over the CANDELS+ERS fields to 19 sources in total, equivalent to 1 candidate per 47 arcmin^2 (1 per 10 WFC3/IR fields). To be comprehensive, we also discuss 28 mostly lower likelihood z~9-10 candidates, including some sources that seem to be reliably at z>8 using the HST+IRAC data alone, but which the ground-based data show are much more likely at z<4. One case example is a bright z~9.4 candidate COS910-8 which seems instead to be at z~2. Based on this expanded sample, we obtain a more robust LF at z~9 and improved constraints on the volume density of bright z~9 and z~10 galaxies. Our improved z~9-10 results again reinforce previous findings for strong evolution in the UV LF at z>8, with a factor of ~10 evolution seen in the luminosity density from z~10 to z~8.

The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$\beta$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H$\beta$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim8$ galaxies, such as the ionizing photon production efficiency with $\log(\xi_\mathrm{ion}/\mathrm{erg}^{-1}\hspace{1mm}\mathrm{Hz})\geq25.77$. We find a relatively tight correlation between the [OIII]+H$\beta$ and UV luminosity, which we use to derive for the first time the [OIII]+H$\beta$ luminosity function (LF) at $z\sim8$. The $z\sim8$ [OIII]+H$\beta$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H$\beta$ luminosity at a given UV luminosity from $z\sim3$ to $z\sim8$. Finally, using the [OIII]+H$\beta$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim30$ mag will be required.

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) provides new constraints for galaxy formation models on the molecular gas properties of galaxies. We compare results from ASPECS to predictions from two cosmological galaxy formation models: the IllustrisTNG hydrodynamical simulations and the Santa Cruz semi-analytic model (SC SAM). We explore several recipes to model the H$_2$ content of galaxies, finding them to be consistent with one another, and take into account the sensitivity limits and survey area of ASPECS. For a canonical CO-to-H$_2$ conversion factor of $\alpha_{\rm CO} = 3.6\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$ the results of our work include: (1) the H$_2$ mass of $z>1$ galaxies predicted by the models as a function of their stellar mass is a factor of 2-3 lower than observed; (2) the models do not reproduce the number of H$_2$-rich ($M_{\rm H2} > 3\times 10^{10}\,\rm{M}_\odot$) galaxies observed by ASPECS; (3) the H$_2$ cosmic density evolution predicted by IllustrisTNG (the SC SAM) is in tension (only just agrees) with the observed cosmic density, even after accounting for the ASPECS selection function and field-to-field variance effects. The tension between models and observations at $z>1$ can be alleviated by adopting a CO-to-H$_2$ conversion factor in the range $\alpha_{\rm CO} = 2.0 - 0.8\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$. Additional work on constraining the CO-to-H$_2$ conversion factor and CO excitation conditions of galaxies through observations and theory will be necessary to more robustly test the success of galaxy formation models.

We use the results from the ALMA large program ASPECS, the spectroscopic survey in the Hubble Ultra Deep Field (HUDF), to constrain CO luminosity functions of galaxies and the resulting redshift evolution of $\rho$(H$_2$). The broad frequency range covered enables us to identify CO emission lines of different rotational transitions in the HUDF at $z>1$. We find strong evidence that the CO luminosity function evolves with redshift, with the knee of the CO luminosity function decreasing in luminosity by an order of magnitude from $\sim$2 to the local universe. Based on Schechter fits, we estimate that our observations recover the majority (up to $\sim$90%, depending on the assumptions on the faint end) of the total cosmic CO luminosity at $z$=1.0-3.1. After correcting for CO excitation, and adopting a Galactic CO-to-H$_2$ conversion factor, we constrain the evolution of the cosmic molecular gas density $\rho$(H$_2$): this cosmic gas density peaks at $z\sim1.5$ and drops by factor of $6.5_{-1.4}^{+1.8}$ to the value measured locally. The observed evolution in $\rho$(H$_2$) therefore closely matches the evolution of the cosmic star formation rate density $\rho_{\rm SFR}$. We verify the robustness of our result with respect to assumptions on source inclusion and/or CO excitation. As the cosmic star formation history can be expressed as the product of the star formation efficiency and the cosmic density of molecular gas, the similar evolution of $\rho$(H$_2$) and $\rho_{\rm SFR}$ leaves only little room for a significant evolution of the average star formation efficiency in galaxies since $z\sim 3$ (85% of cosmic history).

We discuss the nature and physical properties of gas-mass selected galaxies in the ALMA spectroscopic survey (ASPECS) of the Hubble Ultra Deep Field (HUDF). We capitalize on the deep optical integral-field spectroscopy from the MUSE HUDF Survey and multi-wavelength data to uniquely associate all 16 line-emitters, detected in the ALMA data without preselection, with rotational transitions of carbon monoxide (CO). We identify ten as CO(2-1) at $1 < z < 2$, five as CO(3-2) at $2 < z < 3$ and one as CO(4-3) at $z = 3.6$. Using the MUSE data as a prior, we identify two additional CO(2-1)-emitters, increasing the total sample size to 18. We infer metallicities consistent with (super-)solar for the CO-detected galaxies at $z \le 1.5$, motivating our choice of a Galactic conversion factor between CO luminosity and molecular gas mass for these galaxies. Using deep Chandra imaging of the HUDF, we determine an X-ray AGN fraction of 20% and 60% among the CO-emitters at $z \sim 1.4$ and $z \sim 2.6$, respectively. Being a CO-flux limited survey, ASPECS-LP detects molecular gas in galaxies on, above and below the main sequence (MS) at $z \sim 1.4$. For stellar masses $\ge 10^{10} (10^{10.5})$ M$_{\odot}$, we detect about 40% (50%) of all galaxies in the HUDF at $1 < z < 2$ ($2 < z < 3$). The combination of ALMA and MUSE integral-field spectroscopy thus enables an unprecedented view on MS galaxies during the peak of galaxy formation.

The ALMA SPECtroscopic Survey in the {\it Hubble} Ultra Deep Field is an ALMA large program that obtained a frequency scan in the 3\,mm band to detect emission lines from the molecular gas in distant galaxies. We here present our search strategy for emission lines and continuum sources in the HUDF. We compare several line search algorithms used in the literature, and critically account for the line-widths of the emission line candidates when assessing significance. We identify sixteen emission lines at high fidelity in our search. Comparing these sources to multi-wavelength data we find that all sources have optical/infrared counterparts. Our search also recovers candidates that have lower significance that can be used statistically to derive, e.g. the CO luminosity function. We apply the same detection algorithm to obtain a sample of six 3 mm continuum sources. All of these are also detected in the 1.2 mm continuum with optical/near-infrared counterparts. We use the continuum sources to compute 3 mm number counts in the sub-mJy regime, and find them to be higher by an order of magnitude than expected for synchrotron-dominated sources. However, the number counts are consistent with those derived at shorter wavelengths (0.85--1.3\,mm) once extrapolating to 3\,mm with a dust emissivity index of $\beta=1.5$, dust temperature of 35\,K and an average redshift of $z=2.5$. These results represent the best constraints to date on the faint end of the 3 mm number counts.

One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

We take advantage of new 100-hour Spitzer/IRAC observations available for MACS0647-JD, a strongly lensed $z\approx11$ galaxy candidate, to provide improved constraints on its physical properties. Probing the physical properties of galaxies at $z>8$ is challenging due to the inherent faintness of such sources and the limited wavelength coverage available. Thanks to the high $\approx$2-6$\times$ lensing magnification of the multiple images of MACS0647-JD, we can use the sensitive Spitzer/IRAC data to probe the rest-frame optical fluxes of MACS0647-JD and investigate its physical properties including the age and the stellar mass. In deriving Spitzer/IRAC fluxes for MACS0647-JD, great care is taken in coping with the impact of three bright ($\approx$8-16 mag) stars in our field to ensure robust results. Assuming a constant star formation rate, the age, stellar mass, and rest-frame UV slope we estimate for MACS0647-JD based on a stack of the photometry are log$_{10}$(age/yr) = 8.6$^{+0.1}_{-2.1}$, log$_{10}$(M$_{*}$/M$_{\odot}$) = 9.1$^{+0.2}_{-1.4}$, and $\beta = -$1.3$\pm$0.6, respectively. We compare our results with expectations from the EAGLE simulation and find that MACS0647-JD has properties consistent with corresponding to the most massive and rapidly star-forming galaxies in the simulation. We also find that its radius, 105$\pm$28 pc, is a factor of $\approx$2 smaller than the mean size in a separate simulation project DRAGONS. Interestingly enough, the observed size is similar to the small sizes seen in very low-luminosity $z\approx$6-10 galaxies behind lensing clusters.

Deep optical spectroscopic surveys of galaxies provide us a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at ${z \sim 2-4.5}$. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of HeII${\lambda}$1640 emitters at ${z\sim2}$. The deepest areas of our MUSE pointings reach a ${3\sigma}$ line flux limit of 3.1${\times}$10-19 erg/ s/ cm$^2$. After discarding broad line active galactic nuclei we find 13 HeII${\lambda}$1640 detections from MUSE with a median MUV = $-20.1$ and 21 tentative HeII${\lambda}$1640 detections from other public surveys. Excluding Ly${\alpha}$, all except two galaxies in our sample show at least one other rest-UV emission line, with CIII]${\lambda}$1907,${\lambda}$1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample via spectral energy distribution fitting techniques. Taking advantage of the high quality spectra obtained by MUSE (${\sim10 - 30}$h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the HeII${\lambda}$1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame HeII${\lambda}$1640 equivalent widths (${\sim}$ 0.2 - 10 A), thus additional mechanisms are necessary in models to match the observed HeII${\lambda}$1640 properties.

Recent observations of lensed sources have shown that the faintest ($M_{\mathrm{UV}} \approx -15\,\mathrm{mag}$) galaxies observed at z=6-8 appear to be extremely compact. Some of them have inferred sizes of less than 40 pc for stellar masses between $10^6$ and $10^7\,\mathrm{M}_{\odot}$, comparable to individual super star clusters or star cluster complexes at low redshift. High-redshift, low-mass galaxies are expected to show a clumpy, irregular morphology and if star clusters form in each of these well-separated clumps, the observed galaxy size would be much larger than the size of an individual star forming region. As supernova explosions impact the galaxy with a minimum delay time that exceeds the time required to form a massive star cluster, other processes are required to explain the absence of additional massive star forming regions. In this work we investigate whether the radiation of a young massive star cluster can suppress the formation of other detectable clusters within the same galaxy already before supernova feedback can affect the galaxy. We find that in low-mass ($M_{200} \lesssim 10^{10}\,\mathrm{M}_{\odot}$) haloes, the radiation from a compact star forming region with an initial mass of $10^{7}\,\mathrm{M}_{\odot}$ can keep gas clumps with Jeans masses larger than $\approx 10^{7}\,\mathrm{M}_{\odot}$ warm and ionized throughout the galaxy. In this picture, the small intrinsic sizes measured in the faintest $z=6-8$ galaxies are a natural consequence of the strong radiation field that stabilises massive gas clumps. A prediction of this mechanism is that the escape fraction for ionizing radiation is high for the extremely compact, high-z sources.

We investigate the dependence of galaxy clustering at $z \sim 4 - 7$ on UV-luminosity and stellar mass. Our sample consists of $\sim$ 10,000 Lyman-break galaxies (LBGs) in the XDF and CANDELS fields. As part of our analysis, the $M_\star - M_{\rm UV}$ relation is estimated for the sample, which is found to have a nearly linear slope of $d\log_{10} M_\star / d M_{\rm UV} \sim 0.44$. We subsequently measure the angular correlation function and bias in different stellar mass and luminosity bins. We focus on comparing the clustering dependence on these two properties. While UV-luminosity is only related to recent starbursts of a galaxy, stellar mass reflects the integrated build-up of the whole star formation history, which should make it more tightly correlated with halo mass. Hence, the clustering segregation with stellar mass is expected to be larger than with luminosity. However, our measurements suggest that the segregation with luminosity is larger with $\simeq 90\%$ confidence (neglecting contributions from systematic errors). We compare this unexpected result with predictions from the \textsc{Meraxes} semi-analytic galaxy formation model. Interestingly, the model reproduces the observed angular correlation functions, and also suggests stronger clustering segregation with luminosity. The comparison between our observations and the model provides evidence of multiple halo occupation in the small scale clustering.

[abridged] Characterizing the number counts of faint, dusty star-forming galaxies is currently a challenge even for deep, high-resolution observations in the FIR-to-mm regime. They are predicted to account for approximately half of the total extragalactic background light at those wavelengths. Searching for dusty star-forming galaxies behind massive galaxy clusters benefits from strong lensing, enhancing their measured emission while increasing spatial resolution. Derived number counts depend, however, on mass reconstruction models that properly constrain these clusters. We estimate the 1.1 mm number counts along the line of sight of three galaxy clusters, i.e. Abell 2744, MACSJ0416.1-2403 and MACSJ1149.5+2223, which are part of the ALMA Frontier Fields Survey. We perform detailed simulations to correct these counts for lensing effects. We use several publicly available lensing models for the galaxy clusters to derive the intrinsic flux densities of our sources. We perform Monte Carlo simulations of the number counts for a detailed treatment of the uncertainties in the magnifications and adopted source redshifts. We find an overall agreement among the number counts derived for the different lens models, despite their systematic variations regarding source magnifications and effective areas. Our number counts span ~2.5 dex in demagnified flux density, from several mJy down to tens of uJy. Our number counts are consistent with recent estimates from deep ALMA observations at a 3$\sigma$ level. Below $\approx$ 0.1 mJy, however, our cumulative counts are lower by $\approx$ 1 dex, suggesting a flattening in the number counts. In our deepest ALMA mosaic, we estimate number counts for intrinsic flux densities $\approx$ 4 times fainter than the rms level. This highlights the potential of probing the sub-10 uJy population in larger samples of galaxy cluster fields with deeper ALMA observations.

Using an ultra-deep blind survey with the MUSE integral field spectrograph on the ESO Very Large Telescope, we obtain spectroscopic redshifts to a depth never explored before: galaxies with observed magnitudes m > 30 - 32. Specifically, we detect objects via Lyman-alpha emission at 2.9 < z < 6.7 without individual continuum counterparts in areas covered by the deepest optical/near-infrared imaging taken by the Hubble Space Telescope, the Hubble Ultra Deep Field. In total, we find 102 such objects in 9 square arcminutes at these redshifts. Detailed stacking analyses confirm the Lyman-alpha emission as well as the 1216 Angstrom-breaks and faint UV continua (M_UV ~ -15). This makes them the faintest spectroscopically-confirmed objects at these redshifts, similar to the sources believed to reionize the universe. A simple model for the expected fraction of detected/undetected Lyman-alpha emitters as a function of luminosity is consistent with these objects being the high-equivalent width tail of the normal Lyman-alpha-emitter population at these redshifts.

We compare the sizes and luminosities of 307 faint z=6-8 sources revealed by the Hubble Frontier Fields (HFF) program with sources in the nearby universe. Making use of the latest lensing models and data from the first four HFF clusters with an extensive suite of public lens models, we measure both the sizes and luminosities for 153 z~6, 101 z~7, and 53 z~8 galaxies. The sizes range over more than a decade from ~500 to <50 pc. Extremely small sizes are inferred for many of our lowest luminosity sources, reaching individual sizes as small as 10-30 pc (the smallest is 11(-6)(+28) pc). The uncertainty in these measures ranges from 80 pc for the largest sources to typically about 20 pc for the smallest. Such sizes are smaller than extrapolations of the size-luminosity relation, and expectations for the completeness of our faint samples, suggesting a likely break in the size-luminosity relation at ~-17 mag with size proportional to L**(0.50(-0.11)(+0.10)). The sizes and luminosities of the lowest-luminosity sources are similar to those of single star cluster complexes like 30 Doradus in the lower-redshift universe and -- in a few cases -- super star clusters. Remarkably, our identification of these compact, faint star-forming sources in the z~6-8 universe also allow us to set upper limits on the proto-globular cluster LF at z~6. Comparisons with recent models allow us to rule out (with some caveats) some scenarios for proto-globular cluster formation and set useful upper limits on other less extreme ones. Our results suggest we may be very close to discovering a bona-fide population of forming globular clusters at high redshift.

Star-forming galaxies have been found to follow a relatively tight relation between stellar mass ($M_{*}$) and star formation rate (SFR), dubbed the `star formation sequence'. A turnover in the sequence has been observed, where galaxies with $M_{*} < 10^{10} {\rm M}_{\odot}$ follow a steeper relation than their higher mass counterparts, suggesting that the low-mass slope is (nearly) linear. In this paper, we characterise the properties of the low-mass end of the star formation sequence between $7 \leq \log M_{*}[{\rm M}_{\odot}] \leq 10.5$ at redshift $0.11 < z < 0.91$. We use the deepest MUSE observations of the Hubble Ultra Deep Field and the Hubble Deep Field South to construct a sample of 179 star-forming galaxies with high signal-to-noise emission lines. Dust-corrected SFRs are determined from H$\beta$ $\lambda 4861$ and H$\alpha$ $\lambda 6563$. We model the star formation sequence with a Gaussian distribution around a hyperplane between $\log M_{*}$, $\log {\rm SFR}$, and $\log (1+z)$, to simultaneously constrain the slope, redshift evolution, and intrinsic scatter. We find a sub-linear slope for the low-mass regime where $\log {\rm SFR}[{\rm M}_{\odot}/{\rm yr}] = 0.83^{+0.07}_{-0.06} \log M_{*}[{\rm M}_{\odot}] + 1.74^{+0.66}_{-0.68} \log (1+z)$, increasing with redshift. We recover an intrinsic scatter in the relation of $\sigma_{\rm intr} = 0.44^{+0.05}_{-0.04}$ dex, larger than typically found at higher masses. As both hydrodynamical simulations and (semi-)analytical models typically favour a steeper slope in the low-mass regime, our results provide new constraints on the feedback processes which operate preferentially in low-mass halos.

We present the Hubble Deep UV Legacy Survey (HDUV), a 132 orbit imaging program with the WFC3/UVIS camera onboard the Hubble Space Telescope (HST). The HDUV extends and builds on the few previous UV imaging surveys in the two GOODS/CANDELS-Deep fields to provide deep images over a total area of ~100 arcmin2 in the two filters F275W and F336W. Our release also includes all the F275W imaging data taken by the CANDELS survey, which were aligned using a novel approach and combined with the HDUV survey data. By reaching depths of 27.5-28.0 mag (5sigma, in 0.4" apertures), these are the deepest high-resolution UV data over such a large area taken to date. Such unique UV imaging enables a wide range of science by the community. Among the main goals of the HDUV survey are: (1) provide a complete sample of faint star-forming galaxies at z~1-3, (2) constrain the ionizing photon escape fraction from galaxies at z~2-3, and (3) track the build-up of bulges and the disappearance of clumpy disk galaxies through reliable internal stellar population properties at sub-kpc resolution out to z~3. The addition of the HDUV data further enhances the legacy value of the two GOODS/CANDELS-Deep fields, which now include deep 11-band HST imaging as well as very deep ancillary data from X-ray to radio, enabling unique multi-wavelength studies. Here, we provide an overview of the survey design, describe the data reduction, and highlight a few basic analyses on the images which are released to the community as high level science products via the Mikulski Archive for Space Telescopes (MAST).

AIMS: We aim to evaluate the near-infrared colors of brown dwarfs as observed with four major infrared imaging space observatories: the Hubble Space Telescope (HST), the James Webb Space Telescope (JWST), the Euclid mission, and the WFIRST telescope. METHODS: We used the SPLAT SPEX/ISPEX spectroscopic library to map out the colors of the M-, L-, and T-type dwarfs. We have identified which color-color combination is optimal for identifying broad type and which single color is optimal to then identify the subtype (e.g., T0-9). We evaluated each observatory separately as well as the narrow-field (HST and JWST) and wide-field (Euclid and WFIRST) combinations. RESULTS: The Euclid filters perform poorly typing brown dwarfs and WFIRST performs only marginally better, despite a wider selection of filters. WFIRST's W146 and F062 combined with Euclid's Y-band discriminates somewhat better between broad brown dwarf categories. However, subtyping with any combination of Euclid and WFIRST observations remains uncertain due to the lack of medium or narrow-band filters. We argue that a medium band added to the WFIRST filter selection would greatly improve its ability to preselect brown dwarfs its imaging surveys. CONCLUSIONS: The HST filters used in high-redshift searches are close to optimal to identify broad stellar type. However, the addition of F127M to the commonly used broad filter sets would allow for unambiguous subtyping. An improvement over HST is one of two broad and medium filter combinations on JWST: pairing F140M with either F150W or F162M discriminates very well between subtypes.

Massive foreground galaxy clusters magnify and distort the light of objects behind them, permitting a view into both the extremely distant and intrinsically faint galaxy populations. We present here the z ~ 6 - 8 candidate high-redshift galaxies from the Reionization Lensing Cluster Survey (RELICS), a Hubble and Spitzer Space Telescope survey of 41 massive galaxy clusters spanning an area of ~200 arcmin^2. These clusters were selected to be excellent lenses and we find similar high-redshift sample sizes and magnitude distributions as CLASH. We discover 321 candidate galaxies with photometric redshifts between z ~ 6 to z ~ 8, including extremely bright objects with H-band magnitudes of m_AB ~ 23 mag. As a sample, the observed (lensed) magnitudes of these galaxies are among the brightest known at z> 6, comparable to much wider, blank-field surveys. RELICS demonstrates the efficiency of using strong gravitational lenses to produce high-redshift samples in the epoch of reionization. These brightly observed galaxies are excellent targets for follow-up study with current and future observatories, including the James Webb Space Telescope.

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

We present an analysis of all prime HST legacy fields spanning >800 arcmin^2 for the search of z~10 galaxy candidates and the study of their UV luminosity function (LF). In particular, we present new z~10 candidates selected from the full Hubble Frontier Field (HFF) dataset. Despite the addition of these new fields, we find a low abundance of z~10 candidates with only 9 reliable sources identified in all prime HST datasets that include the HUDF09/12, the HUDF/XDF, all the CANDELS fields, and now the HFF survey. Based on this comprehensive search, we find that the UV luminosity function decreases by one order of magnitude from z~8 to z~10 at all luminosities over a four magnitude range. This also implies a decrease of the cosmic star-formation rate density by an order of magnitude within 170 Myr from z~8 to z~10. We show that this accelerated evolution compared to lower redshift can entirely be explained by the fast build-up of the dark matter halo mass function at z>8. Consequently, the predicted UV LFs from several models of galaxy formation are in good agreement with this observed trend, even though the measured UV LF lies at the low end of model predictions. In particular, the number of only 9 observed candidate galaxies is lower, by ~50%, than predicted by galaxy evolution models. The difference is generally still consistent within the Poisson and cosmic variance uncertainties. However, essentially all models predict larger numbers than observed. We discuss the implications of these results in light of the upcoming James Webb Space Telescope mission, which is poised to find much larger samples of z~10 galaxies as well as their progenitors at less than 400 Myr after the Big Bang.

The earliest galaxies are expected to emerge in the first billion years of the Universe during the Epoch of Reionization. However, both the spectroscopic confirmation of photometrically-selected galaxies at this epoch and the characterization of their early dynamical state has been hindered by the lack of bright, accessible lines to probe the velocity structure of their interstellar medium. We present the first ALMA spectroscopic confirmation of such sources at z > 6 using the far-infrared [C II]{\lambda}157.74{\mu}m emission line, and, for the first time, measurement of the velocity structure, for two galaxies at z = 6.8540+/-0.0003 and z = 6.8076+/-0.0002. Remarkably, the [C II] line luminosity from these galaxies is higher than previously found in `normal' star-forming galaxies at z > 6.5. This suggests that we are sampling a part of the galaxy population different from the galaxies found through detection of the Ly{\alpha} line. The luminous and extended [C II] detections reveal clear velocity gradients that, if interpreted as rotation, would suggest these galaxies have turbulent, yet rotation-dominated disks, with similar stellar-to-dynamical mass fractions as observed for H{\alpha} emitting galaxies 2 Gyr later at cosmic noon. Our novel approach for confirming galaxies during Reionization paves the way for larger studies of distant galaxies with spectroscopic redshifts from ALMA. Particularly important, this opens up opportunities for high angular-resolution [C II] dynamics in galaxies less than one billion years after the Big Bang.

We present rest-frame Ly$\alpha$ equivalent widths (EW) of 417 Ly$\alpha$ emitters (LAEs) detected with Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) at $2.9 < z < 6.6$ in the Hubble Ultra Deep Field. Based on the deep MUSE spectroscopy and ancillary Hubble Space Telescope (HST) photometry data, we carefully measured EW values taking into account extended Ly$\alpha$ emission and UV continuum slopes ($\beta$). Our LAEs reach unprecedented depths, both in Ly$\alpha$ luminosities and UV absolute magnitudes, from log($L_{\rm Ly\alpha}$/erg s$^{-1}$) $\sim$41.0 to 43.0 and from Muv $\sim$ -16 to -21 (0.01-1.0 $L^{*}_{\rm z=3}$). The EW values span the range of $\sim$ 5 to 240 \AA\ or larger, and their distribution can be well fitted by an exponential law $N = N_{\rm 0}$ exp($-$EW/$w_{\rm 0}$). Owing to the high dynamic range in Muv, we find that the scale factor, $w_{\rm 0}$, depends on Muv in the sense that including fainter Muv objects increases $w_{\rm 0}$, i.e., the Ando effect. The results indicate that selection functions affect the EW scale factor. Taking these effects into account, we find that our $w_{\rm 0}$ values are consistent with those in the literature within $1\sigma$ uncertainties at $2.9 < z < 6.6$ at a given threshold of Muv and $L_{\rm Ly\alpha}$. Interestingly, we find 12 objects with EW $>200$ \AA\ above $1\sigma$ uncertainties. Two of these 12 LAEs show signatures of merger or AGN activity: the weak CIV $\lambda 1549$ emission line. For the remaining 10 very large EW LAEs, we find that the EW values can be reproduced by young stellar ages ($< 100$ Myr) and low metallicities ($\lesssim 0.02$ $Z_{\rm \odot}$). Otherwise, at least part of the Ly$\alpha$ emission in these LAEs needs to arise from anisotropic radiative transfer effects, fluorescence by hidden AGN or quasi-stellar object activity, or gravitational cooling.

We present a measurement of the fraction of Lyman $\alpha$ (Ly$\alpha$) emitters ($X_{\rm{Ly} \alpha}$) amongst HST continuum-selected galaxies at $3<z<6$ with the Multi-Unit Spectroscopic Explorer (MUSE) on the VLT. Making use of the first 24 MUSE-Wide pointings in GOODS-South, each having an integration time of 1 hour, we detect 100 Ly$\alpha$ emitters and find $X_{\rm{Ly} \alpha}\gtrsim0.5$ for most of the redshift range covered, with 29 per cent of the Ly$\alpha$ sample exhibiting rest equivalent widths (rest-EWs) $\leq$ 15\AA. Adopting a range of rest-EW cuts (0 - 75\AA), we find no evidence of a dependence of $X_{\rm{Ly} \alpha}$ on either redshift or UV luminosity.

The C III] 1907,1909 emission doublet has been proposed as an alternative to Lyman-alpha in redshift confirmations of galaxies at z > 6 since it is not attenuated by the largely neutral intergalactic medium at these redshifts and is believed to be strong in the young, vigorously star-forming galaxies present at these early cosmic times. We present a statistical sample of 17 C III]-emitting galaxies beyond z~1.5 using 30 hour deep VLT/MUSE integral field spectroscopy covering 2 square arcminutes in the Hubble Deep Field South (HDFS) and Ultra Deep Field (UDF), achieving C III] sensitivities of ~2e-17 erg/s/cm^2 in the HDFS and ~7e-18 erg/s/cm^2 in the UDF. The rest-frame equivalent widths range from 2 to 19 Angstroms. These 17 galaxies represent ~3% of the total sample of galaxies found between 1.5 < z < 4. They also show elevated star formation rates, lower dust attenuation, and younger mass-weighted ages than the general population of galaxies at the same redshifts. Combined with deep slitless grism spectroscopy from the HST/WFC3 in the UDF, we can tie the rest-frame ultraviolet C III] emission to rest-frame optical emission lines, namely [O III] 5007, finding a strong correlation between the two. Down to the flux limits that we observe (~1e-18 erg/s/cm^2 with the grism data in the UDF), all objects with a rest-frame [O III] 4959,5007 equivalent width in excess of 250 Angstroms, the so-called Extreme Emission Line Galaxies, have detections of C III] in our MUSE data. More detailed studies of the C III]-emitting population at these intermediate redshifts will be crucial to understand the physical conditions in galaxies at early cosmic times and to determine the utility of C III] as a redshift tracer.

We report on the discovery of three especially bright candidate $z_{phot} \gtrsim 8$ galaxies. Five sources were targeted for follow-up with HST/WFC3, selected from a larger sample of 16 bright ($24.8 \lesssim H\lesssim25.5$~mag) candidate $z\gtrsim 8$ LBGs identified over the 1.6 degrees$^2$ of the COSMOS/UltraVISTA field. These were identified as Y and J dropouts by leveraging the deep (Y-to-$K_{S} \sim 25.3-24.8$~mag, $5\sigma$) NIR data from the UltraVISTA DR3 release, deep ground based optical imaging from the CFHTLS and Subaru Suprime Cam programs and Spitzer/IRAC mosaics combining observations from the SMUVS and SPLASH programs. Through the refined spectral energy distributions, which now also include new HyperSuprime Cam g, r, i, z and Y band data, we confirm that 3/5 galaxies have robust $z_{phot}\sim8.0-8.7$, consistent with the initial selection. The remaining 2/5 galaxies have a nominal $z_{phot}\sim2$. However, if we use the HST data alone, these objects have increased probability of being at $z\sim9$. Furthermore, we measure mean UV continuum slopes $\beta=-1.91\pm0.26$ for the three $z\sim8-9$ galaxies, marginally bluer than similarly luminous $z\sim4-6$ in CANDELS but consistent with previous measurements of similarly luminous galaxies at $z\sim7$. The circularized effective radius for our brightest source is $0.9\pm0.2$ kpc, similar to previous measurements for a bright $z\sim11$ galaxy and bright $z\sim7$ galaxies. Finally, enlarging our sample to include the six brightest $z\sim8$ LBGs identified over UltraVISTA (i.e., including three other sources from Labbe et al. 2017, in prep.) we estimate for the first time the volume density of galaxies at the extreme bright ($M_{UV}\sim-22$~mag) end of the $z\sim8$ UV LF. Despite this exceptional result, the still large statistical uncertainties do not allow us to discriminate between a Schechter and a double power-law form.

We present the MUSE Hubble Ultra Deep Survey, a mosaic of nine MUSE fields covering 90\% of the entire HUDF region with a 10-hour deep exposure time, plus a deeper 31-hour exposure in a single 1.15 arcmin2 field. The improved observing strategy and advanced data reduction results in datacubes with sub-arcsecond spatial resolution (0.65 arcsec at 7000 A) and accurate astrometry (0.07 arcsec rms). We compare the broadband photometric properties of the datacubes to HST photometry, finding a good agreement in zeropoint up to mAB=28 but with an increasing scatter for faint objects. We have investigated the noise properties and developed an empirical way to account for the impact of the correlation introduced by the 3D drizzle interpolation. The achieved 3 sigma emission line detection limit for a point source is 1.5 and 3.1 10-19 erg.s-1.cm-2 for the single ultra-deep datacube and the mosaic, respectively. We extracted 6288 sources using an optimal extraction scheme that takes the published HST source locations as prior. In parallel, we performed a blind search of emission line galaxies using an original method based on advanced test statistics and filter matching. The blind search results in 1251 emission line galaxy candidates in the mosaic and 306 in the ultradeep datacube, including 72 sources without HST counterparts (mAB>31). In addition 88 sources missed in the HST catalog but with clear HST counterparts were identified. This data set is the deepest spectroscopic survey ever performed. In just over 100 hours of integration time, it provides nearly an order of magnitude more spectroscopic redshifts compared to the data that has been accumulated on the UDF over the past decade. The depth and high quality of these datacubes enables new and detailed studies of the physical properties of the galaxy population and their environments over a large redshift range.

We aim at constraining the dust mass in high-redshift ($z\gtrsim 5$) galaxies using the upper limits obtained by ALMA in combination with the rest-frame UV--optical spectral energy distributions (SEDs). For SED fitting, because of degeneracy between dust extinction and stellar age, we focus on two extremes: continuous star formation (Model A) and instantaneous star formation (Model B). We apply these models to Himiko (as a representative UV-bright object) and a composite SED of $z>5$ Lyman break galaxies (LBGs). For Himiko, Model A requires a significant dust extinction, which leads to a high dust temperature $>70$ K for consistency with the ALMA upper limit. This high dust temperature puts a strong upper limit on the total dust mass $M_\mathrm{d}\lesssim 2\times 10^6$ M$_{\odot}$, and the dust mass produced per supernova (SN) $m_\mathrm{d,SN}\lesssim 0.1$ M$_{\odot}$. Such a low $m_\mathrm{d,SN}$ suggests significant loss of dust by reverse shock destruction or outflow, and implies that SNe are not the dominant source of dust at high $z$. Model B allows $M_\mathrm{d}\sim 2\times 10^7$ M$_{\odot}$ and $m_\mathrm{d,SN}\sim 0.3$ M$_{\odot}$.} We could distinguish between Models A and B if we observe Himiko at {wavelength $<$ 1.2 mm by ALMA. For the LBG sample, we obtain $M_\mathrm{d}\lesssim 2\times 10^6$ M$_{\odot}$ for a typical LBG at $z>5$, but this only puts an upper limit for $m_\mathrm{d,SN}$ as $\sim 2$ M$_{\odot}$. This clarifies the importance of observing UV-bright objects (like Himiko) to constrain the dust production by SNe.

We present an analysis of MUSE observations obtained on the massive Frontier Fields cluster Abell 2744. This new dataset covers the entire multiply-imaged region around the cluster core. We measure spectroscopic redshifts for HST-selected continuum sources together with line emitters blindly detected in the datacube. The combined catalog consists of 514 spectroscopic redshifts (with 414 new identifications), including 156 cluster members and 326 magnified background sources. We use this redshift information to perform a strong-lensing analysis of all multiple images previously found in the deep Frontier Field images, and add three new MUSE-detected multiply-imaged systems with no obvious HST counterpart. The combined strong lensing constraints include a total of 60 systems producing 188 images altogether, out of which 29 systems and 83 images are spectroscopically confirmed, making Abell 2744 one of the most well-constrained clusters to date. A parametric mass model including two cluster-scale components in the core and several group-scale substructures at larger radii accurately reproduces all the spectroscopic multiple systems, reaching an rms of 0.67" in the image plane. Overall, the large number of spectroscopic redshifts gives us a robust model and we estimate the systematics on the mass density and magnification within the cluster core to be typically ~9%.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

[abridged] The Hubble and Spitzer Space Telescope surveys of the Frontier Fields (FF) provide extremely deep images around six massive, strong-lensing clusters of galaxies. The ALMA FF survey aims to cover the same fields at 1.1mm, with maps reaching (unlensed) sensitivities of $<$70$\mu$Jy, in order to explore the properties of background dusty star-forming galaxies. We report on the multi-wavelength photometric analysis of all 12 significantly detected ($>$5$\sigma$) sources in the first three FF clusters observed by ALMA, based on data from Hubble and Spitzer, the VLT and Herschel. We measure the total photometry in all available bands and determine the photometric redshifts and the physical properties of the counterparts via SED-fitting. In particular, we carefully estimate the FIR photometry using 1.1mm priors to limit the misidentification of blended FIR counterparts, which strongly affect some flux estimates in previous FIR catalogs. We identify robust near-infrared (NIR) counterparts for all 11 sources with K$_s$ detection, the majority of which are quite red, with eight having $F814W-K_s\gtrsim 4$ and five having $F160W-[4.5]\gtrsim3$. From the FIR point of view, all our objects have $z_{phot}$$\sim$1--3, whereas based on the optical SED one object prefers a high-$z$ solution ($z\geq\ $7). Five objects among our sample have spectroscopic redshifts from the GLASS survey for which we can reproduce their SEDs with existing templates. This verification confirms the validity of our photometric redshift methodology. The mean redshift of our sample is $z_{phot}$=1.99$\pm$0.27. All 1.1mm selected objects are massive (10.0$<\log[M_{\star}(M_{\odot})]<$ 11.5), with high star formation rates ($<\log[SFR(M_{\odot}/yr)]> \approx$1.6) and high dust contents (8.1 $<\log[M_{dust} (M_{\odot})]<$8.8), consistent with previous ALMA surveys.

We present the first determination of the galaxy luminosity function (LF) at z~4, 5, 6 and 7 in the rest-frame optical at lambda_rest~900 nm (z' band). The rest-frame optical light traces the content in low-mass evolved stars (~stellar mass - Mstar), minimizing potential measurement biases for Mstar. Moreover it is less affected by nebular line emission contamination and dust attenuation, is independent of stellar population models, and can be probed up to z~8 through Spitzer/IRAC. Our analysis leverages the unique full depth Spitzer/IRAC 3.6um-to-8.0um data over the CANDELS/GOODS-N, CANDELS/GOODS-S and COSMOS/UltraVISTA fields. We find that at absolute magnitudes M_z' fainter than >-23 mag, M_z' linearly correlates with M_UV,1600. At brighter M_z', M_UV,1600 presents a turnover, suggesting that the stellar mass-to-light ratio Mstar/L_UV,1600 could be characterised by a very broad range of values at high stellar masses. Median-stacking analysis recovers a Mstar/L_z' roughly independent on M_z' for M_z'>-23 mag, but exponentially increasing at brighter magnitudes. We find that the evolution of the LF marginally prefers a pure evolution in luminosity over a pure evolution in density, with the characteristic luminosity decreasing by a factor ~5x between z~4 and z~7. Direct application of the recovered Mstar/L_z' generates stellar mass functions consistent with average measurements from the literature. Measurements of the stellar-to-halo mass ratio at fixed cumulative number density show that it is roughly constant with redshift for Mh>10^12Msun. This is also supported by the fact that the evolution of the LF at 4<z<7 can be accounted for by a rigid displacement in luminosity corresponding to the evolution of the halo mass from abundance matching.

We use the largest sample of z~6 galaxies to date from the first four Hubble Frontier Fields clusters to set constraints on the shape of the z~6 luminosity functions (LFs) to fainter than Muv=-14 mag. We quantify, for the first time, the impact of magnification uncertainties on LF results and thus provide more realistic constraints than other recent work. Our simulations reveal that for the highly-magnified sources the systematic uncertainties can become extremely large fainter than -14 mag, reaching several orders of magnitude at 95% confidence at ~-12 mag. Our new forward-modeling formalism incorporates the impact of magnification uncertainties into the LF results by exploiting the availability of many independent magnification models for the same cluster. One public magnification model is used to construct a mock high-redshift galaxy sample that is then analyzed using the other magnification models to construct a LF. Large systematic errors occur at high magnifications (mu>30) because of differences between the models. The volume densities we derive for faint (>-17 mag) sources are ~3-4x lower than one recent report and give a faint-end slope alpha=-1.92+/-0.04, which is 3.0-3.5sigma shallower (including or not including the size uncertainties, respectively). We introduce a new curvature parameter delta to model the faint end of the LF and demonstrate that the observations permit (at 68% confidence) a turn-over at z~6 in the range -15.3 to -14.2 mag, depending on the assumed lensing model. The present consideration of magnification errors and new size determinations raise doubts about previous reports regarding the form of the LF at >-14 mag. We discuss the implications of our turn-over constraints in the context of recent theoretical predictions.

We provide the first observational constraints on the sizes of the faintest galaxies lensed by the Hubble Frontier Fields (HFF) clusters. Ionizing radiation from faint galaxies likely drives cosmic reionization, and the HFF initiative provides a key opportunity to find such galaxies. Yet, we cannot really assess their ionizing emissivity without a robust measurement of their sizes, since this is key to quantifying both their prevalence and the faint-end slope to the UV luminosity function. Here we provide the first such size constraints with 2 new techniques. The first utilizes the fact that the detectability of highly-magnified galaxies as a function of shear is very dependent on a galaxy's size. Only the most compact galaxies will remain detectable in regions of high shear (vs. a larger detectable size range for low shear), a phenomenon we carefully quantify using simulations. Remarkably, however, no correlation is found between the surface density of faint galaxies and the predicted shear, using 87 faint high-magnification mu>10 z~2-8 galaxies seen behind the first 4 HFF clusters. This can only be the case if such faint (~-15 mag) galaxies have significantly smaller sizes than luminous galaxies. We constrain their half-light radii to be <~30 mas (<160-240 pc). As a 2nd size probe, we rotate and stack 26 faint high-magnification sources along the major shear axis. Less elongation is found than even for objects with an intrinsic half-light radius of 10 mas. Together these results indicate that extremely faint z~2-8 galaxies have near point-source profiles in the HFF dataset (half-light radii conservatively <30 mas and likely 5-10 mas). These results suggest smaller completeness corrections and hence much lower volume densities for faint z~2-8 galaxies and shallower faint-end slopes than have been derived in many recent studies (by factors of ~2-3 and by dalpha>~0.1-0.3).

We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged]

We have submitted to MAST the 1.5 version data release of the Hubble Legacy Fields (HLF) project covering a 25 x 25 arcmin area over the GOODS-S (ECDF-S) region from the HST archival program AR-13252. The release combines exposures from Hubble's two main cameras, the Advanced Camera for Surveys (ACS/WFC) and the Wide Field Camera 3 (WFC3/IR), taken over more than a decade between mid-2002 to the end of 2016. The HLF includes essentially all optical (ACS/WFC F435W, F606W, F775W, F814W and F850LP filters) and infrared (WFC3/ IR F098M, F105W, F125W, F140W and F160W filters) data taken by Hubble over the original CDF-S region including the GOODS-S, ERS, CANDELS and many other programs (31 in total). The data has been released at https://archive.stsci.edu/prepds/hlf/ as images with a common astrometric reference frame, with corresponding inverse variance weight maps. We provide one image per filter of WFC3/IR images at 60 mas per pixel resolution and two ACS/WFC images per filter, at both 30 and 60 mas per pixel. Since this comprehensive dataset combines data from 31 programs on the GOODS-S/CDF-S, the AR proposal identified the MAST products by the global name "Hubble Legacy Field", with this region being identified by "HLF-GOODS-S". This dataset complements that of the Frontier Fields program. The total incorporated in the HLF-GOODS-S is 5.8 Msec in 7211 exposures from 2442 orbits. This is ~70% of a HST full cycle!

The selection of high redshift sources from broad-band photometry using the Lyman-break galaxy (LBG) technique is a well established methodology, but the characterization of its contamination for the faintest sources is still incomplete. We use the optical and near-IR data from four (ultra)deep Hubble Space Telescope legacy fields to investigate the contamination fraction of LBG samples at z~5-8 selected using a colour-colour method. Our approach is based on characterizing the number count distribution of interloper sources, that is galaxies with colors similar to those of LBGs, but showing detection at wavelengths shorter than the spectral break. Without sufficient sensitivity at bluer wavelengths, a subset of interlopers may not be properly classified, and contaminate the LBG selection. The surface density of interlopers in the sky gets steeper with increasing redshift of LBG selections. Since the intrinsic number of dropouts decreases significantly with increasing redshift, this implies increasing contamination from misclassified interlopers with increasing redshift, primarily by intermediate redshift sources with unremarkable properties (intermediate ages, lack of ongoing star formation and low/moderate dust content). Using Monte Carlo simulations, we estimate that the CANDELS deep data have contamination induced by photometric scatter increasing from ~2% at z~5 to ~6% at z~8 for a typical dropout color >1 mag, with contamination naturally decreasing for a more stringent dropout selection. Contaminants are expected to be located preferentially near the detection limit of surveys, ranging from 0.1 to 0.4 contaminants per arcmin2 at J=30, depending on the field considered. This analysis suggests that the impact of contamination in future studies of z>10 galaxies needs to be carefully considered.

We present six galaxies at z~2 that show evidence of Lyman continuum (LyC) emission based on the newly acquired UV imaging of the Hubble Deep UV legacy survey (HDUV) conducted with the WFC3/UVIS camera on the Hubble Space Telescope (HST). At the redshift of these sources, the HDUV F275W images partially probe the ionizing continuum. By exploiting the HST multi-wavelength data available in the HDUV/GOODS fields, models of the UV spectral energy distributions, and detailed Monte Carlo simulations of the intergalactic medium absorption, we estimate the absolute ionizing photon escape fractions of these galaxies to be very high -- typically >60% (>13% for all sources at 90% likelihood). Our findings are in broad agreement with previous studies that found only a small fraction of galaxies to show high escape fraction. These six galaxies comprise the largest sample yet of LyC leaking candidates at z~2 whose inferred LyC flux has been cleanly observed at HST resolution. While three of our six candidates show evidence of hosting an active galactic nucleus (AGN), two of these are heavily obscured and their LyC emission appears to originate from star-forming regions rather than the central nucleus. This suggests an AGN-aided pathway for LyC escape from these sources. Extensive multi-wavelength data in the GOODS fields, especially the near-IR grism spectra from the 3D-HST survey, enable us to study the candidates in detail and tentatively test some recently proposed indirect methods to probe LyC leakage -- namely, the [OIII]/[OII] line ratio and the H$\beta-$UV slope diagram. High-resolution spectroscopic followup of our candidates will help constrain such indirect methods which are our only hope of studying $f_{esc}$ at z~5-9 in the fast-approaching era of the James Webb Space Telescope.

We study the production rate of ionizing photons of a sample of 588 H$\alpha$ emitters (HAEs) and 160 Lyman-$\alpha$ emitters (LAEs) at $z=2.2$ in the COSMOS field in order to assess the implied emissivity from galaxies, based on their UV luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/$NUV$ data, we find f$_{\rm esc} < 2.8\, (6.4)$% through median (mean) stacking. By combining the H$\alpha$ luminosity density with IGM emissivity measurements from absorption studies, we find a globally averaged $\langle$f$_{\rm esc}\rangle$ of $5.9^{+14.5}_{-4.2}$ % at $z=2.2$ if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global $\langle$f$_{\rm esc}\rangle$ at $z\approx3-5$, also ruling out a high $\langle$f$_{\rm esc}\rangle$ at $z<5$. These low escape fractions allow us to measure $\xi_{ion}$, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical $\xi_{ion} \approx 10^{24.77\pm0.04}$ Hz erg$^{-1}$ for HAEs and $\xi_{ion} \approx 10^{25.14\pm0.09}$ Hz erg$^{-1}$ for LAEs. LAEs and low mass HAEs at $z=2.2$ show similar values of $\xi_{ion}$ as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing $\xi_{ion}$ with increasing EW(H$\alpha$), $\xi_{ion}$ likely increases with redshift. This evolution alone is fully in line with the observed evolution of $\xi_{ion}$ between $z\approx2-5$, indicating a typical value of $\xi_{ion} \approx 10^{25.4}$ Hz erg$^{-1}$ in the reionization era.

We derive H{\alpha} fluxes for a large spectroscopic and photometric-redshift-selected sample of sources over GOODS-North and South in the redshift range z = 3.8-5.0 with deep HST, Spitzer/IRAC, and ground-based observations. The H{\alpha} flux is inferred based on the offset between the IRAC 3.6 {\mu}m flux and that predicted from the best-fit SED. We demonstrate that the H{\alpha} flux correlates well with dust- corrected UV star-formation rate (SFR) and therefore can serve as an independent SFR indicator. However, we also find a systematic offset in the SFR_H{\alpha}/SFR_UV ratios for z ~ 4-5 galaxies relative to local relations (assuming the same dust corrections for nebular regions and stellar light). We show that we can resolve the modest tension in the inferred SFRs by assuming bluer intrinsic UV slopes (increasing the dust correction), a rising star-formation history or assuming a low metallicity stellar population with a hard ionizing spectrum (increasing the L_H{\alpha}/SFR ratio). Using H{\alpha} as a SFR indicator, we find a normalization of the star formation main sequence in good agreement with recent SED-based determinations and also derive the SFR functions at z ~ 4-8. In addition, we assess for the first time the burstiness of star formation in z ~ 4 galaxies on <100 Myr time scales by comparing UV and H{\alpha}-based sSFRs; their one-to-one relationship argues against significantly bursty star-formation histories. Further progress will be made on these results, by incorporating new results from ALMA to constrain the dust-obscured star formation in high-redshift UV-selected samples.

Dusty star-forming galaxies are among the most prodigious systems at high redshift (z>1), characterized by high star formation rates and huge dust reservoirs. The bright end of this population has been well characterized in recent years, but considerable uncertainties remain for fainter dusty star-forming galaxies, which are responsible for the bulk of star formation at high redshift and thus play a key role in galaxy growth and evolution. In this first paper of our series, we describe our methods for finding high redshift faint dusty galaxies using millimeter observations with ALMA. We obtained ALMA 1.1 mm mosaic images for three strong-lensing galaxy clusters from the Frontier Fields survey. The 2'x2' mosaics overlap with the deep HST WFC3/IR footprints and encompass the high magnification regions of each cluster. The combination of extremely high ALMA sensitivity and the magnification power of these clusters allows us to systematically probe the sub-mJy population of dusty star-forming galaxies over a large surveyed area. We present a description of the reduction and analysis of the ALMA continuum observations for the galaxy clusters Abell 2744 (z=0.308), MACSJ0416.1-2403 (z=0.396) and MACSJ1149.5+2223 (z=0.543), for which we reach observed rms sensitivities of 55, 59 and 71 $\mu$Jy beam$^{-1}$ respectively. We detect 12 dusty star-forming galaxies at S/N>5.0 across the three clusters, all of them presenting coincidence with NIR detected counterparts in the HST images. None of the sources fall close to the lensing caustics, thus they are not strongly lensed. The observed 1.1 mm flux densities for the total sample of galaxies range from 0.41 to 2.82 mJy, with observed effective radii spanning <0.05" to 0.37"$\pm$0.21". The lensing-corrected sizes of the detected sources appear to be in the same range as those measured in brighter samples, albeit with possibly larger dispersion.

We study the molecular gas properties of high-$z$ galaxies observed in the ALMA Spectroscopic Survey (ASPECS) that targets a $\sim1$ arcmin$^2$ region in the Hubble Ultra Deep Field (UDF), a blind survey of CO emission (tracing molecular gas) in the 3mm and 1mm bands. Of a total of 1302 galaxies in the field, 56 have spectroscopic redshifts and correspondingly well-defined physical properties. Among these, 11 have infrared luminosities $L_{\rm{}IR}>10^{11}$ L$_\odot$, i.e. a detection in CO emission was expected. Out these, 7 are detected at various significance in CO, and 4 are undetected in CO emission. In the CO-detected sources, we find CO excitation conditions that are lower than typically found in starburst/SMG/QSO environments. We use the CO luminosities (including limits for non-detections) to derive molecular gas masses. We discuss our findings in context of previous molecular gas observations at high redshift (star-formation law, gas depletion times, gas fractions): The CO-detected galaxies in the UDF tend to reside on the low-$L_{\rm{}IR}$ envelope of the scatter in the $L_{\rm{}IR}-L'_{\rm{}CO}$ relation, but exceptions exist. For the CO-detected sources, we find an average depletion time of $\sim$ 1 Gyr, with significant scatter. The average molecular-to-stellar mass ratio ($M_{\rm{}H2}$/$M_*$) is consistent with earlier measurements of main sequence galaxies at these redshifts, and again shows large variations among sources. In some cases, we also measure dust continuum emission. On average, the dust-based estimates of the molecular gas are a factor $\sim$2-5$\times$ smaller than those based on CO. Accounting for detections as well as non-detections, we find large diversity in the molecular gas properties of the high-redshift galaxies covered by ASPECS.

In this paper we use ASPECS, the ALMA Spectroscopic Survey in the {\em Hubble} Ultra Deep Field (UDF) in band 3 and band 6, to place blind constraints on the CO luminosity function and the evolution of the cosmic molecular gas density as a function of redshift up to $z\sim 4.5$. This study is based on galaxies that have been solely selected through their CO emission and not through any other property. In all of the redshift bins the ASPECS measurements reach the predicted `knee' of the CO luminosity function (around $5\times10^{9}$ K km/s pc$^2$). We find clear evidence of an evolution in the CO luminosity function with respect to $z\sim 0$, with more CO luminous galaxies present at $z\sim 2$. The observed galaxies at $z\sim 2$ also appear more gas-rich than predicted by recent semi-analytical models. The comoving cosmic molecular gas density within galaxies as a function of redshift shows a factor 3-10 drop from $z \sim 2$ to $z \sim 0$ (with significant error bars), and possibly a decline at $z>3$. This trend is similar to the observed evolution of the cosmic star formation rate density. The latter therefore appears to be at least partly driven by the increased availability of molecular gas reservoirs at the peak of cosmic star formation ($z\sim2$).

We present an analysis of a deep (1$\sigma$=13 $\mu$Jy) cosmological 1.2-mm continuum map based on ASPECS, the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field. In the 1 arcmin$^2$ covered by ASPECS we detect nine sources at $>3.5\sigma$ significance at 1.2-mm. Our ALMA--selected sample has a median redshift of $z=1.6\pm0.4$, with only one galaxy detected at z$>$2 within the survey area. This value is significantly lower than that found in millimeter samples selected at a higher flux density cut-off and similar frequencies. Most galaxies have specific star formation rates similar to that of main sequence galaxies at the same epoch, and we find median values of stellar mass and star formation rates of $4.0\times10^{10}\ M_\odot$ and $\sim40~M_\odot$ yr$^{-1}$, respectively. Using the dust emission as a tracer for the ISM mass, we derive depletion times that are typically longer than 300 Myr, and we find molecular gas fractions ranging from $\sim$0.1 to 1.0. As noted by previous studies, these values are lower than using CO--based ISM estimates by a factor $\sim$2. The 1\,mm number counts (corrected for fidelity and completeness) are in agreement with previous studies that were typically restricted to brighter sources. With our individual detections only, we recover $55\pm4\%$ of the extragalactic background light (EBL) at 1.2 mm measured by the Planck satellite, and we recover $80\pm7\%$ of this EBL if we include the bright end of the number counts and additional detections from stacking. The stacked contribution is dominated by galaxies at $z\sim1-2$, with stellar masses of (1-3)$\times$10$^{10}$ M$_\odot$. For the first time, we are able to characterize the population of galaxies that dominate the EBL at 1.2 mm.

The deep, wide-area (~800-900 arcmin**2) near-infrared/WFC3/IR + Spitzer/IRAC observations over the CANDELS fields have been a remarkable resource for constraining the bright end of high redshift UV luminosity functions (LFs). However, the lack of HST 1.05-micron observations over the CANDELS fields has made it difficult to identify z~9-10 sources robustly, since such data are needed to confirm the presence of an abrupt Lyman break at 1.2 microns. We report here on the successful identification of many such z~9-10 sources from a new HST program (z9-CANDELS) that targets the highest-probability z~9-10 galaxy candidates with observations at 1.05 microns, to search for a robust Lyman-break at 1.2 microns. The potential z~9-10 candidates are preselected from the full HST, Spitzer/IRAC S-CANDELS observations, and the deepest-available ground-based optical+near-infrared observations. We identified 15 credible z~9-10 galaxies over the CANDELS fields. Nine of these galaxies lie at z~9 and 5 are new identifications. Our targeted follow-up strategy has proven to be very efficient in making use of scarce HST time to secure a reliable sample of z~9-10 galaxies. Through extensive simulations, we replicate the selection process for our sample (both the preselection and follow-up) and use it to improve current estimates for the volume density of bright z~9 and z~10 galaxies. The volume densities we find are 5(-2)(+3)x and 8(-3)(+9)x lower, respectively, than found at z~8. When compared with the best-fit evolution (i.e., dlog_{10} rho(UV)/dz=-0.29+/-0.02) in the UV luminosities densities from z~8 to z~4 integrated to 0.3L*(z=3) (-20 mag), these luminosity densities are 2.6(-0.9)(+1.5)x and 2.2(-1.1)(+2.0)x lower, respectively, than the extrapolated trends. Our new results are broadly consistent with the "accelerated evolution" scenario at z>8, as seen in many theoretical models.

We present a search for [CII] line and dust continuum emission from optical dropout galaxies at $z>6$ using ASPECS, our ALMA Spectroscopic Survey in the Hubble Ultra-Deep Field (UDF). Our observations, which cover the frequency range $212-272$ GHz, encompass approximately the range $6<z<8$ for [CII] line emission and reach a limiting luminosity of L$_{\rm [CII]}\sim$(1.6-2.5)$\times$10$^{8}$ L$_{\odot}$. We identify fourteen [CII] line emitting candidates in this redshift range with significances $>$4.5 $\sigma$, two of which correspond to blind detections with no optical counterparts. At this significance level, our statistical analysis shows that about 60\% of our candidates are expected to be spurious. For one of our blindly selected [CII] line candidates, we tentatively detect the CO(6-5) line in our parallel 3-mm line scan. None of the line candidates are individually detected in the 1.2 mm continuum. A stack of all [CII] candidates results in a tentative detection with $S_{1.2mm}=14\pm5\mu$Jy. This implies a dust-obscured star formation rate (SFR) of $(3\pm1)$ M$_\odot$ yr$^{-1}$. We find that the two highest--SFR objects have candidate [CII] lines with luminosities that are consistent with the low-redshift $L_{\rm [CII]}$ vs. SFR relation. The other candidates have significantly higher [CII] luminosities than expected from their UV--based SFR. At the current sensitivity it is unclear whether the majority of these sources are intrinsically bright [CII] emitters, or spurious sources. If only one of our line candidates was real (a scenario greatly favored by our statistical analysis), we find a source density for [CII] emitters at $6<z<8$ that is significantly higher than predicted by current models and some extrapolations from galaxies in the local universe.

We make use of deep 1.2mm-continuum observations (12.7microJy/beam RMS) of a 1 arcmin^2 region in the Hubble Ultra Deep Field to probe dust-enshrouded star formation from 330 Lyman-break galaxies spanning the redshift range z=2-10 (to ~2-3 Msol/yr at 1sigma over the entire range). Given the depth and area of ASPECS, we would expect to tentatively detect 35 galaxies extrapolating the Meurer z~0 IRX-beta relation to z>~2 (assuming T_d~35 K). However, only 6 tentative detections are found at z>~2 in ASPECS, with just three at >3sigma. Subdividing z=2-10 galaxies according to stellar mass, UV luminosity, and UV-continuum slope and stacking the results, we only find a significant detection in the most massive (>10^9.75 Msol) subsample, with an infrared excess (IRX=L_{IR}/L_{UV}) consistent with previous z~2 results. However, the infrared excess we measure from our large selection of sub-L* (<10^9.75 Msol) galaxies is 0.11(-0.42)(+0.32) and 0.14(-0.14)(+0.15) at z=2-3 and z=4-10, respectively, lying below even an SMC IRX-beta relation (95% confidence). These results demonstrate the relevance of stellar mass for predicting the IR luminosity of z>~2 galaxies. We furthermore find that the evolution of the IRX-stellar mass relationship depends on the evolution of the dust temperature. If the dust temperature increases monotonically with redshift (as (1+z)^0.32) such that T_d~44-50 K at z>=4, current results are suggestive of little evolution in this relationship to z~6. We use these results to revisit recent estimates of the z>~3 SFR density. One less obvious implication is in interpreting the high Halpha EWs seen in z~5 galaxies: our results imply that star-forming galaxies produce Lyman-continuum photons at twice the efficiency (per unit UV luminosity) as implied in conventional models. Star-forming galaxies can then reionize the Universe, even if the escape fraction is <10%.

We present the rationale for and the observational description of ASPECS: The ALMA SPECtroscopic Survey in the Hubble Ultra-Deep Field (UDF), the cosmological deep field that has the deepest multi-wavelength data available. Our overarching goal is to obtain an unbiased census of molecular gas and dust continuum emission in high-redshift (z$>$0.5) galaxies. The $\sim$1$'$ region covered within the UDF was chosen to overlap with the deepest available imaging from HST. Our ALMA observations consist of full frequency scans in band 3 (84-115 GHz) and band 6 (212-272 GHz) at approximately uniform line sensitivity ($L'_{\rm CO}\sim$2$\times$10$^{9}$ K km/s pc$^2$), and continuum noise levels of 3.8 $\mu$Jy beam$^{-1}$ and 12.7 $\mu$Jy beam$^{-1}$, respectively. The molecular surveys cover the different rotational transitions of the CO molecule, leading to essentially full redshift coverage. The [CII] emission line is also covered at redshifts $6.0<z<8.0$. We present a customized algorithm to identify line candidates in the molecular line scans, and quantify our ability to recover artificial sources from our data. Based on whether multiple CO lines are detected, and whether optical spectroscopic redshifts as well as optical counterparts exist, we constrain the most likely line identification. We report 10 (11) CO line candidates in the 3mm (1mm) band, and our statistical analysis shows that $<$4 of these (in each band) are likely spurious. Less than 1/3 of the total CO flux in the low-J CO line candidates are from sources that are not associated with an optical/NIR counterpart. We also present continuum maps of both the band 3 and band 6 observations. The data presented here form the basis of a number of dedicated studies that are presented in subsequent papers.

The James Webb Space Telescope (JWST) will be an exquisite new near-infrared observatory with imaging and multi-object spectroscopy through ESA's NIRspec instrument with its unique Micro-Shutter Array (MSA), allowing for slits to be positioned on astronomical targets by opening specific 0.002"-wide micro shutter doors. To ensure proper target acquisition, the on-sky position of the MSA needs to be verified before spectroscopic observations start. An onboard centroiding program registers the position of pre-identified guide stars in a Target Acquisition (TA) image, a short pre-spectroscopy exposure without dispersion (image mode) through the MSA with all shutters open. The outstanding issue is the availability of Galactic stars in the right luminosity range for TA relative to typical high redshift targets. We explore this here using the stars and $z\sim8$ candidate galaxies identified in the source extractor catalogs of Brightest of Reionizing Galaxies survey (BoRG[z8]), a pure-parallel program with Hubble Space Telescope Wide-Field Camera 3. We find that (a) a single WFC3 field contains enough Galactic stars to satisfy the NIRspec astrometry requirement (20 milli-arcseconds), provided its and the NIRspec TA's are $m_{lim}>24.5$ AB in WFC3 F125W, (b) a single WFC3 image can therefore serve as the pre-image if need be, (c) a WFC3 mosaic and accompanying TA image satisfy the astrometry requirement at $\sim23$ AB mag in WFC3 F125W, (d) no specific Galactic latitude requires deeper TA imaging due to a lack of Galactic stars, and (e) a depth of $\sim24$ AB mag in WFC3 F125W is needed if a guide star in the same MSA quadrant as a target is required. We take the example of a BoRG identified $z\sim8$ candidate galaxy and require a Galactic star within 20" of it. In this case, a depth of 25.5 AB in F125W is required (with $\sim$97% confidence).

The Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) enabled the search for the first galaxies observed at z ~ 8 - 11 (500 - 700 Myr after the Big Bang). To continue quantifying the number density of the most luminous galaxies (M_AB ~ -22.0) at the earliest epoch observable with HST, we search for z ~ 10 galaxies (F125W-dropouts) in archival data from the Brightest of Reionizing Galaxies (BoRG[z8]) survey, originally designed for detection of z ~ 8 galaxies (F098M-dropouts). By focusing on the deepest 293 arcmin^2 of the data along 62 independent lines of sight, we identify six z ~ 10 candidates satisfying the color selection criteria, detected at S/N > 8 in F160W with M_AB = -22.8 to -21.1 if at z = 10. Three of the six sources, including the two brightest, are in a single WFC3 pointing (~ 4 arcmin^2), suggestive of significant clustering, which is expected from bright galaxies at z ~ 10. However, the two brightest galaxies are too extended to be likely at z ~ 10, and one additional source is unresolved and possibly a brown dwarf. The remaining three candidates have m_AB ~ 26, and given the area and completeness of our search, our best estimate is a number density of sources that is marginally higher but consistent at 2{\sigma} with searches in legacy fields. Our study highlights that z ~ 10 searches can yield a small number of candidates, making tailored follow-ups of HST pure-parallel observations viable and effective.

We investigate the clustering of Lyman-break galaxies (LBGs) at $z\sim4$. Using the hierarchical galaxy formation model GALFORM, we predict, for the first time using a semi-analytical model with feedback from active galactic nuclei (AGN), the angular correlation function (ACF) of LBGs and find agreement within $3\,\sigma$ with new measurements of the ACF from surveys including the Hubble eXtreme Deep Field (XDF) and CANDELS field. Our simulations confirm the conclusion reached using independent models that although the predicted ACFs reproduce the trend of increased clustering with luminosity, the dependence is less strong than observed. We find that for the detection limits of the XDF field central LBGs at $z\sim 4$ predominantly reside in haloes of mass $\sim 10^{11}-10^{12}h^{-1}M_{\rm \odot}$ and that satellites reside in larger haloes of mass $\sim 10^{12}-10^{13}h^{-1}M_{\rm \odot}$. The model predicts fewer bright satellite LBGs at $z\sim4$ than is inferred from measurements of the ACF at small scales. By analysing the halo occupation distribution (HOD) predicted by the model, we find evidence that AGN feedback affects the HOD of central LBGs in massive haloes. This is a new high-redshift test of this important feedback mechanism. We investigate the effect of photometric errors in the observations on the ACF predictions. We find that the observational uncertainty in the galaxy luminosity reduces the clustering amplitude and that this effect increases towards faint galaxies, particularly on small scales. To compare properties of model with observed LBGs this uncertainty must be considered.

We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z~7-9 selected to have intense optical line emission by Roberts-Borsani et al. (2016). Previous follow-up has revealed Lyman-alpha in two of the four galaxies. Our new MOSFIRE observations confirm that Lyman-alpha is present in the entire sample. We detect Lyman-alpha emission in COS-zs7-1, confirming its redshift as z=7.154, and we detect Lyman-alpha in EGS-zs8-2 at z=7.477, verifying a tentative detection presented in an earlier study. The ubiquity of Lyman-alpha in this sample is puzzling given that the IGM is likely significantly neutral over 7<z<9. To investigate this result in more detail, we have initiated a campaign to target UV metal emission in the four Lyman-alpha emitters as a probe of both the radiation field and the velocity offset of Lyman-alpha. Here we present the detection of intense CIII] emission in EGS-zs8-1, a galaxy from this sample previously shown to have Lyman-alpha at z=7.73. Photoionization models indicate that an intense radiation field and low metallicity are required to reproduce the intense CIII] and optical line emission. We argue that this extreme radiation field is likely to affect the local environment, increasing the transmission of Lyman-alpha through the galaxy. Moreover, the centroid of CIII] indicates that Lyman-alpha is redshifted from the systemic value by 340 km/s. This velocity offset is larger than that seen in less luminous systems, providing an additional explanation for the transmission of Lyman-alpha emission through the IGM. Since the transmission is further enhanced by the likelihood that such systems are also situated in the densest regions with the largest ionized bubbles, the visibility of Lyman-alpha at z>7 is expected to be strongly luminosity-dependent, with the most effective transmission occurring in systems with intense star formation.

Galaxies represent one of the preferred candidate sources to drive the reionization of the universe. Even as gains are made in mapping the galaxy UV luminosity density to z>6, significant uncertainties remain regarding the conversion to the implied ionizing emissivity. The relevant unknowns are the Lyman-continuum (LyC) photon production efficiency xi_{ion} and the escape fraction f_{esc}. As we show here, the first of these unknowns is directly measureable in z=4-5 galaxies, based on the impact the Halpha line has on the observed IRAC fluxes. By computing a LyC photon production rate from the implied Halpha luminosities for a broad selection of z=4-5 galaxies and comparing this against the dust-corrected UV-continuum luminosities, we provide the first-ever direct estimates of the LyC photon production efficiency xi_{ion} for the z>~4 galaxy population. We find log_{10} xi_{ion}/[Hz/ergs] to have a mean value of 25.27_{-0.03}^{+0.03} and 25.34_{-0.02}^{+0.02} for sub-L* z=4-5 galaxies adopting Calzetti and SMC dust laws, respectively. Reassuringly, both values are consistent with standardly assumed xi_{ion}'s in reionization models, with a slight preference for higher xi_{ion}'s (by ~0.1 dex) adopting the SMC dust law. A modest ~0.03-dex increase in these estimates would result if the escape fraction for ionizing photons is non-zero and galaxies dominate the ionizing emissivity at z~4.4. High values of xi_{ion} (~25.5-25.8 dex) are derived for the bluest galaxies (beta<-2.3) in our samples, independent of dust law and consistent with results for a z=7.045 galaxy. Such elevated values of xi_{ion} would have important consequences, indicating that f_{esc} cannot be in excess of 13% unless the galaxy UV luminosity function does not extend down to -13 mag or the clumping factor is greater than 3. A low escape fraction would fit well with the low rate of LyC leakage observed at z~3.

We identify 4 unusually bright (H < 25.5) galaxies from HST and Spitzer CANDELS data with probable redshifts z ~ 7-9. These identifications include the brightest-known galaxies to date at z > 7.5. As Y-band observations are not available over the full CANDELS program to perform a standard Lyman-break selection of z > 7 galaxies, we employ an alternate strategy using deep Spitzer/IRAC data. We identify z ~ 7.1 - 9.1 galaxies by selecting z >~ 6 galaxies from the HST CANDELS data that show quite red IRAC [3.6]-[4.5] colors, indicating strong [OIII]+Hbeta lines in the 4.5 micron band. This selection strategy was validated using a modest sample for which we have deep Y-band coverage, and subsequently used to select the brightest z > 7 sources. Applying the IRAC criteria to all HST-selected optical-dropout galaxies over the full ~900 arcmin**2 of the CANDELS survey revealed four unusually bright z ~ 7.1, 7.6, 7.9 and 8.6 candidates. The median [3.6]-[4.5] color of our selected z ~ 7.1-9.1 sample is consistent with rest-frame [OIII]+Hbeta EWs of ~1500A, in the [4.5] band. Keck/MOSFIRE spectroscopy has been independently reported for two of our selected sources, showing Ly-alpha at redshifts of 7.7302+/-0.0006 and 8.683^+0.001_-0.004, respectively. We present similar Keck/MOSFIRE spectroscopy for a third selected galaxy with a probable 4.7sigma Ly-alpha line at z_spec=7.4770+/-0.0008. All three have H-band magnitudes of ~25 mag and are ~0.5 mag more luminous (M(UV) ~ -22.0) than any previously discovered z ~ 8 galaxy, with important implications for the UV LF. Our 3 brightest, highest redshift z > 7 galaxies all lie within the CANDELS EGS field, providing a dramatic illustration of the potential impact of field-to-field variance.

We present K-band spectra of rest-frame optical emission lines for 24 star-forming galaxies at z~3.2-3.7 using MOSFIRE on the Keck 1 telescope. Strong rest-frame optical [O III] and Hbeta emission lines were detected in 18 LBGs. The median flux ratio of [O III]5007 to Hbeta is 5.1+/-0.5, a factor of 5-10x higher than in local galaxies with similar stellar masses. The observed Hbeta luminosities are in good agreement with expectations from the estimated star-formation rates, and none of our sources are detected in deep X-ray stacks, ruling out significant contamination by active galactic nuclei. Combining our sample with a variety of LBGs from the literature, including 49 galaxies selected in a very similar manner, we find a high median ratio of [OIII]/Hbeta = 4.8+0.8-1.7. This high ratio seems to be an ubiquitous feature of z~3-4 LBGs, very different from typical local star-forming galaxies at similar stellar masses. The only comparable systems at z~0 are those with similarly high specific star-formation rates, though ~5x lower stellar masses. High specific star-formation rates either result in a much higher ionization parameter or other unusual conditions for the interstellar medium, which result in a much higher [OIII]/Hbeta line ratio. This implies a strong relation between a global property of a galaxy, the specific star-formation rate, and the local conditions of ISM in star-forming regions.

We have identified 274 M-type Brown Dwarfs in the Hubble Space Telescope's Wide Field Camera 3 (WFC3) pure parallel fields from the Brightest of Reionizing Galaxies (BoRG) survey for high redshift galaxies. These are near-infrared observations with multiple lines-of-sight out of our Milky Way. Using these observed M-type Brown Dwarfs we fitted a Galactic disk and halo model with a Markov chain Monte Carlo (MCMC) analysis. This model worked best with the scale length of the disk fixed at $h$ = 2.6 kpc. For the scale height of the disk, we found $z_0 = 0.29^{+0.02}_{-0.019}$ kpc and for the central number density $\rho_0 = 0.29^{+0.20}_{-0.13}$ \#/pc$^3$. For the halo we derived a flattening parameter $\kappa$ = 0.45$\pm{0.04}$ and a power-law index $p$ = 2.4$\pm{0.07}$. We found the fraction of M-type brown dwarfs in the local density that belong to the halo to be $f_{h}$ = 0.0075$^{+0.0025}_{-0.0019}$. We found no correlation between subtype of M-dwarf and any model parameters. The total number of M-type Brown Dwarfs in the disk and halo was determined to be $58.2^{+9.81}_{-6.70} \times10^{9}$. We found an upper limit for the fraction of M-type Brown Dwarfs in the halo of 7$^{+5}_{-4}$\%. The upper limit for the total Galactic Disk mass in M-dwarfs is $4.34^{+0.73}_{-0.5}\times10^{9}$ $M_{\odot}$, assuming all M-type Brown Dwarfs have a mass of $80 M_J$.

The Lyman Continuum photon production efficiency ($\xi_{\rm ion}$) is a critical ingredient for inferring the number of photons available to reionise the intergalactic medium. To estimate the theoretical production efficiency in the high-redshift Universe we couple the BlueTides cosmological hydrodynamical simulation with a range of stellar population synthesis models. We find Lyman Continuum photon production efficiencies of $\log_{10}(\xi_{\rm ion}/{\rm erg^{-1}\, Hz})\approx 25.1-25.5$ depending on the choice of stellar population synthesis model. These results are broadly consistent with recent observational constraints at high-redshift though favour a model incorporating the effects of binary evolution

We present Hubble WFC3/IR slitless grism spectra of a remarkably bright $z\gtrsim10$ galaxy candidate, GN-z11, identified initially from CANDELS/GOODS-N imaging data. A significant spectroscopic continuum break is detected at $\lambda=1.47\pm0.01~\mu$m. The new grism data, combined with the photometric data, rule out all plausible lower redshift solutions for this source. The only viable solution is that this continuum break is the Ly$\alpha$ break redshifted to ${z_\mathrm{grism}=11.09^{+0.08}_{-0.12}}$, just $\sim$400 Myr after the Big Bang. This observation extends the current spectroscopic frontier by 150 Myr to well before the Planck (instantaneous) cosmic reionization peak at z~8.8, demonstrating that galaxy build-up was well underway early in the reionization epoch at z>10. GN-z11 is remarkably and unexpectedly luminous for a galaxy at such an early time: its UV luminosity is 3x larger than L* measured at z~6-8. The Spitzer IRAC detections up to 4.5 $\mu$m of this galaxy are consistent with a stellar mass of ${\sim10^{9}~M_\odot}$. This spectroscopic redshift measurement suggests that the James Webb Space Telescope (JWST) will be able to similarly and easily confirm such sources at z>10 and characterize their physical properties through detailed spectroscopy. Furthermore, WFIRST, with its wide-field near-IR imaging, would find large numbers of similar galaxies and contribute greatly to JWST's spectroscopy, if it is launched early enough to overlap with JWST.

We present a tally of Milky Way late-type dwarf stars in 68 WFC3 pure-parallel fields (227 arcmin^2) from the Brightest of Reionizing Galaxies (BoRG) survey for high-redshift galaxies. Using spectroscopically identified M-dwarfs in two public surveys, the CANDELS and the ERS mosaics, we identify a morphological selection criterion using the half-light radius (r50), a near-infrared J-H, G-J color region where M-dwarfs are found, and a V-J relation with M-dwarf subtype. We apply this morphological selection of stellar objects, color-color selection of M-dwarfs and optical-near-infrared color subtyping to compile a catalog of 274 M-dwarfs belonging to the disk of the Milky Way with a limiting magnitude of m_F125W < 24. Based on the M-dwarfs statistics, we conclude that (a) the previously identified North/South discrepancy in M-dwarf numbers persists in our sample; there are more M-dwarfs in the Northern fields on average than in Southern ones, (b) the Milky Way's single disk scale-height for M-dwarfs is 0.3-4 kpc, depending on sub-type, (c) {\bf ERRATUM:} we present corrected coordinates (AstroPy) and distances and find a constant $z_0$=600 pc for all types. (d) a second component is visible in the vertical distribution, with a different, much higher scale-height. We report the M-dwarf component of the Sagittarius stream in one of our fields with 11 confirmed M-dwarfs, 7 of which are at the stream's distance. The dwarf scale-height and the relative low incidence in our fields of L- and T-dwarfs in these fields makes it unlikely that these stars will be interlopers in great numbers in color-selected samples of high-redshift galaxies. The relative ubiquity of M-dwarfs however will make them ideal tracers of Galactic Halo substructure with EUCLID and reference stars for JWST observations.

We present the first results and design from the redshift z~9-10 Brightest of the Reionizing Galaxies {\it Hubble Space Telescope} survey BoRG[z9-10], aimed at searching for intrinsically luminous unlensed galaxies during the first 700 Myr after the Big Bang. BoRG[z9-10] is the continuation of a multi-year pure-parallel near-IR and optical imaging campaign with the Wide Field Camera 3. The ongoing survey uses five filters, optimized for detecting the most distant objects and offering continuous wavelength coverage from {\lambda}=0.35{\mu}m to {\lambda}=1.7{\mu}m. We analyze the initial ~130 arcmin$^2$ of area over 28 independent lines of sight (~25% of the total planned) to search for z>7 galaxies using a combination of Lyman break and photometric redshift selections. From an effective comoving volume of (5-25) $times 10^5$ Mpc$^3$ for magnitudes brighter than $m_{AB}=26.5-24.0$ in the $H_{160}$-band respectively, we find five galaxy candidates at z~8.3-10 detected at high confidence (S/N>8), including a source at z~8.4 with mAB=24.5 (S/N~22), which, if confirmed, would be the brightest galaxy identified at such early times (z>8). In addition, BoRG[z9-10] data yield four galaxies with $7.3 \lesssim z \lesssim 8$. These new Lyman break galaxies with m$\lesssim26.5$ are ideal targets for follow-up observations from ground and space based observatories to help investigate the complex interplay between dark matter growth, galaxy assembly, and reionization.

Star-forming galaxies in the early universe provide us with perhaps the most natural way of explaining the reionization of the universe. Current observational results are sufficiently comprehensive, as to allow us to approximately calculate how the ionizing radiation from galaxies varies as a function of cosmic time. Important uncertainties in modeling reionization by galaxies revolve around the escape fraction and its luminosity and redshift dependence, a possible truncation of the galaxy luminosity function at the faint end, and an evolution in the production efficiency of Lyman-continuum photons with cosmic time. Despite these uncertainties, plausible choices for these parameters naturally predict a cosmic ionizing emissivity at z~6-10 whose evolution and overall normalization is in excellent agreement with that derived from current observational constraints. This strongly suggests that galaxies provide the necessary photons to reionize the universe.

Measurements of the UV-continuum slopes provide valuable information on the physical properties of galaxies forming in the early universe, probing the dust reddening, age, metal content, and even the escape fraction. While constraints on these slopes generally become more challenging at higher redshifts as the UV continuum shifts out of the Hubble Space Telescope bands (particularly at z>7), such a characterisation actually becomes abruptly easier for galaxies in the redshift window z=9.5-10.5 due to the Spitzer/IRAC 3.6um-band probing the rest-UV continuum and the long wavelength baseline between this Spitzer band and the Hubble H-band. Higher S/N constraints on the UV slope are possible at z~10 than at z=8. Here we take advantage of this opportunity and five recently discovered bright z=9.5-10.5 galaxies to present the first measurements of the mean slope for a multi-object sample of galaxy candidates at z~10. We find the measured observed slopes of these candidates are $-2.1\pm0.3\pm0.2$ (random and systematic), only slightly bluer than the measured slopes at 3.5<z<7.5 for galaxies of similar luminosities. Small increases in the stellar ages, metallicities, and dust content of the galaxy population from z~10 to z~7 could easily explain the apparent evolution in slopes.

Recent Spitzer/InfraRed Array Camera (IRAC) photometric observations have revealed that rest-frame optical emission lines contribute signficantly to the broadband fluxes of high-redshift galaxies. Specifically, in the narrow redshift range z~5.1-5.4 the [3.6]-[4.5] color is expected to be very red, due to contamination of the 4.5-micron band by the dominant Halpha line, while the 3.6-micron filter is free of nebular emission lines. We take advantage of new reductions of deep Spitzer/IRAC imaging over the GOODS-North+South fields (Labbe+2015) to obtain a clean measurement of the mean Halpha equivalent width from the [3.6]-[4.5] color in the redshift range z=5.1-5.4. The selected sources either have measured spectroscopic redshifts (13 sources) or lie very confidently in the redshift range z=5.1-5.4 based on the photometric redshift likelihood intervals (11 sources). Our z_{phot}=5.1-5.4 sample and z_{spec}=5.10-5.40 spectroscopic sample have a mean [3.6]-[4.5] color of 0.31+/-0.05 mag and 0.35+/-0.07 mag, implying a rest-frame equivalent width EW(Halpha+[NII]+[SII]) of 665+/-53 Angstroms and 707+/-74 Angstroms, respectively, for sources in these samples. These values are consistent albeit slightly higher than derived by Stark+2013 at z~4, suggesting an evolution to higher values of the Halpha+[NII]+[SII] EW at z>2. Using the 3.6micron band, which is free of emission line contamination, we perform robust SED fitting and find a median specific star formation rate of sSFR = 17_{-5}^{+2} Gyr^{-1}, 7_{-2}^{+1}x higher than at z~2. We find no strong correlation (<2sigma) between the Halpha+[NII]+[SII] EW and the stellar mass of sources. Before the advent of JWST, improvements in these results will come through an expansion of current spectroscopic samples and deeper Spitzer/IRAC measurements.

Thomson optical depth tau measurements from Planck provide new insights into the reionization of the universe. In pursuit of model-independent constraints on the properties of the ionising sources, we determine the empirical evolution of the cosmic ionizing emissivity. We use a simple two-parameter model to map out the evolution in the emissivity at z>~6 from the new Planck optical depth tau measurements, from the constraints provided by quasar absorption spectra and from the prevalence of Ly-alpha emission in z~7-8 galaxies. We find the redshift evolution in the emissivity dot{N}_{ion}(z) required by the observations to be d(log Nion)/dz=-0.15(-0.11)(+0.08), largely independent of the assumed clumping factor C_{HII} and entirely independent of the nature of the ionising sources. The trend in dot{N}_{ion}(z) is well-matched by the evolution of the galaxy UV-luminosity density (dlog_{10} rho_UV/dz=-0.11+/-0.04) to a magnitude limit >~-13 mag, suggesting that galaxies are the sources that drive the reionization of the universe. The role of galaxies is further strengthened by the conversion from the UV luminosity density rho_UV to dot(N)_{ion}(z) being possible for physically-plausible values of the escape fraction f_{esc}, the Lyman-continuum photon production efficiency xi_{ion}, and faint-end cut-off $M_{lim}$ to the luminosity function. Quasars/AGN appear to match neither the redshift evolution nor normalization of the ionizing emissivity. Based on the inferred evolution in the ionizing emissivity, we estimate that the z~10 UV-luminosity density is 8(-4)(+15)x lower than at $z~6, consistent with the observations. The present approach of contrasting the inferred evolution of the ionizing emissivity with that of the galaxy UV luminosity density adds to the growing observational evidence that faint, star-forming galaxies drive the reionization of the universe.

We report the discovery of Lyman-alpha emission (Ly$\alpha$) in the bright galaxy EGSY-2008532660 (hereafter EGSY8p7) using the MOSFIRE spectrograph at the Keck Observatory. First reported by Roberts-Borsani et al. (2015), it was selected for spectroscopic observations because of its photometric redshift ($z_{phot}=8.57^{+0.22}_{-0.43}$), apparent brightness (H$_{160}=25.26\pm0.09$) and red Spitzer/IRAC [3.6]-[4.5] color indicative of contamination by strong oxygen emission in the [4.5] band. With a total integration of $\sim$4.3 hours, our data reveal an emission line at $\simeq$11776 {\AA} which we argue is likely Ly$\alpha$ at a redshift $z_{spec}=8.683^{+0.001}_{-0.004}$, in good agreement with the photometric estimate. The line was detected independently on two nights using different slit orientations and its detection significance is $\sim7.5\sigma$. An overlapping skyline contributes significantly to the uncertainty on the total line flux although the significance of the detected line is robust to a variety of skyline-masking procedures. By direct addition and a Gaussian fit, we estimate a 95\% confidence range of 1.0--2.5$\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$, corresponding to a rest-frame equivalent width of 17--42 {\AA}. EGSY8p7 is the most distant galaxy confirmed spectroscopically to date, and the third luminous source in the EGS field beyond $z_{phot}\gtrsim7.5$ with detectable Ly$\alpha$ emission viewed at a time when the intergalactic medium is believed to be fairly neutral. Although the reionization process was probably patchy, we discuss whether luminous sources with prominent IRAC color excesses may harbor harder ionizing spectra than the dominant fainter population thereby creating earlier ionized bubbles. Further spectroscopic follow-up of such bright sources promises important insight into the early formation of galaxies.

The IRAC ultradeep field (IUDF) and IRAC Legacy over GOODS (IGOODS) programs are two ultradeep imaging surveys at 3.6{\mu}m and 4.5{\mu}m with the Spitzer Infrared Array Camera (IRAC). The primary aim is to directly detect the infrared light of reionization epoch galaxies at z > 7 and to constrain their stellar populations. The observations cover the Hubble Ultra Deep Field (HUDF), including the two HUDF parallel fields, and the CANDELS/GOODS-South, and are combined with archival data from all previous deep programs into one ultradeep dataset. The resulting imaging reaches unprecedented coverage in IRAC 3.6{\mu}m and 4.5{\mu}m ranging from > 50 hour over 150 arcmin^2, > 100 hour over 60 sq arcmin2, to 200 hour over 5 - 10 arcmin$^2$. This paper presents the survey description, data reduction, and public release of reduced mosaics on the same astrometric system as the CANDELS/GOODS-South WFC3 data. To facilitate prior-based WFC3+IRAC photometry, we introduce a new method to create high signal-to-noise PSFs from the IRAC data and reconstruct the complex spatial variation due to survey geometry. The PSF maps are included in the release, as are registered maps of subsets of the data to enable reliability and variability studies. Simulations show that the noise in the ultradeep IRAC images decreases approximately as the square root of integration time over the range 20 - 200 hours, well below the classical confusion limit, reaching 1{\sigma} point source sensitivities as faint as of 15 nJy (28.5 AB) at 3.6{\mu}m and 18 nJy (28.3 AB) at 4.5{\mu}m. The value of such ultradeep IRAC data is illustrated by direct detections of z = 7 - 8 galaxies as faint as HAB = 28.

We search the complete Hubble Frontier Field dataset of Abell 2744 and its parallel field for z~10 sources to further refine the evolution of the cosmic star-formation rate density (SFRD) at z>8. We independently confirm two images of the recently discovered triply-imaged z~9.8 source by Zitrin et al. (2014) and set an upper limit for similar z~10 galaxies with red colors of J_125-H_160>1.2 in the parallel field of Abell 2744. We utilize extensive simulations to derive the effective selection volume of Lyman-break galaxies at z~10, both in the lensed cluster field and in the adjacent parallel field. Particular care is taken to include position-dependent lensing shear to accurately account for the expected sizes and morphologies of highly-magnified sources. We show that both source blending and shear reduce the completeness at a given observed magnitude in the cluster, particularly near the critical curves. These effects have a significant, but largely overlooked, impact on the detectability of high-redshift sources behind clusters, and substantially reduce the expected number of highly-magnified sources. The detections and limits from both pointings result in a SFRD which is higher by 0.4+-0.4 dex than previous estimates at z~10 from blank fields. Nevertheless, the combination of these new results with all other estimates remain consistent with a rapidly declining SFRD in the 170 Myr from z~8 to z~10 as predicted by cosmological simulations and dark-matter halo evolution in LambdaCDM. Once biases introduced by magnification-dependent completeness are accounted for, the full six cluster and parallel Frontier Field program will be an extremely powerful new dataset to probe the evolution of the galaxy population at z>8 before the advent of the JWST.

Recently, a small sample of six $z\sim9-10$ candidates was discovered in CANDELS that are $\sim10-20\times$ more luminous than any of the previous $z\sim9-10$ galaxies identified over the HUDF/XDF and CLASH fields. We measure the sizes of these candidates to map out the size evolution of galaxies from the earliest observable times. Their sizes are also used to provide a valuable constraint on whether these unusual galaxy candidates are at high redshift. Using galfit to derive sizes from the CANDELS F160W images of these candidates, we find a mean size of 0.13$\pm$0.02" (or 0.5$\pm$0.1 kpc at $z\sim9-10$). This handsomely matches the 0.6 kpc size expected extrapolating lower redshift measurements to $z\sim9-10$, while being much smaller than the 0.59" mean size for lower-redshift interlopers to $z\sim9-10$ photometric selections lacking the blue IRAC color criterion. This suggests that source size may be an effective constraint on contaminants from $z\sim9-10$ selections lacking IRAC data. Assuming on the basis of the strong photometric evidence that the Oesch et al. 2014 sample is entirely at $z\sim9-10$, we can use this sample to extend current constraints on the size-luminosity, size-mass relation, and size evolution of galaxies to $z\sim10$. We find that the $z\sim9-10$ candidate galaxies have broadly similar sizes and luminosities as $z\sim6$-8 counterparts with star-formation-rate surface densities in the range of $\rm \Sigma_{SFR}=1-20\, M_\odot~ yr^{-1}\, kpc^{-2}$. The stellar mass-size relation is uncertain, but shallower than those inferred for lower-redshift galaxies. In combination with previous size measurements at z=4-7, we find a size evolution of $(1+z)^{-m}$ with $m=1.0\pm0.1$ for $>0.3L^*_{z=3}$ galaxies, consistent with the evolution previously derived from $2 < z < 8$ galaxies.

We present a spectroscopic redshift measurement of a very bright Lyman break galaxy at z=7.7302+-0.0006 using Keck/MOSFIRE. The source was pre-selected photometrically in the EGS field as a robust z~8 candidate with H=25.0 mag based on optical non-detections and a very red Spitzer/IRAC [3.6]-[4.5] broad-band color driven by high equivalent width [OIII]+Hbeta line emission. The Lyalpha line is reliably detected at 6.1 sigma and shows an asymmetric profile as expected for a galaxy embedded in a relatively neutral inter-galactic medium near the Planck peak of cosmic reionization. The line has a rest-frame equivalent width of EW0=21+-4 A and is extended with V_FWHM=360+90-70 km/s. The source is perhaps the brightest and most massive z~8 Lyman break galaxy in the full CANDELS and BoRG/HIPPIES surveys, having assembled already 10^(9.9+-0.2) M_sol of stars at only 650 Myr after the Big Bang. The spectroscopic redshift measurement sets a new redshift record for galaxies. This enables reliable constraints on the stellar mass, star-formation rate, formation epoch, as well as combined [OIII]+Hbeta line equivalent widths. The redshift confirms that the IRAC [4.5] photometry is very likely dominated by line emission with EW0(OIII+Hbeta)= 720-150+180 A. This detection thus adds to the evidence that extreme rest-frame optical emission lines are a ubiquitous feature of early galaxies promising very efficient spectroscopic follow-up in the future with infrared spectroscopy using JWST and, later, ELTs.

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

We report on our first set of spectroscopic Hubble Space Telescope observations of the z~11 candidate galaxy strongly lensed by the MACSJ0647.7+7015 galaxy cluster. The three lensed images are faint and we show that these early slitless grism observations are of sufficient depth to investigate whether this high-redshift candidate, identified by its strong photometric break at ~1.5 micron, could possibly be an emission line galaxy at a much lower redshift. While such an interloper would imply the existence of a rather peculiar object, we show here that such strong emission lines would clearly have been detected. Comparing realistic, two-dimensional simulations to these new observations we would expect the necessary emission lines to be detected at >5 sigma while we see no evidence for such lines in the dispersed data of any of the three lensed images. We therefore exclude that this object could be a low redshift emission line interloper, which significantly increases the likelihood of this candidate being a bona fide z~11 galaxy.

Detection of Lyman-Break Galaxies (LBGs) at high-redshift can be affected by gravitational lensing induced by foreground deflectors not only in galaxy clusters, but also in blank fields. We quantify the impact of strong magnification in the samples of $B$, $V$, $i$, $z$ $\&$ $Y$ LBGs ($4\lesssim z \lesssim8$) observed in the XDF and GOODS/CANDELS fields, by investigating the proximity of dropouts to foreground objects. We find that $\sim6\%$ of bright LBGs ($m_{H_{160}}<26$) at $z\sim7$ have been strongly lensed ($\mu>2$) by foreground objects. This fraction decreases from $\sim 3.5\%$ at $z\sim6$ to $\sim1.5\%$ at $z\sim4$. Since the observed fraction of strongly lensed galaxies is a function of the shape of the luminosity function (LF), it can be used to derive Schechter parameters, $\alpha$ and $M_{\star}$, independently from galaxy number counts. Our magnification bias analysis yields Schechter-function parameters in close agreement with those determined from galaxy counts albeit with larger uncertainties. Extrapolation of our analysis to $z\gtrsim 8$ suggests that future surveys with JSWT, WFIRST and EUCLID should find excess LBGs at the bright-end, even if there is an intrinsic exponential cutoff of number counts. Finally, we highlight how the magnification bias measurement near the detection limit can be used as probe of the population of galaxies too faint to be detected. Preliminary results using this novel idea suggest that the magnification bias at $M_{UV}\sim -18$ is not as strong as expected if $\alpha\lesssim -1.7$ extends well below the current detection limits in the XDF. At face value this implies a flattening of the LF at $M_{UV}\gtrsim-16.5$. However, selection effects and completeness estimates are difficult to quantify precisely. Thus, we do not rule out a steep LF extending to $M_{UV}\gtrsim -15$.

The remarkable HST datasets from the CANDELS, HUDF09, HUDF12, ERS, and BoRG/HIPPIES programs have allowed us to map out the evolution of the UV LF from z~10 to z~4. We have identified 5859, 3001, 857, 481, 217, and 6 galaxy candidates at z~4, z~5, z~6, z~7, z~8, and z~10, respectively from the ~1000 arcmin**2 area probed. The selection of z~4-8 galaxies over the five CANDELS fields allows us to assess the cosmic variance; the largest variations are apparent at z>=7. Our new LF determinations at z~4 and z~5 span a 6-mag baseline (-22.5 to -16 AB mag). These determinations agree well with previous estimates, but the larger samples and volumes probed here result in a more reliable sampling of >L* galaxies and allow us to reassess the form of the UV LFs. Our new LF results strengthen our earlier findings to 3.4 sigma significance for a steeper faint-end slope to the UV LF at z>4, with alpha evolving from alpha=-1.64+/-0.04 at z~4 to alpha=-2.06+/-0.13 at z~7 (and alpha = -2.02+/-0.23 at z~8), consistent with that expected from the evolution of the halo mass function. With our improved constraints at the bright end, we find less evolution in the characteristic luminosity M* over the redshift range z~4 to z~7; the observed evolution in the LF is now largely represented by changes in phi*. No evidence for a non-Schechter-like form to the z~4-8 LFs is found. A simple conditional LF model based on halo growth and evolution in the M/L ratio of halos ((1+z)**-1.5) provides a good representation of the observed evolution.

One of the most challenging aspects of studying galaxies in the z>~7 universe is the infrequent confirmation of their redshifts through spectroscopy, a phenomenon thought to occur from the increasing opacity of the intergalactic medium to Lya photons at z>6.5. The resulting redshift uncertainties inhibit the efficient search for [C II] in z~7 galaxies with sub-mm instruments such as ALMA, given their limited scan speed for faint lines. One means by which to improve the precision of the inferred redshifts is to exploit the potential impact of strong nebular emission lines on the colors of z~4-8 galaxies as observed by Spitzer/IRAC. At z~6.8, galaxies exhibit IRAC colors as blue as [3.6]-[4.5] ~-1, likely due to the contribution of [O III]+Hb to the 3.6 mum flux combined with the absence of line contamination in the 4.5 mum band. In this paper we explore the use of extremely blue [3.6]-[4.5] colors to identify galaxies in the narrow redshift window z~6.6-6.9. When combined with an I-dropout criterion, we demonstrate that we can plausibly select a relatively clean sample of z~6.8 galaxies. Through a systematic application of this selection technique to our catalogs from all five CANDELS fields, we identify 20 probable z~6.6-6.9 galaxies. We estimate that our criteria select the ~50% strongest line emitters at z~6.8 and from the IRAC colors we estimate a typical [O III]+Hb rest-frame equivalent width of 1085A for this sample. The small redshift uncertainties on our sample make it particularly well suited for follow-up studies with facilities such as ALMA.

We observed the Hubble Deep Field South with the new panoramic integral field spectrograph MUSE that we built and just commissioned at the VLT. The data cube resulting from 27 hours of integration covers one arcmin^2 field of view at an unprecedented depth with a 1 sigma emission line surface brightness limit of 1x$10^{-19}$ erg/s/cm$^2$/arcsec$^2$ and contains ~90,000 spectra. We present the combined and calibrated data cube, and we perform a first-pass analysis of the sources detected in the HDF-S imaging. We measured the redshifts of 189 sources up to a magnitude F814W = 29.5, increasing by more than an order of magnitude the number of known spectroscopic redshifts in this field. We also discovered 26 Lya emitting galaxies which are not detected in the HST WFPC2 deep broad band images. The intermediate spectral resolution of 2.3{\AA} allows us to separate resolved asymmetric Lya emitters, [O II] emitters, and C III] emitters and the large instantaneous wavelength range of 4500{\AA} helps to identify single emission lines. We also show how the three dimensional information of MUSE helps to resolve sources which are confused at ground-based image quality. Overall, secure identifications are provided for 83% of the 227 emission line sources detected in the MUSE data cube and for 32% of the 586 sources identified in the HST catalog of Casertano et al 2000. The overall redshift distribution is fairly flat to z=6.3, with a reduction between z=1.5 to 2.9, in the well-known redshift desert. The field of view of MUSE also allowed us to detect 17 groups within the field. We checked that the number counts of [O II] and Ly-a emitters are roughly consistent with predictions from the literature. Using two examples we demonstrate that MUSE is able to provide exquisite spatially resolved spectroscopic information on intermediate redshift galaxies present in the field.

We utilise a two-color Lyman-Break selection criterion to search for z~9-10 galaxies over the first 19 clusters in the CLASH program. A systematic search yields three z~9-10 candidates. While we have already reported the most robust of these candidates, MACS1149-JD, two additional z~9 candidates are also found and have H_{160}-band magnitudes of ~26.2-26.9. A careful assessment of various sources of contamination suggests <~1 contaminants for our z~9-10 selection. To determine the implications of these search results for the LF and SFR density at z~9, we introduce a new differential approach to deriving these quantities in lensing fields. Our procedure is to derive the evolution by comparing the number of z~9-10 galaxy candidates found in CLASH with the number of galaxies in a slightly lower redshift sample (after correcting for the differences in selection volumes), here taken to be z~8. This procedure takes advantage of the fact that the relative volumes available for the z~8 and z~9-10 selections behind lensing clusters are not greatly dependent on the details of the lensing models. We find that the normalization of the UV LF at z~9 is just 0.28_{-0.20}^{+0.39}\times that at z~8, ~1.4_{-0.8}^{+3.0}x lower than extrapolating z~4-8 LF results. While consistent with the evolution in the UV LF seen at z~4-8, these results marginally favor a more rapid evolution at z>8. Compared to similar evolutionary findings from the HUDF, our result is less insensitive to large-scale structure uncertainties, given our many independent sightlines on the high-redshift universe.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for $z\sim6-8$ galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at $z\sim6$, $z\sim7$, and $z\sim8$, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at $z\sim 6-8$ presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed $z \gtrsim 6$ candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the $z\sim6$, $z\sim7$, and $z\sim8$ samples, respectively, over an average area of 4.5 arcmin$^2$ per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to $\sim27$ mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at $z\sim6$, where we have the best statistics.

We measure the UV-continuum slope beta for over 4000 high-redshift galaxies over a wide range of redshifts z~4-8 and luminosities from the HST HUDF/XDF, HUDF09-1, HUDF09-2, ERS, CANDELS-N, and CANDELS-S data sets. Our new beta results reach very faint levels at z~4 (-15.5 mag: 0.006 L*(z=3)), z~5 (-16.5 mag: 0.014L*(z=3)), and z~6 and z~7 (-17 mag: 0.025 L*(z=3)). Inconsistencies between previous studies led us to conduct a comprehensive review of systematic errors and develop a new technique for measuring beta that is robust against biases that arise from the impact of noise. We demonstrate, by object-by-object comparisons, that all previous studies, including our own and those done on the latest HUDF12 dataset, suffer from small systematic errors in beta. We find that after correcting for the systematic errors (typically d(beta) ~0.1-0.2) all beta results at z~7 from different groups are in excellent agreement. The mean beta we measure for faint (-18 mag: 0.1L*(z=3)) z~4, z~5, z~6, and z~7 galaxies is -2.03+/-0.03+/-0.06 (random and systematic errors), -2.14+/-0.06+/-0.06, -2.24+/-0.11+/-0.08, and -2.30+/-0.18+/-0.13, respectively. Our new beta values are redder than we have reported in the past, but bluer than other recent results. Our previously reported trend of bluer beta's at lower luminosities is confirmed, as is the evolution to bluer beta's at high redshifts. beta appears to show only a mild luminosity dependence faintward of M(UV,AB) ~ -19 mag, suggesting that the mean beta asymptotes to ~ -2.2 to -2.4 for faint z>~4 galaxies. At z~7, the observed beta's suggest non-zero, but low dust extinction, and they agree well with values predicted in cosmological hydrodynamical simulations.

Lyman-Break Galaxy (LBG) samples observed during reionization ($z\gtrsim6$) with Hubble Space Telescope's Wide Field Camera 3 are reaching sizes sufficient to characterize their clustering properties. Using a combined catalog from the Hubble eXtreme Deep Field and CANDELS surveys, containing $N=743$ LBG candidates at z>6.5 at a mean redshift of $z=7.2$, we detect a clear clustering signal in the angular correlation function (ACF) at $\sim4\sigma$, corresponding to a real-space correlation length $r_{0}=6.7^{+0.9}_{-1.0}h^{-1}$cMpc. The derived galaxy bias $b=8.6^{+0.9}_{-1.0}$ is that of dark-matter halos of $M=10^{11.1^{+0.2}_{-0.3}}$M$_{\odot}$ at $z=7.2$, and highlights that galaxies below the current detection limit ($M_{AB}\sim-17.7$) are expected in lower-mass halos ($M\sim10^{8}-10^{10.5}$M$_{\odot}$). We compute the ACF of LBGs at $z\sim3.8-z\sim5.9$ in the same surveys. A trend of increasing bias is found from $z=3.8$ ($b\sim3.0$) to $z=7.2$ ($b\sim8.6$), broadly consistent with galaxies at fixed luminosity being hosted in dark-matter halos of similar mass at $4<z<6$, followed by a slight rise in halo masses at $z\sim7$ ($\sim2\sigma$ confidence). Separating the data at the median luminosity of the $z=7.2$ sample ($M_{UV}=-19.4$) shows higher clustering at $z=5.9$ for bright galaxies ($r_{0}=5.5^{+1.4}_{-1.5}h^{-1}$cMpc, $b=6.2^{+1.2}_{-1.5}$) compared to faint galaxies ($r_{0}=1.9^{+1.1}_{-1.0}h^{-1}$cMpc, $b=2.7\pm1.2$) implying a constant mass-to-light ratio $\frac{dlogM}{dlogL}\sim1.2^{+1.8}_{-0.8}$. A similar trend is present in the $z=7.2$ sample with larger uncertainty. Finally, our bias measurements allow us to investigate the fraction of dark-matter halos hosting UV-bright galaxies (the duty-cycle, $\epsilon_{DC}$). At $z=7.2$ values near unity are preferred, which may be explained by the shortened halo assembly time at high-redshift.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

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.

[abridged] We present the discovery of four surprisingly bright (H_160 ~ 26 - 27 mag AB) galaxy candidates at z~9-10 in the complete HST CANDELS WFC3/IR GOODS-N imaging data, doubling the number of z~10 galaxy candidates that are known, just ~500 Myr after the Big Bang. Two similarly bright sources are also detected in a systematic re-analysis of the GOODS-S data set. Three of the four galaxies in GOODS-N are significantly detected at 4.5-6.2sigma in the very deep Spitzer/IRAC 4.5 micron data, as is one of the GOODS-S candidates. Furthermore, the brightest of our candidates (at z=10.2+-0.4) is robustly detected also at 3.6 micron (6.9sigma), revealing a flat UV spectral energy distribution with a slope beta=-2.0+-0.2, consistent with demonstrated trends with luminosity at high redshift. The abundance of such luminous candidates suggests that the luminosity function evolves more significantly in phi_* than in L_* at z>~8 with a higher number density of bright sources than previously expected. Despite the discovery of these luminous candidates, the cosmic star formation rate density for galaxies with SFR >0.7 M_sun/yr shows an order-of-magnitude increase in only 170 Myr from z ~ 10 to z ~ 8, consistent with previous results. Based on the IRAC detections, we derive galaxy stellar masses at z~10, finding that these luminous objects are typically 10^9 M_sun. The cosmic stellar mass density at z~10 is log10 rho_* = 4.7^+0.5_-0.8 M_sun Mpc^-3 for galaxies brighter than M_UV~-18. The remarkable brightness, and hence luminosity, of these z~9-10 candidates highlights the opportunity for deep spectroscopy to determine their redshift and nature, demonstrates the value of additional search fields covering a wider area to understand star-formation in the very early universe, and highlights the opportunities for JWST to map the buildup of galaxies at redshifts much earlier than z~10.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

We present an analysis of $\sim$1500 H160-selected photometric galaxies detected to a limiting magnitude of 27.8 in the HUDF, using imaging from the HST WFC3/IR camera in combination with archival UV, optical, and NIR imaging. We fit photometric redshifts and stellar population estimates for all galaxies with well-determined Spitzer IRAC fluxes, allowing for the determination of the cumulative mass function within the range $1<z<6$. By selecting samples of galaxies at a constant cumulative number density, we explore the co-evolution of stellar masses and star formation rates from z$\sim$6. We find a steady increase in the SFRs of galaxies at constant number density from z$\sim$6 to z$\sim$3. The peak epoch of star formation is found to shift to later times for galaxies with increasing number densities, in agreement with the expectations from cosmic downsizing. The observed SFRs can fully account for the mass growth to z$\sim$2 amongst galaxies with cumulative number densities greater than 10$^{-3.5}$ Mpc$^{-3}$. For galaxies with a lower constant number density we find the observed stellar masses are $\sim$3 times greater than that which may be accounted for by the observed star formation alone at late times, implying that growth from mergers plays an important role at $z<2$. We additionally observe a decreasing sSFR, equivalent to approximately one order of magnitude, from z$\sim$6 to z$\sim$2 amongst galaxies with number densities less than 10$^{-3.5}$ Mpc$^{-3}$ along with significant evidence that at any redshift the sSFR is higher for galaxies at higher number density. The combination of these findings can qualitatively explain the previous findings of a sSFR plateau at high redshift. Tracing the evolution of the fraction of quiescent galaxies for samples matched in cumulative number density over this redshift range, we find no unambiguous examples of quiescent galaxies at $z>4$.

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

We measure the evolution of the specific star formation rate (sSFR = SFR / Mstellar) between redshift 4 and 6 to investigate the previous reports of "constant" sSFR at z>2. We obtain photometry on a large sample of galaxies at z~4-6 located in the GOODS-S field that have high quality imaging from HST and Spitzer. We have derived stellar masses and star formation rates (SFRs) through stellar population modeling of their spectral energy distributions (SEDs). We estimate the dust extinction from the observed UV colors. In the SED fitting process we have studied the effects of assuming a star formation history (SFH) both with constant SFR and one where the SFR rises exponentially with time. The latter SFH is chosen to match the observed evolution of the UV luminosity function. We find that neither the mean SFRs nor the mean stellar masses change significantly when the rising SFR (RSF) model is assumed instead of the constant SFR model. When focusing on galaxies with Mstar ~ 5x10^9 Msun, we find that the sSFR evolves weakly with redshift (sSFR(z) \propto (1+z)^(0.6+/-0.1) Gyr^-1), consistent with previous results and with recent estimates of the sSFR at z~2-3 using similar assumptions. We have also investigated the impact of optical emission lines on our results. We estimate that the contribution of emission lines to the rest-frame optical fluxes is only modest at z~4 and 5 but it could reach ~50% at z~6. When emission lines of this strength are taken into account, the sSFR shows somewhat higher values at high redshifts, according to the relation sSFR(z) \propto (1+z)^(1.0+/-0.1) Gyr^-1, i.e., ~2.3x higher at z~6 than at z~2. However, the observed evolution is substantially weaker than that found at z<2 or that expected from current models (which corresponds to sSFR(z) \propto (1+z)^(2.5) Gyr^-1). -abridged-

Using new ultradeep Spitzer/IRAC photometry from the IRAC Ultradeep Field program (IUDF), we investigate the stellar populations of a sample of 63 Y-dropout galaxy candidates at z~8, only 650Myr after the Big Bang. The sources are selected from HST/ACS+WFC3/IR data over the Hubble Ultra Deep Field (HUDF), two HUDF parallel fields, and wide area data over the CANDELS/GOODS-South. The new Spitzer/IRAC data increase the coverage at 3.6 micron and 4.5 micron to ~120h over the HUDF reaching depths of ~28 (AB,1 sigma). The improved depth and inclusion of brighter candidates result in direct >3 sigma IRAC detections of 20/63 sources, of which 11/63 are detected at > 5 sigma. The average [3.6]-[4.5] colors of IRAC detected galaxies at z~8 are markedly redder than those at z~7, observed only 130Myr later. The simplest explanation is that we witness strong rest-frame optical emission lines (in particular [OIII]4959,5007+Hbeta) moving through the IRAC bandpasses with redshift. Assuming that the average rest-frame spectrum is the same at both z~7 and z~8 we estimate a rest-frame equivalent width of W([OIII]4959,5007+Hbeta) = 670 (+260,-170) Angstrom contributing 0.56 (+0.16,-0.11) mag to the 4.5 micron filter at z~8. The corresponding W(Halpha) = 430 (+160,-110) Angstrom implies an average specific star formation rate of sSFR = 11 (+11,-5) Gyr^-1 and a stellar population age of 100 (+100,-50) Myr. Correcting the spectral energy distribution for the contribution of emission lines lowers the average best-fit stellar masses and mass-to-light ratios by x3, decreasing the integrated stellar mass density to rho*(z=8,MUV<-18)=0.6 (+0.4,-0.3) x 10^6 Msun Mpc^-3.

The eXtreme Deep Field (XDF) combines data from ten years of observations with the HST Advanced Camera for Surveys (ACS) and the Wide-Field Camera 3 Infra-Red (WFC3/IR) into the deepest image of the sky ever in the optical/near-IR. Since the initial observations on the Hubble Ultra-Deep Field (HUDF) in 2003, numerous surveys and programs, including supernova followup, HUDF09, CANDELS, and HUDF12 have contributed additional imaging data across the HUDF region. Yet these have never been combined and made available as one complete ultra-deep optical and near-infrared image dataset. We do so now for the eXtreme Deep Field (XDF) program. Our new and improved processing techniques provide higher quality reductions of the total dataset. All WFC3 near-IR and optical ACS data sets have been fully combined and accurately matched, resulting in the deepest imaging ever taken at these wavelengths ranging from 29.1 to 30.3 AB mag (5sigma in a 0.35" diameter aperture) in 9 filters. The gains in the optical for the four filters done in the original ACS HUDF correspond to a typical improvement of 0.15 mag, with gains of 0.25 mag in the deepest areas. Such gains are equivalent to adding ~130 to ~240 orbits of ACS data to the HUDF. Improved processing alone results in a typical gain of ~0.1 mag. Our 5sigma (optical+near-IR) SExtractor catalogs reveal about 14140 sources in the full field and about 7121 galaxies in the deepest part of the XDF (the HUDF09 region). The XDF is the deepest image of the universe ever taken, reaching, in the combined image for a flat f_nu source, to 31.2 AB mag 5sigma (32.9 at 1sigma) in a 0.35" diameter aperture.

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5\sigma_\kappa. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2{\sigma} departure above the maximum simulated halo mass at this redshift.

We provide a new observational test for a key prediction of the \Lambda CDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $\sigma$ agreement with \Lambda CDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged)

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $\sigma=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

Using stellar kinematics measurements, we investigate the growth of massive, quiescent galaxies from z~2 to today. We present X-Shooter spectra from the UV to NIR and dynamical mass measurements of 5 quiescent massive (>10^11 Msun) galaxies at z~2. This triples the sample of z>1.5 galaxies with well constrained (dsigma <100 km/s) velocity dispersion measurements. From spectral population synthesis modeling we find that these galaxies have stellar ages that range from 0.5-2 Gyr, with no signs of ongoing star formation. We measure velocity dispersions (290-450 km/s) from absorption lines and find that they are 1.6-2.1 times higher than those of galaxies in the SDSS at the same mass. Sizes are measured using GALFIT from HST-WFC3 H160 and UDS K-band images. The dynamical masses correspond well to the SED based stellar masses, with dynamical masses that are ~15% higher. We find that M_*/M_dyn may decrease slightly with time, which could reflect the increase of the dark matter fraction within an increasing effective radius. We combine different stellar kinematic studies from the literature, and examine the structural evolution from z~2 to z~0: we confirm that at fixed dynamical mass, the effective radius increases by a factor of ~2.8, and the velocity dispersion decreases by a factor of ~1.7. The mass density within one effective radius decreases by a factor of ~20, while within a fixed physical radius (1 kpc) it decreases only mildly (~2). When we allow for an evolving mass limit by selecting a population of galaxies at fixed number density, a stronger size growth with time is found (~4), velocity dispersion decreases by a factor of ~1.4, and interestingly, the mass density within 1 kpc is consistent with no evolution. This finding suggests that massive quiescent galaxies at ~2 grow inside-out, consistent with the expectations from minor mergers.

Ultra-deep WFC3/IR observations on the HUDF from the HUDF09 program revealed just one plausible z~10 candidate UDFj-39546284. UDFj-39546284 had all the properties expected of a galaxy at z~10 showing (1) no detection in the deep ACS+WFC3 imaging data blueward of the F160W band, exhibiting (2) a blue spectral slope redward of the break, and showing (3) no prominent detection in deep IRAC observations. The new, similarly deep WFC3/IR HUDF12 F160W observations over the HUDF09/XDF allow us to further assess this candidate. These observations show that this candidate, previously only detected at ~5.9 sigma in a single band, clearly corresponds to a real source. It is detected at ~5.3 sigma in the new H-band data and at ~7.8 sigma in the full 85-orbit H-band stack. Interestingly, the non-detection of the source (<1 sigma) in the new F140W observations suggests a higher redshift. Formally, the best-fit redshift of the source utilizing all the WFC3+ACS (and IRAC+K-band) observations is 11.8+/-0.3. However, we consider the z~12 interpretation somewhat unlikely, since the source would either need to be ~20x more luminous than expected or show very high-EW Ly-alpha emission (which seems improbable given the extensive neutral gas prevalent early in the reionization epoch). Lower-redshift solutions fail if only continuum models are allowed. Plausible lower-redshift solutions require that the H-band flux be dominated by line emission such as Halpha or [OIII] with extreme EWs. The tentative detection of line emission at 1.6 microns in UDFj-39546284 in a companion paper suggests that such emission may have already been found.

We present a comprehensive analysis of z>8 galaxies based on ultra-deep WFC3/IR data. We constrain the evolution of the UV luminosity function (LF) and luminosity densities from z~11 to z~8 by exploiting all the WFC3/IR data over the Hubble Ultra-Deep Field from the HUDF09 and the new HUDF12 program, in addition to the HUDF09 parallel field data, as well as wider area WFC3/IR imaging over GOODS-South. Galaxies are selected based on the Lyman Break Technique in three samples centered around z~9, z~10 and z~11, with seven z~9 galaxy candidates, and one each at z~10 and z~11. We confirm a new z~10 candidate (with z=9.8+-0.6) that was not convincingly identified in our first z~10 sample. The deeper data over the HUDF confirms all our previous z>~7.5 candidates as genuine high-redshift candidates, and extends our samples to higher redshift and fainter limits (H_160~29.8 mag). We perform one of the first estimates of the z~9 UV LF and improve our previous constraints at z~10. Extrapolating the lower redshift UV LF evolution should have revealed 17 z~9 and 9 z~10 sources, i.e., a factor ~3x and 9x larger than observed. The inferred star-formation rate density (SFRD) in galaxies above 0.7 M_sun/yr decreases by 0.6+-0.2 dex from z~8 to z~9, in good agreement with previous estimates. The low number of sources found at z>8 is consistent with a very rapid build-up of galaxies across z~10 to z~8. From a combination of all current measurements, we find a best estimate of a factor 10x decrease in the SFRD from z~8 to z~10, following (1+z)^(-11.4+-3.1). Our measurements thus confirm our previous finding of an accelerated evolution beyond z~8, and signify a rapid build-up of galaxies with M_UV<-17.7 within only ~200 Myr from z~10 to z~8, in the heart of cosmic reionization.

We present deep WFC3 grism observations of the candidate z~12 galaxy UDFj-39546284 in the HST Ultra Deep Field (UDF), by combining spectroscopic data from the 3D-HST and CANDELS surveys. The total exposure time is 40.5 ks and the spectrum covers 1.10 < {\lambda} < 1.65{\mu}m. We search for faint emission lines by cross-correlating the 2D G141 spectrum with the observed H160 morphology, a technique that is unique to slitless spectroscopy at HST resolution. We find a 2.7{\sigma} detection of an emission line at 1.599 {\mu}m---just redward of the JH140 filter---with flux 3.5 +/- 1.3x10-18 erg/s/cm2. Assuming the line is real, it contributes 110 +/- 40% of the observed H160 flux and has an observed equivalent width > 7300 {\AA}. If the line is confirmed, it could be Ly-{\alpha} at z = 12.12. However, a more plausible interpretation, given current results, could be a lower redshift feature such as [O III]{\lambda}4959,5007 at z = 2.19. We find two other 3D-HST [O III] emitters within 1000 km/s of that redshift in the GOODS-South field. Additional support for this interpretation comes from the discovery of a bright "[O III] blob" with a secure G141 grism redshift of z = 1.605. This object has a strikingly large observed equivalent width of nearly 9000 {\AA} that results in similar "dropout" colors as UDFj-39546284.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

We present a candidate for the most distant galaxy known to date with a photometric redshift z = 10.7 +0.6 / -0.4 (95% confidence limits; with z < 9.5 galaxies of known types ruled out at 7.2-sigma). This J-dropout Lyman Break Galaxy, named MACS0647-JD, was discovered as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). We observe three magnified images of this galaxy due to strong gravitational lensing by the galaxy cluster MACSJ0647.7+7015 at z = 0.591. The images are magnified by factors of ~8, 7, and 2, with the brighter two observed at ~26th magnitude AB (~0.15 uJy) in the WFC3/IR F160W filter (~1.4 - 1.7 um) where they are detected at >~ 12-sigma. All three images are also confidently detected at >~ 6-sigma in F140W (~1.2 - 1.6 um), dropping out of detection from 15 lower wavelength HST filters (~0.2 - 1.4 um), and lacking bright detections in Spitzer/IRAC 3.6um and 4.5um imaging (~3.2 - 5.0 um). We rule out a broad range of possible lower redshift interlopers, including some previously published as high redshift candidates. Our high redshift conclusion is more conservative than if we had neglected a Bayesian photometric redshift prior. Given CLASH observations of 17 high mass clusters to date, our discoveries of MACS0647-JD at z ~ 10.8 and MACS1149-JD1 at z ~ 9.6 are consistent with a lensed luminosity function extrapolated from lower redshifts. This would suggest that low luminosity galaxies could have reionized the universe. However given the significant uncertainties based on only two galaxies, we cannot yet rule out the sharp drop off in number counts at z >~ 10 suggested by field searches.

We present a study of rest-frame UV-to-optical color distributions for z~4 galaxies based on the combination of deep HST/ACS+WFC3/IR data with Spitzer/IRAC imaging. In particular, we use new, ultra-deep data from the IRAC Ultradeep Field program (IUDF10). Our sample contains a total of ~2600 galaxies selected as B-dropout Lyman Break Galaxies (LBGs) in the HUDF and one of its deep parallel fields, the HUDF09-2, as well as GOODS-North and South. This sample is used to investigate the UV continuum slopes beta and Balmer break colors (J_125-[4.5]) as a function of rest-frame optical luminosity. The [4.5] filter is chosen to avoid potential contamination by strong rest-frame optical emission lines. We find that galaxies at M_z<-21.5 (roughly corresponding to L*[z~4]) are significantly redder than their lower luminosity counterparts. The UV continuum slopes and the J_125-[4.5] colors are well correlated. The most simple explanation for this correlation is that the dust reddening at these redshifts is better described by an SMC-like extinction curve, rather than the typically assumed Calzetti reddening. After correcting for dust, we find that the galaxy population shows mean stellar population ages in the range 10^8.5 to 10^9 yr, with a dispersion of ~0.5 dex, and only weak trends as a function of luminosity. In contrast to some results from the literature, we find that only a small fraction of galaxies shows Balmer break colors which are consistent with extremely young ages, younger than 100 Myr. Under the assumption of smooth star-formation histories, this fraction is only 12-19% for galaxies at M_z<-19.75. Our results are consistent with a gradual build-up of stars and dust in galaxies at z>4, with only a small fraction of stars being formed in short, intense bursts of star-formation.

We present new z~8 galaxy candidates from a search over ~95 arcmin^2 of WFC3/IR data, tripling the previous search area for bright z~8 galaxies. Our analysis uses newly acquired WFC3/IR imaging data from the CANDELS Multi-Cycle Treasury program over the GOODS South field. These new data are combined with existing deep optical ACS imaging to search for relatively bright (M_UV < -19.5 mag) z~8 galaxy candidates using the Lyman Break technique. These new candidates are used to determine the bright end of the UV luminosity function (LF) of star-forming galaxies at z~8. To minimize contamination from lower redshift galaxies, we make full use of all optical ACS data and impose strict non-detection criteria based on an optical chi^2_opt flux measurement. In the whole search area we identify 16 candidate z~8 galaxies, spanning a magnitude range H_160 = 25.7-27.9 mag. The new data show that the UV LF is a factor ~1.7x lower at M_UV < -19.5 mag than determined from the HUDF09 and ERS data alone. Combining this new sample with the previous candidates from the HUDF09 and ERS data allows us to perform the most accurate measurement of the z~8 UV LF yet. Schechter function fits to the combined data result in a best-fit characteristic magnitude of M*(z=8) = -20.04+-0.46 mag. The faint-end slope is very steep, though quite uncertain, with alpha = -2.06+-0.32. A combination of wide area data with additional ultra-deep imaging will be required to significantly reduce the uncertainties on these parameters in the future.

We present results from deep X-ray stacking of >4000 high redshift galaxies from z~1 to 8 using the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky to date. The galaxy samples were selected using the Lyman break technique based primarily on recent HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log SFR/M* > -8.7, we expect that the observed properties of these LBGs are dominated by young stellar populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z<4 LBGs (covering ~90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with redshift (z) and SFR as log LX = 0.93 log (1+z) + 0.65 log SFR + 39.80. By comparing our observations with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution of LX/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 2-10 keV luminosity per Mpc^3) and find excellent agreement. While there are no significant stacked detections at z>5, we use our upper limits from 5<z<8 LBGs to constrain the SMBH accretion history of the Universe around the epoch of reionization.

We report the discovery of 33 Lyman break galaxy (LBG) candidates at z~8 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 arcmin^2) with Y_098 (or Y_105), J_125, and H_160 band coverage needed to search for z~8 galaxies, about three times the current CANDELS area, and slightly larger than what will be the final CANDELS wide component with Y_105 data. Our sample of 33 relatively bright Y_098-dropout galaxies have J_125 band magnitudes between 25.5 and 27.4 mag. This is the largest sample of bright (J_125 <~ 27.4) z~8 galaxy candidates presented to date. Combining our dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet galaxy luminosity function at z~8 over the widest dynamic range currently available. The combined datasets are well fitted by a Schechter function, i.e. \phi(L) = \phi_* (L/L_*)^{\alpha}\ e^{-(L/L_*)}, without evidence for an excess of sources at the bright end. At 68% confidence, we derive \phi_* = (4.3^{+3.5}_{-2.1}) \times 10^{-4} Mpc^{-3}, M_* = -20.26^{+0.29}_{-0.34}, and a very steep faint-end slope \alpha = -1.98^{+0.23}_{-0.22}. While the best-fit parameters still have a strong degeneracy, especially between \phi_* and M_*, our improved coverage at the bright end has reduced the uncertainty of the faint-end power-law slope at z~8 compared to the best previous determination at +/-0.4. With a future expansion of the BoRG survey, combined with planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, thereby confirming the key role played by small galaxies in the reionization of the universe.

For the first time, we study the evolution of the stellar mass-size relation for star-forming galaxies from z ~ 4 to z ~ 7 from Hubble-WFC3/IR camera observations of the HUDF and Early Release Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images in the rest-frame 2100 \AA \ with the stellar masses estimated from SED fitting to rest-frame optical (from Spitzer/IRAC) and UV fluxes. We show that the stellar mass-size relation of Lyman-break galaxies (LBGs) persists, at least to z ~ 5, and the median size of LBGs at a given stellar mass increases towards lower redshifts. For galaxies with stellar masses of 9.5<Log(M*/Msun)<10.4 sizes evolve as $(1+z)^{-1.20\pm0.11}$. This evolution is very similar for galaxies with lower stellar masses of 8.6<Log(M*/Msun)<9.5 which is $r_{e} \propto (1+z)^{-1.18\pm0.10}$, in agreement with simple theoretical galaxy formation models at high z. Our results are consistent with previous measurements of the LBGs mass-size relation at lower redshifts (z ~ 1-3).

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

We combine recent estimates of dust extinction at z~4-7 with UV luminosity function (LF) determinations to derive star formation rate (SFR) functions at z~4, 5, 6 and 7. SFR functions provide a more physical description of galaxy build-up at high redshift and allow for direct comparisons to SFRs at lower redshifts determined by a variety of techniques. Our SFR functions are derived from well-established z~4-7 UV LFs, UV-continuum slope trends with redshift and luminosity, and IRX-beta relations. They are well-described by Schechter relations. We extend the comparison baseline for SFR functions to z~2 by considering recent determinations of the H{\alpha} and mid-IR luminosity functions. The low-end slopes of the SFR functions are flatter than for the UV LFs, \Delta\alpha\sim+0.13, and show no clear evolution with cosmic time (z~0-7). In addition, we find that the characteristic value SFR* from the Schechter fit to SFR function exhibits consistent, and substantial, linear growth as a function of redshift from ~5 M_sun/yr at z~8, 650 Myr after the Big Bang, to ~100 M_sun/yr at z~2, ~2.5 Gyr later. Recent results at z~10, close to the onset of galaxy formation, are consistent with this trend. The uniformity of this evolution is even greater than seen in the UV LF over the redshift range z~2-8, providing validation for our dust corrections. These results provide strong evidence that galaxies build up uniformly over the first 3 Gyr of cosmic time.

Ultra-deep ACS and WFC3/IR HUDF+HUDF09 data, along with the wide-area GOODS+ERS+CANDELS data over the CDF-S GOODS field, are used to measure UV colors, expressed as the UV-continuum slope beta, of star-forming galaxies over a wide range in luminosity (0.1L*(z=3) to 2L*(z=3)) at high redshift (z~7 to z~4). Beta is measured using all ACS and WFC3/IR passbands uncontaminated by Ly_alpha and spectral breaks. Extensive tests show that our beta measurements are only subject to minimal biases. Using a different selection procedure, Dunlop et al. recently found large biases in their beta measurements. To reconcile these different results, we simulated both approaches and found that beta measurements for faint sources are subject to large biases if the same passbands are used both to select the sources and to measure beta. High-redshift galaxies show a well-defined rest-frame UV color-magnitude (CM) relationship that becomes systematically bluer towards fainter UV luminosities. No evolution is seen in the slope of the UV CM relationship in the first 1.5 Gyr, though there is a small evolution in the zero-point to redder colors from z~7 to z~4. This suggests that galaxies are evolving along a well-defined sequence in the L(UV)-color (beta) plane (a "star-forming sequence"?). Dust appears to be the principal factor driving changes in the UV color (beta) with luminosity. These new larger beta samples lead to improved dust extinction estimates at z~4-7 and confirm that the extinction is still essentially zero at low luminosities and high redshifts. Inclusion of the new dust extinction results leads to (i) excellent agreement between the SFR density at z~4-8 and that inferred from the stellar mass density, and (ii) to higher SSFRs at z>~4, suggesting the SSFR may evolve modestly (by factors of ~2) from z~4-7 to z~2.

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

The HUDF09 data are the deepest near-IR observations ever, reaching to 29.5 mag. Luminosity functions (LF) from these new HUDF09 data for 132 z\sim7 and z\sim8 galaxies are combined with new LFs for z\sim5-6 galaxies and the earlier z\sim4 LF to reach to very faint limits (<0.05 L*(z=3)). The faint-end slopes alpha are steep: -1.79+/-0.12 (z\sim5), -1.73+/-0.20 (z\sim6), -2.01+/-0.21 (z\sim7), and -1.91+/-0.32 (z\sim8). Slopes alpha\lesssim-2 lead to formally divergent UV fluxes, though galaxies are not expected to form below \sim-10 AB mag. These results have important implications for reionization. The weighted mean slope at z\sim6-8 is -1.87+/-0.13. For such steep slopes, and a faint-end limit of -10 AB mag, galaxies provide a very large UV ionizing photon flux. While current results show that galaxies can reionize the universe by z\sim6, matching the Thomson optical depths is more challenging. Extrapolating the current LF evolution to z>8, taking alpha to be -1.87+/-0.13 (the mean value at z\sim6-8), and adopting typical parameters, we derive Thomson optical depths of 0.061_{-0.006}^{+0.009}. However, this result will change if the faint-end slope alpha is not constant with redshift. We test this hypothesis and find a weak, though uncertain, trend to steeper slopes at earlier times (dalpha/dz\sim-0.05+/-0.04), that would increase the Thomson optical depths to 0.079_{-0.017}^{+0.063}, consistent with recent WMAP estimates (tau=0.088+/-0.015). It may thus not be necessary to resort to extreme assumptions about the escape fraction or clumping factor. Nevertheless, the uncertainties remain large. Deeper WFC3/IR+ACS observations can further constrain the ionizing flux from galaxies.

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 {\mu}m, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

We report the discovery of a z_{phot}=6.18^{+0.05}_{-0.07} (95% confidence level) dwarf galaxy, lensed into four images by the galaxy cluster MACS J0329.6-0211 (z_{l}=0.45). The galaxy is observed as a high-redshift dropout in HST/ACS/WFC3 CLASH and Spitzer/IRAC imaging. Its redshift is securely determined due to a clear detection of the Lyman-break in the 18-band photometry, making this galaxy one of the highest-redshift multiply-lensed objects known to date with an observed magnitude of F125W=24.00\pm0.04 AB mag for its highest-magnified image. We also present the first strong-lensing analysis of this cluster uncovering 15 additional multiply-imaged candidates of five lower-redshift sources spanning the range z_{s}~2-4. The mass model independently supports the high photometric redshift and reveals magnifications of 11.6^{+8.9}_{-4.1}, 17.6^{+6.2}_{-3.9}, 3.9^{+3.0}_{-1.7}, and 3.7^{+1.3}_{-0.2}, respectively, for the four images of the high-redshift galaxy. With this we construct a source image with a physical resolution of ~200 pc when the universe was ~0.9 Gyr old, where the z~6.2 galaxy occupies a source-plane area of approximately 2.2 kpc^{2}. Modeling the observed spectral energy distribution using population synthesis models, we find a demagnified stellar mass of ~10^{9} {M}_{sun}, subsolar metallicity (Z/Z_{sun}~0.5), low dust content (A_{V}~0.1 mag), a demagnified star formation rate (SFR) of ~3.2 {M}_{sun} yr^{-1}, and a specific SFR of ~3.4 Gyr^{-1}, all consistent with the properties of local dwarf galaxies.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

We report the discovery of seven strongly lensed Lyman break galaxy (LBG) candidates at z~7 detected in Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging of Abell 1703. The brightest candidate, called A1703-zD1, has an observed (lensed) magnitude of 24.0 AB (26 sigma) in the WFC3/IR F160W band, making it 0.2 magnitudes brighter than the z_850-dropout candidate recently reported behind the Bullet Cluster and 0.7 magnitudes brighter than the previously brightest known z~7.6 galaxy, A1689-zD1. With a cluster magnification of ~9, this source has an intrinsic magnitude of H_160 = 26.4 AB, a strong z_850 - J_125 break of 1.7 magnitudes, and a photometric redshift of z~6.7. Additionally, we find six other bright LBG candidates with H_160 band magnitudes of 24.9-26.4, photometric redshifts z~6.4 - 8.8, and magnifications mu~3-40. Stellar population fits to the ACS, WFC3/IR, and \Spitzer/IRAC data for A1703-zD1 and A1703-zD4 yield stellar masses (0.7 - 3.0) x 10^{9} M_sun, stellar ages 5-180 Myr, and star-formation rates ~7.8 M_sun/yr, and low reddening with A_V <= 0.7. The source-plane reconstruction of the exceptionally bright candidate A1703-zD1 exhibits an extended structure, spanning ~4 kpc in the z~6.7 source plane, and shows three resolved star-forming knots of radius r~0.4 kpc.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.

We use the ultra-deep HUDF09 and the deep ERS data from the HST WFC3/IR camera, along with the wide area Spitzer/IRAC data from GOODS-S to derive SEDs of star-forming galaxies from the rest-frame UV to the optical over a wide luminosity range (M_1500 ~ -21 to M_1500 ~ -18) from z ~ 7 to z ~ 4. The sample contains ~ 400 z ~ 4, ~ 120 z ~ 5, ~ 60 z ~ 6, and 36 prior z ~ 7 galaxies. Median stacking enables the first comprehensive study of very faint high-z galaxies at multiple redshifts (e.g., [3.6] = 27.4 +/- 0.1 AB mag for the M_1500 ~ -18 sources at z ~ 4). At z ~ 4 our faint median-stacked SEDs reach to ~ 0.06 L*(z=4) and are combined with recently published results at high luminosity L > L* that extend to M_1500 ~ -23. We use the observed SEDs and template fits to derive rest frame UV-to-optical colors (U - V) at all redshifts and luminosities. We find that this color does not vary significantly with redshift at a fixed luminosity. The UV-to-optical color does show a weak trend with luminosity, becoming redder at higher luminosities. This is most likely due to dust. At z >~ 5 we find blue colors [3.6]-[4.5] ~ -0.3 mag that are most likely due to rest-frame optical emission lines contributing to the flux in the IRAC filter bandpasses. The scatter across our derived SEDs remains substantial, but the results are most consistent with a lack of any evolution in the SEDs with redshift at a given luminosity. The similarity of the SEDs suggests a self-similar mode of evolution over a timespan from 0.7 Gyr to 1.5 Gyr that encompasses very substantial growth in the stellar mass density in the universe (from ~ 4x10^6 to ~ 2x10^7 Msun Mpc^-3).

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We search for z~10 galaxies over ~160 arcmin^2 of WFC3/IR data in the Chandra Deep Field South, using the public HUDF09, ERS, and CANDELS surveys, that reach to 5sigma depths ranging from 26.9 to 29.4 in H_160 AB mag. z>~9.5 galaxy candidates are identified via J_125-H_160>1.2 colors and non-detections in any band blueward of J_125. Spitzer IRAC photometry is key for separating the genuine high-z candidates from intermediate redshift (z~2-4) galaxies with evolved or heavily dust obscured stellar populations. After removing 16 sources of intermediate brightness (H_160~24-26 mag) with strong IRAC detections, we only find one plausible z~10 galaxy candidate in the whole data set, previously reported in Bouwens et al. (2011). The newer data cover a 3x larger area and provide much stronger constraints on the evolution of the UV luminosity function (LF). If the evolution of the z~4-8 LFs is extrapolated to z~10, six z~10 galaxies are expected in our data. The detection of only one source suggests that the UV LF evolves at an accelerated rate before z~8. The luminosity density is found to increase by more than an order of magnitude in only 170 Myr from z~10 to z~8. This increase is >=4x larger than expected from the lower redshift extrapolation of the UV LF. We are thus likely witnessing the first rapid build-up of galaxies in the heart of cosmic reionization. Future deep HST WFC3/IR data, reaching to well beyond 29 mag, can enable a more robust quantification of the accelerated evolution around z~10.

We use the high angular resolution in the near-infrared of the WFC3 on HST to determine YHVz color-color selection criteria to identify and characterize 1.5<z<3.5 galaxies in the HUDF09 and ERS (GOODS-South) fields. The WFC3 NIR images reveal galaxies at these redshifts that were undetected in the rest-frame UV HUDF/GOODS images, as well as true centers and regular disks in galaxies classified as highly irregular in rest-frame UV light. Across the 1.5<z<2.15 redshift range, regular disks are unveiled in the WFC3 images of ~25% of both intermediate and high mass galaxies, i.e., above 10^10 Msun. Meanwhile, galaxies maintaining diffuse and/or irregular morphologies in the rest-frame optical light---i.e., not yet dynamically settled---at these epochs are almost entirely restricted to masses below 10^11 Msun. In contrast at 2.25 < z < 3.5 these diffuse and/or irregular structures overwhelmingly dominate the morphological mix in both the intermediate and high mass regimes, while no regular disks, and only a small fraction (25%) of smooth spheroids, are evident above 10^11 Msun. Strikingly, by 1.5 < z < 2.25 roughly 2 out of every 3 galaxies at the highest masses are spheroids. In our small sample, the fraction of star-forming galaxies at these mass scales decreases concurrently from ~60% to ~5%. If confirmed, this indicates that z~2 is the epoch of both the morphological transformation and quenching of star-formation which assemble the first substantial population of massive ellipticals.

We present uncontaminated rest-frame u - R colors of 78 X-ray-selected AGN hosts at 0.5 < z < 1.5 in the Chandra Deep Fields measured with HST/ACS/NICMOS and VLT/ISAAC imaging. We also present spatially-resolved NUV - R color gradients for a subsample of AGN hosts imaged by HST/WFC3. Integrated, uncorrected photometry is not reliable for comparing the mean properties of soft and hard AGN host galaxies at z ~ 1 due to color contamination from point-source AGN emission. We use a cloning simulation to develop a calibration between concentration and this color contamination and use this to correct host galaxy colors. The mean u - R color of the unobscured/soft hosts beyond ~6 kpc is statistically equivalent to that of the obscured/hard hosts (the soft sources are 0.09 +/- 0.16 magnitudes bluer). Furthermore, the rest-frame V - J colors of the obscured and unobscured hosts beyond ~6 kpc are statistically equivalent, suggesting that the two populations have similar distributions of dust extinction. For the WFC3/IR sample, the mean NUV - R color gradients of unobscured and obscured sources differ by less than ~0.5 magnitudes for r > 1.1 kpc. These three observations imply that AGN obscuration is uncorrelated with the star formation rate beyond ~1 kpc. These observations favor a unification scenario for intermediate-luminosity AGNs in which obscuration is determined geometrically. Scenarios in which the majority of intermediate-luminosity AGN at z ~ 1 are undergoing rapid, galaxy-wide quenching due to AGN-driven feedback processes are disfavored.

We identify 73 z~7 and 59 z~8 candidate galaxies in the reionization epoch, and use this large 26-29.4 AB mag sample of galaxies to derive very deep luminosity functions to <-18 AB mag and the star formation rate density at z~7 and z~8. The galaxy sample is derived using a sophisticated Lyman-Break technique on the full two-year WFC3/IR and ACS data available over the HUDF09 (~29.4 AB mag, 5 sigma), two nearby HUDF09 fields (~29 AB mag, 14 arcmin) and the wider area ERS (~27.5 AB mag) ~40 arcmin**2). The application of strict optical non-detection criteria ensures the contamination fraction is kept low (just ~7% in the HUDF). This very low value includes a full assessment of the contamination from lower redshift sources, photometric scatter, AGN, spurious sources, low mass stars, and transients (e.g., SNe). From careful modelling of the selection volumes for each of our search fields we derive luminosity functions for galaxies at z~7 and z~8 to <-18 AB mag. The faint-end slopes alpha at z~7 and z~8 are uncertain but very steep at alpha = -2.01+/-0.21 and alpha=-1.91+/-0.32, respectively. Such steep slopes contrast to the local alpha<~-1.4 and may even be steeper than that at z~4 where alpha=-1.73+/-0.05. With such steep slopes (alpha<~-1.7) lower luminosity galaxies dominate the galaxy luminosity density during the epoch of reionization. The star formation rate densities derived from these new z~7 and z~8 luminosity functions are consistent with the trends found at later times (lower redshifts). We find reasonable consistency, with the SFR densities implied from reported stellar mass densities, being only ~40% higher at z<7. This suggests that (1) the stellar mass densities inferred from the Spitzer IRAC photometry are reasonably accurate and (2) that the IMF at very high redshift may not be very different from that at later times.

We use ultradeep HST WFC3/IR imaging of the HUDF to investigate the rest-frame optical morphologies of a mass-selected sample of galaxies at z~2. We find a large variety of galaxy morphologies, ranging from large, blue, disk-like galaxies to compact, red, early-type galaxies. We derive rest-frame u-g color profiles for these galaxies and show that most z~2 galaxies in our sample have negative color gradients such that their cores are red. Although these color gradients may partly be caused by radial variations in dust content, they point to the existence of older stellar populations in the centers of z~2 galaxies. This result is consistent with an "inside-out" scenario of galaxy growth. We find that the median color gradient is fairly constant with redshift: [Delta(u-g_{rest})/Delta(log r)]_{median} = -0.46, -0.44 and -0.49 for z~2, z~1 and z=0, respectively. Using structural parameters derived from surface brightness profiles we confirm that at z~2 galaxy morphology correlates well with specific star formation rate. At the same mass, star forming galaxies have larger effective radii, bluer rest-frame u-g colors and lower Sersic indices than quiescent galaxies. These correlations are very similar to those at lower redshift, suggesting that the relations that give rise to the Hubble sequence at z=0 are already in place for massive galaxies at this early epoch.

We derive stellar masses from SED fitting to rest-frame optical and UV fluxes for 401 star-forming galaxies at z 4, 5, and 6 from Hubble-WFC3/IR observations of the ERS combined with the deep GOODS-S Spitzer/IRAC data (and include a previously-published z 7 sample). A mass-luminosity relation with strongly luminosity-dependent M/Luv ratios is found for the largest sample (299 galaxies) at z 4. The relation M \propto L_{UV,1500}^(1.7+/-0.2) has a well-determined intrinsic sample variance of 0.5 dex. This relation is also consistent with the more limited samples at z 5-7. This z 4 mass-luminosity relation, and the well-established faint UV luminosity functions at z 4-7, are used to derive galaxy mass functions (MF) to masses M~10^8 at z 4-7. A bootstap approach is used to derive the MFs to account for the large scatter in the M--Luv relation and the luminosity function uncertainties, along with an analytical crosscheck. The MFs are also corrected for the effects of incompleteness. The incompleteness-corrected MFs are steeper than previously found, with slopes \alpha_M-1.4 to -1.6 at low masses. These slopes are, however, still substantially flatter than the MFs obtained from recent hydrodynamical simulations. We use these MFs to estimate the stellar mass density (SMD) of the universe to a fixed M_{UV,AB}<-18 as a function of redshift and find a SMD growth \propto(1+z)^{-3.4 +/-0.8} from z 7 to z 4. We also derive the SMD from the completeness-corrected MFs to a mass limit M~10^{8} Msun. Such completeness-corrected MFs and the derived SMDs will be particularly important for model comparisons as future MFs reach to lower masses.

Searches for very-high-redshift galaxies over the past decade have yielded a large sample of more than 6,000 galaxies existing just 900-2,000 million years (Myr) after the Big Bang (redshifts 6 > z > 3; ref. 1). The Hubble Ultra Deep Field (HUDF09) data have yielded the first reliable detections of z ~ 8 galaxies that, together with reports of a gamma-ray burst at z ~ 8.2 (refs 10, 11), constitute the earliest objects reliably reported to date. Observations of z ~ 7-8 galaxies suggest substantial star formation at z > 9-10. Here we use the full two-year HUDF09 data to conduct an ultra-deep search for z ~ 10 galaxies in the heart of the reionization epoch, only 500 Myr after the Big Bang. Not only do we find one possible z ~ 10 galaxy candidate, but we show that, regardless of source detections, the star formation rate density is much smaller (~10%) at this time than it is just ~200 Myr later at z ~ 8. This demonstrates how rapid galaxy build-up was at z ~ 10, as galaxies increased in both luminosity density and volume density from z ~ 8 to z ~ 10. The 100-200 Myr before z ~ 10 is clearly a crucial phase in the assembly of the earliest galaxies.

We present the first results on the search for very bright (M_AB -21) galaxies at redshift z~8 from the Brightest of Reionizing Galaxies (BoRG) survey. BoRG is a Hubble Space Telescope Wide Field Camera 3 pure-parallel survey that is obtaining images on random lines of sight at high Galactic latitudes in four filters (F606W, F098M, F125W, F160W), with integration times optimized to identify galaxies at z>7.5 as F098M-dropouts. We discuss here results from a search area of approximately 130 arcmin^2 over 23 BoRG fields, complemented by six other pure-parallel WFC3 fields with similar filters. This new search area is more than two times wider than previous WFC3 observations at z~8. We identify four F098M-dropout candidates with high statistical confidence (detected at greater than 8sigma confidence in F125W). These sources are among the brightest candidates currently known at z~8 and approximately ten times brighter than the z=8.56 galaxy UDFy-38135539. They thus represent ideal targets for spectroscopic followup observations and could potentially lead to a redshift record, as our color selection includes objects up to z~9. However, the expected contamination rate of our sample is about 30% higher than typical searches for dropout galaxies in legacy fields, such as the GOODS and HUDF, where deeper data and additional optical filters are available to reject contaminants.

We report evidence of a fully established galaxy cluster at z=2.07, consisting of a ~20sigma overdensity of red, compact spheroidal galaxies spatially coinciding with extended X-ray emission detected with XMM-Newton. We use VLT VIMOS and FORS2 spectra and deep Subaru, VLT and Spitzer imaging to estimate the redshift of the structure from a prominent z=2.07 spectroscopic redshift spike of emission-line galaxies, concordant with the accurate 12-band photometric redshifts of the red galaxies. Using NICMOS and Keck AO observations, we find that the red galaxies have elliptical morphologies and compact cores. While they do not form a tight red sequence, their colours are consistent with that of a >1.3$~Gyr population observed at z~2.1. From an X-ray luminosity of .2*10^43 erg s^-1 and the stellar mass content of the red galaxy population, we estimate a halo mass of 5.3-8*10^13 Msun, comparable to the nearby Virgo cluster. These properties imply that this structure could be the most distant, mature cluster known to date and that X-ray luminous, elliptical-dominated clusters are already forming at substantially earlier epochs than previously known.

We present UV luminosity functions (LFs) at 1500 Angstrom derived from the HST Early Release Science WFC3/UVIS data acquired over ~50 arcmin^2 of the GOODS-South field. The LFs are determined over the entire redshift range z=0.75-2.5 using two methods, similar to those used at higher redshifts for Lyman Break Galaxies (LBGs): (1) 13-band UV+optical+NIR photometric redshifts to study galaxies in the range z=0.5-2 in three bins of dz=0.5, and (2) dropout samples in three redshift windows centered at z~1.5, z~1.9, and z~2.5. The characteristic luminosity dims by ~1.5 mag from z=2.5 to z=0.75, consistent with earlier work. However, the other Schechter function parameters, the faint-end slope and the number density, are found to be remarkably constant over the range z=0.75-2.5. Using these LF determinations we find the UV luminosity density to increase by ~1.4 dex according to (1+z)^{2.58+-0.15} from z~0 to its peak at z~2.5. Strikingly, the inferred faint-end slopes for our LFs are all steeper than alpha=-1.5, in agreement with higher-redshift LBG studies. Since the faint-end slope in the local universe is found to be much flatter with alpha~=-1.2, this poses the question as to when and how the expected flattening occurs. Despite relatively large uncertainties, our data suggest alpha~=-1.7 at least down to z~1. These new results from such a shallow early dataset demonstrate very clearly the remarkable potential of WFC3/UVIS for the thorough characterization of galaxy evolution over the full redshift range z~0.5 to z~3.

We use ~88 arcmin**2 of deep (>~26.5 mag at 5 sigma) NICMOS data over the two GOODS fields and the HDF South to conduct a search for bright z>~7 galaxy candidates. This search takes advantage of an efficient preselection over 58 arcmin**2 of NICMOS H-band data where only plausible z>~7 candidates are followed up with NICMOS J-band observations. ~248 arcmin**2 of deep ground-based near-infrared data (>~25.5 mag, 5 sigma) is also considered in the search. In total, we report 15 z-dropout candidates over this area -- 7 of which are new to these search fields. Two possible z~9 J-dropout candidates are also found, but seem unlikely to correspond to z~9 galaxies. The present z~9 search is used to set upper limits on the prevalence of such sources. Rigorous testing is undertaken to establish the level of contamination of our selections by photometric scatter, low mass stars, supernovae (SNe), and spurious sources. The estimated contamination rate of our z~7 selection is ~24%. Through careful simulations, the effective volume available to our z>~7 selections is estimated and used to establish constraints on the volume density of luminous (L*(z=3), or -21 mag) galaxies from these searches. We find that the volume density of luminous star-forming galaxies at z~7 is 13_{-5}^{+8}x lower than at z~4 and >25x lower (1 sigma) at z~9 than at z~4. This is the most stringent constraint yet available on the volume density of >~L* galaxies at z~9. The present wide-area, multi-field search limits cosmic variance to <20%. The evolution we find at the bright end of the UV LF is similar to that found from recent Subaru Suprime-Cam, HAWK-I or ERS WFC3/IR searches. The present paper also includes a complete summary of our final z~7 z-dropout sample (18 candidates) identified from all NICMOS observations to date (over the two GOODS fields, the HUDF, galaxy clusters).

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 multi-wavelength study of galaxy populations in the core of the massive, X-ray luminous cluster XMMU J2235 at z=1.39, based on VLT and HST optical and near-infrared photometry. Luminosity functions in the z, H, and Ks bands show a faint-end slope consistent with being flat, and a characteristic magnitude M* close to passive evolution predictions of M* of local massive clusters, with a formation redshift z>2. The color-magnitude and color-mass diagrams show evidence of a tight red sequence of massive galaxies, with overall old stellar populations, generally early-type morphology, typically showing early-type spectral features and rest-frame far-UV emission consistent with very low star formation rates (SFR<0.2Msun/yr). Star forming spectroscopic members, with SFRs of up to ~100Msun/yr, are all located at clustercentric distances >~250kpc, with the central cluster region already appearing effectively quenched. Massive galaxies in the core of this cluster appear to be in an advanced evolutionary stage in terms of both star formation and mass assembly. The high-mass end of the galaxy stellar mass function is essentially already in place, and the stellar mass fraction estimated within r500 (~1%, Kroupa IMF) is already similar to that of local massive clusters. On the other hand, morphological analysis of the massive red sequence galaxies suggests that they are smaller than similarly massive local early-types. While possibly affected by systematics and biases, this result might imply that, in spite of the overall early assembly of these sources, their evolution is not complete, and processes like minor (and likely dry) merging might still shape their structural properties to resemble those of their local counterparts, without substantially affecting their stellar mass or host stellar populations.[abridged]

(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 investigate the Spitzer/IRAC properties of 36 z~7 z-dropout galaxies and 3 z~8 Y-dropout galaxies derived from deep/wide-area WFC3/IR data of the Early Release Science, the ultradeep HUDF09, and wide-area NICMOS data. We fit stellar population synthesis models to the SEDs to derive mean redshifts, stellar masses, and ages. The z~7 galaxies are best characterized by substantial ages (>100 Myr) and M/L_V~0.2. The main trend with decreasing luminosity is that of bluing of the far-UV slope from beta~-2.0 to beta~-3.0. This can be explained by decreasing metallicity, except for the lowest luminosity galaxies (0.1 L*{z=3}), where low metallicity and smooth SFHs fail to match the blue far-UV and moderately red H-[3.6] color. Such colors may require episodic SFHs with short periods of activity and quiescence ("on-off" cycles) and/or a contribution from emission lines. The stellar mass of our sample of z~7 star forming galaxies correlates with SFR according to log M* = 8.70 (+-0.09) + 1.06(+-0.10) log SFR, implying star formation may have commenced at z>10. No galaxies are found with SFRs much higher or lower than the past averaged SFR suggesting that the typical star formation timescales are probably a substantial fraction of the Hubble time. We report the first IRAC detection of Y-dropout galaxies at z~8. The average rest-frame U-V=0.3 (AB) of the 3 galaxies are similar to faint z~7 galaxies, implying similar M/L. The stellar mass density to M_{UV,AB}<-18 is rho*(z=8) = 1.8^{+0.7}_{-1.0} x 10^6 M_sun Mpc^-3, following log rho*(z)= 10.6(+-0.6) - 4.4(+-0.7) log(1+z) [M_sun Mpc^-3] over 3<z< 8.

We present very deep Wide Field Camera 3 (WFC3) photometry of a massive, compact galaxy located in the Hubble Ultra Deep Field. This quiescent galaxy has a spectroscopic redshift z=1.91 and has been identified as an extremely compact galaxy by Daddi et al. 2005. We use new H-F160W imaging data obtained with Hubble Space Telescope/WFC3 to measure the deconvolved surface brightness profile to H = 28 mag arcsec**-2. We find that the surface brightness profile is well approximated by an n=3.7 Sersic profile. Our deconvolved profile is constructed by a new technique which corrects the best-fit Sersic profile with the residual of the fit to the observed image. This allows for galaxy profiles which deviate from a Sersic profile. We determine the effective radius of this galaxy: r_e=0.42 +- 0.14 kpc in the observed H-F160W-band. We show that this result is robust to deviations from the Sersic model used in the fit. We test the sensitivity of our analysis to faint "wings" in the profile using simulated galaxy images consisting of a bright compact component and a faint extended component. We find that due to the combination of the WFC3 imaging depth and our method's sensitivity to extended faint emission we can accurately trace the intrinsic surface brightness profile, and that we can therefore confidently rule out the existence of a faint extended envelope around the observed galaxy down to our surface brightness limit. These results confirm that the galaxy lies a factor of 10 off from the local mass-size relation.

The new Wide Field Camera 3/IR observations on the Hubble Ultra-Deep Field started investigating the properties of galaxies during the reionization epoch. To interpret these observations, we present a novel approach inspired by the conditional luminosity function method. We calibrate our model to observations at z=6 and assume a non-evolving galaxy luminosity versus halo mass relation. We first compare model predictions against the luminosity function measured at z=5 and z=4. We then predict the luminosity function at z>=7 under the sole assumption of evolution in the underlying dark-matter halo mass function. Our model is consistent with the observed z>6.5 galaxy number counts in the HUDF survey and suggests a possible steepening of the faint-end slope of the luminosity function: alpha(z>8)< -1.9 compared to alpha=-1.74 at z=6. Although we currently see only the brightest galaxies, a hidden population of lower luminosity objects (L/L_{*}> 10^{-4}) might provide >75% of the total reionizing flux. Assuming escape fraction f_{esc}~0.2, clumping factor C~5, top heavy-IMF and low metallicity, galaxies below the detection limit produce complete reionization at z>8. For solar metallicity and normal stellar IMF, reionization finishes at z>6, but a smaller C/f_{esc} is required for an optical depth consistent with the WMAP measurement. Our model highlights that the star formation rate in sub-L_* galaxies has a quasi-linear relation to dark-matter halo mass, suggesting that radiative and mechanical feedback were less effective at z>6 than today.

We use a robust sample of 11 z~7 galaxies (z-dropouts) to estimate the stellar mass density of the universe when it was only ~750 Myr old. We combine the very deep optical to near-Infrared photometry from the HST ACS and NICMOS cameras with mid-Infrared Spitzer IRAC imaging available through the GOODS program. After carefully removing the flux from contaminating foreground sources we have obtained reliable photometry in the 3.6 and 4.5 micron IRAC channels. The spectral shapes of these sources, including their rest frame optical colors, strongly support their being at z~7 with a mean photometric redshift of <z>=7.2+/-0.5. We use Bruzual & Charlot (2003) synthetic stellar population models to constrain their stellar masses and star formation histories. We find stellar masses that range over 0.1 -12x10^9 M_sol and average ages from 20 Myr to up to 425 Myr with a mean of ~300 Myr, suggesting that in some of these galaxies most of the stars were formed at z>8 (and probably at z>~10). The best fits to the observed SEDs are consistent with little or no dust extinction, in agreement with recent results at z~4-8. The star formation rates (SFR) are in the range from 5-20 M_sol/yr. From this sample we measure a stellar mass density of 6.6_{-3.3}^{+5.4}x10^5 M_sol/Mpc^3 to a limit of M_{UV,AB}<-20 (or 0.4 L*(z=3)). Combined with a fiducial lower limit for their ages (80 Myr) this implies a maximum SFR density of 0.008 M_sol/yr/Mpc^3. This is well below the critical level needed to reionize the universe at z~8 using standard assumptions. However, this result is based on luminous sources (>L*) and does not include the dominant contribution of the fainter galaxies. Strikingly, we find that the specific SFR is constant from z~7 to z~2 but drops substantially at more recent times.

We use the ultra-deep WFC3/IR data over the HUDF and the Early Release Science WFC3/IR data over the CDF-South GOODS field to quantify the broadband spectral properties of candidate star-forming galaxies at z~7. We determine the UV-continuum slope beta in these galaxies, and compare the slopes with galaxies at later times to measure the evolution in beta. For luminous L*(z=3) galaxies, we measure a mean UV-continuum slope beta of -2.0+/-0.2, which is comparable to the beta~-2 derived at similar luminosities at z~5-6. However, for the lower luminosity 0.1L*(z=3) galaxies, we measure a mean beta of -3.0+/-0.2. This is substantially bluer than is found for similar luminosity galaxies at z~4, just 800 Myr later, and even at z~5-6. In principle, the observed beta of -3.0 can be matched by a very young, dust-free stellar population, but when nebular emission is included the expected beta becomes >~-2.7. To produce these very blue beta's (i.e., beta~-3), extremely low metallicities and mechanisms to reduce the red nebular emission are likely required. For example, a large escape fraction (i.e., f_{esc}>~0.3) could minimize the contribution from this red nebular emission. If this is correct and the escape fraction in faint z~7 galaxies is >~0.3, it may help to explain how galaxies reionize the universe.

We utilize the newly-acquired, ultra-deep WFC3/IR observations over the HUDF to search for star-forming galaxies at z~8-8.5, only 600 million years from recombination, using a Y_{105}-dropout selection. The new 4.7 arcmin**2 WFC3/IR observations reach to ~28.8 AB mag (5 sigma) in the Y_{105}J_{125}H_{160} bands. These remarkable data reach ~1.5 AB mag deeper than the previous data over the HUDF, and now are an excellent match to the HUDF optical ACS data. For our search criteria, we use a two-color Lyman-Break selection technique to identify z~8-8.5 Y_{105}-dropouts. We find 5 likely z~8-8.5 candidates. The sources have H_{160}-band magnitudes of ~28.3 AB mag and very blue UV-continuum slopes, with a median estimated beta of <~-2.5 (where f_{\lambda}\propto \lambda^{\beta}). This suggests that z~8 galaxies are not only essentially dust free but also may have very young ages or low metallicities. The observed number of Y_{105}-dropout candidates is smaller than the 20+/-6 sources expected assuming no evolution from z~6, but is consistent with the 5 expected extrapolating the Bouwens et al. 2008 LF results to z~8. These results provide evidence that the evolution in the LF seen from z~7 to z~3 continues to z~8. The remarkable improvement in the sensitivity of WFC3/IR has enabled HST to cross a threshold, revealing star-forming galaxies at z~8-9.

We study the Spitzer Infrared Array Camera (IRAC) mid-infrared (rest-frame optical) fluxes of 14 newly WFC3/IR-detected z=7 z_{850}-dropout galaxies and 5 z=8 Y_{105}-dropout galaxies. The WFC3/IR depth and spatial resolution allow accurate removal of contaminating foreground light, enabling reliable flux measurements at 3.6 micron and 4.5 micron. None of the galaxies are detected to [3.6]=26.9 (AB, 2 sigma), but a stacking analysis reveals a robust detection for the z_{850}-dropouts and an upper limit for the Y_{105}-dropouts. We construct average broadband SEDs using the stacked ACS, WFC3, and IRAC fluxes and fit stellar population synthesis models to derive mean redshifts, stellar masses, and ages. For the z_{850}-dropouts, we find z=6.9^{+0.1}_{-0.1}, (U-V)_{rest}=0.4, reddening A_V=0, stellar mass M*=1.2^{+0.3}_{-0.6} x 10^9 M_sun (Salpeter IMF). The best-fit ages ~300Myr, M/L_V=0.2, and SSFR=1.7Gyr^{-1} are similar to values reported for luminous z=7 galaxies, indicating the galaxies are smaller but not younger. The sub-L* galaxies observed here contribute significantly to the stellar mass density and under favorable conditions may have provided enough photons for sustained reionization at 7<z<11. In contrast, the z=8.3^{+0.1}_{-0.2} Y_{105}-dropouts have stellar masses that are uncertain by 1.5 dex due to the near-complete reliance on far-UV data. Adopting the 2 sigma upper limit on the M/L(z=8), the stellar mass density to M_{UV,AB} < -18 declines from rho*(z=7)=3.7^{+1.0}_{-1.8} x 10^6 M_sun Mpc^{-3} to rho*(z=8) < 8 x 10^5 M_sun Mpc^{-3}, following (1+z)^{-6} over 3<z<8. Lower masses at z=8 would signify more dramatic evolution, which can be established with deeper IRAC observations, long before the arrival of the James Webb Space Telescope.

We present a first morphological study of z~7-8 Lyman Break galaxies (LBGs) from Oesch et al. 2009 and Bouwens et al. 2009 detected in ultra-deep near-infrared imaging of the Hubble Ultra Deep field (HUDF) by the HUDF09 program. With an average intrinsic size of 0.7+-0.3 kpc these galaxies are found to be extremely compact having an average observed surface brightness of mu_J ~= 26 mag arcsec^(-2), and only two out of the full sample of 16 z~7 galaxies show extended features with resolved double cores. By comparison to lower redshift LBGs it is found that only little size evolution takes place from z~7 to z~6, while galaxies between z~4-5 show more extended wings in their apparent profiles. The average size scales as (1+z)^(-m) with m=1.12+-0.17 for galaxies with luminosities in the range (0.3-1)L*_{z=3} and with m=1.32+-0.52 for (0.12-0.3)L*_{z=3}, consistent with galaxies having constant comoving sizes. The peak of the size distribution changes only slowly from z~7 to z~4. However, a tail of larger galaxies (>~ 1.2 kpc) is gradually built up towards later cosmic times, possibly via hierarchical build-up or via enhanced accretion of cold gas. Additionally, the average star-formation surface density of LBGs with luminosities (0.3-1)L*_{z=3} is nearly constant at Sigma_{SFR}=1.9 Msun/yr/kpc^2 over the entire redshift range z~4-7 suggesting similar star-formation efficiencies at these early epochs. The above evolutionary trends seem to hold out to z~8 though the sample is still small and possibly incomplete.

We present a sample of 16 robust z~7 z-drop galaxies detected by the newly installed WFC3/IR camera on the Hubble Space Telescope. Our analysis is based on the first epoch data of the HUDF09 program covering the Hubble Ultra Deep Field with 60 orbits of Y, J, and H observations. These remarkable data cover 4.7 arcmin^2 and are the deepest NIR images ever taken, reaching to ~29 mag AB (5 sigma). The 16 z~6.5-7.5 galaxies have been identified based on the Lyman Break technique utilizing (z-Y) vs. (Y-J) colors. They have magnitudes J = 26.0-29.0 (AB), an average apparent half-light radius of ~0.16 arcsec (<~1 kpc), and show very blue colors (some even beta<~-2.5), in particular at low luminosities. The WFC3/IR data confirms previous NICMOS detections indicating that the dropout selection at z~7 is very reliable. Our data allow a first determination of the faint end slope of the z~7 luminosity function, reaching down to M_UV ~ -18, a full magnitude fainter than previous measurements. When fixing phi*=1.4e-3/Mpc^3/mag to the value previously measured at z~6, we find a best-fit value of alpha=-1.77+-0.20, with a characteristic luminosity of M_*=-19.91+-0.09. This steep slope is similar to what is seen at z~2-6 and indicates that low luminosity galaxies could potentially provide adequate flux to reionize the universe. The remarkable depth and resolution of these new images provide insights into the coming power of JWST.

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 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 new information on galaxies in the vicinity of luminous radio galaxies and quasars at z=4,5,6. These fields were previously found to contain overdensities of Lyman Break Galaxies (LBGs) or spectroscopic Lyman alpha emitters. We use HST and Spitzer data to infer stellar masses, and contrast our results with large samples of LBGs in more average environments as probed by the Great Observatories Origins Deep Survey (GOODS). The following results were obtained. First, LBGs in both overdense regions and in the field at z=4-5 lie on a very similar sequence in a z'-[3.6] versus [3.6] color-magnitude diagram. This is interpreted as a sequence in stellar mass (log[M*/Msun] = 9-11) in which galaxies become increasingly red due to dust and age as their star formation rate (SFR) increases. Second, the two radio galaxies are among the most massive objects (log[M*/Msun]~11) known to exist at z~4-5, and are extremely rare based on the low number density of such objects as estimated from the ~25x larger area GOODS survey. We suggest that the presence of these massive galaxies and supermassive black holes has been boosted through rapid accretion of gas or merging inside overdense regions. Third, the total stellar mass found in the z=4 ``proto-cluster'' TN1338 accounts for <30% of the stellar mass on the cluster red sequence expected to have formed at z>4, based on a comparison with the massive X-ray cluster Cl1252 at z=1.2. Although future near-infrared observations should determine whether any massive galaxies are currently being missed, one possible explanation for this mass difference is that TN1338 evolves into a smaller cluster than Cl1252. This raises the interesting question of whether the most massive protocluster regions at z>4 remain yet to be discovered.

PILOT (the Pathfinder for an International Large Optical Telescope) is a proposed 2.5 m optical/infrared telescope to be located at Dome C on the Antarctic plateau. The atmospheric conditions at Dome C deliver a high sensitivity, high photometric precision, wide-field, high spatial resolution, and high-cadence imaging capability to the PILOT telescope. These capabilities enable a unique scientific potential for PILOT, which is addressed in this series of papers. The current paper presents a series of projects dealing with the distant (redshift >) Universe, that have been identified as key science drivers for the PILOT facility. The potential for PILOT to detect the first populations of stars to form in the early Universe, via infrared projects searching for pair-instability supernovae and gamma-ray burst afterglows, is investigated. Two projects are proposed to examine the assembly and evolution of structure in the Universe: an infrared survey searching for the first evolved galaxies at high redshift, and an optical survey aimed at characterising moderate-redshift galaxy clusters. Finally, a large-area weak-lensing survey and a program to obtain supernovae infrared light-curves are proposed to examine the nature and evolution of dark energy and dark matter.

PILOT (the Pathfinder for an International Large Optical Telescope) is a proposed 2.5 m optical/infrared telescope to be located at Dome C on the Antarctic plateau. Conditions at Dome C are known to be exceptional for astronomy. The seeing (above ~30 m height), coherence time, and isoplanatic angle are all twice s good as at typical mid-latitude sites, while the water-vapour column, and the atmosphere and telescope thermal emission are all an order of magnitude better. These conditions enable a unique scientific capability for PILOT, which is addressed in this series of papers. The current paper presents an overview of the optical and instrumentation suite for PILO and its expected performance, a summary of the key science goals and observational approach for the facility, a discussion of the synergies between the science goals for PILOT and other telescopes, and a discussion of the future of Antarctic astronomy. Paper II and Paper III present details of the science projects divided, respectively, between the distant Universe (i.e., studies of first light, and the assembly and evolution of structure) and the nearby Universe (i.e., studies of Local Group galaxies, the Milky Way, and the Solar System).

We report on the discovery of three bright, strongly-lensed objects behind Abell 1703 and CL0024+16 from a dropout search over 25 square arcminutes of deep NICMOS data, with deep ACS optical coverage. They are undetected in the deep ACS images below 8500 A and have clear detections in the J and H bands. Fits to the ACS, NICMOS and IRAC data yield robust photometric redshifts in the range z~6-7 and largely rule out the possibility that they are low-redshift interlopers. All three objects are extended, and resolved into a pair of bright knots. The bright i-band dropout in Abell 1703 has an H-band AB magnitude of 23.9, which makes it one of the brightest known galaxy candidates at z>5.5. Our model fits suggest a young, massive galaxy only ~ 60 million years old with a mass of ~ 1E10 solar mass. The dropout galaxy candidates behind CL0024+16 are separated by 2.5" (~ 2 kpc in the source plane), and have H-band AB magnitudes of 25.0 and 25.6. Lensing models of CL0024+16 suggest that the objects have comparable intrinsic magnitudes of AB ~ 27.3, approximately one magnitude fainter than L* at z~6.5. Their similar redshifts, spectral energy distribution, and luminosities, coupled with their very close proximity on the sky, suggest that they are spatially associated, and plausibly are physically bound. Combining this sample with two previously reported, similarly magnified galaxy candidates at z~6-8, we find that complex systems with dual nuclei may be a common feature of high-redshift galaxies.

We conduct a search for z>~7 dropout galaxies behind 11 massive lensing clusters using 21 arcmin**2 of deep HST NICMOS, ACS, and WFPC2 image data. In total, over this entire area, we find only one robust z~7 z-dropout candidate (previously reported around Abell 1689). Four less robust z-dropout and J-dropout candidates are also found. The nature of the four weaker candidates could not be precisely determined due to the limited depth of the available optical data, but detailed simulations suggest that all four could be low-redshift interlopers. We compare these numbers with what we might expect using the z~7 UV luminosity function (LF) determined from field searches. We predict 2.7 z~7 z-dropouts and 0.3 z~9 J-dropouts over our cluster search area, in reasonable agreement with our observational results, given the small numbers. The number of z>~7 candidates we find in the present search are much lower than has been reported in several previous studies of the prevalence of z>~7 galaxies behind lensing clusters. To understand these differences, we examined z>~7 candidates in other studies and conclude that only a small fraction are likely to be z>~7 galaxies. Our findings support models that show that gravitational lensing from clusters is of the most value for detecting galaxies at magnitudes brighter than L* (H<~ 27) where the LF is expected to be very steep. Use of these clusters to constrain the faint-end slope or determine the full LF is likely of less value due to the shallower effective slope measured for the LF at fainter magnitudes, as well as significant uncertainties introduced from modelling both the gravitational lensing and incompleteness.

We calculate the major pair fraction and derive the major merger fraction and rate for 82 massive ($M_{*}>10^{11}M_{\odot}$) galaxies at $1.7 < z < 3.0$ utilising deep HST NICMOS data taken in the GOODS North and South fields. For the first time, our NICMOS data provides imaging with sufficient angular resolution and depth to collate a sufficiently large sample of massive galaxies at z $>$ 1.5 to reliably measure their pair fraction history. We find strong evidence that the pair fraction of massive galaxies evolves with redshift. We calculate a pair fraction of $f_{m}$ = 0.29 +/- 0.06 for our whole sample at $1.7 < z < 3.0$. Specifically, we fit a power law function of the form $f_{m}=f_{0}(1+z)^{m}$ to a combined sample of low redshift data from Conselice et al. (2007) and recently acquired high redshift data from the GOODS NICMOS Survey. We find a best fit to the free parameters of $f_{0}$ = 0.008 +/- 0.003 and $m$ = 3.0 +/- 0.4. We go on to fit a theoretically motivated Press-Schechter curve to this data. This Press-Schechter fit, and the data, show no sign of levelling off or turning over, implying that the merger fraction of massive galaxies continues to rise with redshift out to z $\sim$ 3. Since previous work has established that the merger fraction for lower mass galaxies turns over at z $\sim$ 1.5 - 2.0, this is evidence that higher mass galaxies experience more mergers earlier than their lower mass counterparts, i.e. a galaxy assembly downsizing. Finally, we calculate a merger rate at z = 2.6 of $\Re$ $<$ 5 $\times$ 10$^{5}$ Gpc$^{-3}$ Gyr$^{-1}$, which experiences no significant change to $\Re$ $<$ 1.2 $\times$ 10$^{5}$ Gpc$^{-3}$ Gyr$^{-1}$ at z = 0.5. This corresponds to an average $M_{*}>10^{11}M_{\odot}$ galaxy experiencing 1.7 +/- 0.5 mergers between z = 3 and z = 0.

Motivated by recent observational studies of the environment of z~6 QSOs, we have used the Millennium Run (MR) simulations to construct a very large (~20 deg^2) mock redshift survey of star-forming galaxies at z~6. We use this simulated survey to study the relation between density enhancements in the distribution of i-dropouts and Lya-emitters, and their relation to the most massive halos and protocluster regions at z~6. Our simulation predicts significant variations in surface density across the sky with some voids and filaments extending over scales of 1 degree, much larger than probed by current surveys. Approximately one third of all z~6 halos hosting i-dropouts brighter than z=26.5 mag (~M* at z=6) become part of z=0 galaxy clusters. i-dropouts associated with protocluster regions are found in regions where the surface density is enhanced on scales ranging from a few to several tens of arcminutes on the sky. We analyze two structures of i-dropouts and Lya-emitters observed with the Subaru Telescope and show that these structures must be the seeds of massive clusters-in-formation. In striking contrast, six z~6 QSO fields observed with HST show no significant enhancements in their i-dropout number counts. With the present data, we cannot rule out the QSOs being hosted by the most massive halos. However, neither can we confirm this widely used assumption. We conclude by giving detailed recommendations for the interpretation and planning of observations by current and future ground- and space based instruments that will shed new light on questions related to the large-scale structure at z~6.

We measure the sizes of 82 massive (M>10^11 M_sun) galaxies at 1.7<z<3 utilizing deep HST NICMOS data taken in the GOODS North and South fields. Our sample is almost an order of magnitude larger than previous studies at these redshifts, providing the first statistical study of massive galaxy sizes at z>2, confirming the extreme compactness of these galaxies. We split our sample into disk-like (n<2) and spheroid-like (n>2) galaxies based on their Sersic indices, and find that at a given stellar mass disk-like galaxies at z~2.3 are a factor of 2.6+/-0.3 smaller than present day equal mass systems, and spheroid-like galaxies at the same redshifts are 4.3+/-0.7 smaller than comparatively massive elliptical galaxies today. At z>2 our results are compatible with both a leveling off, or a mild evolution in size. Furthermore, the high density (~2x10^10 M_sun kpc^-3) of massive galaxies at these redshifts, which are similar to present day globular clusters, possibly makes any further evolution in sizes beyond z=3 unlikely.

Using HST and Spitzer IRAC imaging, we report the discovery of a very bright strongly lensed Lyman break galaxy (LBG) candidate at z~7.6 in the field of the massive galaxy cluster Abell 1689. The galaxy candidate, which we refer to as A1689-zD1, shows a strong z-J break of at least 2.2 mag and is completely undetected (<1 sigma) in HST/ACS g, r, i, and z-band data. These properties, combined with the very blue J-H and H-[4.5] colors, are exactly the properties of an z~7.6 LBG and can only be reasonably fit by a star-forming galaxy at z=7.6 +/- 0.4. Attempts to reproduce these properties with a model galaxy at z<4 yield particularly poor fits. A1689-zD1 has an observed (lensed) magnitude of 24.7 AB (8 sigma) in the NICMOS H band and is ~1.3 mag brighter than the brightest-known z-dropout galaxy. When corrected for the cluster magnification of 9.3 at z~7.6, the candidate has an intrinsic magnitude of H=27.1 AB, or about an L* galaxy at z~7.6. The source-plane deprojection shows that the star formation is occurring in compact knots of size ~<300 pc. The best-fit stellar population synthesis models yield a median redshift of 7.6, stellar masses (1.6-3.9) x 10^9 M_sun, stellar ages 45-320 Myr, star-formation rates ~<7.6 M_sun/yr, and low reddening with A_V <= 0.3. These properties are generally similar to those of LBGs found at z~5-6. The inferred stellar ages suggest a formation redshift of z~8-10 (t~<0.63 Gyr). A1689-zD1 is the brightest observed, highly reliable z>7.0 galaxy candidate found to date.

We use all available deep optical ACS and near-IR data over both the HUDF and the two GOODS fields to search for star-forming galaxies at z>~7 and constrain the UV LF within the first 700 Myrs. Our data set includes ~23 arcmin^2 of deep NICMOS J+H data and ~248 arcmin^2 of ground-based (ISAAC+MOIRCS) data, coincident with ACS optical data of greater or equal depths. In total, we find 8 <z>~7.3 z-dropouts in our search fields, but no z~9 J-dropout candidates. A careful consideration of a wide variety of different contaminants suggest an overall contamination level of just ~12% for our z-dropout selection. After performing detailed simulations to accurately estimate the selection volumes, we derive constraints on the UV LFs at z~7 and z~9. For a faint-end slope alpha=-1.74, our most likely values for M*(UV) and phi* at z~7 are -19.8+/-0.4 mag and 1.1_{-0.7}^{+1.7} x 10^{-3} Mpc^{-3}, respectively. Our search results for z~9 J-dropouts set a 1 sigma lower limit on M*(UV) of -19.6 mag assuming that phi* and alpha are the same as their values at slightly later times. This lower limit on M*(UV) is 1.4 mag fainter than our best-fit value at z~4, suggesting that the UV LF has undergone substantial evolution over this time period. No evolution is ruled out at 99% confidence from z~7 to z~6 and at 80% confidence from z~9 to z~7. The inferred brightening in M*(UV) with redshift (i.e., M*(UV) = (-21.02+/-0.09) + (0.36+/-0.08)(z - 3.8)) matches the evolution expected in the halo mass function, if the mass-to-light ratio of halos evolves as ~(1+z)**{-1}. Finally, we consider the shape of the UV LF at z>~5 and discuss the implications of the Schechter-like form of the observed LFs, particularly the unexpected abrupt cut-off at the bright end.

Using deep near-infrared spectroscopy Kriek et al. (2006) found that ~45% of massive galaxies at z~2.3 have evolved stellar populations and little or no ongoing star formation. Here we determine the sizes of these quiescent galaxies using deep, high-resolution images obtained with HST/NIC2 and laser guide star-assisted Keck/AO. Considering that their median stellar mass is 1.7x10^11 Solar masses the galaxies are remarkably small, with a median effective radius of 0.9 kpc. Galaxies of similar mass in the nearby Universe have sizes of ~5 kpc and average stellar densities which are two orders of magnitude lower than the z~2.3 galaxies. These results extend earlier work at z~1.5 and confirm previous studies at z>2 which lacked spectroscopic redshifts and imaging of sufficient resolution to resolve the galaxies. Our findings demonstrate that fully assembled early-type galaxies make up at most ~10% of the population of K-selected quiescent galaxies at z~2.3, effectively ruling out simple monolithic models for their formation. The galaxies must evolve significantly after z~2.3, through dry mergers or other processes, consistent with predictions from hierarchical models.

We use the ACS BViz data from the HUDF and all other deep HST ACS fields (including the GOODS fields) to find large samples of star-forming galaxies at z~4 and z~5 and to extend our previous z~6 sample. These samples contain 4671, 1416, and 627 B, V, and i dropouts, respectively, and reach to extremely low luminosities (0.01-0.04 L* or M(UV)~-16 to -17), allowing us to determine the rest-frame UV luminosity function (LF) and faint-end slope alpha at z~4-6 to high accuracy. We find faint-end slopes alpha of -1.73+/-0.05 at z~4, -1.66+/-0.09 at z~5, and -1.74+/-0.16 at z~6 -- suggesting that the faint-end slope is very steep and shows little evolution with cosmic time. We find that M*(UV) brightens considerably in the 0.7 Gyr from z~6 to z~4 (by ~0.7 mag from M*=-20.24+/-0.19 to M*=-20.98+/-0.10). The observed increase in the characteristic luminosity over this range is almost identical to that expected for the halo mass function -- suggesting that the observed evolution is likely due to the hierarchical coalescence and merging of galaxies. The evolution in phi* is not significant. The UV luminosity density at z~6 is modestly lower (0.45+/-0.09 times) than that at z~4 (integrated to -17.5 AB mag) though a larger change is seen in the dust-corrected star-formation rate density. We thoroughly examine published LF results and assess the reasons for their wide dispersion. We argue that the results reported here are the most robust available. The extremely steep faint-end slopes alpha found here suggest that lower luminosity galaxies play a significant role in reionizing the universe. Finally, we consider recent search results for galaxies at z~7-8 and use them to extend our estimates of the evolution in M* from z~7-8 to z~4.

We have detected 506 i-dropouts (z~6 galaxies) in deep, wide-area HST ACS fields: HUDF, enhanced GOODS, and HUDF-Parallel ACS fields (HUDF-Ps). The contamination levels are <=8% (i.e., >=92% are at z~6). With these samples, we present the most comprehensive, quantitative analyses of z~6 galaxies yet and provide optimal measures of the UV luminosity function (LF) and luminosity density at z~6, and their evolution to z~3. We redetermine the size and color evolution from z~6 to z~3. Field-to-field variations (cosmic variance), completeness, flux, and contamination corrections are modelled systematically and quantitatively. After corrections, we derive a rest-frame continuum UV (~1350 A) LF at z~6 that extends to M_{1350,AB} ~ -17.5 (0.04L*(z=3)). There is strong evidence for evolution of the LF between z~6 and z~3, most likely through a brightening (0.6+/-0.2 mag) of M* (at 99.7% confidence) though the degree depends upon the faint-end slope. As expected from hierarchical models, the most luminous galaxies are deficient at z~6. Density evolution (phi*) is ruled out at >99.99% confidence. Despite large changes in the LF, the luminosity density at z~6 is similar (0.82+/-0.21x) to that at z~3. Changes in the mean UV color of galaxies from z~6 to z~3 suggest an evolution in dust content, indicating the true evolution is substantially larger: at z~6 the star formation rate density is just ~30% of the z~3 value. Our UV luminosity function is consistent with z~6 galaxies providing the necessary UV flux to reionize the universe.

M85 optical transient 2006-1 (M85 OT 2006-1) is the most luminous member of the small family of V838 Mon-like objects, whose nature is still a mystery. This event took place in the Virgo cluster of galaxies and peaked at an absolute magnitude of I~-13. Here we present Hubble Space Telescope images of M85 OT 2006-1 and its environment, taken before and after the eruption, along with a spectrum of the host galaxy at the transient location. We find that the progenitor of M85 OT 2006-1 was not associated with any star forming region. The g and z-band absolute magnitudes of the progenitor were fainter than about -4 and -6 mag, respectively. Therefore, we can set a lower limit of ~50 Myr on the age of the youngest stars at the location of the progenitor that corresponds to a mass of <7 solar mass. Previously published line indices suggest that M85 has a mean stellar age of 1.6+/-0.3 Gyr. If this mean age is representative of the progenitor of M85 OT 2006-1, then we can further constrain its mass to be less than 2 solar mass. We compare the energetics and mass limit derived for the M85 OT 2006-1 progenitor with those expected from a simple model of violent stellar mergers. Combined with further modeling, these new clues may ultimately reveal the true nature of these puzzling events.

We present deep HST/ACS observations in g,r,i,z towards the z=4.1 radio galaxy TN J1338-1942 and its overdensity of >30 spectroscopically confirmed Lya emitters (LAEs). We select 66 g-band dropouts to z=27, 6 of which are also a LAE. Although our color-color selection results in a relatively broad redshift range centered on z=4.1, the field of TN J1338-1942 is richer than the average field at the >5 sigma significance, based on a comparison with GOODS. The angular distribution is filamentary with about half of the objects clustered near the radio galaxy, and a small, excess signal (2 sigma) in the projected pair counts at separations of <10" is interpreted as being due to physical pairs. The LAEs are young (a few x 10^7 yr), small (<r_50> = 0.13") galaxies, and we derive a mean stellar mass of ~10^8-9 Msun based on a stacked K-band image. We determine star formation rates, sizes, morphologies, and color-magnitude relations of the g-dropouts and find no evidence for a difference between galaxies near TN J1338-1942 and in the field. We conclude that environmental trends as observed in clusters at much lower redshift are either not yet present, or are washed out by the relatively broad selection in redshift. The large galaxy overdensity, its corresponding mass overdensity and the sub-clustering at the approximate redshift of TN J1338-1942 suggest the assemblage of a >10^14 Msun structure, confirming that it is possible to find and study cluster progenitors in the linear regime at z>4.

We present a comprehensive mass reconstruction of the rich galaxy cluster Cl 0024+17 at z~0.4 from ACS data, unifying both strong- and weak-lensing constraints. The weak-lensing signal from a dense distribution of background galaxies (~120 per square arcmin) across the cluster enables the derivation of a high-resolution parameter-free mass map. The strongly-lensed objects tightly constrain the mass structure of the cluster inner region on an absolute scale, breaking the mass-sheet degeneracy. The mass reconstruction of Cl 0024+17 obtained in such a way is remarkable. It reveals a ringlike dark matter substructure at r~75" surrounding a soft, dense core at r~50". We interpret this peculiar sub-structure as the result of a high-speed line-of-sight collision of two massive clusters 1-2 Gyr ago. Such an event is also indicated by the cluster velocity distribution. Our numerical simulation with purely collisionless particles demonstrates that such density ripples can arise by radially expanding, decelerating particles that originally comprised the pre-collision cores. Cl 0024+17 can be likened to the bullet cluster 1E0657-56, but viewed $along$ the collision axis at a much later epoch. In addition, we show that the long-standing mass discrepancy for Cl 0024+17 between X-ray and lensing can be resolved by treating the cluster X-ray emission as coming from a superposition of two X-ray systems. The cluster's unusual X-ray surface brightness profile that requires a two isothermal sphere description supports this hypothesis.

We present complimentary techniques to find emission-line targets and measure their properties in a semi-automated fashion from grism observations obtained with the Advanced Camera for Surveys aboard the Hubble Space Telescope. The first technique is to find all likely sources in a direct image, extract their spectra and search them for emission lines. The second method is to look for emission-line sources as compact structures in an unsharp masked version of the grism image. Using these methods we identify 46 emission-line targets in the Hubble Deep Field North using a modest (3 orbit) expenditure of HST observing time. Grism spectroscopy is a powerful tool for efficiently identifying interesting low luminosity, moderate redshift emission-line field galaxies. The sources found here have a median i band flux 1.5 mag fainter than the spectroscopic redshift catalog of Cohen et al. They have redshift z <= 1.42, high equivalent widths (typically EW > 100{\AA}), and are usually less luminous than the characteristic luminosity at the same redshift. The chief obstacle in interpreting the results is line identification, since the majority of sources have a single emission line and the spectral resolution is low. Photometric redshifts are useful for providing a first guess redshift. However, even at the depth of the state-of-the-art data used here, photometric errors can result in uncertainties in line identifications, especially for sources with i > ~24.5 ABmag. Reliable line identification for the faintest emission-line galaxies requires additional ground-based spectroscopy for confirmation. Of particular concern are the faint high EW [OII] emitters which could represent a strongly evolving galaxy population if the possibility that they are mis-identified lower redshift interlopers can be ruled out. (Slightly abridged)

We report on deep spectroscopy using LRIS on Keck I and FORS2 on the VLT of a sample of 22 candidate z~6 Lyman Break galaxies (LBGs) selected by the i-z> 1.3 dropout criterion. Redshifts could be measured for eight objects. These redshifts are all in the range z = 5.5 - 6.1, confirming the efficiency of the i-z color selection technique. Six of the confirmed galaxies show Ly-alpha emission. Assuming that the 14 objects without redshifts are z~6 LBGs, but lack detectable Ly-alpha emission lines, we infer that the fraction of Ly-alpha emitting LBGs with Ly-alpha equivalent widths greater than 20 Angstroms among z~6 LBGs is ~30%, similar to that found at z~3. Every Ly-alpha emitting object in our sample is compact with r <= 0.14". Furthermore, all the Ly-alpha emitting objects in our sample are more compact than average relative to the observed size-magnitude relation of a large i-dropout sample (332 candidate z~6 objects). We can reject the hypothesis that the Ly-alpha emitting population is a subset of the rest of the z~6 LBG population at >97% confidence. We speculate the small sizes of Ly-alpha emitting LBGs are due to these objects being less massive than other LBGs at z~6.

We measured the angular clustering at z~6 from a large sample of i-dropout galaxies (293 with z<27.5 from GOODS and 95 with z<29.0 from the UDF). Our largest and most complete subsample (having L>0.5L*) shows the presence of clustering at 94% significance. For this sample we derive a (co-moving) correlation length of r_0=4.5^{+2.1}_{-3.2} h_{72}^{-1} Mpc and bias b=4.1^{+1.5}_{-2.6}, using an accurate model for the redshift distribution. No clustering could be detected in the much deeper but significantly smaller UDF, yielding b<4.4 (1 sigma). We compare our findings to Lyman break galaxies at z=3-5 at a fixed luminosity. Our best estimate of the bias parameter implies that i-dropouts are hosted by dark matter halos having masses of ~10^11 M_sun, similar to that of V-dropouts at z~5. We evaluate a recent claim that at z>5 star formation might have occurred more efficiently compared to that at z=3-4. This may provide an explanation for the very mild evolution observed in the UV luminosity density between z=6 and z=3. Although our results are consistent with such a scenario, the errors are too large to find conclusive evidence for this.

The first 900 million years (Myr) to redshift z~6 (the first seven per cent of the age of the Universe) remains largely unexplored for the formation of galaxies. Large samples of galaxies have been found at z~6, but detections at earlier times are uncertain and unreliable. It is not at all clear how galaxies built up from the first stars when the Universe was ~300 Myr old (z~12-15) to z~6, just 600 Myr later. Here we report the results of a search for galaxies at z~7-8, about 700 Myr after the Big Bang, using the deepest near-infrared and optical images ever taken. Under conservative selection criteria we find only one candidate galaxy at z~7-8, where ten would be expected if there were no evolution in the galaxy population between z~7-8 and z~6. Using less conservative criteria, there are four candidates, where 17 would be expected with no evolution. This demonstrates that very luminous galaxies are quite rare 700 Myr after the Big Bang. The simplest explanation is that the Universe is just too young to have built up many luminous galaxies at z~7-8 by the hierarchical merging of small galaxies.

Using Spitzer IRAC mid-infrared imaging from the Great Observatories Origins Deep Survey, we study z_850-dropout sources in the Hubble Ultra Deep Field. After carefully removing contaminating flux from foreground sources, we clearly detect two z_850-dropouts at 3.6 micron and 4.5 micron, while two others are marginally detected. The mid-infrared fluxes strongly support their interpretation as galaxies at z~7, seen when the Universe was only 750 Myr old. The IRAC observations allow us for the first time to constrain the rest-frame optical colors, stellar masses, and ages of the highest redshift galaxies. Fitting stellar population models to the spectral energy distributions, we find photometric redshifts in the range 6.7-7.4, rest-frame colors U-V=0.2-0.4, V-band luminosities L_V=0.6-3 x 10^10 L_sun, stellar masses 1-10 x 10^9 M_sun, stellar ages 50-200 Myr, star formation rates up to ~25 M_sun/yr, and low reddening A_V<0.4. Overall, the z=7 galaxies appear substantially less massive and evolved than Lyman break galaxies or Distant Red Galaxies at z=2-3, but fairly similar to recently identified systems at z=5-6. The stellar mass density inferred from our z=7 sample is rho* = 1.6^{+1.6}_{-0.8} x 10^6 M_sun Mpc^-3 (to 0.3 L*(z=3)), in apparent agreement with recent cosmological hydrodynamic simulations, but we note that incompleteness and sample variance may introduce larger uncertainties. The ages of the two most massive galaxies suggest they formed at z>8, during the era of cosmic reionization, but the star formation rate density derived from their stellar masses and ages is not nearly sufficient to reionize the universe. The simplest explanation for this deficiency is that lower-mass galaxies beyond our detection limit reionized the universe.

Today's deep high resolution multiwavelength surveys contain a wealth of information about galaxies at different epochs. To fully exploit this information, it is useful to be able to produce highly realistic simulations to compare with the observations. Here, we describe one such simulator for producing imaging data for deep galaxy fields. Based upon a pixel-by-pixel modelling of object SEDs and their selection volumes, this simulator allows users to make realistic multicolor simulations of galaxy fields from galaxy samples at all redshifts. Input samples for this simulator range from z~0 samples selected from the Sloan Digital Sky Survey to z~1-6 samples selected from the Great Observatories Origins Deep Survey. Users can specify the passbands, noise, and PSFs, or equivalently the exposure times on well-known instruments like HST or ground-based telescopes.

We present aperture-matched PSF-corrected BVi'z'JH photometry and Bayesian photometric redshifts (BPZ) for objects detected in the Hubble Ultra Deep Field (UDF), 8,042 of which are detected at the 10-sigma level (e.g., i'<29.01 or z'<28.43). Most of our objects are defined identically to those in the public STScI catalogs, enabling straightforward object-by-object comparison. We have combined detections from i', z', J+H, and B+V+i'+z' images into a single comprehensive segmentation map. Using a new program called SExSeg we are able to force this segmentation map into SExtractor for photometric analysis. The resulting photometry is corrected for the wider NIC3 PSFs using our ColorPro software. We also correct for the ACS z'-band PSF halo. The NIC3 magnitudes are found to be too faint relative to the ACS fluxes. Based on BPZ SED fits to objects of know spectroscopic redshift, we derived corrections of -0.30 +/- 0.03 mag in J and -0.18 +/- 0.04 mag in H. The offsets appear to be supported by a recent recalibration of the UDF NIC3 images combined with non-linearity measured in NICMOS itself. The UDF reveals a large population of faint blue galaxies (presumably young starbursts), bluer than those observed in the original Hubble Deep Fields (HDF). To accommodate these galaxies, we have added two new starburst templates to the SED library used in previous BPZ papers. The resulting photometric redshifts are accurate to within 0.04 * 1+z_spec out to z < 6. Finally, we measure galaxy morphology, including Sersic index and asymmetry. Our full catalog, software packages, and more are available at http://adcam.pha.jhu.edu/~coe/UDF/ (abridged)

We study the photometric and structural properties of spectroscopically confirmed members in the two massive X-ray--selected z=0.83 galaxy clusters MS1054-03 and RXJ0152-1357 using three-band mosaic imaging with the Hubble Space Telescope Advanced Camera for Surveys. The samples include 105 and 140 members of MS1054-03 and RXJ0152-1357, respectively, with ACS F775W magnitude < 24.0. We develop a promising new structural classification method, based on a combination of the best-fit Sersic indices and the normalized root-mean-square residuals from the fits; the resulting classes agree well with the visual ones, but are less affected by galaxy orientation. We examine the color--magnitude relations in detail and find that the color residuals correlate with the local mass density measured from our weak lensing maps; we identify a threshold density of $\Sigma \approx 0.1$, in units of the critical density, above which the star formation appears to cease. For RXJ0152-1357, we also find a trend in the color residuals with velocity, resulting from an offset of about 980 km/s in the mean redshifts of the early- and late-type galaxies. Analysis of the color--color diagrams indicates that a range of star formation time-scales are needed to reproduce the loci of the galaxy colors. We also identify some cluster galaxies whose colors can only be explained by large amounts, $A_V \approx 1$ mag, of internal dust extinction. [Abstract shortened]

(Abridged.) We present F435W (B), F606W (Broad V), and F814W (Broad I) coronagraphic images of the debris disk around Beta Pictoris obtained with HST's Advanced Camera for Surveys. We confirm that the previously reported warp in the inner disk is a distinct secondary disk inclined by ~5 deg from the main disk. The main disk's northeast extension is linear from 80 to 250 AU, but the southwest extension is distinctly bowed with an amplitude of ~1 AU over the same region. Both extensions of the secondary disk appear linear, but not collinear, from 80 to 150 AU. Within ~120 AU of the star, the main disk is ~50% thinner than previously reported. The surface-brightness profiles along the spine of the main disk are fitted with four distinct radial power laws between 40 and 250 AU, while those of the secondary disk between 80 and 150 AU are fitted with single power laws. These discrepancies suggest that the two disks have different grain compositions or size distributions. The F606W/F435W and F814W/F435W flux ratios of the composite disk are nonuniform and asymmetric about both projected axes of the disk. Within ~120 AU, the m_F435W-m_F606W and m_F435W-m_F814W colors along the spine of the main disk are ~10% and ~20% redder, respectively, than those of Beta Pic. These colors increasingly redden beyond ~120 AU, becoming 25% and 40% redder, respectively, than the star at 250 AU. We compare the observed red colors within ~120 AU with the simulated colors of non-icy grains having a radial number density ~r^-3 and different compositions, porosities, and minimum grain sizes. The observed colors are consistent with those of compact or moderately porous grains of astronomical silicate and/or graphite with sizes >0.15-0.20 um, but the colors are inconsistent with the blue colors expected from grains with porosities >90%.

HST ACS and NICMOS data are now of sufficient depth and areal coverage to place strong constraints on the formation and evolution of galaxies during the first 1-2 Gyrs of the universe. Of particular interest are galaxies at z~6 since they represent the earliest epoch accessible by current high-efficiency optical instrumentation. To this end, we have constructed a sample of 506 z~6 objects from all the deepest wide-area HST data (UDF, UDF-Parallel, and GOODS fields). They have been used to construct an optimal determination of the rest-frame continuum UV LF at z~6. Our LF extends to over 3 magnitudes below L*, fainter than has been done at z~3. Over the interval z~6 to z~3, we find strong evidence for evolution in the UV LF. Though constraints on the faint-end slope remain modest (and are consistent with no-evolution), the characteristic luminosity appears to have approximately doubled over the interval z~6 to z~3, consistent with hierarchical expectations. Remarkably, this shift to lower luminosities extends to even higher redshifts. Using all deep J+H NICMOS observations (800 orbits in total), we have been able to demonstrate that the bright end of the LF (>0.3L*) is at least 5 times lower at z~10 than at z~4, with a similar deficit being established from our recent detections and first statistical sample of z~7-8 galaxies using our UDF NICMOS data. In these proceedings, we discuss what is known about the UV LF and UV luminosity density at z~6-10 from current data and its evolution relative to z~3. We also describe several exciting prospects for advance in this area over the next year.

The color-magnitude relation has been determined for the RDCS J0910+5422 cluster of galaxies at redshift z = 1.106. Cluster members were selected from HST ACS images, combined with ground--based near--IR imaging and optical spectroscopy. The observed early--type color--magnitude relation (CMR) in (i_775 -z_850) versus z_850 shows intrinsic scatters in color of 0.042 +/- 0.010 mag and 0.044 +/- 0.020 mag for ellipticals and S0s, respectively. From the scatter about the CMR, a mean luminosity--weighted age t > 3.3 Gyr (z > 3) is derived for the elliptical galaxies. Strikingly, the S0 galaxies in RDCS J0910+5422 are systematically bluer in (i_775 - z_850) by 0.07 +/- 0.02 mag, with respect to the ellipticals. The ellipticity distribution as a function of color indicates that the face-on S0s in this particular cluster have likely been classified as elliptical. Thus, if anything, the offset in color between the elliptical and S0 populations may be even more significant. The color offset between S0 and E corresponds to an age difference of ~1 Gyr, for a single-burst solar metallicity model. A solar metallicity model with an exponential decay in star formation will reproduce the offset for an age of 3.5 Gyr, i.e. the S0s have evolved gradually from star forming progenitors. The early--type population in this cluster appears to be still forming. The blue early-type disk galaxies in RDCS J0910+5422 likely represent the direct progenitors of the more evolved S0s that follow the same red sequence as ellipticals in other clusters. Thirteen red galaxy pairs are observed and the galaxies associated in pairs constitute ~40% of the CMR galaxies in this cluster.

A five square arcminute region around the luminous radio-loud quasar SDSS J0836+0054 (z=5.8) hosts a wealth of associated galaxies, characterized by very red (1.3 < i_775 - z_{850} < 2.0) color. The surface density of these z~5.8 candidates is approximately six times higher than the number expected from deep ACS fields. This is one of the highest galaxy overdensities at high redshifts, which may develop into a group or cluster. We also find evidence for a substructure associated with one of the candidates. It has two very faint companion objects within two arcseconds, which are likely to merge. The finding supports the results of a recent simulation that luminous quasars at high redshifts lie on the most prominent dark-matter filaments and are surrounded by many fainter galaxies. The quasar activity from these regions may signal the buildup of a massive system.

We study the impact of outflows driven by active galactic nuclei (AGN) on galaxy formation. Outflows move into the surrounding intergalactic medium (IGM) and heat it sufficiently to prevent it from condensing onto galaxies. In the dense, high-redshift IGM, such feedback requires highly energetic outflows, driven by a large AGN. However, in the more tenuous low-redshift IGM, equivalently strong feedback can be achieved by less energetic winds (and thus smaller galaxies). Using a simple analytic model, we show that this leads to the anti-hierarchical quenching of star-formation in large galaxies, consistent with current observations. At redshifts prior to the formation of large AGN, galaxy formation is hierarchical and follows the growth of dark-matter halos. The transition between the two regimes lies at the z ~ 2 peak of AGN activity.

We present HST/ACS observations of the most distant radio galaxy known, TN J0924-2201 at z=5.2. This radio galaxy has 6 spectroscopically confirmed Lya emitting companion galaxies, and appears to lie within an overdense region. The radio galaxy is marginally resolved in i_775 and z_850 showing continuum emission aligned with the radio axis, similar to what is observed for lower redshift radio galaxies. Both the half-light radius and the UV star formation rate are comparable to the typical values found for Lyman break galaxies at z~4-5. The Lya emitters are sub-L* galaxies, with deduced star formation rates of 1-10 Msun/yr. One of the Lya emitters is only detected in Lya. Based on the star formation rate of ~3 Msun/yr calculated from Lya, the lack of continuum emission could be explained if the galaxy is younger than ~2 Myr and is producing its first stars. Observations in V_606, i_775, and z_850 were used to identify additional Lyman break galaxies associated with this structure. In addition to the radio galaxy, there are 22 V-break (z~5) galaxies with z_850<26.5 (5sigma), two of which are also in the spectroscopic sample. We compare the surface density of 2/arcmin^2 to that of similarly selected V-dropouts extracted from GOODS and the UDF Parallel fields. We find evidence for an overdensity to very high confidence (>99%), based on a counts-in-cells analysis applied to the control field. The excess is suggestive of the V-break objects being associated with a forming cluster around the radio galaxy.

We have obtained a high S/N (22.3 hr integration) UV continuum VLT FORS2 spectrum of an extremely bright (z_850 = 24.3) z = 5.515 +/- 0.003 starforming galaxy (BD38) in the field of the z = 1.24 cluster RDCS 1252.9-2927. This object shows substantial continuum (0.41 +/- 0.02 \muJy at \lambda1300) and low-ionization interstellar absorption features typical of LBGs at lower redshift (z ~ 3); this is the highest redshift LBG confirmed via metal absorption spectral features. The equivalent widths of the absorption features are similar to z ~ 3 strong Ly\alpha absorbers. No noticeable Ly\alpha emission was detected (F <= 1.4 * 10^-18 ergs cm^-2 s^-1, 3\sigma). The half-light radius of this object is 1.6 kpc (0\farcs25) and the star formation rate derived from the rest-frame UV luminosity is SFR_UV = 38 h^-2_0.7 M_sun yr^-1 (142 h^-2_0.7 M_sun yr^-1 corrected for dust extinction). In terms of recent determinations of the z ~ 6 UV luminosity function, this object appears to be 6L*. The Spitzer IRAC fluxes for this object are 23.3 and 23.2 AB mag (corrected for 0.3 mag of cluster lensing) in the 3.6\mu and 4.5\mu channels, respectively, implying a mass of 1-6 * 10^10 M_sun from population synthesis models. This galaxy is brighter than any confirmed z ~ 6 i-dropout to date in the z_850 band, and both the 3.6\mu and 4.5\mu channels, and is the most massive starbursting galaxy known at z > 5. -- Abstract Abridged

This paper describes the observations and data reduction techniques for the version 2.0 images and catalog of the NICMOS Ultra Deep Field Treasury program. All sources discussed in this paper are based on detections in the combined NICMOS F110W and F160W bands only. The NICMOS images are drizzled to 0.09 arc second pixels and aligned to the ACS UDF F850LP image which was rebinned to the same pixel scale. These form the NICMOS version 2.0 UDF images. The catalog sources are chosen with a conservative detection limit to avoid the inclusion of numerous spurious sources. The catalog contains 1293 objects in the 144 x 144 arc sececonds NICMOS subfield of the UDF. The 5 sigma signal to noise level is an average 0.6 arc second diameter aperture AB magnitude of ~27.7 at 1.1 and 1.6 microns. The catalog sources, listed in order of right ascension, satisfy a minimum signal to noise criterion of 1.4 sigma in at least 7 contiguous pixels of the combined F110W and F160W image

We study the dynamical origin of the structures observed in the scattered-light images of the resolved debris disk around HD 141569A. We explore the roles of radiation pressure from the central star, gas drag from the gas disk, and the tidal forces from two nearby stars in creating and maintaining these structures. We use a simple one-dimensional axisymmetric model to show that the presence of the gas helps confine the dust and that a broad ring of dust is produced if a central hole exists in the disk. This model also suggests that the disk is in a transient, excited dynamical state, as the observed dust creation rate applied over the age of the star is inconsistent with submillimeter mass measurements. We model in two dimensions the effects of a fly-by encounter between the disk and a binary star in a prograde, parabolic, coplanar orbit. We track the spatial distribution of the disk's gas, planetesimals, and dust. We conclude that the surface density distribution reflects the planetesimal distribution for a wide range of parameters. Our most viable model features a disk of initial radius 400 AU, a gas mass of 50 M_earth, and beta = 4 and suggests that the system is being observed within 4000 yr of the fly-by periastron. The model reproduces some features of HD 141569A's disk, such as a broad single ring and large spiral arms, but it does not reproduce the observed multiple spiral rings or disk asymmetries nor the observed clearing in the inner disk. For the latter, we consider the effect of a 5 M_Jup planet in an eccentric orbit on the planetesimal distribution of HD 141569A.

We take advantage of the Hubble Ultra Deep Field (UDF) data to study the restframe optical and ultra violet (UV) morphologies of the novel population of Distant Red Galaxies (DRGs). Six galaxies with J-Ks > 2.3 are found to Ks=21.5, five of which have photometric redshifts z_phot > 2, corresponding to a surface density of 0.9/arcmin^2. The surface brightness distributions of the z_phot > 2 galaxies are better represented by exponential disks than R^{1/4}-laws. Two of the z_phot > 2 galaxies are extended, while three have compact morphologies. The restframe optical morphology of the z_phot > 2 galaxies is quite different from the restframe UV morphology: all the galaxies have red central components which dominate in the NICMOS H_{160}-band images, and distinct off-center blue features which show up in (and often dominate) the ACS images. The mean measured effective radius of the z_phot > 2 galaxies is <r_e> =1.9+/-1.4 kpc, similar (within the errors) to the mean size of LBGs at similar redshifts. All the DRGs are resolved in the ACS images, while four are resolved in the NICMOS images. Two of the z_phot > 2 galaxies are bright X-ray sources and hence host AGN. The diverse restframe optical and UV morphological properties of DRGs derived here suggest that they have complex stellar populations, consisting of both evolved populations that dominate the mass and the restframe optical light, and younger populations, which show up as patches of star formation in the restframe UV light; in many ways resembling the properties of normal local galaxies. This interpretation is supported by fits to the broadband SEDs, which for all five z_phot > 2 are best represented by models with extended star formation histories and substantial amounts of dust.

We use all available fields with deep NICMOS imaging to search for J dropouts (H<28) at z~10. Our primary data set for this search were the two J+H NICMOS parallel fields taken with the ACS HUDF. The 5 sigma limiting mags were 28.6 in J and 28.5 in H. Several shallower fields were also used: J+H NICMOS frames available over the HDF North, the HDF South NICMOS parallel, and the ACS HUDF. The primary selection criterion was (J-H)>1.8. 11 such sources were found in all search fields using this criterion. 8 of these were clearly ruled out as credible z~10 sources, either as a result of detections (>2 sigma) blueward of J or their colors redward of the break (H-K~1.5). The nature of the 3 remaining sources could not be determined from the data. The number appears consistent with the expected contamination from low-z interlopers. Analysis of the stacked images for the 3 candidates also suggests contamination. Regardless of their true redshifts, the actual number of z~10 sources must be <=3. To assess the significance of these results, two lower redshift samples (a z~3.8 B-dropout and z~6 i-dropout sample) were projected to z~8-12 using a (1+z)^{-1} size scaling. They were added to the image frames, and the selection repeated, giving 15.6 and 4.8 J-dropouts, respectively. This suggests that to the limit of this probe (0.3 L*) there has been evolution from z~3.8 and possibly from z~6. This is consistent with the strong evolution already noted at z~6 and z~7.5 relative to z~3-4. Even assuming that 3 sources from this probe are at z~10, the rest-frame continuum UV (~1500 A) luminosity density at z~10 (integrated down to 0.3 L*) is just 0.19_{-0.09}^{+0.13}x that at z~3.8 (or 0.19_{-0.10}^{+0.15}x including cosmic variance). However, if none of our sources is at z~10, this ratio has a 1 sigma upper limit of 0.07. (abridged)

We investigate the evolution in the distribution of surface brightness, as a function of size, for elliptical and S0 galaxies in the two clusters RDCS J1252.9-2927, z=1.237 and RX J0152.7-1357, z=0.837. We use multi-color imaging with the Advanced Camera for Surveys on the Hubble Space Telescope to determine these sizes and surface brightnesses. Using three different estimates of the surface brightnesses, we find that we reliably estimate the surface brightness for the galaxies in our sample with a scatter of < 0.2 mag and with systematic shifts of \lesssim 0.05 mag. We construct samples of galaxies with early-type morphologies in both clusters. For each cluster, we use a magnitude limit in a band which closely corresponds to the rest-frame B, to magnitude limit of M_B = -18.8 at z=0, and select only those galaxies within the color-magnitude sequence of the cluster or by using our spectroscopic redshifts. We measure evolution in the rest-frame B surface brightness, and find -1.41 \+/- 0.14 mag from the Coma cluster of galaxies for RDCS J1252.9-2927 and -0.90 \+/- 0.12 mag of evolution for RX J0152.7-1357, or an average evolution of (-1.13 \+/- 0.15) z mag. Our statistical errors are dominated by the observed scatter in the size-surface brightness relation, sigma = 0.42 \+/- 0.05 mag for RX J0152.7-1357 and sigma = 0.76 \+/- 0.10 mag for RDCS J1252.9-2927. We find no statistically significant evolution in this scatter, though an increase in the scatter could be expected. Overall, the pace of luminosity evolution we measure agrees with that of the Fundamental Plane of early-type galaxies, implying that the majority of massive early-type galaxies observed at z =~ 1 formed at high redshifts.

We present new measurements of the galaxy luminosity function (LF) and its dependence on local galaxy density, color, morphology, and clustocentric radius for the massive z=0.83 cluster MS1054-0321. Our analyses are based on imaging performed with the ACS onboard the HST in the F606W, F775W and F850LP passbands and extensive spectroscopic data obtained with the Keck LRIS. Our main results are based on a spectroscopically selected sample of 143 cluster members with morphological classifications derived from the ACS observations. Our three primary findings are (1) the faint-end slope of the LF is steepest in the bluest filter, (2) the LF in the inner part of the cluster (or highest density regions) has a flatter faint-end slope, and (3) the fraction of early-type galaxies is higher at the bright end of the LF, and gradually decreases toward fainter magnitudes. These characteristics are consistent with those in local galaxy clusters, indicating that, at least in massive clusters, the common characteristics of cluster LFs are established at z=0.83. We also find a 2sigma deficit of intrinsically faint, red galaxies (i-z>0.5, Mi>-19) in this cluster. This trend may suggest that faint, red galaxies (which are common in z<0.1 rich clusters) have not yet been created in this cluster at z=0.83. The giant-to-dwarf ratio in MS1054-0321 starts to increase inwards of the virial radius or when Sigma>30 Mpc^-2, coinciding with the environment where the galaxy star formation rate and the morphology-density relation start to appear. (abridged)

We combine imaging data from the Advanced Camera for Surveys (ACS) with VLT/FORS optical spectroscopy to study the properties of star-forming galaxies in the z=0.837 cluster CL0152-1357. We have morphological information for 24 star-forming cluster galaxies, which range in morphology from late-type and irregular to compact early-type galaxies. We find that while most star-forming galaxies have $r_{625}-i_{775}$ colors bluer than 1.0, eight are in the red cluster sequence. Among the star-forming cluster population we find five compact early-type galaxies which have properties consistent with their identification as progenitors of dwarf elliptical galaxies. The spatial distribution of the star-forming cluster members is nonuniform. We find none within $R\sim 500$ Mpc of the cluster center, which is highly suggestive of an intracluster medium interaction. We derive star formation rates from [OII] $\lambda\lambda 3727$ line fluxes, and use these to compare the global star formation rate of CL0152-1357 to other clusters at low and intermediate redshifts. We find a tentative correlation between integrated star formation rates and $T_{X}$, in the sense that hotter clusters have lower integrated star formation rates. Additional data from clusters with low X-ray temperatures is needed to confirm this trend. We do not find a significant correlation with redshift, suggesting that evolution is either weak or absent between z=0.2-0.8.

We present a sample of over 50 luminous field bulges (including ellipticals) found in the Groth Strip Survey (GSS), with 0.73< z < 1.04 and with bulge magnitudes I <= 23. The exponential disk light is removed via decomposition of HST images using GIM2D. We find that 85% of these bulges are nearly as red as local E/S0's and have a shallow slope and a small color dispersion in the color-luminosity relation, suggesting roughly coeval formation. The surface brightnesses of these bulges are about 1 mag higher than local bulges. These results are explained adopting a "drizzling" scenario where a metal-rich early formation is later polluted by small amounts of additional star formation. Almost all disks have the same or bluer colors than their accompanying bulges, regardless of the bulge-disk ratio and bulge luminosity, as expected from semi-analytic hierarchical galaxy formation models. We present evidence that the few blue bulge candidates are not likely to be genuine blue ellipticals or bulges. Our deeper, more extensive, and less disk-contaminated observations challenge prior claims that 30% to 50% of field bulges or ellipticals are in a blue, star-forming phase at z < 1. We conclude that field bulges and ellipticals at z ~ 1, like luminous early- type cluster galaxies at the same redshift, are already dominated by metal-rich, old stellar populations that have been fading from a formation epoch earlier than z ~ 1.5. (abridged)

Galaxies undergoing formation and evolution can now be observed over a time baseline of some 12 Gyr. An inherent difficulty with high-redshift observations is that the objects are very faint and the best resolution (HST) is only ~0.5 kpc. Such studies thereby combine in a highly synergistic way with the great detail that can be obtained for nearby galaxies. 3 new developments are highlighted. First is the derivation of stellar masses for galaxies from SEDs using HST and now Spitzer data, and dynamical masses from both sub-mm observations of CO lines and near-IR observations of optical lines like Halpha. A major step has been taken with evidence that points to the z~2-3 LBGs having masses that are a few x 10^10 Msolar. Second is the discovery of a population of evolved red galaxies at z~2-3 which appear to be the progenitors of the more massive early-type galaxies of today, with dynamical masses around a few x 10^11 Msolar. Third are the remarkable advances that have occurred in characterizing dropout galaxies to z~6 and beyond, < 1 Gyr from recombination. The HST ACS has played a key role here, with the dropout technique being applied to i & z images in several deep ACS fields, yielding large samples of these objects. This has allowed a detailed determination of their properties and meaningful comparisons against lower-z samples. The use of cloning techniques has overcome many of the strong selection biases affecting the study of these objects. A clear trend of size with redshift has been identified, and its impact on the luminosity density and SFR estimated. There is a significant though modest decrease in the SFR from z~2.5 to z~6. The latest data also allow for the first robust determination of the LF at z~6. Finally, the latest UDF ACS and NICMOS data has resulted in the detection of some galaxies at z~7-8.

We have detected likely z~7-8 galaxies in the 144"x144" NICMOS observations of the Hubble Ultra Deep Field. Objects are required to be >3 sigma detections in both NICMOS bands, J and H. The selection criteria for this sample are (z-J)>0.8, (z-J) > 0.66 (J-H)+0.8, (J-H)<1.2, and no detection at <8500 Angstroms. The 5 selected sources have total magnitudes H~27. Four of the five sources are quite blue compared to typical lower-redshift dropout galaxies and are clustered within a 1 arcmin**2 region. Because all 5 sources are near the limit of the NICMOS data, we have carefully evaluated their reality. Each of the candidates is visible in different splits of the data and a median stack. We analyzed several noise images and estimate the number of spurious sources to be 1+/-1. A search using an independent reduction of this same data set clearly revealed 3 of the 5 candidates and weakly detected a 4th candidate, suggesting the contamination could be higher. For comparison with predictions from lower redshift samples we take a conservative approach and adopt four z~7-8 galaxies as our sample. With the same detection criteria on simulated datasets, assuming no-evolution from z~3.8, we predict 10 sources at z~7-8, or 14 if we use a more realistic (1+z)^{-1} size scaling. We estimate that the rest-frame continuum UV (~1800 Angstroms) luminosity density at z~7.5 (integrated down to 0.3L*) is just 0.20_{-0.08}^{+0.12}x that found at z~3.8. Effectively this sets an upper limit on the luminosity density down to 0.3 L*. This result is consistent with significant evolution at the bright end of the luminosity function from z~7.5 to z~3.8. Even with the lower UV luminosity density at z~7.5, it appears that galaxies could still play an important role in reionization at these redshifts, though definitive measurements remain to be made. [abridged]

We analyse deep multi-colour Advanced Camera images of the largest known gravitational lens, A1689. Radial and tangential arcs delineate the critical curves in unprecedented detail and many small counter-images are found near the center of mass. We construct a flexible light deflection field to predict the appearance and positions of counter-images. The model is refined as new counter-images are identified and incorporated to improve the model, yielding a total of 106 images of 30 multiply lensed background galaxies, spanning a wide redshift range, 1.0$<$z$<$5.5. The resulting mass map is more circular in projection than the clumpy distribution of cluster galaxies and the light is more concentrated than the mass within $r<50kpc/h$. The projected mass profile flattens steadily towards the center with a shallow mean slope of $d\log\Sigma/d\log r \simeq -0.55\pm0.1$, over the observed range, r$<250kpc/h$, matching well an NFW profile, but with a relatively high concentration, $C_{vir}=8.2^{+2.1}_{-1.8}$. A softened isothermal profile ($r_{core}=20\pm2$\arcs) is not conclusively excluded, illustrating that lensing constrains only projected quantities. Regarding cosmology, we clearly detect the purely geometric increase of bend-angles with redshift. The dependence on the cosmological parameters is weak due to the proximity of A1689, $z=0.18$, constraining the locus, $\Omega_M+\Omega_{\Lambda} \leq 1.2$. This consistency with standard cosmology provides independent support for our model, because the redshift information is not required to derive an accurate mass map. Similarly, the relative fluxes of the multiple images are reproduced well by our best fitting lens model.

We measure the luminosity function of morphologically selected E/S0 galaxies from $z=0.5$ to $z=1.0$ using deep high resolution Advanced Camera for Surveys imaging data. Our analysis covers an area of $48\Box\arcmin$ (8$\times$ the area of the HDF-N) and extends 2 magnitudes deeper ($I\sim24$ mag) than was possible in the Deep Groth Strip Survey (DGSS). At $0.5<z<0.75$, we find $M_B^*-5\log h_{0.7}=-21.1\pm0.3$ and $\alpha=-0.53\pm0.2$, and at $0.75<z<1.0$, we find $M_B^*-5\log h_{0.7}=-21.4\pm0.2$. These luminosity functions are similar in both shape and number density to the luminosity function using morphological selection (e.g., DGSS), but are much steeper than the luminosity functions of samples selected using morphological proxies like the color or spectral energy distribution (e.g., CFRS, CADIS, or COMBO-17). The difference is due to the `blue', $(U-V)_0<1.7$, E/S0 galaxies, which make up to $\sim30%$ of the sample at all magnitudes and an increasing proportion of faint galaxies. We thereby demonstrate the need for {\it both morphological and structural information} to constrain the evolution of galaxies. We find that the `blue' E/S0 galaxies have the same average sizes and Sersic parameters as the `red', $(U-V)_0>1.7$, E/S0 galaxies at brighter luminosities ($M_B<-20.1$), but are increasingly different at fainter magnitudes where `blue' galaxies are both smaller and have lower Sersic parameters. Fits of the colors to stellar population models suggest that most E/S0 galaxies have short star-formation time scales ($\tau<1$ Gyr), and that galaxies have formed at an increasing rate from $z\sim8$ until $z\sim2$ after which there has been a gradual decline.

We use the exceptional depth of the Ultra Deep Field (UDF) and UDF-Parallel ACS fields to study the sizes of high redshift (z~2-6) galaxies and address long-standing questions about possible biases in the cosmic star formation rate due to surface brightness dimming. Contrasting B, V, and i-dropout samples culled from the deeper data with those obtained from the shallower GOODS fields, we demonstrate that the shallower data are essentially complete at bright magnitudes to z<5.5 and that the principal effect of depth is to add objects at the magnitude limit. This indicates that high redshift galaxies are compact in size (~0.1-0.3") and that large (>0.4", >3 kpc) low surface brightness galaxies are rare. A simple comparison of the half-light radii of the HDF-N + HDF-S U-dropouts with B, V, and i-dropouts from the UDF shows that the sizes follow a (1+z)^{-1.05+/-0.21} scaling towards high redshift. A more rigorous measurement compares different scalings of our U-dropout sample with the mean profiles for a set of intermediate magnitude (26.0<z_{850,AB}<27.5) i-dropouts from the UDF. The best-fit is found with a (1+z)^{-0.94_{-0.25} ^{+0.19}} size scaling (for fixed luminosity). This result is then verified by repeating this experiment with different size measures, low redshift samples, and magnitude ranges. Very similar scalings are found for all comparisons. A robust measurement of size evolution is thereby demonstrated for galaxies from z~6 to z~2.5 using data from the UDF.

The properties of Ultra Compact Dwarf (UCD) galaxy candidates in Abell 1689 (z=0.183) are investigated, based on deep high resolution ACS images. A UCD candidate has to be unresolved, have i<28 (M_V<-11.5) mag and satisfy color limits derived from Bayesian photometric redshifts. We find 160 UCD candidates with 22<i<28 mag. It is estimated that about 100 of these are cluster members, based on their spatial distribution and photometric redshifts. For i>26.8 mag, the radial and luminosity distribution of the UCD candidates can be explained well by Abell 1689's globular cluster (GC) system. For i<26.8 mag, there is an overpopulation of 15 +/- 5 UCD candidates with respect to the GC luminosity function. For i<26 mag, the radial distribution of UCD candidates is more consistent with the dwarf galaxy population than with the GC system of Abell 1689. The UCD candidates follow a color-magnitude trend with a slope similar to that of Abell 1689's genuine dwarf galaxy population, but shifted fainter by about 2-3 mag. Two of the three brightest UCD candidates (M_V ~ -17 mag) are slightly resolved. At the distance of Abell 1689, these two objects would have King-profile core radii of ~35 pc and r_eff ~300 pc, implying luminosities and sizes 2-3 times those of M32's bulge. Additional photometric redshifts obtained with late type stellar and elliptical galaxy templates support the assignment of these two resolved sources to Abell 1689. Our findings imply that in Abell 1689 there are at least 10 UCDs with M_V<-12.7 mag. Compared to the UCDs in the Fornax cluster they are brighter, larger and have colors closer to normal dwarf galaxies. This suggests that they may be in an intermediate stage of the stripping process. Spectroscopy is needed to definitely confirm the existence of UCDs in Abell 1689.

(Abridged) We study the internal color properties of a morphologically selected sample of spheroidal galaxies taken from HST/ACS ERO program of UGC 10214 (``The Tadpole''). By taking advantage of the unprecedented high resolution of the ACS in this very deep dataset we are able to characterize spheroids at sub-arcseconds scales. Using the V_606W and I_814W bands, we construct V-I color maps and extract color gradients for a sample of spheroids at I_814W < 24 mag. We investigate the existence of a population of morphologically classified spheroids which show extreme variation in their internal color properties similar to the ones reported in the HDFs. These are displayed as blue cores and inverse color gradients with respect to those accounted from metallicity variations. Following the same analysis we find a similar fraction of early-type systems (~30%-40%) that show non-homologous internal colors, suggestive of recent star formation activity. We present two statistics to quantify the internal color variation in galaxies and for tracing blue cores, from which we estimate the fraction of non-homogeneous to homogeneous internal colors as a function of redshift up to z<1.2. We find that it can be described as about constant as a function of redshift, with a small increase with redshift for the fraction of spheroids that present strong color dispersions. The implications of a constant fraction at all redshifts suggests the existence of a relatively permanent population of evolving spheroids up to z~1. We discuss the implications of this in the context of spheroidal formation.

We report on the i-dropouts detected in two exceptionally deep ACS fields (B_{435}, V_{606}, i_{775}, and z_{850} with 10 sigma limits of 28.8, 29.0, 28.5, and 27.8, respectively) taken in parallel with the UDF NICMOS observations. Using an i-z>1.4 cut, we find 30 i-dropouts over 21 arcmin^2 down to z_AB=28.1, or 1.4 i-dropouts arcmin^{-2}, with significant field-to-field variation (as expected from cosmic variance). This extends i-dropout searches some ~0.9^m further down the luminosity function than was possible in the GOODS field, netting a ~7x increase in surface density. An estimate of the size evolution for UV bright objects is obtained by comparing the composite radial flux profile of the bright i-dropouts (z<27.2) with scaled versions of the HDF-N + HDF-S U-dropouts. The best-fit is found with a (1+z)^{-1.57_{-0.53} ^{+0.50}} scaling in size (for fixed luminosity), extending lower redshift (1<z<5) trends to z~6. Adopting this scaling and the brighter i-dropouts from both GOODS fields, we make incompleteness estimates and construct a z~6 LF in the rest-frame continuum UV (~1350 A) over a 3.5 magnitude baseline, finding a shape consistent with that found at lower redshift. To evaluate the evolution in the LF from z~3.8, we make comparisons against different scalings of a lower redshift B-dropout sample. Though a strong degeneracy is found between luminosity and density evolution, our best-fit model scales as (1+z)^{-2.8} in number and (1+z)^0.1 in luminosity, suggesting a rest-frame continuum UV luminosity density at z~6 which is just 0.38_{-0.07} ^{+0.09}x that at z~3.8. Our inclusion of size evolution makes the present estimate lower than previous z~6 estimates.

The unique high-resolution wide-field imaging capabilities of HST with ACS have allowed the characterization of galaxies at redshift 6, less than 1 Gyr from recombination. The dropout technique, applied to deep ACS i, z images in the RCDS 1252-2927, GOODS and UDF-Parallel fields has yielded large samples of these objects, allowing a detailed determination of their properties (e.g., size, color) and meaningful comparisons against lower redshift dropout samples. The use of cloning techniques has enabled us to control for many of the strong selection biases that affect the study of high redshift populations. A clear trend of size with redshift has been identified, and its impact on the luminosity density and star formation rate can be estimated. There is a significant, though modest, decrease in the star formation rate from redshifts z~2.5 out through z~6. The latest data also allow for the first robust determination of the luminosity function at z~6.

We report the discovery of a strongly lensed arc system around a field elliptical galaxy in Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) images of a parallel field observed during NICMOS observations of the HST Ultra-Deep Field. The ACS parallel data comprise deep imaging in the F435W, F606W, F775W, and F850LP bandpasses. The main arc is at a radius of 1.6 arcsec from the galaxy center and subtends about 120 deg. Spectroscopic follow-up at Magellan Observatory yields a redshift z=0.6174 for the lensing galaxy, and we photometrically estimate z_phot = 2.4\pm0.3 for the arc. We also identify a likely counter-arc at a radius of 0.6 arcsec, which shows structure similar to that seen in the main arc. We model this system and find a good fit to an elliptical isothermal potential of velocity dispersion $\sigma \approx 300$ \kms, the value expected from the fundamental plane, and some external shear. Several other galaxies in the field have colors similar to the lensing galaxy and likely make up a small group.

We have discovered three globular clusters beyond the Holmberg radius in Hubble Space Telescope Advanced Camera for Surveys images of the gas-rich dark matter dominated blue compact dwarf galaxy NGC2915. The clusters, all of which start to resolve into stars, have M_{V606} = -8.9 to -9.8 mag, significantly brighter than the peak of the luminosity function of Milky Way globular clusters. Their colors suggest a metallicity [Fe/H] ~ -1.9 dex, typical of metal-poor Galactic globular clusters. The specific frequency of clusters is at a minimum normal, compared to spiral galaxies. However, since only a small portion of the system has been surveyed it is more likely that the luminosity and mass normalized cluster content is higher, like that seen in elliptical galaxies and galaxy clusters. This suggests that NGC2915 resembles a key phase in the early hierarchical assembly of galaxies - the epoch when much of the old stellar population has formed, but little of the stellar disk. Depending on the subsequent interaction history, such systems could go on to build-up larger elliptical galaxies, evolve into normal spirals, or in rare circumstances remain suspended in their development to become systems like NGC2915.

We present the analysis of the faint galaxy population in the Advanced Camera for Surveys (ACS) Early Release Observation fields VV 29 (UGC 10214) and NGC 4676. Here we attempt to thoroughly consider all aspects relevant for faint galaxy counting and photometry, developing methods which are based on public software and that are easily reproducible by other astronomers. Using simulations we determine the best SExtractor parameters for the detection of faint galaxies in deep HST observations, paying special attention to the issue of deblending, which significantly affects the normalization and shape of the number count distribution. We confirm, as claimed by Bernstein, Freedman and Madore (2002), that Kron-like magnitudes, such as the ones generated by SExtractor, can miss more than half of the light of faint galaxies, what dramatically affects the slope of the number counts. We present catalogs for the VV 29 and NGC 4676 fields with photometry in the g,V and I bands. We also show that combining the bayesian software BPZ with superb ACS data and new spectral templates enables us to estimate reliable photometric redshifts for a significant fraction of galaxies with as few as three filters. After correcting for selection effects, we measure slopes of 0.32+- 0.01 for 22 < g < 28, 0.34+-0.01 for 22< V <27.5 and 0.33+-0.01 for 22 < I < 27. The counts do not flatten (except perhaps in the g-band), up to the depth of our observations. We find that the faint counts m_{AB}> 25.5 can be well approximated in all our filters by a passive luminosity evolution model based on the COMBO-17 luminosity function (\alpha=-1.5), with a strong merging rate following the prescription of Glazebrook et al. (1994), \phi^*\propto (1+Qz), with Q=4.

We investigate the color-magnitude (CM) relation of galaxies in the distant X-ray selected cluster RDCS1252.9-2927 at z=1.24 using images obtained with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope in the F775W and F850LP bandpasses. We select galaxies based on morphological classifications extending about 3.5 mag down the galaxy luminosity function, augmented by spectroscopic membership information. At the core of the cluster is an extensive early-type galaxy population surrounding a central pair of galaxies that show signs of dynamical interaction. The early-type population defines a tight sequence in the CM diagram, with an intrinsic scatter in observed (i-z) of $0.029\pm0.007$ mag based on 52 galaxies, or $0.024\pm0.008$ mag for $\sim30$ ellipticals. Simulations using the latest stellar population models indicate an age scatter for the ellipticals of about 34%, with a mean age $\taul\gta2.6$ Gyr (corresponding to $z_L\gta2.7$), and the last star formation occurring at $z_{end}\gta1.5$. Transforming to rest-frame $(U-B)$, we conclude that the slope and scatter in the CM relation for morphologically selected early-type galaxies show little or no evidence for evolution out to $z\approx1.2$. Thus, elliptical galaxies were already well established in X-ray luminous clusters when the universe was a third of its present age.

Using an i-z dropout criterion, we determine the space density of z~6 galaxies from two deep ACS GTO fields with deep optical-IR imaging. A total of 23 objects are found over 46 arcmin^2, or ~0.5 objects/arcmin^2 down to z~27.3 (6 sigma; all AB mag) (including one probable z~6 AGN). Combining deep ISAAC data for our RDCS1252-2927 field (J~25.7 and Ks~25.0 (5 sigma)) and NICMOS data for the HDF North (JH~27.3 (5 sigma)), we verify that these dropouts have flat spectral slopes. i-dropouts in our sample range in luminosity from ~1.5 L* (z~25.6) to ~0.3 L* (z~27.3) with the exception of one very bright candidate at z~24.2. The half-light radii vary from 0.09" to 0.29", or 0.5 kpc to 1.7 kpc. We derive the z~6 rest-frame UV luminosity density using three different procedures, each utilizing simulations based on a CDF South V dropout sample. First, we compare our findings with a no-evolution projection of this V-dropout sample. We find 23+/-25% more i-dropouts than we predict. Adopting previous results to z~5, this works out to a 20+/-29% drop in the luminosity density from z~3 to z~6. Second, we use these same V-dropout simulations to derive a selection function for our i-dropout sample and compute the UV-luminosity density (7.2+/-2.5 x 10^25 ergs/s/Hz/Mpc^3 down to z~27). We find a 39+/-21% drop over the same redshift range. This is our preferred value and suggests a star formation rate of 0.0090+/-0.0031 M_sol/yr/Mpc^3 to z~27, or ~0.036+/- 0.012 M_sol/yr/Mpc^3 extrapolating the LF to the faint limit. Third, we follow a very similar procedure, but assume no incompleteness, finding a luminosity density which is ~2-3X lower. This final estimate constitutes a lower limit. All three estimates are within the canonical range of luminosity densities necessary for reionization of the universe at this epoch. (abridged)

The evolution of high redshift galaxies in the two Hubble Deep Fields, HDF-N and HDF-S, is investigated using a cloning technique that replicates z~ 2-3 U dropouts to higher redshifts, allowing a comparison with the observed B and V dropouts at higher redshifts (z ~ 4-5). We treat each galaxy selected for replication as a set of pixels that are k-corrected to higher redshift, accounting for resampling, shot-noise, surface-brightness dimming, and the cosmological model. We find evidence for size evolution (a 1.7x increase) from z ~ 5 to z ~ 2.7 for flat geometries (Omega_M+Omega_LAMBDA=1.0). Simple scaling laws for this cosmology predict that size evolution goes as (1+z)^{-1}, consistent with our result. The UV luminosity density shows a similar increase (1.85x) from z ~ 5 to z ~ 2.7, with minimal evolution in the distribution of intrinsic colors for the dropout population. In general, these results indicate less evolution than was previously reported, and therefore a higher luminosity density at z ~ 4-5 (~ 50% higher) than other estimates. We argue the present technique is the preferred way to understand evolution across samples with differing selection functions, the most relevant differences here being the color cuts and surface brightness thresholds (e.g., due to the (1+z)^4 cosmic surface brightness dimming effect).

Multicolor coronagraphic images of the circumstellar disk around HD141569A have been obtained with the Hubble Space Telescope (HST) Advanced Camera for Surveys. B, V, and I images show that the disk's previously-described multiple-ring structure is actually a continuous distribution of dust with a tightly-wound spiral structure. Two, more open spiral arms extend from the disk, one of which appears to reach the nearby binary star HD141569BC. Diffuse dust is seen up to 1200AU from HD141569A. Although planets may exist in the inner region of the disk, tidal interaction with HD141569BC seems more likely to be the cause of these phenomena. The disk appears redder than the star (B-V = 0.21 and V-I = 0.25), and its color is spatially uniform. A scattering asymmetry factor of g = 0.25-0.35 is derived. The azimuthal density distribution is asymmetric, varying by a factor of ~3 at some radii.

We present observations of the first two supernovae discovered with the recently installed Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. The supernovae were found in Wide Field Camera images of the Hubble Deep Field North taken with the F775W, F850LP, and G800L optical elements as part of the ACS guaranteed time observation program. Spectra extracted from the ACS G800L grism exposures confirm that the objects are Type Ia supernovae (SNe Ia) at redshifts z=0.47 and z=0.95. Follow-up HST observations have been conducted with ACS in F775W and F850LP and with NICMOS in the near-infrared F110W bandpass, yielding a total of 9 flux measurements in the 3 bandpasses over a period of 50 days in the observed frame. We discuss many of the important issues in doing accurate photometry with the ACS. We analyze the multi-band light curves using two different fitting methods to calibrate the supernovae luminosities and place them on the SNe Ia Hubble diagram. The resulting distances are consistent with the redshift-distance relation of the accelerating universe model, although evolving intergalactic grey dust remains as a less likely possibility. The relative ease with which these SNe Ia were found, confirmed, and monitored demonstrates the potential ACS holds for revolutionizing the field of high-redshift SNe Ia, and therefore of testing the accelerating universe cosmology and constraining the "epoch of deceleration".

The off-axis location of the Advanced Camera for Surveys (ACS) is the chief (but not sole) cause of strong geometric distortion in all detectors: the Wide Field Camera (WFC), High Resolution Camera (HRC), and Solar Blind Camera (SBC). Dithered observations of rich star cluster fields are used to calibrate the distortion. We describe the observations obtained, the algorithms used to perform the calibrations and the accuracy achieved.

We present the first Advanced Camera for Surveys (ACS) observations of young star clusters in the colliding/merging galaxy UGC 10214. The observations were made as part of the Early Release Observation (ERO) program for the newly installed ACS during service mission SM3B for the Hubble Space Telescope (HST). Many young star clusters can be identified in the tails of UGC 10214, with ages ranging from ~3 Myr to 10 Myr. The extreme blue V-I (F606W-F814W) colors of the star clusters found in the tail of UGC 10214 can only be explained if strong emission lines are included with a young stellar population. This has been confirmed by our Keck spectroscopy of some of these bright blue stellar knots. The most luminous and largest of these blue knots has an absolute magnitude of M_V = -14.45, with a half-light radius of 161 pc, and if it is a single star cluster, would qualify as a super star cluster (SSC). Alternatively, it could be a superposition of multiple scaled OB associations or clusters. With an estimated age of ~ 4-5 Myr, its derived mass is < 1.3 x 10^6 solar masses. Thus the young stellar knot is unbound and will not evolve into a normal globular cluster. The bright blue clusters and associations are much younger than the dynamical age of the tail, providing strong evidence that star formation occurs in the tail long after it was ejected. UGC 10214 provides a nearby example of processes that contributed to the formation of halos and intra-cluster media in the distant and younger Universe.

We present simple models for disk evolution based on two different approaches: a forward approach based on predictions generic to hierarchical models for structure formation (e.g., Mo, Mao, & White 1998) and a backwards approach based on detailed modeling of the Milky Way galaxy (e.g., Bouwens, Cayon, & Silk 1997). We normalize these models to local observations and predict high-redshift luminosities, sizes, circular velocities, and surface brightnesses. Both approaches yield somewhat similar predictions for size, surface brightness, and luminosity evolution though they clearly differ in the amount of number evolution. These predictions seem to be broadly consistent with the high-redshift observations of Simard et al. (1999), suggesting that the B-band surface brightness of disks has indeed evolved by ~1.5 mag from z~0 to z~1 similar to the models and is not an artifact of selection effects as previously claimed. We also find a lack of low surface brightness galaxies in several high redshift samples relative to model predictions based on local samples (de Jong & van der Kruit 1994; Mathewson, Ford, & Buchhorn 1992).

The optical-IR images of the Northern and Southern Hubble Deep Fields are used to measure the spectral and density evolution of early-type galaxies. The mean optical SED is found to evolve passively towards a mid F-star dominated spectrum by z ~ 2. We demonstrate with realistic simulations that hotter ellipticals would be readily visible if evolution progressed blueward and brightward at z > 2, following a standard IMF. The colour distributions are best fitted by a `red' IMF, deficient above ~2 M_solar and with a spread of formation in the range 1.5 < z_f < 2.5. Traditional age dating is spurious in this context, a distant elliptical can be young but appear red, with an apparent age >3 Gyrs independent of its formation redshift. Regarding density evolution, we demonstrate that the sharp decline in numbers claimed at z > 1 results from a selection bias against distant red galaxies in the optical, where the flux is too weak for morphological classification, but is remedied with relatively modest IR exposures revealing a roughly constant space density to z ~ 2. We point out that the lack of high mass star-formation inferred here and the requirement of metals implicates cooling-flows of pre-enriched gas in the creation of the stellar content of spheroidal galaxies. Deep-field X-ray images will be very helpful to examine this possibility.

A substantial population of high redshift early-type galaxies is detected in very deep UBVRIJHK images towards the HDF-South. Four elliptical profile galaxies are identified in the redshift range z=1-2, all with very red SEDs, implying ages of >2 Gyrs for standard passive evolution. We also find later type IR-luminous galaxies at similarly high redshift, (10 objects with z>1, H<25), whose weak UV emission implies single burst ages of >1 Gyr. The number and luminosity-densities of these galaxies are comparable with the local E/SO-Sbc populations for \Omega_m>0.2, and in the absence of a significant cosmological constant, we infer that the major fraction of luminous Hubble-sequence galaxies have evolved little since z~2. A highly complete photometric redshift distribution is constructed to H=25 (69 galaxies) showing a broad spread of redshift, peaking at z~1.5, in reasonable agreement with some analyses of the HDF. Five `dropout' galaxies are detected at z~3.8, which are compact in the IR, ~0.5 kpc/h at rest 3500\AA. No example of a blue IR luminous elliptical is found, restricting the star-formation epoch of ellipticals to z>10 for a standard IMF and modest extinction.

The forwards approach to galaxy formation and evolution is extremely powerful but leaves several questions unanswered. Foremost among these is the origin of disks. A backwards approach is able to provide a more realistic treatment of star formation and feedback and provides a practical guide to eventually complement galaxy formation ab initio.

In this paper, we present three qualitatively different scenarios for bulge formation: a secular evolution model in which bulges form after disks and undergo several central starbursts, a primordial collapse model in which bulges and disks form simultaneously, and an early bulge formation model in which bulges form prior to disks. We normalize our models to the local z=0 observations of de Jong & van der Kruit (1994) and Peletier & Balcells (1996) and make comparisons with high redshift observations. We consider model predictions relating directly to bulge-to-disk properties. As expected, smaller bulge-to-disk ratios and bluer bulge colors are predicted by the secular evolution model at all redshifts, although uncertainties in the data are currently too large to differentiate strongly between the models.

Bulges, often identified with the spheroidal component of a galaxy, have a complex pedigree. Massive bulges are generally red and old, but lower mass bulges have broader dispersions in color that may be correlated with disk colors. This suggests different formation scenarios. I will review possible formation sequences for bulges, describe the various signatures that distinguish these scenarios, and discuss implications for the high redshift universe.

In a companion paper we outlined a methodology for generating parameter-free, model-independent ``no-evolution'' fields of faint galaxy images, demonstrating the need for significant evolution in the HDF at faint magnitudes. Here we incorporate evolution into our procedure, by transforming the input bright galaxy images with redshift, for comparison with the HDF at faint magnitudes. Pure luminosity evolution is explored assuming that galaxy surface brightness evolves uniformly, at a rate chosen to reproduce the I-band counts. This form of evolution exacerbates the size discrepancy identified by our no-evolution simulations, by increasing the area of a galaxy visible to a fixed isophote. Reasonable dwarf-augmented models are unable to generate the count excess invoking moderate rates of stellar evolution. A plausible fit to the counts and sizes is provided by `mass-conserving' density-evolution, consistent with small-scale hierarchical growth, where the product of disk area and space density is conserved with redshift. Here the increased surface brightness generated by stellar evolution is accomodated by the reduced average galaxy size, for a wide range of geometries. These models are useful for assessing the limitations of the HDF images, by calculating their rates of incompleteness and the degree of over-counting. Finally we demonstrate the potential for improvement in quantifying evolution at fainter magnitudes using the HST Advanced Camera, with its superior UV and optical performance.

We present a model-independent method of quantifying galaxy evolution in high-resolution images, which we apply to the Hubble Deep Field (HDF). Our procedure is to k-correct all pixels belonging to the images of a complete set of bright galaxies and then to replicate each galaxy image to higher redshift by the product of its space density, 1/V_{max}, and the cosmological volume. The set of bright galaxies is itself selected from the HDF, because presently the HDF provides the highest quality UV images of a redshift-complete sample of galaxies (31 galaxies with I<21.9, \bar{z}=0.5, and for which V/V_{max} is spread fairly). These galaxies are bright enough to permit accurate pixel-by-pixel k-corrections into the restframe UV (\sim 2000 A). We match the shot noise, spatial sampling and PSF smoothing of the HDF data, resulting in entirely empirical and parameter-free ``no-evolution'' deep fields of galaxies for direct comparison with the HDF. In addition, the overcounting rate and the level of incompleteness can be accurately quantified by this procedure. We obtain the following results. Faint HDF galaxies (I>24) are much smaller, more numerous, and less regular than our ``no-evolution'' extrapolation, for any interesting geometry. A higher proportion of HDF galaxies ``dropout'' in both U and B, indicating that some galaxies were brighter at higher redshifts than our ``cloned'' z\sim0.5 population.

We present a model-independent method of quantifying galaxy evolution in high- resolution images, which we apply to the Hubble Deep Field (HDF). Our procedure is to k-correct the pixels belonging to the images of a complete set of bright galaxies and then to replicate each galaxy image to higher redshift by the product of its space density, 1/V_{max}, and the cosmological volume. The set of bright galaxies is itself selected from the HDF because presently the HDF provides the highest quality UV images of a redshift-complete sample of galaxies (31 galaxies with I<21.9, \bar{z}=0.5, and for which V/V_{max} is spread fairly). These galaxies are bright enough to permit accurate pixel-by-pixel k-corrections into the restframe UV (\sim 2000 A). We match the shot noise, spatial sampling and PSF smoothing of the HDF, resulting in entirely empirical and parameter free ``no-evolution'' deep fields of galaxies for direct comparison with the HDF. We obtain the following results. Faint HDF galaxies (I>24) are much smaller, more numerous, and less regular than our ``no-evolution'' extrapolation, for any relevant geometry. A higher proportion of HDF galaxies ``dropout'' in both U and B, indicating that some galaxies were brighter at higher redshifts than our ``cloned'' z\sim0.5 population. By simple image transformations we demonstrate that bolometric luminosity evolution generates galaxies which are too large and the contribution of any evolving dwarf population is uninterestingly small. A plausible fit is provided by `mass-conserving' density-evolution, consistent with hierarchical growth of small-scale structure. Finally, we show the potential for improvement using the Advanced Camera, with its superior UV and optical performance.

This paper presents the results of simulated NGST observations. The goal of this work is to study the detectability of very high redshift galaxies and to try to characterize Wide Deep-Field observations. The main points made in this paper are 1) high-z galaxies observable by NGST would be similar to local starburst galaxies 2) even if they exist, large spiral galaxies similar to local ones would not be observed by a NGST-like telescope and 3) the density of high-z galaxies detected so far is apparently moderate.

We develop an idealized inside-out formation model for disk galaxies to include a realistic mix of galaxy types and luminosities that provides a fair match to the traditional observables. The predictions of our infall models are compared against identical models with no-size evolution by generating fully realistic simulations of the HDF, from which we recover the angular size distributions. We find that our infall models produce nearly identical angular size distributions to those of our no-size evolution models in the case of a Omega = 0 geometry but produce slightly smaller sizes in the case of a Omega = 1 geometry, a difference we associate with the fact that there is a different amount of cosmic time in our two models for evolving to relatively low redshifts (z \approx 1-2). Our infall models also predict a slightly smaller (11% - 29%) number of large (disk scale lengths > 4 h_{50} ^{-1} kpc) galaxies at z \approx 0.7 for the CFRS as well as different increases in the central surface brightness of the disks for early-type spirals, the infall model predicting an increase by 1.2 magnitudes out to z \approx 2 (Omega = 0), 1 (Omega = 1), while our no-size evolution models predict an increase of only 0.5 magnitude. This result suggests that infall models could be important for explaining the 1.2-1.6 magnitude increase in surface brightness reported by Schade et al. (1995, 1996a, 1996b).

A constant age population of blue galaxies, postulated in the model of Gronwall & Koo (1995), seems to provide an attractive explanation of the excess of very blue galaxies in the deep galaxy counts. Such a population may be generated by a set of galaxies with cycling star formation rates, or at the other extreme, be maintained by the continual formation of new galaxies which fade after they reach the age specified in the Gronwall and Koo model. For both of these hypotheses, we have calculated the luminosity functions including the respective selection criteria, the redshift distributions, and the number counts in the B_J and K bands. We find a substantial excess in the number of galaxies at low redshift (0 < z < 0.05) over that observed in the CFH redshift survey (Lilly et al. 1995) and at the faint end of the Las Campanas luminosity function (Lin et al. 1996). Passive or mild evolution fails to account for the deep galaxy counts because of the implications for low redshift determinations of the I-selected redshift distribution and the r-selected luminosity function in samples where the faded counterparts of the star-forming galaxies would be detectable.