Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level' galaxy formation physics to `top-level' cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

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

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

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

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

The James Webb Space Telescope (JWST) is expected to enable transformational progress in studying galaxy populations in the very early Universe, during the Epoch of Reionization (EoR). A critical parameter for understanding the sources that reionized the Universe is the Lyman-continuum production efficiency, $\xi_{\rm ion}$, defined as the rate of production of ionizing photons divided by the intrinsic UV luminosity. In this work, we combine self-consistent star formation and chemical enrichment histories predicted by semi-analytic models of galaxy formation with stellar population synthesis (SPS) models to predict the expected dependence of $\xi_{\rm ion}$ on galaxy properties and cosmic epoch from $z=4-10$. We then explore the sensitivity of the production rate of ionizing photons, $\dot{N}_{\rm ion}$, to the choice of SPS model and the treatment of stellar feedback in our galaxy formation model. We compare our results to those of other simulations, constraints from empirical models, and observations. We find that adopting SPS models that include binary stars predict about a factor of two more ionizing radiation than models that only assume single stellar populations. We find that UV-faint, low-mass galaxies have values of $\xi_{\rm ion}$ about 0.25 dex higher than those of more massive galaxies, but find weak evolution with cosmic time, about 0.2 dex from $z\sim 12-4$ at fixed rest-UV luminosity. We provide predictions of $\dot{N}_{\rm ion}$ as a function of $M_{\rm h}$ and a number of other galaxy properties. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

Looking for bright galaxies born in the early universe is fundamental to investigating the Epoch of Reionization, the era when the first stars and galaxies ionized the intergalactic medium. We utilize Hubble Space Telescope pure parallel imaging to select galaxy candidates at a time 500 to 650 million years after the Big Bang, which corresponds to redshifts z ~ 8-10. These data come from the Brightest of Reionizing Galaxies Survey (BoRG) Cycle 22 dataset, which consists of pure-parallel imaging in ~ 90 different lines of sight that sum up to an area of ~ 420 arcmin^2. This survey uses five filters and has the advantage (compared to the Cycle 21 BoRG program) of including imaging in the JH140 band, covering continuous wavelengths from the visible to near-infrared (lambda = 0.35um - 1.7um). This allows us to perform reliable selection of galaxies at z>8 using the photometric redshift technique. We use these galaxy candidates to constrain the bright end of the rest-frame ultraviolet luminosity function in this epoch. These candidates are excellent targets for follow-up observations, particularly with the James Webb Space Telescope.

We study the evolution in the number density of the highest mass galaxies over $0.4<z<1.5$ (covering 9 Gyr). We use the Spitzer/HETDEX Exploratory Large-Area (SHELA) Survey, which covers 17.5 $\mathrm{deg}^2$ with eight photometric bands spanning 0.3-4.5 $\mu$m within the SDSS Stripe 82 field. This size produces the lowest counting uncertainties and cosmic variance yet for massive galaxies at $z\sim1.0$. We study the stellar mass function (SMF) for galaxies with $\log(M_\ast/M_\odot)>10.3$ using a forward-modeling method that fully accounts for statistical and systematic uncertainties on stellar mass. From $z$=0.4 to 1.5 the massive end of the SMF shows minimal evolution in its shape: the characteristic mass ($M^\ast$) evolves by less than 0.1 dex ($\pm$0.05 dex); the number density of galaxies with $\log (M_\ast/M_\odot) >11$ stays roughly constant at $\log (n/\mathrm{Mpc}^{-3})$ $\simeq$ $-$3.4 ($\pm$0.05), then declines to $\log (n/\mathrm{Mpc}^{-3})$=$-$3.7 ($\pm$0.05) at $z$=1.5. We discuss the uncertainties in the SMF, which are dominated by assumptions in the star formation history and details of stellar population synthesis models for stellar mass estimations. For quiescent galaxies, the data are consistent with no (or slight) evolution ($\lesssim0.1$ dex) in the characteristic mass nor number density from $z\sim 1.5$ to the present. This implies that any mass growth (presumably through "dry" mergers) of the quiescent massive galaxy population must balance the rate of mass losses from late-stage stellar evolution and the formation of quenching galaxies from the star-forming population. We provide a limit on this mass growth from $z=1.0$ to 0.4 of $\Delta M_\ast/M_\ast\leq$ 45% (i.e., $\simeq0.16$ dex) for quiescent galaxies more massive than $10^{11}$ $M_\odot$.

We present new estimates of the luminosity function (LF) and star formation rate density (SFRD) for an H$\alpha$ selected sample at $z\sim0.62$ from the Deep And Wide Narrow-band (DAWN) survey. Our results are based on a new H$\alpha$ sample in the extended COSMOS region (compared to Coughlin et al. 2018) with the inclusion of flanking fields, resulting in a total area coverage of $\sim$1.5 deg$^2$. A total of 241 H$\alpha$ emitters were selected based on robust selection criteria using spectro-photometric redshifts and broadband color-color classification. We explore the effect of different dust correction prescriptions by calculating the LF and SFRD using a constant dust extinction correction, A{$_{\textrm{H}\alpha}=1$} mag, a luminosity-dependent correction, and a stellar-mass dependent correction. The resulting H$\alpha$ LFs are well fitted using Schechter functions with best-fit parameters: L$^*=10^{42.24}$ erg s$^{-1}$, $\phi^*=10^{-2.85}$ Mpc$^{-3}$, $\alpha = -1.62$ for constant dust correction, L$^*=10^{42.31}$ erg s$^{-1}$, $\phi^*=10^{-2.8}$ Mpc$^{-3}$, $\alpha=-1.39$ for luminosity-dependent dust correction, and L$^*=10^{42.36}$ erg s$^{-1}$, $\phi^*=10^{-2.91}$ Mpc$^{-3}$, $\alpha = -1.48$, for stellar mass-dependent dust correction. The deep and wide nature of the DAWN survey effectively samples H$\alpha$ emitters over a wide range of luminosities, thereby providing better constraints on both the faint and bright end of the LF. Also, the SFRD estimates $\rho_{\textrm{SFR}}=10^{-1.39}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (constant dust correction), $\rho_{\textrm{SFR}}=10^{-1.47}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (luminosity-dependent dust correction), and $\rho_{\textrm{SFR}}=10^{-1.49}$ M$_{\odot}$yr$^{-1}$Mpc$^{-3}$ (stellar mass-dependent dust correction) are in good agreement with the evolution of SFRD across redshifts ($0 < z < 2$) seen from previous H$\alpha$ surveys.

We examine the stacked thermal Sunyaev-Zel\text{'}dovich (SZ) signals for a sample of galaxy cluster candidates from the Spitzer-HETDEX Exploratory Large Area (SHELA) Survey, which are identified in combined optical and infrared SHELA data using the redMaPPer algorithm. We separate the clusters into three richness bins, with average photometric redshifts ranging from 0.70 to 0.80. The richest bin shows a clear temperature decrement at 148 GHz in the Atacama Cosmology Telescope data, which we attribute to the SZ effect. All richness bins show an increment at 220 GHz, which we attribute to dust emission from cluster galaxies. We correct for dust emission using stacked profiles from Herschel Stripe 82 data, and allow for synchrotron emission using stacked profiles created by binning source fluxes from NVSS data. We see dust emission in all three richness bins, but can only confidently detect the SZ decrement in the highest richness bin, finding $M_{500}$ = $8.7^{+1.7}_{-1.3} \times 10^{13} M_\odot$. Neglecting the correction for dust depresses the inferred mass by 26 percent, indicating a partial fill-in of the SZ decrement from thermal dust and synchrotron emission by the cluster member galaxies. We compare our corrected SZ masses to two redMaPPer mass--richness scaling relations and find that the SZ mass is lower than predicted by the richness. We discuss possible explanations for this discrepancy, and note that the SHELA richnesses may differ from previous richness measurements due to the inclusion of IR data in redMaPPer.

While most of the inter-galactic medium (IGM) today is permeated by ionized hydrogen, it was largely filled with neutral hydrogen for the first 700 million years after the Big Bang. The process that ionized the IGM (cosmic reionization) is expected to be spatially inhomogeneous, with fainter galaxies playing a significant role. However, we still have only a few direct constraints on the reionization process. Here we report the first spectroscopic confirmation of two galaxies and very likely a third galaxy in a group (hereafter EGS77) at redshift z = 7.7, merely 680 Myrs after the Big Bang. The physical separation among the three members is < 0.7 Mpc. We estimate the radius of ionized bubble of the brightest galaxy to be about 1.02 Mpc, and show that the individual ionized bubbles formed by all three galaxies likely overlap significantly, forming a large yet localized ionized region, which leads to the spatial inhomogeneity in the reionization process. It is striking that two of three galaxies in EGS77 are quite faint in the continuum, thanks to our selection of reionizing sources using their Lyman-alpha line emission. Indeed, one is the faintest spectroscopically confirmed galaxy yet discovered at such high redshifts. Our observations provide direct constraints in the process of cosmic reionization, and allow us to investigate the properties of sources responsible for reionizing the universe.

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

We present a catalog of 208 $0.3 < z < 2.1$ Emission Line Galaxies (ELG) selected from 1D slitless spectroscopy obtained using Hubble's WFC3 G102 grism, as part of the Faint Infrared Grism Survey (FIGS). We identify ELG candidates by searching for significant peaks in all continuum-subtracted G102 spectra, and, where possible, confirm candidates by identifying consistent emission lines in other available spectra or with published spectroscopic redshifts. We provide derived emission line fluxes and errors, redshifts, and equivalent widths (EW) for H$\alpha$ $\lambda6563$, [OIII]$\lambda\lambda4959,5007$, and [OII]$\lambda\lambda3727$ emission lines, for emission line galaxies down to AB(F105W) $ > 28$ and $> 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ line flux. We use the resulting line catalog to investigate a possible relationship between line emission and a galaxy's environment. We use 7th-nearest-neighbor distances to investigate the typical surroundings of ELGs compared to non-ELGs, and we find that [OIII] emitters are preferentially found at intermediate galaxy densities near galaxy groups. We characterize these ELGs in terms of the galaxy specific star formation rate (SSFR) versus stellar mass, and find no significant influence of environment on that relation. We calculate star formation rates (SFR), and find no dependence of SFR on local galaxy surface density for $0.3 < z < 0.8$ H$\alpha$ emitters and for $0.8<z<1.3$ [OIII] emitters. We find similar rates of close-pair interaction between ELGs and non-ELGs. For galaxy surface densities $\Sigma \leq 30$ Mpc$^{-2}$, we find no consistent effect of environment on star formation.

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

We present the high-mass end of the galaxy stellar mass function using the largest sample to date (5,352) of star-forming galaxies with $M_{\star} > 10^{11} M_{\odot}$ at cosmic noon, $1.5 < z < 3.5$. This sample is uniformly selected across 17.2 deg$^2$ ($\sim$0.44 Gpc$^3$ comoving volume from $1.5 < z < 3.5$), mitigating the effects of cosmic variance and encompassing a wide range of environments. This area, a factor of 10 larger than previous studies, provides robust statistics at the high-mass end. Using multi-wavelength data in the Spitzer/HETDEX Exploratory Large Area (SHELA) footprint we find that the SHELA footprint star-forming galaxy stellar mass function is steeply declining at the high-mass end probing values as high as $\sim$$10^{-4}$ Mpc$^3$/dex and as low as $\sim$5$\times$$10^{-8}$ Mpc$^3$/dex across a stellar mass range of log($M_\star$/$M_\odot$) $\sim$ 11 - 12. We compare our empirical star-forming galaxy stellar mass function at the high mass end to three types of numerical models: hydrodynamical models from IllustrisTNG, abundance matching from the UniverseMachine, and three different semi-analytic models (SAMs; SAG, SAGE, GALACTICUS). At redshifts $1.5 < z < 3.5$ we find that results from IllustrisTNG and abundance matching models agree within a factor of $\sim$2 to 10, however the three SAMs strongly underestimate (up to a factor of 1,000) the number density of massive galaxies. We discuss the implications of these results for our understanding of galaxy evolution.

We present a physical characterization of MMJ100026.36+021527.9 (a.k.a. ``MAMBO-9''), a dusty star-forming galaxy (DSFG) at $z=5.850\pm0.001$. This is the highest redshift unlensed DSFG (and fourth most distant overall) found to-date, and is the first source identified in a new 2mm blank-field map in the COSMOS field. Though identified in prior samples of DSFGs at 850$\mu$m-1.2mm with unknown redshift, the detection at 2mm prompted further follow-up as it indicated a much higher probability that the source was likely to sit at $z>4$. Deep observations from the Atacama Large Millimeter and submillimeter Array (ALMA) presented here confirm the redshift through the secure detection of $^{12}$CO($J\!=$6$\rightarrow$5) and p-H$_{2}$O(2$_{1,1}\!\rightarrow$2$_{0,2}$). MAMBO-9 is comprised of a pair of galaxies separated by 6kpc with corresponding star-formation rates of 590M$_\odot$yr$^{-1}$ and 220M$_\odot$yr$^{-1}$ total molecular hydrogen gas mass of (1.7$\pm$0.4)$\times10^{11}$M$_\odot$, dust mass of (1.3$\pm$0.3)$\times10^{9}$M$_\odot$ and stellar mass of (3.2$^{+1.0}_{-1.5}$)$\times10^{9}$M$_\odot$. The total halo mass, (3.3$\pm$0.8)$\times10^{12}$M$_\odot$, is predicted to exceed $>10^{15}$M$_\odot$ by $z=0$. The system is undergoing a merger-driven starburst which will increase the stellar mass of the system tenfold in $\tau_{\rm depl}=40-80$Myr, converting its large molecular gas reservoir (gas fraction of 96$^{+1}_{-2}$%) into stars. MAMBO-9 evaded firm spectroscopic identification for a decade, following a pattern that has emerged for some of the highest redshift DSFGs found. And yet, the systematic identification of unlensed DSFGs like MAMBO-9 is key to measuring the global contribution of obscured star-formation to the star-formation rate density at $z>4$, the formation of the first massive galaxies, and the formation of interstellar dust at early times ($<$1Gyr).

The long-anticipated James Webb Space Telescope (JWST) will be able to directly detect large samples of galaxies at very high redshift. Using the well-established, computationally efficient Santa Cruz semi-analytic model, with recently implemented multiphase gas partitioning and H2-based star formation recipes, we make predictions for a wide variety of galaxy properties for galaxy populations at $z$ = 4-10. In this work, we provide forecasts for the physical properties of high-redshift galaxies and links to their photometric properties. With physical parameters calibrated only to $z\sim0$ observations, our model predictions are in good agreement with current observational constraints on stellar mass and star formation rate distribution functions up to $z \sim 8$. We also provide predictions representing wide, deep, and lensed JWST survey configurations. We study the redshift evolution of key galaxy properties and the scaling relations among them. Taking advantage of our models' high computational efficiency, we study the impact of systematically varying the model parameters. All distribution functions and scaling relations presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We assemble an unbiased sample of 29 galaxies with [O II] $\lambda 3727$ and/or [O III] $\lambda 5007$ detections at $z < 0.15$ from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Pilot Survey (HPS). HPS finds galaxies without pre-selection based on their detected emission lines via integral field spectroscopy. Sixteen of these objects were followed up with the second-generation, low resolution spectrograph (LRS2) on the upgraded Hobby-Eberly Telescope. Oxygen abundances were then derived via strong emission lines using a Bayesian approach. We find most of the galaxies fall along the mass-metallicity relation derived from photometrically selected star forming galaxies in the Sloan Digital Sky Survey (SDSS). However, two of these galaxies have low metallicity (similar to the very rare green pea galaxies in mass-metallicity space). The star formation rates of this sample fall in an intermediate space between the SDSS star forming main sequence and the extreme green pea galaxies. We conclude that spectroscopic selection fills part of the mass-metallicity-SFR phase space that is missed in photometric surveys with pre-selection like SDSS, i.e., we find galaxies that are actively forming stars but are faint in continuum. We use the results of this pilot investigation to make predictions for the upcoming unbiased, large spectroscopic sample of local line emitters from HETDEX. With the larger HETDEX survey we will determine if galaxies selected spectroscopically without continuum brightness pre-selection have metallicities that fall on a continuum that bridges typical star forming and rarer, more extreme systems like green peas.

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

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

We utilize GIZMO, coupled with newly developed sub-grid models for Population~III (Pop~III) and Population~II (Pop~II), to study the legacy of star formation in the pre-reionization Universe. We find that the Pop~II star formation rate density (SFRD), produced in our simulation ($\sim 10^{-2}\ \Msun{\rm yr^{-1}\, Mpc^{-3}}$ at $z\simeq 10$), matches the total SFRD inferred from observations within a factor of $<2$ at $7\lesssim z \lesssim10$. The Pop~III SFRD, however, reaches a plateau at $\sim10^{-3}\ \Msun{\rm yr^{-1}\, Mpc^{-3}}$ by $z\approx10$, remaining largely unaffected by the presence of Pop~II feedback. At $z$=7.5, $\sim 20\%$ of Pop~III star formation occurs in isolated haloes which have never experienced any Pop~II star formation (i.e. primordial haloes). We predict that Pop~III-only galaxies exist at magnitudes $M_{\rm UV}\gtrsim -11$, beyond the limits for direct detection with the {\it James Webb Space Telescope (JWST)}. We assess that our stellar mass function (SMF) and UV luminosity function (UVLF) agree well with the observed low mass/faint-end behavior at $z=8$ and $10$. However, beyond the current limiting magnitudes, we find that both our SMF and UVLF demonstrate a deviation/turnover from the expected power-law slope ($M_{\rm UV,turn}= -13.4\pm1.1$ at $z$=10). This could impact observational estimates of the true SFRD by a factor of $2 (10)$ when integrating to $M_{\rm UV} = -$12 ($-$8) at $z\sim 10$, depending on integration limits. Our turnover correlates well with the transition from dark matter haloes dominated by molecular cooling to those dominated by atomic cooling, for a mass $M_{\rm halo}\approx 10^{8} \Msun$ at $z\simeq 10$.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

We present Keck/MOSFIRE $H$-band spectroscopy targeting C$\,$III] $\lambda$1907,1909 in a $z=7.5056$ galaxy previously identified via Ly$\alpha$ emission. We detect strong line emission at $1.621\pm0.002\,\mu$m with a line flux of ($2.63\pm0.52$)$\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. We tentatively identify this line as [C$\,$III] $\lambda$1907, but we are unable to detect C$\,$III] $\lambda$1909 owing to sky emission at the expected location. This gives a galaxy systemic redshift, $z_{sys}=7.5032\pm0.0003$, with a velocity offset to Ly$\alpha$ of $\Delta$v$_{Ly\alpha}$ = $88\pm27$ km s$^{-1}$. The ratio of combined C$\,$III]/Ly$\alpha$ is 0.30-0.45, one of the highest values measured for any $z>2$ galaxy. We do not detect Si$\,$III] $\lambda\lambda$1883, 1892, and place an upper limit on Si$\,$III]/C$\,$III] $<$ 0.35 ($2\sigma$). Comparing our results to photoionization models, the C$\,$III] equivalent width (W$_{\text{CIII]}} = 16.23\pm2.32\,$\AA), low Si$\,$III]/C$\,$III] ratio, and high implied [O$\,$III] equivalent width (from the $Spitzer$/IRAC [3.6]$-$[4.5]$\simeq$0.8 mag color) require sub-Solar metallicities ($Z\simeq0.1-0.2 Z_{\odot}$) and a high ionization parameter, log$\,$U $\gtrsim -1.5$. These results favor models that produce higher ionization, such as the BPASS models for the photospheres of high-mass stars, and that include both binary stellar populations and/or an IMF that extends to 300 M$_{\odot}$. The combined C$\,$III] equivalent width and [3.6]$-$[4.5] color are more consistent with ionization from young stars than AGN, however we cannot rule out ionization from a combination of an AGN and young stars. We make predictions for $James~Webb~Space~Telescope$ spectroscopy using these different models, which will ultimately test the nature of the ionizing radiation in this source.

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

One of the key questions in understanding the formation and evolution of galaxies is how starbursts affect the assembly of stellar populations in galaxies over time. We define a burst indicator ($\eta$), which compares a galaxy's star formation rates on short ($\sim10$ Myr) and long ($\sim100$ Myr) timescales. To estimate $\eta$, we apply the detailed time-luminosity relationship for H$\alpha$ and near-ultraviolet emission to simulated star formation histories (SFHs) from semi-analytic models and the Mufasa hydrodynamical cosmological simulations. The average of $\eta$ is not a good indicator of star formation stochasticity (burstiness); indeed, we show that this average should be close to zero unless the population has an average SFH which is rising or falling rapidly. Instead, the width of the $\eta$ distribution characterizes the burstiness of a galaxy population's recent star formation. We find this width to be robust to variations in stellar initial mass function and metallicity. We apply realistic noise and selection effects to the models to generate mock HST and JWST galaxy catalogs and compare these catalogs with 3D-HST observations of 956 galaxies at $0.65<z<1.5$ detected in H$\alpha$. Measurements of $\eta$ are unaffected by dust measurement errors under the assumption that $E(B-V)_\mathrm{stars}=0.44\,E(B-V)_\mathrm{gas}$ (i.e., $Q_\mathrm{sg}=0.44$). However, setting $Q_\mathrm{sg}=0.8^{+0.1}_{-0.2}$ removes an unexpected dependence of the average value of $\eta$ upon dust attenuation and stellar mass in the 3D-HST sample while also resolving disagreements in the distribution of star formation rates. However, even varying the dust law cannot resolve all discrepancies between the simulated and the observed galaxies.

We use cosmological hydrodynamic zoom-in simulations to study early structure formation in two dark matter (DM) cosmologies, the standard CDM model, and a thermal warm DM (WDM) model with a particle mass of $m_{\chi}c^{2}=3\ \mathrm{keV}$. We focus on DM haloes with virial masses $M\sim 10^{10}\ M_{\odot}$. We find that the first star formation activity is delayed by $\sim 200\ \mathrm{Myr}$ in the WDM model, with similar delays for metal enrichment and the formation of the second generation of stars. However, the differences between the two models in globally-averaged properties, such as star formation rate density and mean metallicity, decrease towards lower redshifts ($z\lesssim 10$). Metal enrichment in the WDM cosmology is restricted to dense environments, while low-density gas can also be significantly enriched in the CDM case. The free-free contribution from early structure formation at redshifts $z>6$ to the cosmic radio background (CRB) is $3_{-1.5}^{+13}\%$ ($8_{-3.5}^{+33}\%$) of the total signal inferred from radio experiments such as ARCADE 2, in the WDM (CDM) model. The direct detection of the $\mathrm{H_{2}}$ emission from early structure formation ($z\gtrsim 7.2$), originating from the low-mass haloes explored here, will be challenging even with the next generation of far-infrared space telescopes, unless the signal is magnified by at least a factor of 10 via gravitational lensing or shocks. However, more massive haloes with $M\gtrsim 10^{12}\ M_{\odot}$ may be observable for $z\gtrsim 10$, even without magnification, provided that our extrapolation from the scale of our simulated haloes is valid.

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

We constrain the stellar population properties of a sample of 52 massive galaxies, with stellar mass log Ms>10.5, over the redshift range 0.5<z<2 by use of observer-frame optical and near-infrared slitless spectra from HST's ACS and WFC3 grisms. The deep exposures (~100 ks) allow us to target individual spectra of massive galaxies to F160W=22.5AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via an MCMC method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Dt) that is found to dominate the variance of the retrieved parameters according to Principal Component Analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Dt appears to weakly anti-correlate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

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.)

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The Cosmic Dawn Intensity Mapper (CDIM) will transform our understanding of the era of reionization when the Universe formed the first stars and galaxies, and UV photons ionized the neutral medium. CDIM goes beyond the capabilities of upcoming facilities by carrying out wide area spectro-imaging surveys, providing redshifts of galaxies and quasars during reionization as well as spectral lines that carry crucial information on their physical properties. CDIM will make use of unprecedented sensitivity to surface brightness to measure the intensity fluctuations of reionization on large-scales to provide a valuable and complementary dataset to 21-cm experiments. The baseline mission concept is an 83-cm infrared telescope equipped with a focal plane of 24 \times 20482 detectors capable of R = 300 spectro-imaging observations over the wavelength range of 0.75 to 7.5 {\mu}m using Linear Variable Filters (LVFs). CDIM provides a large field of view of 7.8 deg2 allowing efficient wide area surveys, and instead of moving instrumental components, spectroscopic mapping is obtained through a shift-and-stare strategy through spacecraft operations. CDIM design and capabilities focus on the needs of detecting faint galaxies and quasars during reionization and intensity fluctuation measurements of key spectral lines, including Lyman-{\alpha} and H{\alpha} radiation from the first stars and galaxies. The design is low risk, carries significant science and engineering margins, and makes use of technologies with high technical readiness level for space observations.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role.

The epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from today's 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than today's deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.

The rest-frame UV emission from massive stars contains a wealth of information about the physical nature and conditions of star formation in galaxies. Using studies of the rest-frame UV, the past decade has witnessed the beginning of knowledge about the existence and properties of galaxies during the first few billion years after the Big Bang. This period of history corresponds to the formation of the first stars, the rapid formation of galaxy stellar populations, the reionization of the IGM, the production and dissemination of heavy elements, and the formation of the first black holes. Massive stars in these galaxies drive all of these events, and their light dominates the spectral energy distributions of galaxies. As we look to the 2020s, fundamental questions remain about the nature of these stellar populations and their evolution, from just before the peak of the cosmic star formation density (z~3), up to the epoch of reionization (z > 6). This next decade will provide transformative gains both in our ability to identify star-forming galaxies and accreting supermassive black holes at these early epochs with imaging surveys in the rest-frame UV (e.g., LSST, WFIRST). Ground-based, rest-frame UV spectroscopy on >20 m-class telescopes (e.g., GMT/TMT) offers the ability to investigate the astrophysical conditions in galaxies at the earliest cosmic times. This includes studies of the evolution in galaxy stellar populations, gas ionization (temperature, pressure), metallicity, and interstellar (and circumgalactic) gas kinematics and covering fractions. In this white paper, we describe the scientific prospects and the requirements for research in this area.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

We derive direct measurement gas-phase metallicities of $7.4 < 12 + \log(O/H) < 8.4$ for 14 low-mass Emission Line Galaxies (ELGs) at $0.3 < z < 0.8$ identified in the Faint Infrared Grism Survey (FIGS). We use deep slitless G102 grism spectroscopy of the Hubble Ultra Deep Field (HUDF), dispersing light from all objects in the field at wavelengths between 0.85 and 1.15 microns. We run an automatic search routine on these spectra to robustly identify 71 emission line sources, using archival data from VLT/MUSE to measure additional lines and confirm redshifts. We identify 14 objects with $0.3 < z < 0.8$ with measurable O[III]$\lambda$4363 \AA\ emission lines in matching VLT/MUSE spectra. For these galaxies, we derive direct electron-temperature gas-phase metallicities with a range of $7.4 < 12 + \log(O/H) < 8.4$. With matching stellar masses in the range of $10^{7.9} M_{\odot} < M_{\star} < 10^{10.4} M_{\odot}$, we construct a mass-metallicity (MZ) relation and find that the relation is offset to lower metallicities compared to metallicities derived from alternative methods (e.g.,$R_{23}$, O3N2, N2O2) and continuum selected samples. Using star formation rates (SFR) derived from the $H\alpha$ emission line, we calculate our galaxies' position on the Fundamental Metallicity Relation (FMR), where we also find an offset toward lower metallicities. This demonstrates that this emission-line-selected sample probes objects of low stellar masses but even lower metallicities than many comparable surveys. We detect a trend suggesting galaxies with higher Specific Star Formation (SSFR) are more likely to have lower metallicity. This could be due to cold accretion of metal-poor gas that drives star formation, or could be because outflows of metal-rich stellar winds and SNe ejecta are more common in galaxies with higher SSFR.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

Realizing the utility of Ly$\alpha$ emission to trace the evolution of the intergalactic medium (IGM) during the epoch of reionization requires deep spectroscopy across the boundary of optical and near-infrared (NIR) spectrographs at $z\sim7.2$ when Ly$\alpha$ emission is at $\sim$1$\mu$m. Our Texas Spectroscopic Search for Ly$\alpha$ Emission at the End of Reionization includes 18 nights of deep spectroscopic observations using the Keck DEIMOS (optical) and MOSFIRE (NIR) spectrographs. Within this dataset we observe Ly$\alpha$ emission from 183 photometric-redshift selected galaxies at $z =$ 5.5 - 8.3 from the Cosmic Assembly Near infrared Deep Extragalactic Legacy Survey (CANDELS). Our overlapping MOSFIRE observations, over 84 galaxies total, provide the deepest NIR spectroscopic data yet obtained for Ly$\alpha$ from galaxies $z > 7$, with $>16$ hr integration time for four observed galaxies. Here we analyze these four targets, and we report the discovery of a new $z = 7.60$ Ly$\alpha$ detection as well as provide an updated observation of the previously confirmed $z=7.51$ Ly$\alpha$ emission from Finkelstein et al. (2013) with a $\sim$3$\times$ longer exposure time. Our analysis of these Ly$\alpha$ emission line profiles reveal a significant asymmetric shape. The two detected Ly$\alpha$ emission lines from bright sources ($M_{\text{UV}}<-20.25$) could imply that these bright galaxies inhabit ionized bubbles in a partially neutral IGM, although deeper exposures may yet reveal Ly$\alpha$ emission in the fainter sources.

We present the $ugriz$-band Dark Energy Camera (DECam) plus 3.6 and 4.5 $\mu$m IRAC catalogs for the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim24$ deg$^{2}$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 region, with seven bandpasses spanning a wavelength range of 0.35 to 4.5 $\mu$m. SHELA falls within the footprint of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which will provide spectroscopic redshifts for $\sim200{,}000$ Ly$\alpha$ emitters at $1.9<z<3.5$ and also for $\sim200{,}000$ [OII] emitters at $z<0.5$. SHELA's deep, wide-area multiwavelength images combined with HETDEX's spectroscopic information, will facilitate many extragalactic studies, including measuring the evolution of galaxy stellar mass, halo mass, and environment from $1.5<z<3.5$. Here we present $riz$-band selected $ugriz$-band DECam catalogs that reach a $5\sigma$ depth of $\sim24.5$ AB mag (for point sources with an aperture that encloses $70\%$ of the total flux) and cover $17.5$ deg$^{2}$ of the overall SHELA field. We validate our DECam catalog by comparison to the DECam Legacy Survey (DECaLS) DR5 and the Dark Energy Survey (DES) DR1. We perform IRAC forced photometry with The Tractor image modeling code to measure 3.6 and 4.5 $\mu$m fluxes for all objects within our DECam catalog. We demonstrate the utility of our catalog by computing galaxy number counts and estimating photometric redshifts. Our photometric redshifts recover the available $\left\langle z \right\rangle = 0.33 $ SDSS spectroscopic redshifts with a $1\sigma$ scatter in $\Delta z/(1 +z)$ of 0.04.

In anticipation of the upcoming deployment of the James Webb Space Telescope (JWST), we present high-redshift predictions by the well-established Santa Cruz semi-analytic model. We update the models by re-calibrating them after adopting cosmological parameters consistent with recent constraints from Planck. We provide predictions for rest-frame UV luminosity functions for galaxy populations over a wide range of $M_\text{UV}$ from $\sim-6$ to $\sim-24$ between $z = 4 - 10$. In addition, we present the corresponding predictions for observed-frame galaxy number counts in different redshift bins in the full set of NIRCam filters. We provide predictions of the quantitative effect on these observables of varying the physical recipes implemented in the models, such as the molecular gas depletion time (star formation efficiency) scalings or the scalings of outflow rates driven by stars and supernovae with galaxy circular velocity. Based on these results, we discuss what may be learned about the physical processes that shape galaxy formation from JWST observations of galaxy number densities at different intrinsic luminosities. All data tables for the results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

We present the results from the application of a two-dimensional emission line detection method, EMission-line two-Dimensional (EM2D), to the near-infrared G102 grism observations obtained with the Wide-Field Camera 3 (WFC3) as part of the Cycle 22 {\em Hubble Space Telescope} Treasury Program: the Faint Infrared Grism Survey (FIGS). Using the EM2D method, we have assembled a catalog of emission line galaxies (ELGs) with resolved star formation from each of the four FIGS fields. Not only can one better assess the global properties of ELGs, but the EM2D method allows for the analysis and an improved study of the individual emission-line region {\it within} each galaxy. This paper includes a description of the methodology, advantages, and the first results of the EM2D method applied to ELGs in FIGS. The advantage of 2D emission line measurements includes significant improvement of galaxy redshift measurements, approaching the level of accuracy seen in high-spectral-resolution data, but with greater efficiency; and the ability to identify and measure the properties of multiple sites of star-formation and over scales of $\sim$ 1 kpc within individual galaxies out to z $\sim$ 4. The EM2D method also significantly improves the reliability of high-redshift ($z\sim7$) Lyman-$\alpha$ detections. Coupled with the wide field of view and high efficiency of space-based grism observations, EM2D provides a noteworthy improvement on the physical parameters that can be extracted from grism observations.

Lyman-alpha (Ly{\alpha}) photons from ionizing sources and cooling radiation undergo a complex resonant scattering process that generates unique spectral signatures in high-redshift galaxies. We present a detailed Ly{\alpha} radiative transfer study of a cosmological zoom-in simulation from the Feedback In Realistic Environments (FIRE) project. We focus on the time, spatial, and angular properties of the Ly{\alpha} emission over a redshift range of z = 5-7, after escaping the galaxy and being transmitted through the intergalactic medium (IGM). Over this epoch, our target galaxy has an average stellar mass of $M_{\rm star} \approx 5 \times 10^8 {\rm M}_\odot$. We find that many of the interesting features of the Ly{\alpha} line can be understood in terms of the galaxy's star formation history. The time variability, spatial morphology, and anisotropy of Ly{\alpha} properties are consistent with current observations. For example, the rest frame equivalent width has a ${\rm EW}_{{\rm Ly}\alpha,0} > 20 {\rm \AA}$ duty cycle of 62% with a non-negligible number of sightlines with $> 100 {\rm \AA}$, associated with outflowing regions of a starburst with greater coincident UV continuum absorption, as these conditions generate redder, narrower (or single peaked) line profiles. The lowest equivalent widths correspond to cosmological filaments, which have little impact on UV continuum photons but efficiently trap Ly{\alpha} and produce bluer, broader lines with less transmission through the IGM. We also show that in dense self-shielding, low-metallicity filaments and satellites Ly{\alpha} radiation pressure can be dynamically important. Finally, despite a significant reduction in surface brightness with increasing redshift, Ly{\alpha} detections and spectroscopy of high-$z$ galaxies with the upcoming James Webb Space Telescope is feasible.

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

To achieve a fuller understanding of galaxy evolution, SED fitting can be used to recover quantities beyond stellar masses (M$_*$) and star formation rates (SFRs). We use Star Formation Histories (SFHs) reconstructed via the Dense Basis method of Iyer \& Gawiser (2017) for a sample of $17,873$ galaxies at $0.5<z<6$ in the CANDELS GOODS-S field to study the nature and evolution of the SFR-M$_*$ correlation. The reconstructed SFHs represent trajectories in SFR-M$_*$ space, enabling us to study galaxies at epochs earlier than observed by propagating them backwards in time along these trajectories. We study the SFR-M$_*$ correlation at $z=1,2,3,4,5,6$ using both direct fits to galaxies observed at those epochs and SFR-M$_*$ trajectories of galaxies observed at lower redshifts. The SFR-M$_*$ correlations obtained using the two approaches are found to be consistent with each other through a KS test. Validation tests using SFHs from semi-analytic models and cosmological hydrodynamical simulations confirm the sensitivity of the method to changes in the slope, normalization and shape of the SFR-M$_*$ correlation. This technique allows us to further probe the low-mass regime of the correlation at high-z by $\sim 1$ dex and over an effective volume of $\sim 10\times$ larger than possible with just direct fits. We find that the SFR-M$_*$ correlation is consistent with being linear down to M$_*\sim 10^7 M_\odot$ at $z>4$. The evolution of the correlation is well described by $\log SFR= (0.80\pm 0.029 - 0.017\pm 0.010\times t_{univ})\log M_*$ $- (6.487\pm 0.282-0.039\pm 0.008\times t_{univ})$, where $t_{univ}$ is the age of the universe in Gyr.

We present the results of CANDELSz7, an ESO large program aimed at confirming spectroscopically a homogeneous sample of z~6 and z~7 star forming galaxies. The candidates were selected in the GOODS-South, UDS and COSMOS fields using the official CANDELS catalogs based on H160-band detections. Standard color criteria, which were tailored depending on the ancillary multi-wavelength data available for each field, were applied to select more than 160 candidate galaxies at z~6 and z~7. Deep medium resolution FORS2 spectroscopic observations were then conducted with integration times ranging from 12 to 20 hours, to reach a Lyalpha flux limit of approximately 1-3x 10-18 erg/s/cm^2 at 3sigma. For about 40% of the galaxies we could determine a spectroscopic redshift, mainly through the detection of a single emission line that we interpret as Lyalpha emission, or for some of the brightest objects (H160< 25.5) from the presence of faint continuum and sharp drop that we interpret as a Lyman break. In this paper we present the redshifts and main properties of 65 newly confirmed high redshift galaxies. Adding previous proprietary and archival data we assemble a sample of ~260 galaxies that we use to explore the evolution of the Lyalpha fraction in Lyman break galaxies and the change in the shape of the emission line between z~6 and z~7. We also discuss the accuracy of the CANDELS photometric redshifts in this redshift range.

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

Deep, pencil-beam surveys from ALMA at 1.1-1.3mm have uncovered an apparent absence of high-redshift dusty galaxies, with existing redshift distributions peaking around $z\sim1.5-2.5$. This has led to a perceived dearth of dusty systems at $z>4$, and the conclusion, according to some models, that the early Universe was relatively dust-poor. In this paper, we extend the backward evolution galaxy model described by Casey et al. (2018) to the ALMA regime (in depth and area) and determine that the measured number counts and redshift distributions from ALMA deep field surveys are fully consistent with constraints of the infrared luminosity function (IRLF) at $z<2.5$ determined by single-dish submillimeter and millimeter surveys conducted on much larger angular scales ($\sim1-10$deg$^{2}$). We find that measured 1.1-1.3mm number counts are most constraining for the measurement of the faint-end slope of the IRLF at $z<2.5$ instead of the prevalence of dusty galaxies at $z>4$. Recent studies have suggested that UV-selected galaxies at $z>4$ may be particularly dust-poor, but we find their millimeter-wave emission cannot rule out consistency with the Calzetti dust attenuation law even by assuming relatively typical, cold-dust ($T_{\rm dust}\approx30\,$K) SEDs. Our models suggest that the design of ALMA deep fields requires substantial revision to constrain the prevalence of $z>4$ early Universe obscured starbursts. The most promising avenue for detection and characterization of such early dusty galaxies will come from future ALMA 2mm blank field surveys covering a few hundred arcmin$^{2}$ and the combination of existing and future dual-purpose 3mm datasets.

We present a joint analysis of the rest-frame ultraviolet (UV) luminosity functions of continuum-selected star-forming galaxies and galaxies dominated by active galactic nuclei (AGNs) at $z \sim$ 4. These 3,740 $z \sim$ 4 galaxies are selected from broad-band imaging in nine photometric bands over 18 deg$^2$ in the \textit{Spitzer}/HETDEX Exploratory Large Area Survey (SHELA) field. The large area and moderate depth of our survey provide a unique view of the intersection between the bright end of the galaxy UV luminosity function (M$_{AB}<-$22) and the faint end of the AGN UV luminosity function. We do not separate AGN-dominated galaxies from star-formation-dominated galaxies, but rather fit both luminosity functions simultaneously. These functions are best fit with a double power-law (DPL) for both the galaxy and AGN components, where the galaxy bright-end slope has a power-law index of $-3.80\pm0.10$, and the corresponding AGN faint-end slope is $\alpha_{AGN} = -1.49^{+0.30}_{-0.21}$. We cannot rule out a Schechter-like exponential decline for the galaxy UV luminosity function, and in this scenario the AGN luminosity function has a steeper faint-end slope of $-2.08^{+0.18}_{-0.11}$. Comparison of our galaxy luminosity function results with a representative cosmological model of galaxy formation suggests that the molecular gas depletion time must be shorter, implying that star formation is more efficient in bright galaxies at $z=4$ than at the present day. If the galaxy luminosity function does indeed have a power-law shape at the bright end, the implied ionizing emissivity from AGNs is not inconsistent with previous observations. However, if the underlying galaxy distribution is Schechter, it implies a significantly higher ionizing emissivity from AGNs at this epoch.

Though half of cosmic starlight is absorbed by dust and reradiated at long wavelengths (3$\mu$m-3mm), constraints on the infrared through millimeter galaxy luminosity function (the `IRLF') are poor in comparison to the rest-frame ultraviolet and optical galaxy luminosity function, particularly at z>2.5. Here we present a backward evolution model for interpreting number counts, redshift distributions, and cross-band flux density correlations in the infrared and submillimeter sky, from 70$\mu$m-2mm, using a model for the IRLF out to the epoch of reionization. Mock submillimeter maps are generated by injecting sources according to the prescribed IRLF and flux densities drawn from model spectral energy distributions that mirror the distribution of SEDs observed in $0<z<5$ dusty star-forming galaxies (DSFGs). We explore two extreme hypothetical case-studies: a dust-poor early Universe model, where DSFGs contribute negligibly ($<$10%) to the integrated star-formation rate density at $z>4$, and an alternate dust-rich early Universe model, where DSFGs dominate $\sim$90% of $z>4$ star-formation. We find that current submm/mm datasets do not clearly rule out either of these extreme models. We suggest that future surveys at 2mm will be crucial to measuring the IRLF beyond $z\sim4$. The model framework developed in this paper serves as a unique tool for the interpretation of multiwavelength IR/submm extragalactic datasets and will enable more refined constraints on the IRLF than can be made from direct measurements of individual galaxies' integrated dust emission.

We report a new analysis of the Hubble Frontier Fields clusters Abell 2744 and MACS 0416 using wavelet decomposition to remove the cluster light, enabling the detection of highly magnified (>50x) galaxies a factor of 10x fainter in luminosity than previous studies. We find 167 galaxies at z > 6, and with this sample we are able to characterize the UV luminosity function to M_UV = -12.5 at z ~ 6, -14 at z ~ 7 and -15 at z ~ 8. We find a steep faint-end slope (alpha <-2), and with our improved statistics at the faint end we reduce the fractional uncertainty on alpha to <2% at z ~ 6 - 7 and 4% at z ~ 8. We also investigate the systematic uncertainty due to the lens modelling by using every available lens model individually and comparing the results; this systematic fractional uncertainty on alpha is <4% at all redshifts. We now directly observe galaxies in the luminosity regime where simulations predict a change in the faint-end slope of the luminosity function, yet our results provide statistically very strong evidence against any turnover in the luminosity range probed, more consistent with simulations in which stars form in lower-mass halos. Thus we find strong support for the extension of the steep luminosity function to M_UV = -13 at z > 6, consistent with the number of faint galaxies needed to reionize the Universe under standard assumptions.

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

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

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

We improve the accuracy of photometric redshifts by including low-resolution spectral data from the G102 grism on the Hubble Space Telescope, which assists in redshift determination by further constraining the shape of the broadband Spectral Energy Disribution (SED) and identifying spectral features. The photometry used in the redshift fits includes near-IR photometry from FIGS+CANDELS, as well as optical data from ground-based surveys and HST ACS, and mid-IR data from Spitzer. We calculated the redshifts through the comparison of measured photometry with template galaxy models, using the EAZY photometric redshift code. For objects with F105W $< 26.5$ AB mag with a redshift range of $0 < z < 6$, we find a typical error of $\Delta z = 0.03 * (1+z)$ for the purely photometric redshifts; with the addition of FIGS spectra, these become $\Delta z = 0.02 * (1+z)$, an improvement of 50\%. Addition of grism data also reduces the outlier rate from 8\% to 7\% across all fields. With the more-accurate spectrophotometric redshifts (SPZs), we searched the FIGS fields for galaxy overdensities. We identified 24 overdensities across the 4 fields. The strongest overdensity, matching a spectroscopically identified cluster at $z=0.85$, has 28 potential member galaxies, of which 8 have previous spectroscopic confirmation, and features a corresponding X-ray signal. Another corresponding to a cluster at $z=1.84$ has 22 members, 18 of which are spectroscopically confirmed. Additionally, we find 4 overdensities that are detected at an equal or higher significance in at least one metric to the two confirmed clusters.

The distribution of Ly$\alpha$ emission is an presently accessible method for studying the state of the intergalactic medium (IGM) into the reionization era. We carried out deep spectroscopic observations in order to search for Ly$\alpha$ emission from galaxies with photometric redshifts $z$ = 5.5 - 8.3 selected from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). Utilizing data from the Keck/DEIMOS spectrograph, we explore a wavelength coverage of Ly$\alpha$ emission at $z$ ~ 5 - 7 with four nights of spectroscopic observations for 118 galaxies, detecting five emission lines with ~ 5$\sigma$ significance: three in the GOODS-N and two in the GOODS-S field. We constrain the equivalent width (EW) distribution of Ly$\alpha$ emission by comparing the number of detected objects with the expected number constructed from detailed simulations of mock emission lines that account for the observational conditions (e.g., exposure time, wavelength coverage, and sky emission) and galaxy photometric redshift probability distribution functions. The Ly$\alpha$ EW distribution is well described by an exponential form, $\text{dN/dEW}\propto \text{exp(-EW/}W_0)$, characterized by the $e$-folding scale ($W_0$) of ~ 60 - 100$\AA$ at 0.3 < $z$ < 6. By contrast, our measure of the Ly$\alpha$ EW distribution at 6.0 < $z$ < 7.0 rejects a Ly$\alpha$ EW distribution with $W_0$ > 36.4$\AA$ (125.3$\AA$) at 1$\sigma$ (2$\sigma$) significance. This provides additional evidence that the EW distribution of Ly$\alpha$ declines at $z$ > 6, suggesting an increasing fraction of neutral hydrogen in the IGM at that epoch.

Using cosmological volume simulations and a custom built sub-grid model for Pop~III star formation, we examine the baseline dust extinction in the first galaxies due to Pop~III metal enrichment in the first billion years of cosmic history. We find that while the most enriched, high-density lines of sight in primordial galaxies can experience a measurable amount of extinction from Pop~III dust ($E(B-V)_{\rm max}=0.07,\ A_{\rm V,max}\approx0.28$), the average extinction is very low with $\left< E(B-V) \right> \lesssim 10^{-3}$. We derive a power-law relationship between dark matter halo mass and extinction of $E(B-V)\propto M_{\rm halo}^{0.80}$. Performing a Monte Carlo parameter study, we establish the baseline reddening of the UV spectra of dwarf galaxies at high redshift due to Pop~III enrichment only. With this method, we find $\left<\beta_{\rm UV}\right>-2.51\pm0.07$, which is both nearly halo mass and redshift independent.

We utilize the hydrodynamic and N-body code {\small GIZMO} coupled with our newly developed sub-grid Population~III (Pop~III) Legacy model, designed specifically for cosmological volume simulations, to study the baseline metal enrichment from Pop~III star formation at $z>7$. In this idealized numerical experiment, we only consider Pop~III star formation. We find that our model Pop~III star formation rate density (SFRD), which peaks at $\sim 10^{-3}\ {\rm M_\odot yr^{-1} Mpc^{-1}}$ near $z\sim10$, agrees well with previous numerical studies and is consistent with the observed estimates for Pop~II SFRDs. The mean Pop~III metallicity rises smoothly from $z=25-7$, but does not reach the critical metallicity value, $Z_{\rm crit}=10^{-4}\ Z_\odot$, required for the Pop~III to Pop~II transition in star formation mode until $z\simeq7$. This suggests that, while individual halos can suppress in-situ Pop~III star formation, the external enrichment is insufficient to globally terminate Pop~III star formation. The maximum enrichment from Pop~III star formation in star forming dark matter halos is $Z\sim10^{-2}\ Z_\odot$, whereas the minimum found in externally enriched haloes is $Z\gtrsim10^{-7}\ Z_\odot$. Finally, mock observations of our simulated IGM enriched with Pop~III metals produce equivalent widths similar to observations of an extremely metal poor damped Lyman alpha (DLA) system at $z=7.04$, which is thought to be enriched by Pop~III star formation only.

The rate of major galaxy-galaxy merging is theoretically predicted to steadily increase with redshift during the peak epoch of massive galaxy development ($1{\leq}z{\leq}3$). We use close-pair statistics to objectively study the incidence of massive galaxies (stellar $M_{1}{\geq}2{\times}10^{10}M_{\odot}$) hosting major companions ($1{\leq}M_{1}/M_{2}{\leq}4$; i.e., $<$4:1) at six epochs spanning $0{<}z{<}3$. We select companions from a nearly complete, mass-limited ($\geq5{\times}10^{9}M_{\odot}$) sample of 23,696 galaxies in the five CANDELS fields and the SDSS. Using $5-50$ kpc projected separation and close redshift proximity criteria, we find that the major companion fraction $f_{\mathrm{mc}}(z)$ based on stellar mass-ratio (MR) selection increases from 6% ($z{\sim}0$) to 16% ($z{\sim}0.8$), then turns over at $z{\sim}1$ and decreases to 7% ($z{\sim}3$). Instead, if we use a major F160W flux ratio (FR) selection, we find that $f_{\mathrm{mc}}(z)$ increases steadily until $z=3$ owing to increasing contamination from minor (MR$>$4:1) companions at $z>1$. We show that these evolutionary trends are statistically robust to changes in companion proximity. We find disagreements between published results are resolved when selection criteria are closely matched. If we compute merger rates using constant fraction-to-rate conversion factors ($C_{\mathrm{merg,pair}}{=}0.6$ and $T_{\mathrm{obs,pair}}{=}0.65\mathrm{Gyr}$), we find that MR rates disagree with theoretical predictions at $z{>}1.5$. Instead, if we use an evolving $T_{\mathrm{obs,pair}}(z){\propto}(1+z)^{-2}$ from Snyder et al., our MR-based rates agree with theory at $0{<}z{<}3$. Our analysis underscores the need for detailed calibration of $C_{\mathrm{merg,pair}}$ and $T_{\mathrm{obs,pair}}$ as a function of redshift, mass and companion selection criteria to better constrain the empirical major merger history.

We present the results of an unbiased search for Ly{\alpha} emission from continuum-selected 6 < z < 8 galaxies. Our dataset consists of 160 orbits of G102 slitless grism spectroscopy obtained with the Hubble Space Telescope (HST) Wide Field Camera 3 as part of the Faint Infrared Grism Survey (FIGS; PI: Malhotra), which obtains deep slitless spectra of all sources in four fields, and was designed to minimize contamination in observations of previously-identified high-redshift galaxy candidates. The FIGS data can potentially spectroscopically confirm the redshifts of galaxies, and as Ly{\alpha} emission is resonantly scattered by neutral gas, FIGS can also constrain the ionization state of the intergalactic medium (IGM) during the epoch of reionization. These data have sufficient depth to detect Ly{\alpha} emission in this epoch, as Tilvi et al. (2016) have published the FIGS detection of previously known (Finkelstein et al. 2013) Ly{\alpha} emission at z = 7.51. The FIGS data use five separate roll-angles of HST to mitigate the contamination by nearby galaxies. We created a method that accounts for and removes the contamination from surrounding galaxies, and also removes any dispersed continuum light from each individual spectrum (Pirzkal et al. 2017). We searched for significant (> 4{\sigma}) emission lines using two different automated detection methods, free of any visual inspection biases. Applying these methods on photometrically-selected high-redshift candidates between 6 < z < 8 we find two emission lines, one previously published by Tilvi et al. (2016), and a new line at 1.028{\mu}m. We identify this lines as Ly{\alpha} at z = 7.452 +/- 0.003. This newly spectroscopically confirmed galaxy has the highest Ly{\alpha} rest-frame equivalent width (EW Ly{\alpha}) yet published at z > 7 (140.3 +/- 19.0{\AA}).

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

The fraction of Lyman-$\alpha$ emitters among the galaxy population has been found to increase from $z\sim0$ to $z\sim6$ and drop dramatically at $z>6$. This drop has been interpreted as an effect of an increasingly neutral intergalactic medium with increasing redshift, while a LyC escape fraction evolving with redshift. We report the result of a large VLT/FORS2 program aiming to confirm spectroscopically a large galaxy sample at $z\geq6$ that has been selected in several independent fields through the Lyman Break technique. Combining those data with archival data, we create a large and homogeneous sample of $z\sim6$ galaxies ($N=127$), complete in terms of Ly$\alpha$ detection at $>95\%$ for EW(Ly$\alpha)\geq25\AA$. We use this sample to derive a new measurement of the LAE fraction at $z\sim6$ and derive the physical properties of these galaxies through spectral energy distribution fitting. We find a median LAE fraction at $z\sim6$ lower than in previous studies. The main difference between LAEs and non-LAEs is that the latter are significantly dustier. Using predictions of our SED fitting code accounting for nebular emission, we find an effective Ly$\alpha$ escape fraction $f^{eff}_{esc}(Ly\alpha)=0.23^{+0.36}_{-0.17}$ remarkably consistent with the value derived by comparing UV luminosity function with Ly$\alpha$ luminosity function. We conclude that the drop in the LAE fraction from $z\sim6$ to $z>6$ is less dramatic than previously found and the effect of an increasing IGM neutral fraction is possibly observed at $5<z<6$. Based on our derived $f^{eff}_{esc}(Ly\alpha)$, we find that the IGM has a relatively small impact on Ly$\alpha$ photon visibility at $z\sim6$, with a lower limit for the IGM transmission to \lya\ photons, $T_{IGM}\gtrsim0.20$, likely due to the presence of outflows. [abdridged]

We present a flux-limited sample of $z\sim0.3$ Ly$\alpha$ emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism spectroscopic data. The published GALEX $z\sim0.3$ LAE sample is pre-selected from continuum-bright objects and thus is biased against high equivalent width (EW) LAEs. We remove this continuum pre-selection and compute the EW distribution and the luminosity function of the Ly$\alpha$ emission line directly from our sample. We examine the evolution of these quantities from $z\sim0.3$ to $2.2$ and find that the EW distribution shows little evidence for evolution over this redshift range. As shown by previous studies, the Ly$\alpha$ luminosity density from star-forming galaxies declines rapidly with declining redshift. However, we find that the decline in Ly$\alpha$ luminosity density from $z=2.2$ to $z=0.3$ may simply mirror the decline seen in the H$\alpha$ luminosity density from $z=2.2$ to $z=0.4$, implying little change in the volumetric Ly$\alpha$ escape fraction. Finally, we show that the observed Ly$\alpha$ luminosity density from AGNs is comparable to the observed Ly$\alpha$ luminosity density from star-forming galaxies at $z=0.3$. We suggest that this significant contribution from AGNs to the total observed Ly$\alpha$ luminosity density persists out to $z\sim2.2$.

Narrowband imaging is a highly successful approach for finding large numbers of high redshift Lya emitting galaxies (LAEs) up to z~6.6. However, at z>~7 there are as yet only 3 narrowband selected LAEs with spectroscopic confirmations (two at z~6.9-7.0, one at z~7.3), which hinders extensive studies on cosmic reionization and galaxy evolution at this key epoch. We have selected 23 candidate z~6.9 LAEs in COSMOS field with the large area narrowband survey LAGER (Lyman-Alpha Galaxies at the End of Reionization). In this work we present spectroscopic followup observations of 12 candidates using IMACS on Magellan. For 9 of these, the observations are sufficiently deep to detect the expected lines. Lya emission lines are identified in six sources (yielding a success rate of 2/3), including 3 luminous LAEs with Lya luminosities of L(Lya) ~ 10^{43.5} erg/s, the highest among known spectroscopically confirmed galaxies at >~7.0. This triples the sample size of spectroscopically confirmed narrowband selected LAEs at z>~7, and confirms the bright end bump in the Lya luminosity function we previously derived based on the photometric sample, supporting a patchy reionization scenario. Two luminous LAEs appear physically linked with projected distance of 1.1 pMpc and velocity difference of ~ 170 km/s. They likely sit in a common ionized bubble produced by themselves or with close neighbors, which reduces the IGM attenuation of Lya. A tentative narrow NV${\lambda}$1240 line is seen in one source, hinting at activity of a central massive black hole with metal rich line emitting gas.

The Faint Infrared Grism Survey (FIGS) is a deep Hubble Space Telescope (HST) WFC3/IR (Wide Field Camera 3 Infrared) slitless spectroscopic survey of four deep fields. Two fields are located in the Great Observatories Origins Deep Survey-North (GOODS-N) area and two fields are located in the Great Observatories Origins Deep Survey-South (GOODS-S) area. One of the southern fields selected is the Hubble Ultra Deep Field. Each of these four fields were observed using the WFC3/G102 grism (0.8$\mu m$-1.15$\mu m$ continuous coverage) with a total exposure time of 40 orbits (~ 100 kilo-seconds) per field. This reaches a 3 sigma continuum depth of ~26 AB magnitudes and probes emission lines to $\approx 10^{-17}\ erg\ s^{-1} \ cm^{-2}$. This paper details the four FIGS fields and the overall observational strategy of the project. A detailed description of the Simulation Based Extraction (SBE) method used to extract and combine over 10000 spectra of over 2000 distinct sources brighter than m_F105W=26.5 mag is provided. High fidelity simulations of the observations is shown to significantly improve the background subtraction process, the spectral contamination estimates, and the final flux calibration. This allows for the combination of multiple spectra to produce a final high quality, deep, 1D-spectra for each object in the survey.

We analyze a sample of $z$-dropout galaxies in the CANDELS GOODS South and UDS fields that have been targeted by a dedicated spectroscopic campaign aimed at detecting their Ly$\alpha$ line. Deep IRAC observations at 3.6 and 4.5 $\mu$m are used to determine the strength of optical emission lines affecting these bands at z$\sim$6.5-6.9 in order to i) investigate possible physical differences between Ly$\alpha$ emitting and non-emitting sources; ii) constrain the escape fraction of ionizing photons; iii) provide an estimate of the specific star-formation rate at high redshifts. We find evidence of strong [OIII]+H$\beta$ emission in the average (stacked) SEDs of galaxies both with and without Ly$\alpha$ emission. The blue IRAC [3.6]-[4.5] color of the stack with detected Ly$\alpha$ line can be converted into a rest-frame equivalent width EW([OIII]+H$\beta$)=1500$^{+530}_{-440}\AA$ assuming a flat intrinsic stellar continuum. This strong optical line emission enables a first estimate of f$_{esc}\lesssim$20% on the escape fraction of ionizing photons from Ly$\alpha$ detected objects. The objects with no Ly$\alpha$ line show less extreme EW([OIII]+H$\beta$)=520$^{+170}_{-150}\AA$ suggesting different physical conditions of the HII regions with respect to Ly$\alpha$-emitting ones, or a larger f$_{esc}$. The latter case is consistent with a combined evolution of f$_{esc}$ and the neutral hydrogen fraction as an explanation of the lack of bright Ly$\alpha$ emission at z$>$6. A lower limit on the specific star formation rate, SSFR$>$9.1$Gyr^{-1}$ for $M_{star}=2 \times 10^9 M_{\odot}$ galaxies at these redshifts can be derived from the spectroscopically confirmed sample.

We present the rest-1500\AA\ UV luminosity functions (LF) for star-forming galaxies during the cosmic \textit{high noon} -- the peak of cosmic star formation rate at $1.5<z<3$. We use deep NUV imaging data obtained as part of the \textit{Hubble} Ultra-Violet Ultra Deep Field (UVUDF) program, along with existing deep optical and NIR coverage on the HUDF. We select F225W, F275W and F336W dropout samples using the Lyman break technique, along with samples in the corresponding redshift ranges selected using photometric redshifts and measure the rest-frame UV LF at $z\sim1.7,2.2,3.0$ respectively, using the modified maximum likelihood estimator. We perform simulations to quantify the survey and sample incompleteness for the UVUDF samples to correct the effective volume calculations for the LF. We select galaxies down to $M_{UV}=-15.9,-16.3,-16.8$ and fit a faint-end slope of $\alpha=-1.20^{+0.10}_{-0.13}, -1.32^{+0.10}_{-0.14}, -1.39^{+0.08}_{-0.12}$ at $1.4<z<1.9$, $1.8<z<2.6$, and $2.4<z<3.6$, respectively. We compare the star formation properties of $z\sim2$ galaxies from these UV observations with results from H\alpha\ and UV$+$IR observations. We find a lack of high SFR sources in the UV LF compared to the H\alpha\ and UV$+$IR, likely due to dusty SFGs not being properly accounted for by the generic $IRX-\beta$ relation used to correct for dust. We compute a volume-averaged UV-to-H\alpha\ ratio by \textit{abundance matching} the rest-frame UV LF and H\alpha\ LF. We find an increasing UV-to-H\alpha\ ratio towards low mass galaxies ($M_\star \lesssim 5\times10^9$ M$_\odot$). We conclude that this could be due to a larger contribution from starbursting galaxies compared to the high-mass end.

We present a multi-wavelength photometric catalog in the COSMOS field as part of the observations by the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The catalog is based on Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) and Advanced Camera for Surveys (ACS) observations of the COSMOS field (centered at RA: $10^h00^m28^s$, Dec:$+02^{\circ}12^{\prime}21^{\prime\prime}$). The final catalog has 38671 sources with photometric data in forty two bands from UV to the infrared ($\rm \sim 0.3-8\,\mu m$). This includes broad-band photometry from the HST, CFHT, Subaru, VISTA and Spitzer Space Telescope in the visible, near infrared and infrared bands along with intermediate and narrow-band photometry from Subaru and medium band data from Mayall NEWFIRM. Source detection was conducted in the WFC3 F160W band (at $\rm 1.6\,\mu m$) and photometry is generated using the Template FITting algorithm. We further present a catalog of the physical properties of sources as identified in the HST F160W band and measured from the multi-band photometry by fitting the observed spectral energy distributions of sources against templates.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We perform the first spatially-resolved stellar population study of galaxies in the early universe (z = 3.5 - 6.5), utilizing the Hubble Space Telescope Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) imaging dataset over the GOODS-S field. We select a sample of 418 bright and extended galaxies at z = 3.5 - 6.5 from a parent sample of ~ 8000 photometric-redshift selected galaxies from Finkelstein et al. (2015). We first examine galaxies at 3.5< z < 4.0 using additional deep K-band survey data from the HAWK-I UDS and GOODS Survey (HUGS) which covers the 4000A break at these redshifts. We measure the stellar mass, star formation rate, and dust extinction for galaxy inner and outer regions via spatially-resolved spectral energy distribution fitting based on a Markov Chain Monte Carlo algorithm. By comparing specific star formation rates (sSFRs) between inner and outer parts of the galaxies we find that the majority of galaxies with the high central mass densities show evidence for a preferentially lower sSFR in their centers than in their outer regions, indicative of reduced sSFRs in their central regions. We also study galaxies at z ~ 5 and 6 (here limited to high spatial resolution in the rest-frame ultraviolet only), finding that they show sSFRs which are generally independent of radial distance from the center of the galaxies. This indicates that stars are formed uniformly at all radii in massive galaxies at z ~ 5 - 6, contrary to massive galaxies at z < 4.

Galaxies at high redshifts provide a valuable tool to study cosmic dawn, and therefore it is crucial to reliably identify these galaxies. Here, we present an unambiguous and first simultaneous detection of both the Lyman-alpha emission and the Lyman break from a z = 7.512+/- 0.004 galaxy, observed in the Faint Infrared Grism Survey (FIGS). These spectra, taken with G102 grism on Hubble Space Telescope (HST), show a significant emission line detection (6 sigma) in multiple observational position angles (PA), with total integrated Ly{\alpha} line flux of 1.06+/- 0.12 e10-17erg s-1cm-2. The line flux is nearly a factor of four higher than the previous MOSFIRE spectroscopic observations of faint Ly{\alpha} emission at {\lambda} = 1.0347{\mu}m, yielding z = 7.5078+/- 0.0004. This is consistent with other recent observations implying that ground-based near-infrared spectroscopy underestimates total emission line fluxes, and if confirmed, can have strong implications for reionization studies that are based on ground-based Lyman-{\alpha} measurements. A 4-{\sigma} detection of the NV line in one PA also suggests a weak Active Galactic Nucleus (AGN), potentially making this source the highest-redshift AGN yet found. Thus, this observation from the Hubble Space Telescope clearly demonstrates the sensitivity of the FIGS survey, and the capability of grism spectroscopy to study the epoch of reionization.

The gas accretion and star-formation histories of galaxies like the Milky Way remain an outstanding problem in astrophysics. Observations show that 8 billion years ago, the progenitors to Milky Way-mass galaxies were forming stars 30 times faster than today and predicted to be rich in molecular gas, in contrast with low present-day gas fractions ($<$10%). Here we show detections of molecular gas from the CO(J=3-2) emission (rest-frame 345.8 GHz) in galaxies at redshifts z=1.2-1.3, selected to have the stellar mass and star-formation rate of the progenitors of today's Milky Way-mass galaxies. The CO emission reveals large molecular gas masses, comparable to or exceeding the galaxy stellar masses, and implying most of the baryons are in cold gas, not stars. The galaxies' total luminosities from star formation and CO luminosities yield long gas-consumption timescales. Compared to local spiral galaxies, the star-formation efficiency, estimated from the ratio of total IR luminosity to CO emission,} has remained nearly constant since redshift z=1.2, despite the order of magnitude decrease in gas fraction, consistent with results for other galaxies at this epoch. Therefore the physical processes that determine the rate at which gas cools to form stars in distant galaxies appear to be similar to that in local galaxies.

We present quantified visual morphologies of approximately 48,000 galaxies observed in three Hubble Space Telescope legacy fields by the Cosmic And Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and classified by participants in the Galaxy Zoo project. 90% of galaxies have z < 3 and are observed in rest-frame optical wavelengths by CANDELS. Each galaxy received an average of 40 independent classifications, which we combine into detailed morphological information on galaxy features such as clumpiness, bar instabilities, spiral structure, and merger and tidal signatures. We apply a consensus-based classifier weighting method that preserves classifier independence while effectively down-weighting significantly outlying classifications. After analysing the effect of varying image depth on reported classifications, we also provide depth-corrected classifications which both preserve the information in the deepest observations and also enable the use of classifications at comparable depths across the full survey. Comparing the Galaxy Zoo classifications to previous classifications of the same galaxies shows very good agreement; for some applications the high number of independent classifications provided by Galaxy Zoo provides an advantage in selecting galaxies with a particular morphological profile, while in others the combination of Galaxy Zoo with other classifications is a more promising approach than using any one method alone. We combine the Galaxy Zoo classifications of "smooth" galaxies with parametric morphologies to select a sample of featureless disks at 1 < z < 3, which may represent a dynamically warmer progenitor population to the settled disk galaxies seen at later epochs.

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection.

Although the universe at redshifts greater than six represents only the first one billion years (<10%) of cosmic time, the dense nature of the early universe led to vigorous galaxy formation and evolution activity which we are only now starting to piece together. Technological improvements have, over only the past decade, allowed large samples of galaxies at such high redshifts to be collected, providing a glimpse into the epoch of formation of the first stars and galaxies. A wide variety of observational techniques have led to the discovery of thousands of galaxy candidates at z > 6, with spectroscopically confirmed galaxies out to nearly z = 9. Using these large samples, we have begun to gain a physical insight into the processes inherent in galaxy evolution at early times. In this review, I will discuss i) the selection techniques for finding distant galaxies, including a summary of previous and ongoing ground and space-based searches, and spectroscopic followup efforts, ii) insights into galaxy evolution gleaned from measures such as the rest-frame ultraviolet luminosity function, the stellar mass function, and galaxy star-formation rates, and iii) the effect of galaxies on their surrounding environment, including the chemical enrichment of the universe, and the reionization of the intergalactic medium. Finally, I conclude with prospects for future observational study of the distant universe, using a bevy of new state-of-the-art facilities coming online over the next decade and beyond.

We present a narrowband survey with three adjacent filters for z=2.8--2.9 Lyman Alpha Emitter (LAE) galaxies in the Extended Chandra Deep Field South (ECDFS), along with spectroscopic followup. With a complete sample of 96 LAEs in the narrowband NB466, we confirm a large-scale structure at z~ 2.8. Compared to the blank field in NB470 and NB475, the LAE density excess in the NB466 field is ~6.0+/-0.8 times the standard deviation expected at z~2.8, assuming a linear bias of 2. The overdense large scale structure in NB466 can be decomposed into 4 protoclusters, whose overdensities are 4.6 - 6.6. These 4 protoclusters are expected to evolve into a Coma-like cluster at z~ 0. In the meanwhile, we investigate the average star-formation rates derived from Ly{\alpha}, rest-frame UV and X-ray, the Ly{\alpha} luminosity functions, the Ly{\alpha} photon densities and their dependence on the environment. We find that the Ly{\alpha} photon density in the overdense field (NB466) is ~50\% higher than that in the blank field (NB470+NB475). The 3 brightest LAEs, including a quasar at z=2.81, are all detected in X-ray and in NB466. These three LAE-AGNs contribute an extra 20--30\% Ly{\alpha} photon density. Furthermore, we find that LAEs in overdense regions are younger and less dusty. We conclude that the structure we found is a significant and rare density peak, and narrowband imaging is an efficient method to detect and study such structures in the high-z universe.

We report the results from some of the deepest Keck/Multi-Object Spectrometer For Infra-Red Exploration data yet obtained for candidate $z \gtrsim 7$ galaxies. Our data show one significant line detection with 6.5$\sigma$ significance in our combined 10 hr of integration which is independently detected on more than one night, thus ruling out the possibility that the detection is spurious. The asymmetric line profile and non-detection in the optical bands strongly imply that the detected line is Ly$\alpha$ emission from a galaxy at $z$(Ly$\alpha)=7.6637 \pm 0.0011$, making it the fourth spectroscopically confirmed galaxy via Ly$\alpha$ at $z>7.5$. This galaxy is bright in the rest-frame ultraviolet (UV; $M_{\rm UV} \sim -21.2$) with a moderately blue UV slope ($\beta=-2.2^{+0.3}_{-0.2}$), and exhibits a rest-frame Ly$\alpha$ equivalent width of EW(Ly$\alpha$) $\sim 15.6^{+5.9}_{-3.6}$ \AA. The non-detection of the 11 other $z \sim$ 7-8 galaxies in our long 10 hr integration, reaching a median 5$\sigma$ sensitivity of 28 \AA\ in the rest-frame EW(Ly$\alpha$), implies a 1.3$\sigma$ deviation from the null hypothesis of a non-evolving distribution in the rest-frame EW(Ly$\alpha$) between $3<z<6$ and $z=$ 7-8. Our results are consistent with previous studies finding a decline in Ly$\alpha$ emission at $z>6.5$, which may signal the evolving neutral fraction in the intergalactic medium at the end of the reionization epoch, although our weak evidence suggests the need for a larger statistical sample to allow for a more robust conclusion.

We present galaxy stellar mass functions (GSMFs) at $z=$ 4-8 from a rest-frame ultraviolet (UV) selected sample of $\sim$4500 galaxies, found via photometric redshifts over an area of $\sim$280 arcmin$^2$ in the CANDELS/GOODS fields and the Hubble Ultra Deep Field. The deepest Spitzer/IRAC data yet-to-date and the relatively large volume allow us to place a better constraint at both the low- and high-mass ends of the GSMFs compared to previous space-based studies from pre-CANDELS observations. Supplemented by a stacking analysis, we find a linear correlation between the rest-frame UV absolute magnitude at 1500 \AA\ ($M_{\rm UV}$) and logarithmic stellar mass ($\log M_*$) that holds for galaxies with $\log(M_*/M_{\odot}) \lesssim 10$. We use simulations to validate our method of measuring the slope of the $\log M_*$-$M_{\rm UV}$ relation, finding that the bias is minimized with a hybrid technique combining photometry of individual bright galaxies with stacked photometry for faint galaxies. The resultant measured slopes do not significantly evolve over $z=$ 4-8, while the normalization of the trend exhibits a weak evolution toward lower masses at higher redshift. We combine the $\log M_*$-$M_{\rm UV}$ distribution with observed rest-frame UV luminosity functions at each redshift to derive the GSMFs, finding that the low-mass-end slope becomes steeper with increasing redshift from $\alpha=-1.55^{+0.08}_{-0.07}$ at $z=4$ to $\alpha=-2.25^{+0.72}_{-0.35}$ at $z=8$. The inferred stellar mass density, when integrated over $M_*=10^8$-$10^{13} M_{\odot}$, increases by a factor of $10^{+30}_{-2}$ between $z=7$ and $z=4$ and is in good agreement with the time integral of the cosmic star formation rate density.

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant Universe. Here, we model the spectral energy distributions (SEDs) of galaxies at z = 1.5--3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, $L_{\text{TIR}}/L_{\text{UV}}$) and UV slope ($\beta$). We generalize the shape of the dust law with an empirical model, $A_{\lambda,\delta}=E(B-V)\ k_\lambda\ (\lambda/\lambda_V)^\delta$ where $k_\lambda$ is the dust law of Calzetti et al. (2000), and show that there exists a correlation between the color excess ${E(B-V)}$ and tilt $\delta$ with ${\delta=(0.62\pm0.05)\log(E(B-V))}$+ ${(0.26~\pm~0.02)}$. Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star-formation rate, or $\beta$. The change in the dust law with color excess is consistent with a model where attenuation is caused by by scattering, a mixed star-dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher ($z>3$) redshifts.

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events.

We present a Bayesian approach to the redshift classification of emission-line galaxies when only a single emission line is detected spectroscopically. We consider the case of surveys for high-redshift Lyman-alpha-emitting galaxies (LAEs), which have traditionally been classified via an inferred rest-frame equivalent width (EW) greater than 20 angstrom. Our Bayesian method relies on known prior probabilities in measured emission-line luminosity functions and equivalent width distributions for the galaxy populations, and returns the probability that an object in question is an LAE given the characteristics observed. This approach will be directly relevant for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which seeks to classify ~10^6 emission-line galaxies into LAEs and low-redshift [O II] emitters. For a simulated HETDEX catalog with realistic measurement noise, our Bayesian method recovers 86% of LAEs missed by the traditional EW > 20 angstrom cutoff over 2 < z < 3, outperforming the EW cut in both contamination and incompleteness. This is due to the method's ability to trade off between the two types of binary classification error by adjusting the stringency of the probability requirement for classifying an observed object as an LAE. In our simulations of HETDEX, this method reduces the uncertainty in cosmological distance measurements by 14% with respect to the EW cut, equivalent to recovering 29% more cosmological information. Rather than using binary object labels, this method enables the use of classification probabilities in large-scale structure analyses. It can be applied to narrowband emission-line surveys as well as upcoming large spectroscopic surveys including Euclid and WFIRST.

In this work, we show how the stellar mass (M) of galaxies affects the 3<z<4.6 Ly-alpha equivalent width (EW) distribution. To this end, we design a sample of 629 galaxies in the M range 7.6 < logM/Msun < 10.6 from the 3D-HST/CANDELS survey. We perform spectroscopic observations of this sample using the Michigan/Magellan Fiber System, allowing us to measure Ly-alpha fluxes and use 3D-HST/CANDELS ancillary data. In order to study the Ly-alpha EW distribution dependence on M, we split the whole sample in three stellar mass bins. We find that, in all bins, the distribution is best represented by an exponential profile of the form dN(M)/dEW= A(M)exp(-EW/W0(M))/W0(M). Through a Bayesian analysis, we confirm that lower M galaxies have higher Ly-alpha EWs. We also find that the fraction A of galaxies featuring emission and the e-folding scale W0 of the distribution anti- correlate with M, recovering expressions of the forms A(M)= -0.26(.13) logM/Msun+3.01(1.2) and W0(M)= -15.6(3.5) logM/Msun +166(34). These results are crucial for proper interpretation of Ly-alpha emission trends reported in the literature that may be affected by strong M selection biases.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

Hot Dust-Obscured Galaxies (Hot DOGs) are a population of hyper-luminous infrared galaxies identified by the WISE mission from their very red mid-IR colors, and characterized by hot dust temperatures ($T>60~\rm K$). Several studies have shown clear evidence that the IR emission in these objects is powered by a highly dust-obscured AGN that shows close to Compton-thick absorption at X-ray wavelengths. Thanks to the high AGN obscuration, the host galaxy is easily observable, and has UV/optical colors usually consistent with those of a normal galaxy. Here we discuss a sub-population of 8 Hot DOGs that show enhanced rest-frame UV/optical emission. We discuss three scenarios that might explain the excess UV emission: (i) unobscured light leaked from the AGN by reflection over the dust or by partial coverage of the accretion disk; (ii) a second unobscured AGN in the system; or (iii) a luminous young starburst. X-ray observations can help discriminate between these scenarios. We study in detail the blue excess Hot DOG WISE J020446.13-050640.8, which was serendipitously observed by Chandra/ACIS-I for 174.5 ks. The X-ray spectrum is consistent with a single, hyper-luminous, highly absorbed AGN, and is strongly inconsistent with the presence of a secondary unobscured AGN. Based on this, we argue that the excess blue emission in this object is most likely either due to reflection or a co-eval starburst. We favor the reflection scenario as the unobscured star-formation rate needed to power the UV/optical emission would be $\gtrsim 1000~\rm M_{\odot}~\rm yr^{-1}$. Deep polarimetry observations could confirm the reflection hypothesis.

Using the cosmological smoothed particle hydrodynamical code GADGET-3 we make a realistic assessment of the technique of using constant cumulative number density as a tracer of galaxy evolution at high redshift. We find that over a redshift range of $3\leq z \leq7$ one can on average track the growth of the stellar mass of a population of galaxies selected from the same cumulative number density bin to within $\sim 0.20$ dex. Over the stellar mass range we probe ($10^{10.39}\leq M_s/M_\odot \leq 10^{10.75}$ at $z =$ 3 and $10^{8.48}\leq M_s/M_\odot \leq 10^{9.55}$ at $z =$ 7) one can reduce this bias by selecting galaxies based on an evolving cumulative number density. We find the cumulative number density evolution exhibits a trend towards higher values which can be quantified by simple linear formulations going as $-0.10\Delta z$ for descendants and $0.12\Delta z$ for progenitors. Utilizing such an evolving cumulative number density increases the accuracy of descendant/progenitor tracking by a factor of $\sim2$. This result is in excellent agreement, within $0.10$ dex, with abundance matching results over the same redshift range. However, we find that our more realistic cosmological hydrodynamic simulations produce a much larger scatter in descendant/progenitor stellar masses than previous studies, particularly when tracking progenitors. This large scatter makes the application of either the constant cumulative number density or evolving cumulative number density technique limited to average stellar masses of populations only, as the diverse mass assembly histories caused by stochastic physical processes such as gas accretion, mergers, and star formation of individual galaxies will lead to a larger scatter in other physical properties such as metallicity and star-formation rate.

The upcoming launch of the James Webb Space Telescope (JWST) in less than three years is certain to bring a revolution in our understanding of many area of astrophysics, with one of the key goals being galaxy evolution. As the first proposals will be due in a little over two years, the time is ripe to take a holistic look at the science goals which the community would wish to accomplish with this observatory. Contrary to our experiences with the Hubble Space Telescope, which has now operated successfully for over two decades due to several timely servicing missions, the lifetime of JWST is finite and relatively short, with a lifetime requirement of five years, and a ten-year goal. Following the discussion session at the "Exploring the Universe with JWST" conference at ESA-ESTEC in October 2015, we highlight in this document the (non-local) extragalactic science goals for JWST. We describe how a concerted community effort could best address these, ensuring that the desired survey can be completed during the JWST mission.

We compare the physical and morphological properties of z ~ 2 Lyman-alpha emitting galaxies (LAEs) identified in the HETDEX Pilot Survey and narrow band studies with those of z ~ 2 optical emission line galaxies (oELGs) identified via HST WFC3 infrared grism spectroscopy. Both sets of galaxies extend over the same range in stellar mass (7.5 < logM < 10.5), size (0.5 < R < 3.0 kpc), and star-formation rate (~1 < SFR < 100). Remarkably, a comparison of the most commonly used physical and morphological parameters -- stellar mass, half-light radius, UV slope, star formation rate, ellipticity, nearest neighbor distance, star formation surface density, specific star formation rate, [O III] luminosity, and [O III] equivalent width -- reveals no statistically significant differences between the populations. This suggests that the processes and conditions which regulate the escape of Ly-alpha from a z ~ 2 star-forming galaxy do not depend on these quantities. In particular, the lack of dependence on the UV slope suggests that Ly-alpha emission is not being significantly modulated by diffuse dust in the interstellar medium. We develop a simple model of Ly-alpha emission that connects LAEs to all high-redshift star forming galaxies where the escape of Ly-alpha depends on the sightline through the galaxy. Using this model, we find that mean solid angle for Ly-alpha escape is 2.4+/-0.8 steradians; this value is consistent with those calculated from other studies.

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

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

We present the results from a stellar population modeling analysis of a sample of 162 z=4.5, and 14 z=5.7 Lyman alpha emitting galaxies (LAEs) in the Bootes field, using deep Spitzer/IRAC data at 3.6 and 4.5 um from the Spitzer Lyman Alpha Survey, along with Hubble Space Telescope NICMOS and WFC3 imaging at 1.1 and 1.6 um for a subset of the LAEs. This represents one of the largest samples of high-redshift LAEs imaged with Spitzer IRAC. We find that 30/162 (19%) of the z=4.5 LAEs and 9/14 (64%) of the z=5.7 LAEs are detected at >3-sigma in at least one IRAC band. Individual z=4.5 IRAC-detected LAEs have a large range of stellar mass, from 5x10^8 to 10^11 Msol. One-third of the IRAC-detected LAEs have older stellar population ages of 100 Myr - 1 Gyr, while the remainder have ages < 100 Myr. A stacking analysis of IRAC-undetected LAEs shows this population to be primarily low mass (8 -- 20 x 10^8 Msol) and young (64 - 570 Myr). We find a correlation between stellar mass and the dust-corrected ultraviolet-based star-formation rate (SFR) similar to that at lower redshifts, in that higher mass galaxies exhibit higher SFRs. However, the z=4.5 LAE correlation is elevated 4-5 times in SFR compared to continuum-selected galaxies at similar redshifts. The exception is the most massive LAEs which have SFRs similar to galaxies at lower redshifts suggesting that they may represent a different population of galaxies than the traditional lower-mass LAEs, perhaps with a different mechanism promoting Lyman alpha photon escape.

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

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

Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.

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

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

We present the discovery of two galaxy overdensities in the Hubble Space Telescope UDF: a proto-cluster, HUDFJ0332.4-2746.6 at $z = 1.84 \pm 0.01$, and a group, HUDFJ0332.5-2747.3 at $z =1.90 \pm 0.01$. Assuming viralization, the velocity dispersion of HUDFJ0332.4-2746.6 implies a mass of $M_{200}= (2.2 \pm 1.8) \times 10^{14} M_{\odot}$, consistent with the lack of extended X-ray emission. Neither overdensity shows evidence of a red sequence. About $50\%$ of their members show interactions and/or disturbed morphologies, which are signatures of merger remnants or disk instability. Most of their ETGs have blue colors and show recent star-formation. These observations reveal for the first time large fractions of spectroscopically confirmed star-forming blue ETGs in proto-clusters at $z\approx 2$. These star-forming ETGs are most likely among the progenitors of the quiescent population in clusters at more recent epochs. Their mass-size relation is consistent with that of passive ETGs in clusters at $z\sim0.7-1.5$. If these galaxies are the progenitors of cluster ETGs at these lower redshifts, their size would evolve according to a similar mass-size relation. It is noteworthy that quiescent ETGs in clusters at $z=1.8-2$ also do not show any significant size evolution over this redshift range, contrary to field ETGs. The ETG fraction is $\lesssim 50\%$, compared to the typical quiescent ETG fraction of $\approx 80\%$ in cluster cores at $z< 1$. The fraction, masses, and colors of the newly discovered ETGs imply that other cluster ETGs will be formed/accreted at a later time.

We present Keck-I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy (SFG) at a redshift $z\sim1.7$. Its spectrum reveals both H$_{\alpha}$ and [NII] emission lines and strong Balmer absorption lines. The H$_{\alpha}$ and Spitzer MIPS 24 $\mu$m fluxes are both weak, thus indicating a low star formation rate of SFR $\lesssim5-10$ M$_{\odot}$ yr$^{-1}$. This, added to a relatively young age of $\sim700$ Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to 3 other galaxies in our sample, by $z\sim1.5$. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, $\sigma^{\rm{gas}}_{\!_{\rm LOS}}=127\pm32$ km s$^{-1}$, is nearly 40% smaller than that of its stars, $\sigma^{\star}_{\!_{\rm LOS}}=215\pm35$ km s$^{-1}$. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of $\sim1.5$ with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.

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

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

In order to derive the AGN contribution to the cosmological ionizing emissivity we have selected faint AGN candidates at $z>4$ in the CANDELS GOODS-South field which is one of the deepest fields with extensive multiwavelength coverage from Chandra, HST, Spitzer and various groundbased telescopes. We have adopted a relatively novel criterion. As a first step high redshift galaxies are selected in the NIR $H$ band down to very faint levels ($H\leq27$) using reliable photometric redshifts. This corresponds at $z>4$ to a selection criterion based on the galaxy rest-frame UV flux. AGN candidates are then picked up from this parent sample if they show X-ray fluxes above a threshold of $F_X\sim 1.5\times 10^{-17}$ cgs (0.5-2 keV). We have found 22 AGN candidates at $z>4$ and we have derived the first estimate of the UV luminosity function in the redshift interval $4<z<6.5$ and absolute magnitude interval $-22.5\lesssim M_{1450} \lesssim -18.5$ typical of local Seyfert galaxies. The faint end of the derived luminosity function is about two/four magnitudes fainter at $z\sim 4-6$ than that derived from previous UV surveys. We have then estimated ionizing emissivities and hydrogen photoionization rates in the same redshift interval under reasonable assumptions and after discussion of possible caveats, the most important being the large uncertainties involved in the estimate of photometric redshift for sources with featureless, almost power-law SEDs and/or low average escape fraction of ionizing photons from the AGN host galaxies. We argue that, under reasonable evaluations of possible biases, the probed AGN population can produce at $z=4-6.5$ photoionization rates consistent with that required to keep highly ionized the intergalactic medium observed in the Lyman-$\alpha$ forest of high redshift QSO spectra, providing an important contribution to the cosmic reionization.

We combine HST/WFC3 imaging and G141 grism observations from the CANDELS and 3D-HST surveys to produce a catalog of grism spectroscopic redshifts for galaxies in the CANDELS/GOODS-South field. The WFC3/G141 grism spectra cover a wavelength range of 1.1<lambda<1.7 microns with a resolving power of R~130 for point sources, thus providing rest-frame optical spectra for galaxies out to z~3.5. The catalog is selected in the H-band (F160W) and includes both galaxies with and without previously published spectroscopic redshifts. Grism spectra are extracted for all H-band detected galaxies with H<24 and a CANDELS photometric redshift z_phot > 0.6. The resulting spectra are visually inspected to identify emission lines and redshifts are determined using cross-correlation with empirical spectral templates. To establish the accuracy of our redshifts, we compare our results against high-quality spectroscopic redshifts from the literature. Using a sample of 411 control galaxies, this analysis yields a precision of sigma_NMAD=0.0028 for the grism-derived redshifts, which is consistent with the accuracy reported by the 3D-HST team. Our final catalog covers an area of 153 square arcmin and contains 1019 redshifts for galaxies in GOODS-S. Roughly 60% (608/1019) of these redshifts are for galaxies with no previously published spectroscopic redshift. These new redshifts span a range of 0.677 < z < 3.456 and have a median redshift of z=1.282. The catalog contains a total of 234 new redshifts for galaxies at z>1.5. In addition, we present 20 galaxy pair candidates identified for the first time using the grism redshifts in our catalog, including four new galaxy pairs at z~2, nearly doubling the number of such pairs previously identified.

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

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

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

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

We have targeted two recently discovered Lyman break galaxies (LBGs) to search for dust continuum and [CII] 158 micron line emission. The strongly lensed z~6.8 LBG A1703-zD1 behind the galaxy cluster Abell 1703, and the spectroscopically confirmed z=7.508 LBG z8-GND-5296 in the GOODS-N field have been observed with the Plateau de Bure interferometer (PdBI) at 1.2mm. These observations have been combined with those of three z>6.5 Lya emitters (named HCM6A, Himiko, and IOK-1), for which deep measurements were recently obtained with the PdBI and ALMA. [CII] is undetected in both galaxies, providing a deep upper limit for Abell1703-zD1, comparable to recent ALMA non-detections. Dust continuum emission from Abell1703-zD1 and z8-GND-5296 is not detected with an rms of 0.12 and 0.16 mJy/beam. From these non-detections we derive upper limits on their IR luminosity and star formation rate, dust mass, and UV attenuation. Thanks to strong gravitational lensing the limit for Abell1703-zD1 is probing the sub-LIRG regime ($L_{IR} <8.1 \times 10^{10}$ Lsun) and very low dust masses ($M_d<1.6 \times 10^7$ Msun). We find that all five galaxies are compatible with the Calzetti IRX-$\beta$ relation, their UV attenuation is compatible with several indirect estimates from other methods (the UV slope, extrapolation of the attenuation measured from the IR/UV ratio at lower redshift, and SED fits), and the dust-to-stellar mass ratio is not incompatible with that of galaxies from z=0 to 3. For their stellar mass the high-z galaxies studied here have an attenuation below the one expected from the mean relation of low redshift (z<1.5) galaxies. More and deeper (sub)-mm data are clearly needed to directly determine the UV attenuation and dust content of the dominant population of high-z star-forming galaxies and to establish more firmly their dependence on stellar mass, redshift, and other properties.

Galaxies' rest-frame ultraviolet (UV) properties are often used to directly infer the degree to which dust obscuration affects the measurement of star formation rates. While much recent work has focused on calibrating dust attenuation in galaxies selected at rest-frame ultraviolet wavelengths, locally and at high-$z$, here we investigate attenuation in dusty, star-forming galaxies (DSFGs) selected at far-infrared wavelengths. By combining multiwavelength coverage across 0.15--500\,$\mu$m in the COSMOS field, in particular making use of {\it Herschel} imaging, and a rich dataset on local galaxies, we find a empirical variation in the relationship between rest-frame UV slope ($\beta$) and ratio of infrared-to-ultraviolet emission ($L_{\rm IR}/L_{\rm UV}\equiv\,IRX$) as a function of infrared luminosity, or total star formation rate, SFR. Both locally and at high-$z$, galaxies above SFR$\gt$50\,M$_\odot$\,yr$^{-1}$ deviate from the nominal $IRX-\beta$ relation towards bluer colors by a factor proportional to their increasing IR luminosity. We also estimate contamination rates of DSFGs on high-$z$ dropout searches of $\ll1$\%\ at $z\lt4-10$, providing independent verification that contamination from very dusty foreground galaxies is low in LBG searches. Overall, our results are consistent with the physical interpretation that DSFGs, e.g. galaxies with $>50$\,M$_\odot$\,yr$^{-1}$, are dominated at all epochs by short-lived, extreme burst events, producing many young O and B stars that are primarily, yet not entirely, enshrouded in thick dust cocoons. The blue rest-frame UV slopes of DSFGs are inconsistent with the suggestion that most DSFGs at $z\sim2$ exhibit steady-state star formation in secular disks.

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

The formation of bars in disk galaxies is a tracer of the dynamical maturity of the population. Previous studies have found that the incidence of bars in disks decreases from the local Universe to z ~ 1, and by z > 1 simulations predict that bar features in dynamically mature disks should be extremely rare. Here we report the discovery of strong barred structures in massive disk galaxies at z ~ 1.5 in deep rest-frame optical images from CANDELS. From within a sample of 876 disk galaxies identified by visual classification in Galaxy Zoo, we identify 123 barred galaxies. Selecting a sub-sample within the same region of the evolving galaxy luminosity function (brighter than L*), we find that the bar fraction across the redshift range 0.5< z < 2 (f_bar = 10.7 +6.3 -3.5% after correcting for incompleteness) does not significantly evolve. We discuss the implications of this discovery in the context of existing simulations and our current understanding of the way disk galaxies have evolved over the last 11 billion years.

We estimate the UV continuum slope, beta, for 923 galaxies in the range 1 < z < 8 in the Hubble Ultradeep Field (HUDF). These data include 460 galaxies at 1 < z < 2 down to an absolute magnitude M_{UV} = -14 (~0.006 L*_{z=1}; 0.02 L*_{z=0}), comparable to dwarf galaxies in the local universe. We combine deep HST/UVIS photometry in F225W, F275W, F336W wavebands (UVUDF) with recent data from HST/WFC3/IR (HUDF12). Galaxies in the range 1 < z < 2 are significantly bluer than local dwarf galaxies. We find their mean (median) values <beta> = -1.382 (-1.830) +/- 0.002 (random) +/- 0.1 (systematic). We find comparable scatter in beta (standard deviation = 0.43) to local dwarf galaxies and 30% larger scatter than z > 2 galaxies. We study the trends of beta with redshift and absolute magnitude for binned sub-samples and find a modest color-magnitude relation, dbeta/dM = -0.11 +/- 0.01 and no evolution in dbeta/dM with redshift. A modest increase in dust reddening with redshift and luminosity, Delta E(B-V) ~ 0.1, and a comparable increase in the dispersion of dust reddening at z < 2, appears likely to explain the observed trends. At z > 2, we find trends that are consistent with previous works; combining our data with the literature in the range 1 < z < 8, we find a color evolution with redshift, dbeta/dz = -0.09 +/-0.01 for low luminosity (0.05 L*_{z=3}), and dbeta/dz = -0.06 +/-0.01 for medium luminosity (0.25 L*_{z=3}) galaxies.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

We present Keck-I MOSFIRE near-infrared spectroscopy for a sample of 13 compact star-forming galaxies (SFGs) at redshift $2\leq z \leq2.5$ with star formation rates of SFR$\sim$100M$_{\odot}$ y$^{-1}$ and masses of log(M/M$_{\odot}$)$\sim10.8$. Their high integrated gas velocity dispersions of $\sigma_{\rm{int}}$=230$^{+40}_{-30}$ km s$^{-1}$, as measured from emission lines of H$_{\alpha}$ and [OIII], and the resultant M$_{\star}-\sigma_{\rm{int}}$ relation and M$_{\star}$$-$M$_{\rm{dyn}}$ all match well to those of compact quiescent galaxies at $z\sim2$, as measured from stellar absorption lines. Since log(M$_{\star}$/M$_{\rm{dyn}}$)$=-0.06\pm0.2$ dex, these compact SFGs appear to be dynamically relaxed and more evolved, i.e., more depleted in gas and dark matter ($<$13$^{+17}_{-13}$\%) than their non-compact SFG counterparts at the same epoch. Without infusion of external gas, depletion timescales are short, less than $\sim$300 Myr. This discovery adds another link to our new dynamical chain of evidence that compact SFGs at $z\gtrsim2$ are already losing gas to become the immediate progenitors of compact quiescent galaxies by $z\sim2$.

We report on two regularly rotating galaxies at redshift z=2, using high resolution spectra of the bright [CII] 158 micron emission line from the HIFI instrument on the Herschel Space Observatory. Both SDSS090122.37+181432.3 ("S0901") and SDSS J120602.09+514229.5 ("the Clone") are strongly lensed and show the double-horned line profile that is typical of rotating gas disks. Using a parametric disk model to fit the emission line profiles, we find that S0901 has a rotation speed v sin(i) = 120 +/- 7 km/s and gas velocity dispersion sigma < 23 km/s. The best fitting model for the Clone is a rotationally supported disk having v sin(i) = 79 +/- 11 km/s and sigma < 4km/s. However the Clone is also consistent with a family of dispersion-dominated models having sigma = 92 +/- 20 km/s. Our results showcase the potential of the [CII] line as a kinematic probe of high redshift galaxy dynamics: [CII] is bright; accessible to heterodyne receivers with exquisite velocity resolution; and traces dense star-forming interstellar gas. Future [CII] line observations with ALMA would offer the further advantage of spatial resolution, allowing a clearer separation between rotation and velocity dispersion.

We present the results of a study investigating the rest-frame ultra-violet (UV) spectral slopes of redshift z~5 Lyman-break galaxies (LBGs). By combining deep Hubble Space Telescope imaging of the CANDELS and HUDF fields with ground-based imaging from the UKIDSS Ultra Deep Survey (UDS), we have produced a large sample of z~5 LBGs spanning an unprecedented factor of >100 in UV luminosity. Based on this sample we find a clear colour-magnitude relation (CMR) at z~5, such that the rest-frame UV slopes (beta) of brighter galaxies are notably redder than their fainter counterparts. We determine that the z~5 CMR is well described by a linear relationship of the form: d beta = (-0.12 +/- 0.02) d Muv, with no clear evidence for a change in CMR slope at faint magnitudes (i.e. Muv > -18.9). Using the results of detailed simulations we are able, for the first time, to infer the intrinsic (i.e. free from noise) variation of galaxy colours around the CMR at z~5. We find significant (12 sigma) evidence for intrinsic colour variation in the sample as a whole. Our results also demonstrate that the width of the intrinsic UV slope distribution of z~5 galaxies increases from Delta(beta)=0.1 at Muv=-18 to Delta(beta)=0.4 at Muv=-21. We suggest that the increasing width of the intrinsic galaxy colour distribution and the CMR itself are both plausibly explained by a luminosity independent lower limit of beta=-2.1, combined with an increase in the fraction of red galaxies in brighter UV-luminosity bins.

We study the far-infrared properties of 498 Lyman Alpha Emitters (LAEs) at z=2.8, 3.1 and 4.5 in the Extended Chandra Deep Field-South, using 250, 350 and 500 micron data from the Herschel Multi-tiered Extragalactic Survey (HerMES) and 870 micron data from the LABOCA ECDFS Submillimeter Survey (LESS). None of the 126, 280 or 92 LAEs at z=2.8, 3.1 and 4.5, respectively, are individually detected in the far-infrared data. We use stacking to probe the average emission to deeper flux limits, reaching $1\sigma$ depths of ~0.1 to 0.4 mJy. The LAEs are also undetected at $\ge3\sigma$ in the stacks, although a $2.5\sigma$ signal is observed at 870 micron for the z=2.8 sources. We consider a wide range of far-infrared spectral energy distributions (SEDs), including a M82 and an Sd galaxy template, to determine upper limits on the far-infrared luminosities and far-infrared-derived star-formation rates of the LAEs. These star-formation rates are then combined with those inferred from the Ly$\alpha$ and UV emission to determine lower limits on the LAEs Ly$\alpha$ escape fraction ($f_{\rm esc}($Ly$\alpha$)). For the Sd SED template, the inferred LAEs $f_{\rm esc}($Ly$\alpha$) are $\gtrsim30%$ ($1\sigma$) at z=2.8, 3.1 and 4.5, which are all significantly higher than the global $f_{\rm esc}($Ly$\alpha$) at these redshifts. Thus, if the LAEs $f_{\rm esc}($Ly$\alpha$) follows the global evolution then they have warmer far-infrared SEDs than the Sd galaxy template. The average and M82 SEDs produce lower limits on the LAE $f_{\rm esc}($Ly$\alpha$) of ~10 to 20% ($1\sigma$), all of which are slightly higher than the global evolution of $f_{\rm esc}($Ly$\alpha$) but consistent with it at the 2 to 3$\sigma$ level.

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.

We present a sample of 33 spectroscopically confirmed z ~ 3.1 Ly$\alpha$-emitting galaxies (LAEs) in the Cosmological Evolution Survey (COSMOS) field. This paper details the narrow-band survey we conducted to detect the LAE sample, the optical spectroscopy we performed to confirm the nature of these LAEs, and a new near-infrared spectroscopic detection of the [O III] 5007 \AA\ line in one of these LAEs. This detection is in addition to two [O III] detections in two z ~ 3.1 LAEs we have reported on previously (McLinden et al 2011). The bulk of the paper then presents detailed constraints on the physical characteristics of the entire LAE sample from spectral energy distribution (SED) fitting. These characteristics include mass, age, star-formation history, dust content, and metallicity. We also detail an approach to account for nebular emission lines in the SED fitting process - wherein our models predict the strength of the [O III] line in an LAE spectrum. We are able to study the success of this prediction because we can compare the model predictions to our actual near-infrared observations both in galaxies that have [O III] detections and those that yielded non-detections. We find a median stellar mass of 6.9 $\times$ 10$^8$ M$_{\odot}$ and a median star formation rate weighted stellar population age of 4.5 $\times$ 10$^6$ yr. In addition to SED fitting, we quantify the velocity offset between the [O III] and Ly$\alpha$ lines in the galaxy with the new [O III] detection, finding that the Ly$\alpha$ line is shifted 52 km s$^{-1}$ redwards of the [O III] line, which defines the systemic velocity of the galaxy.

Ly{\alpha} line EWs provide important clues to the physical nature of high redshift LAEs. However, measuring the Ly{\alpha} EW distribution of high-z narrowband selected LAEs can be hard because many sources do not have broadband photometry. We investigate the possible biases in measuring the intrinsic Ly{\alpha} EW distribution for a LAE sample at z $\sim$ 4.5 in the Extended CDFS. Only weak Malmquist-type bias in both the intrinsic Ly{\alpha} luminosity function and the Ly{\alpha} EW distribution were found. However, the observed EW distribution is severely biased if one only considers LAEs with detections in the continuum. Taking the broadband non-detections into account requires fitting the distribution of the broadband-to-narrowband ratio, which then gives a larger EW distribution scale length. Assuming an exponential form of the intrinsic Ly{\alpha} EW distribution dN/dEW = N exp(-EW/W$_0$), we obtain W$_0$ = 167+/-44{\AA} (uncorrected for IGM absorption of Ly{\alpha}). We discuss the likely range of IGM absorption effects in light of recent measurements of Ly{\alpha} line profiles and velocity offsets. Our data are consistent with Ly{\alpha} EW being independent of UV luminosity (i.e., we do not see evidence for the "Ando" effect). Our simulations also imply that broad-band images should be 0.5-1 magnitude deeper than narrowband images for an effective and reasonably complete LAE survey. Comparing with consistent measurements at other redshifts, we see a strong evolution in Ly{\alpha} EW distribution with redshift which goes as a power-law form of W$_0$ $\prop$(1+z)$^${\xi}, with {\xi} = 1.1$\pm$0.1 (0.6$\pm$0.1) if no IGM corrections are applied to the Ly{\alpha} line; or {\xi} = 1.7$\pm$0.1 (1.2$\pm$0.1) after applying a maximal IGM-absorption correction to Ly{\alpha} line) for an exponential (a gaussian) EW distribution from z = 0.3 to 6.5.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

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

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

We use the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) GOODS-S multi-wavelength catalog to identify counterparts for 20 Lya Emitting (LAE) galaxies at z=2.1. We build several types of stacked Spectral Energy Distributions (SEDs) of these objects. We combine photometry to form average and median flux-stacked SEDs, and postage stamp images to form average and median image-stacked SEDs. We also introduce scaled flux stacks that eliminate the influence of variation in overall brightness. We use the SED fitting code SpeedyMC to constrain the physical properties of individual objects and stacks. Our LAEs at z = 2.1 have stellar masses ranging from 2x10^7 Msun - 8x10^9 Msun (median = 3x10^8 Msun), ages ranging from 4 Myr to 500 Myr (median =100 Myr), and E(B-V) between 0.02 and 0.24 (median = 0.12). We do not observe strong correlations between Lya equivalent width (EW) and stellar mass, age, or E(B-V). The Lya radiative transfer (q) factors of our sample are predominantly close to one and do not correlate strongly with EW or E(B-V), implying that Lya radiative transfer prevents Lya photons from resonantly scattering in dusty regions. The SED parameters of the flux stacks match the average and median values of the individual objects, with the flux-scaled median SED performing best with reduced uncertainties. Median image-stacked SEDs provide a poor representation of the median individual object, and none of the stacking methods captures the large dispersion of LAE properties.

Using spectroscopy from the Large Binocular Telescope and imaging from the Hubble Space Telescope we discovered the first strong galaxy lens at z(lens)>1. The lens has a secure photometric redshift of z=1.53+/-0.09 and the source is spectroscopically confirmed at z=3.417. The Einstein radius (0.35"; 3.0 kpc) encloses 7.6 x 10^10 Msol, with an upper limit on the dark matter fraction of 60%. The highly magnified (40x) source galaxy has a very small stellar mass (~10^8 Msol) and shows an extremely strong [OIII]_5007A emission line (EW_0 ~ 1000A) bolstering the evidence that intense starbursts among very low-mass galaxies are common at high redshift.

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

We present new results on the spectrally resolved Lyman-alpha (LyA) emission of three LyA emitting field galaxies at z~2.4 with high LyA equivalent width (>100 Angstroms) and LyA luminosity (~10^43 erg/s). At 120 km/s (FWHM) spectral resolution, the prominent double-peaked LyA profile straddles the systemic velocity, where the velocity zero-point is determined from spectroscopy of the galaxies' rest-frame optical nebular emission lines. The average velocity offset from systemic of the stronger redshifted emission component for our sample is 176 km/s while the average total separation between the redshifted and main blueshifted emission components is 380 km/s. These measurements are a factor of ~2 smaller than for UV continuum-selected galaxies that show LyA in emission with lower LyA equivalent width. We compare our LyA spectra to the predicted line profiles of a spherical "expanding shell" LyA radiative transfer grid that models large-scale galaxy outflows. Specifically blueward of the systemic velocity where two galaxies show a weak, highly blueshifted (by ~1000 km/s) tertiary emission peak, the model line profiles are a relatively poor representation of the observed spectra. Since the neutral gas column density has a dominant influence over the shape of the LyA line profile, we caution against equating the observed LyA velocity offset with a physical outflow velocity, especially at lower spectral resolution where the unresolved LyA velocity offset is a convoluted function of several degenerate parameters. Referring to rest-frame ultraviolet and optical Hubble Space Telescope imaging, we find that galaxy-galaxy interactions may play an important role in inducing a starburst that results in copious LyA emission, as well as perturbing the gas distribution and velocity field which have strong influence over the LyA emission line profile.

We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and S\'ersic index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (n > 2) reside on the red sequence of the color-magnitude diagram, and most disk-dominated galaxies (n < 2) have colors expected for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U - B colors and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster have larger median effective radii and smaller S\'ersic indices compared to the field with a significance of 2?. These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < Rproj < 1.5 Mpc, which have low S\'ersic indices and possibly larger effective radii, more consistent with star-forming galaxies at this epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are consistent with models where gas accretion onto these galaxies is suppressed from processes associated with the cluster environment.

We present an analysis of the luminosities and equivalent widths of the 284 z < 0.56 [O II]-emitting galaxies found in the 169 square arcmin pilot survey for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). By combining emission-line fluxes obtained from the Mitchell spectrograph on the McDonald 2.7-m telescope with deep broadband photometry from archival data, we derive each galaxy's de-reddened [O II] 3727 luminosity and calculate its total star formation rate. We show that over the last ~5 Gyr of cosmic time there has been substantial evolution in the [O II] emission-line luminosity function, with L* decreasing by ~0.6 +/-0.2 dex in the observed function, and by ~0.9 +/-0.2 dex in the de-reddened relation. Accompanying this decline is a significant shift in the distribution of [O II] equivalent widths, with the fraction of high equivalent-width emitters declining dramatically with time. Overall, the data imply that the relative intensity of star formation within galaxies has decreased over the past ~5 Gyr, and that the star formation rate density of the universe has declined by a factor of ~2.5 between z ~ 0.5 and z ~ 0. These observations represent the first [O II]-based star formation rate density measurements in this redshift range, and foreshadow the advancements which will be generated by the main HETDEX survey.

Star-forming galaxies at redshifts z>6 are likely responsible for the reionization of the universe, and it is important to study the nature of these galaxies. We present three candidates for z~7 Lyman-break galaxies (LBGs) from a 155 arcmin^2 area in the CANDELS/COSMOS field imaged by the deep FourStar Galaxy Evolution (zFourGE) survey. The FourStar medium-band filters provide the equivalent of R~10 spectroscopy, which cleanly distinguishes between z~7 LBGs and brown dwarf stars. The distinction between stars and galaxies based on an object's angular size can become unreliable even when using HST imaging; there exists at least one very compact z~7 candidate (FWHM~0.5-1 kpc) that is indistinguishable from a point source. The medium-band filters provide narrower redshift distributions compared with broad-band-derived redshifts. The UV luminosity function derived using the three z~7 candidates is consistent with previous studies, suggesting an evolution at the bright end (MUV -21.6 mag) from z~7 to z~5. Fitting the galaxies' spectral energy distributions, we predict Lyman-alpha equivalent widths for the two brightest LBGs, and find that the presence of a Lyman-alpha line affects the medium-band flux thereby changing the constraints on stellar masses and UV spectral slopes. This illustrates the limitations of deriving LBG properties using only broad-band photometry. The derived specific star-formation rates for the bright LBGs are ~13 per Gyr, slightly higher than the lower-luminosity LBGs, implying that the star-formation rate increases with stellar mass for these galaxies.

We present a spectroscopically confirmed sample of Lyman alpha emitting galaxies (LAEs) at z ~ 4.5 in the Extended Chandra Deep Field South (ECDFS), which we combine with a sample of z ~ 4.5 LAEs from the Large Area Lyman Alpha (LALA) survey to build a unified Lya luminosity function (LF). We spectroscopically observed 64 candidate LAEs in the ECDFS, confirming 46 objects as z~4.5 LAEs. We did not detect significant flux from neither C_iv 1549\AA\ nor the He_ii 1640\AA\ emission in individual LAE spectra, even with a coadded spectrum. With the coadded line ratio of He_ii to Lya constraining the Population III star formation rate (SFR) to be <0.3% of the total SFR, and <1.25% of the observed SFR (both at the 2-$\sigma$ level). Only one LAE was detected in both the X-ray and radio, while the other objects remained undetected, even when stacked. The Lya LF in our two deepest narrowband filters in the ECDFS differ at >2$\sigma$ significance, and the product $L^*\Phi^*$ differs by a factor of >3. Similar LF differences have been used to infer evolution in the neutral gas fraction in the intergalactic medium at z>6, yet here the difference is likely due to cosmic variance, given that the two samples are from adjoining line-of-sight volumes. Combining our new sample of LAEs with those from previous LALA narrowband surveys at z = 4.5, we obtain one of the best measured Lya LFs to date of L* = 42.83 $\pm$ 0.06 and $\Phi^*$ = -3.48 $\pm$ 0.09. We compare our new LF to others from the literature to study the evolution of the Lya luminosity density at 0 < z < 7. We find tentative evidence for evolution in the product $L^* \Phi^*$, which approximately tracks the cosmic SFR density, but since field-to-field and survey-to-survey variations are in some cases as large as the possible evolution, some caution is needed in interpreting this trend.

We analyze the spectral energy distributions (SEDs) of Lyman break galaxies (LBGs) at z=1-3 selected using the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) UVIS channel filters. These HST/WFC3 observations cover about 50 sq. arcmin in the GOODS-South field as a part of the WFC3 Early Release Science program. These LBGs at z=1-3 are selected using dropout selection criteria similar to high redshift LBGs. The deep multi-band photometry in this field is used to identify best-fit SED models, from which we infer the following results: (1) the photometric redshift estimate of these dropout selected LBGs is accurate to within few percent; (2) the UV spectral slope (beta) is redder than at high redshift (z>3), where LBGs are less dusty; (3) on average, LBGs at z=1-3 are massive, dustier and more highly star-forming, compared to LBGs at higher redshifts with similar luminosities (0.1L*<~L<~2.5L*), though their median values are similar within 1-sigma uncertainties. This could imply that identical dropout selection technique, at all redshifts, find physically similar galaxies; and (4) stellar masses of these LBGs are directly proportional to their UV luminosities with a logarithmic slope of ~0.46, and star-formation rates are proportional to their stellar masses with a logarithmic slope of ~0.90. These relations hold true --- within luminosities probed in this study --- for LBGs from z~1.5 to 5. The star-forming galaxies selected using other color-based techniques show similar correlations at z~2, but to avoid any selection biases, and for direct comparison with LBGs at z>3, a true Lyman break selection at z~2 is essential. The future HST UV surveys, both wider and deeper, covering a large luminosity range are important to better understand LBG properties, and their evolution.

We present the first dynamical mass measurements for Lyman alpha galaxies at high redshift, based on velocity dispersion measurements from rest-frame optical emission lines and size measurements from HST imaging, for a sample of nine galaxies drawn from four surveys. These measurements enable us to study the nature of Lyman alpha galaxies in the context of galaxy scaling relations. The resulting dynamical masses range from 1e9 to 1e10 solar masses. We also fit stellar population models to our sample, and use them to plot the Lyman alpha sample on a stellar mass vs. line width relation. Overall, the Lyman alpha galaxies follow well the scaling relation established by observing star forming galaxies at lower redshift (and without regard for Lyman alpha emission), though in 1/3 of the Lyman alpha galaxies, lower-mass fits are also acceptable. In all cases, the dynamical masses agree with established stellarmass-linewidth relation. Using the dynamical masses as an upper limit on gas mass, we show that Lyman alpha galaxies resemble starbursts (rather than "normal" galaxies) in the relation between gas mass surface density and star formation activity, in spite of relatively modest star formation rates. Finally, we examine the mass densities of these galaxies, and show that their future evolution likely requires dissipational ("wet") merging. In short, we find that Lyman alpha galaxies are low mass cousins of larger starbursts.

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

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

We report the results of a study exploring the stellar populations of 13 luminous (L>L*), spectroscopically confirmed, galaxies in the redshift interval 5.5<z<6.5, all with WFC3/IR and IRAC imaging from the HST/CANDELS and Spitzer/SEDS surveys. Based on fitting the observed photometry with SED templates covering a wide range of different star-formation histories, and a self consistent treatment of Lyman-alpha emission, we find that the derived stellar masses lie within the range 10^9 Msun < M*< 10^10 Msun and are robust to within a factor of two. In contrast, we confirm previous reports that the ages of the stellar populations are poorly constrained. Although the best-fitting models for three objects have ages >= 300 Myr, the degeneracies introduced by dust extinction mean that only two of these objects actually require a >300 Myr old stellar population to reproduce the observed photometry. Moreover, when considering only smoothly-varying star-formation histories, we observe a clear tension between the data and models such that a galaxy SED template with an old age is often chosen in order to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break this tension we explore SED fitting with two-component models (burst plus on-going star-formation) and allow for nebular emission. On average, the inclusion of nebular emission leads to lower stellar-mass estimates (median offset 0.18 dex), moderately higher specific star-formation rates, and allows for a wider range of plausible stellar ages. However, based on our SED modelling, we find no strong evidence for extremely young ages in our sample (<50 Myr). Finally, considering all of the different star-formation histories explored, we find that the median best-fitting ages are of the order 200-300 Myr and that the objects with the tightest constraints indicate ages in the range 50-200 Myr (Abridged).

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

Galaxies selected on the basis of their emission line strength show low metallicities, regardless of their redshifts. We conclude this from a sample of faint galaxies at redshifts between 0.6 < z < 2.4, selected by their prominent emission lines in low-resolution grism spectra in the optical with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) and in the near-infrared using Wide-Field Camera 3 (WFC3). Using a sample of 11 emission line galaxies (ELGs) at 0.6 < z < 2.4 with luminosities of -22 < M_B < -19, which have [OII], H\beta, and [OIII] line flux measurements from the combination of two grism spectral surveys, we use the R23 method to derive the gas-phase oxygen abundances: 7.5 < 12+log(O/H) < 8.5. The galaxy stellar masses are derived using Bayesian based Markov Chain Monte Carlo (\piMC^2) fitting of their Spectral Energy Distribution (SED), and span the mass range 8.1 < log(M_*/M_\sun) < 10.1. These galaxies show a mass-metallicity (M-L) and Luminosity-Metallicity (L-Z) relation, which is offset by --0.6 dex in metallicity at given absolute magnitude and stellar mass relative to the local SDSS galaxies, as well as continuum selected DEEP2 samples at similar redshifts. The emission-line selected galaxies most resemble the local "green peas" galaxies and Lyman-alpha galaxies at z~0.3 and z~2.3 in the M-Z and L-Z relations and their morphologies. The G-M_{20} morphology analysis shows that 10 out of 11 show disturbed morphology, even as the star-forming regions are compact. These galaxies may be intrinsically metal poor, being at early stages of formation, or the low metallicities may be due to gas infall and accretion due to mergers.

We discuss the structural and morphological properties of galaxies in a z=1.62 proto-cluster using near-IR imaging data from Hubble Space Telescope Wide Field Camera 3 data of the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The cluster galaxies exhibit a clear color-morphology relation: galaxies with colors of quiescent stellar populations generally have morphologies consistent with spheroids, and galaxies with colors consistent with ongoing star formation have disk-like and irregular morphologies. The size distribution of the quiescent cluster galaxies shows a deficit of compact (< 1kpc), massive galaxies compared to CANDELS field galaxies at z=1.6. As a result the cluster quiescent galaxies have larger average effective sizes compared to field galaxies at fixed mass at greater than 90% significance. Combined with data from the literature, the size evolution of quiescent cluster galaxies is relatively slow from z~1.6 to the present, growing as (1+z)^(-0.6+/-0.1). If this result is generalizable, then it implies that physical processes associated with the denser cluster region seems to have caused accelerated size growth in quiescent galaxies prior to z=1.6 and slower subsequent growth at z<1.6 compared to galaxies in the lower density field. The quiescent cluster galaxies at z=1.6 have higher ellipticities compared to lower redshift samples at fixed mass, and their surface-brightness profiles suggest that they contain extended stellar disks. We argue the cluster galaxies require dissipationless (i.e., gas-poor or "dry") mergers to reorganize the disk material and to match the relations for ellipticity, stellar mass, size, and color of early-type galaxies in z<1 clusters.

We present Hubble Space Telescope Wide Field Camera 3 slitless grism spectroscopy of 28 emission-line galaxies at z~2, in the GOODS-S region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS). The high sensitivity of these grism observations, with 1-sigma detections of emission lines to f > 2.5x10^{-18} erg/s/cm^2, means that the galaxies in the sample are typically ~7 times less massive (median M_* = 10^{9.5} M_sun) than previously studied z~2 emission-line galaxies. Despite their lower mass, the galaxies have OIII/Hb ratios which are very similar to previously studied z~2 galaxies and much higher than the typical emission-line ratios of local galaxies. The WFC3 grism allows for unique studies of spatial gradients in emission lines, and we stack the two-dimensional spectra of the galaxies for this purpose. In the stacked data the OIII emission line is more spatially concentrated than the Hb emission line with 98.1 confidence. We additionally stack the X-ray data (all sources are individually undetected), and find that the average L(OIII)/L(0.5-10 keV) ratio is intermediate between typical z~0 obscured active galaxies and star-forming galaxies. Together the compactness of the stacked OIII spatial profile and the stacked X-ray data suggest that at least some of these low-mass, low-metallicity galaxies harbor weak active galactic nuclei.

We have coadded X-ray flux of all known Lyman Alpha Emitters (LAEs) in the 4 Msec Chandra Deep Field South (CDF-S) region, to place sensitive upper limits on the average unobscured star-formation rate (SFR_X) in these galaxies. A very small fraction of Lyman-Alpha galaxies in the field are individually detected in the X-rays, implying a low fraction of AGN activity. After excluding the few X-ray detected LAEs, we stack the undetected LAEs located in the 4 Ms CDF-S data and 250 ks Extended CDF-S (ECDFS) data, and compute a 1-\sigma upper limit on SFR_X < 14, 28, 28, 140, 440, 880 M$_{\sun}$ yr$^{-1}$ for LAEs located at z = 2.1, 3.1, 3.2, 4.5, 5.7 and 6.5, respectively. The upper limit of SFR_X in LAEs can be then be compared to SFR$_{Ly\alpha}$ derived from Lyman-Alpha line and thus can constrain on the Lyman-Alpha escape fraction ($f^{Esc}_{Ly\alpha}$). We derive a lower limit on f(Lyman-Alpha Escape) > 14% (84 % confidence level, 1-\sigma lower limit) for LAEs at redshift z ~ 2.1 and z ~ 3.1-3.2. At z > 4, the current LAE samples are not of sufficient size to constrain SFR_X well. By averaging all the LAEs at z> 2, the X-ray non-detection constrains f(Lyman-Alpha Escape) > 17% (84 % confidence level, 1-\sigma lower limit), and rejects f(Lyman-Alpha Escape) < 5.7% at the 99.87% confidence level from 2.1 < z < 6.5.

We have a developed a new method for fitting spectral energy distributions (SEDs) to identify and constrain the physical properties of high-redshift (4 < z < 8) galaxies. Our approach uses an implementation of Bayesian based Markov Chain Monte Carlo (PiMC^2) that allows us to compare observations to arbitrarily complex models and to compute 95% credible intervals that provide robust constraints for the model parameters. The work is presented in 2 sections. In the first, we test PiMC^2 using simulated SEDs to not only confirm the recovery of the known inputs but to assess the limitations of the method and identify potential hazards of SED fitting when applied specifically to high redshift (z>4) galaxies. Our tests reveal five critical results: 1) the ability to confidently constrain metallicity, population ages, and Av all require photometric accuracy better than what is currently achievable (i.e. less than a few percent); 2) the ability to confidently constrain stellar masses (within a factor of two) can be achieved without the need for high-precision photometry; 3) the addition of IRAC photometry does not guarantee that tighter constraints of the stellar masses and ages can be defined; 4) different assumptions about the star formation history can lead to significant biases in mass and age estimates; and 5) we are able to constrain stellar age and Av of objects that are both young and relatively dust free. In the second part of the paper we apply PiMC^2 to 17 4<z<8 objects, including the GRAPES Ly alpha sample (4<z<6), supplemented by HST/WFC3 near-IR observations, and several broad band selected z>6 galaxies. Using PiMC^2, we are able to constrain the stellar mass of these objects and in some cases their stellar age and find no evidence that any of these sources formed at a redshift much larger than z_f=8, a time when the Universe was ~ 0.6 Gyr old.

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

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

We present the analysis of Herschel SPIRE far-infrared (FIR) observations of the z = 2.515 lensed galaxy SMM J163554.2+661225. Combining new 250, 350, and 500 micron observations with existing data, we make an improved fit to the FIR spectral energy distribution (SED) of this galaxy. We find a total infrared (IR) luminosity of L(8--1000 micron) = 6.9 +/- 0.6x10^11 Lsol; a factor of 3 more precise over previous L_IR estimates for this galaxy, and one of the most accurate measurements for any galaxy at these redshifts. This FIR luminosity implies an unlensed star formation rate (SFR) for this galaxy of 119 +/- 10 Msol per yr, which is a factor of 1.9 +/- 0.35 lower than the SFR derived from the nebular Pa-alpha emission line (a 2.5-sigma discrepancy). Both SFR indicators assume identical Salpeter initial mass functions (IMF) with slope Gamma=2.35 over a mass range of 0.1 - 100 Msol, thus this discrepancy suggests that more ionizing photons may be necessary to account for the higher Pa-alpha-derived SFR. We examine a number of scenarios and find that the observations can be explained with a varying star formation history (SFH) due to an increasing star formation rate (SFR), paired with a slight flattening of the IMF. If the SFR is constant in time, then larger changes need to be made to the IMF by either increasing the upper-mass cutoff to ~ 200 Msol, or a flattening of the IMF slope to 1.9 +/- 0.15, or a combination of the two. These scenarios result in up to double the number of stars with masses above 20 Msol, which produce the requisite increase in ionizing photons over a Salpeter IMF with a constant SFH.

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z~1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. Supported by spectroscopic confirmation of strong [OIII] emission lines -- with rest-frame equivalent widths ~1000\AA -- in the four candidates that have HST/WFC3 grism observations, we conclude that these objects are galaxies with 10^8 Msol in stellar mass, undergoing an enormous starburst phase with M_*/(dM_*/dt) of only ~15 Myr. These bursts may cause outflows that are strong enough to produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the co-moving number density (3.7x10^-4 Mpc^-3) can produce in ~4 Gyr much of the stellar mass density that is presently contained in 10^8-10^9 Msol dwarf galaxies. Therefore, our observations provide a strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, short-lived bursts, mostly at z>1.

We show that populations of high redshift galaxies grow more luminous as they grow in linear size. This is because the brightness per unit area has a distinct upper limit due to the self-regulating nature of star-formation. As a corollary, we show that the observed increase in characteristic luminosity of Lyman Break Galaxies (L*) with time can be explained by their increase in size, which scales as 1/H(z). In contrast, Lyman-alpha selected galaxies have a characteristic, constant, small size between redshift z=2.25 and 6.5. Coupled with a characteristic surface brightness, this can explain their non-evolving ultraviolet continuum luminosity function. This compact physical size seems to be a critical determining factor in whether a galaxy will show Lyman-alpha emission. We base these conclusions on new size measurements for a sample of about 170 Lyman-alpha selected galaxies with Hubble Space Telescope broad-band imaging, over the redshift range 2.25 < z < 6. We combine these with a similar number of Lyman-alpha selected galaxies with half-light radii from the literature, and derive surface brightnesses for the entire combined sample.

We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]

We present the first spectroscopic measurements of the [OIII] 5007 A line in two z ~ 3.1 Lyman-alpha emitting galaxies (LAEs) using the new near-infrared instrument LUCIFER1 on the 8.4m Large Binocular Telescope (LBT). We also describe the optical imaging and spectroscopic observations used to identify these Lya emitting galaxies. Using the [OIII] line we have measured accurate systemic redshifts for these two galaxies, and discovered a velocity offset between the [OIII] and Ly-alpha lines in both, with the Lya line peaking 342 and 125 km/s redward of the systemic velocity. These velocity offsets imply that there are powerful outflows in high-redshift LAEs. They also ease the transmission of Lya photons through the interstellar medium and intergalactic medium around the galaxies. By measuring these offsets directly, we can refine both Lya-based tests for reionization, and Lya luminosity function measurements where the Lya forest affects the blue wing of the line. Our work also provides the first direct constraints on the strength of the [OIII] line in high-redshift LAEs. We find [OIII] fluxes of 7 and 36 x 10^-17 erg s^-1 cm^-2 in two z ~ 3.1 LAEs. These lines are strong enough to dominate broad-band flux measurements that include the line (in thiscase, K_s band photometry). Spectral energy distribution fits that do not account for the lines would therefore overestimate the 4000 A (and/or Balmer) break strength in such galaxies, and hence also the ages and stellar masses of such high-z galaxies.

We present Keck/NIRSPEC spectroscopic observations of three Lyman alpha emitting galaxies (LAEs) at z ~ 2.3 discovered with the HETDEX pilot survey. We detect Halpha, [OIII], and Hbeta emission from two galaxies at z = 2.29 and 2.49, designated HPS194 and HPS256, respectively, representing the first detection of multiple rest-frame optical emission lines in galaxies at high-redshift selected on the basis of their Lyman alpha emission. The redshifts of the Lyman alpha emission from these galaxies are offset redward of the systemic redshifts by Delta_v = 162 +/- 37 (photometric) +/- 42 (systematic) km/s for HPS194, and Delta_v = 36 +/- 35 +/- 18 km/s for HPS256. An interpretation for HPS194 is that a large-scale outflow may be occurring in its interstellar medium. The emission line ratios imply that neither LAE hosts an active galactic nucleus. Using the upper limits on the [NII] emission we place meaningful constraints on the gas-phase metallicities in these two LAEs of Z < 0.17 and < 0.28 Zsol (1 sigma). Measuring the stellar masses of these objects via spectral energy distribution (SED) fitting (~ 10^10 and 6 x 10^8 Msol, respectively), we study the nature of LAEs in a mass-metallicity plane. At least one of these two LAEs appears to be more metal poor than continuum-selected star-forming galaxies at the same redshift and stellar mass, implying that objects exhibiting Lyman alpha emission may be systematically less chemically enriched than the general galaxy population. We use the SEDs of these two galaxies to show that neglecting the emission lines when fitting stellar population models to the observed photometry can result in overestimates of the population age by orders of magnitude, and the stellar mass by a factor of ~ 2. This effect is particularly important at z > 7, where similarly strong emission lines may masquerade in the photometry as a 4000 A break (abridged).

We present the results of an optical spectroscopic study of 12 GALEX-discovered star-forming Lyman alpha emitting galaxies (LAEs) at z ~ 0.3. We measure the emission line fluxes from these galaxies by fitting their observed spectra to stellar population models in order to correct for underlying stellar absorption. We revisit earlier stellar population model fitting results, finding that excluding now-known AGNs lowers the typical stellar population age and stellar mass of this sample to ~ 300 Myr and ~ 4 x 10^9 Msol, respectively. We calculate their dust extinction using the Balmer decrement, and find a typical visual attenuation of Av ~ 0.3 mag, similar to that seen in many high-redshift LAEs. Comparing the ratio of Lyalpha/Halpha and the Lyman alpha equivalent widths to the measured dust extinction, we find that the ISMs in these objects appear to be neither enhancing nor seriously attenuating the Lyman alpha equivalent widths, as would be the case in a quasi-clumpy ISM. Lastly, we perform a detailed analysis of the gas-phase metallicities of these galaxies, and we find that most galaxies in our sample have Z < 0.4 Zsol. We find that at a fixed stellar mass, these low-redshift LAE analogs are offset by ~ 0.6 dex lower in metallicity from the general galaxy population at similar redshifts based on the local mass-metallicity relationship. This implies that galaxies with Lyman alpha in emission may be systematically more metal poor than star-forming galaxies at the same stellar mass and redshift, similar to preliminary results at z ~ 2.

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

We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.

Lyman alpha (Lya) emission lines should be attenuated in a neutral intergalactic medium (IGM). Therefore the visibility of Lya emitters at high redshifts can serve as a valuable probe of reionization at about the 50% level. We present an imaging search for z=7.7 Lya emitting galaxies using an ultra-narrowband filter (filter width= 9A) on the NEWFIRM imager at the Kitt Peak National Observatory. We found four candidate Lya emitters in a survey volume of 1.4 x 10^4 Mpc^3, with a line flux brighter than 6x10^-18 erg/cm^2/s (5 sigma in 2" aperture). We also performed a detailed Monte-Carlo simulation incorporating the instrumental effects to estimate the expected number of Lya emitters in our survey, and found that we should expect to detect one Lya emitter, assuming a non-evolving Lya luminosity function (LF) between z=6.5 and z=7.7. Even if one of the present candidates is spectroscopically confirmed as a z~8 Lya emitter, it would indicate that there is no significant evolution of the Lya LF from z=3.1 to z~8. While firm conclusions would need both spectroscopic confirmations and larger surveys to boost the number counts of galaxies, we successfully demonstrate the feasibility of sensitive near-infrared (1.06 um) narrow-band searches using custom filters designed to avoid the OH emission lines that make up most of the sky background.

In this proceeding we present the results from a study of very high-redshift galaxies with the newly commissioned Wide Field Camera 3 on the Hubble Space Telescope. With the deepest near-infrared data ever taken, we discovered 31 galaxies at 6.3 < z < 8.6. The rest-frame ultraviolet (UV) colors of these galaxies are extremely blue, showing significant (> 4 sigma) evolution from z ~ 3, over only 1 Gyr of cosmic time. While we cannot yet diagnose the exact cause of the bluer colors, it appears a low dust content is the primary factor. The stellar masses of these galaxies are less than comparably selected galaxies at 3 < z < 6, highlighting evolution in the stellar mass of characteristic (L*) galaxies with redshift. Lastly, the measured rest-UV luminosity density of galaxies in our sample seems sufficient to sustain reionization at z ~ 7 when we account for the likely contribution from galaxies below our magnitude limit.

We measure the rest-frame colors (dust-corrected), infrared luminosities, star formation rates, and stellar masses of 92 galaxies in a Spitzer-selected cluster at z=1.62. By fitting spectral energy distributions (SEDs) to 10-band photometry (0.4 micron < lambda(obs) <8 micron) and measuring 24 micron fluxes for the 12 spectroscopically confirmed and 80 photometrically selected members, we discover an exceptionally high level of star formation in the cluster core of ~1700 Msun/yr per Mpc^2. The cluster galaxies define a strong blue sequence in (U-V) color and span a range in color. We identify 17 members with L(IR)>10^(11) Lsun, and these IR luminous members follow the same trend of increasing star formation with stellar mass that is observed in the field at z~2. Using rates derived from both the 24 micron imaging and SED fitting, we find that the relative fraction of star-forming members triples from the lowest to highest galaxy density regions, e.g. the IR luminous fraction increases from ~8% at Sigma~10 gal per Mpc^2 to ~25% at Sigma>100 gal per Mpc^2. The observed increase is a reversal of the well-documented trend at z<1 and signals that we have reached the epoch when massive cluster galaxies are still forming a substantial fraction of their stars.

We present the highest redshift detections of resolved Lyman alpha emission, using Hubble Space Telescope/ACS F658N narrowband-imaging data taken in parallel with the Wide Field Camera 3 Early Release Science program in the GOODS CDF-S. We detect Lyman alpha emission from three spectroscopically confirmed z = 4.4 Lyman alpha emitting galaxies (LAEs), more than doubling the sample of LAEs with resolved Lyman alpha emission. Comparing the light distribution between the rest-frame ultraviolet continuum and narrowband images, we investigate the escape of Lyman alpha photons at high redshift. While our data do not support a positional offset between the Lyman alpha and rest-frame ultraviolet (UV) continuum emission, the half-light radii in two out of the three galaxies are significantly larger in Lyman alpha than in the rest-frame UV continuum. This result is confirmed when comparing object sizes in a stack of all objects in both bands. Additionally, the narrowband flux detected with HST is significantly less than observed in similar filters from the ground. These results together imply that the Lyman alpha emission is not strictly confined to its indigenous star-forming regions. Rather, the Lyman alpha emission is more extended, with the missing HST flux likely existing in a diffuse outer halo. This suggests that the radiative transfer of Lyman alpha photons in high-redshift LAEs is complicated, with the interstellar-medium geometry and/or outflows playing a significant role in galaxies at these redshifts.

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

We report the first X-ray detection of 113 Lyman-alpha emitters at redshift z ~ 4.5. Only one source (J033127.2-274247) is detected in the Extended Chandra Deep Field South (ECDF-S) X-ray data, and has been spectroscopically confirmed as a z = 4.48 quasar with $L_X = 4.2\times 10^{44}$ erg/s. The single detection gives a Lyman-alpha quasar density consistent with the X-ray luminosity function of quasars. The coadded counts of 22 Lyman-alpha emitters (LAEs) in the central Chandra Deep Field South (CDF-S) region yields a S/N=2.4 (p=99.83%) detection at soft band, with an effective exposure time of ~36 Ms. Further analysis of the equivalent width (EW) distribution shows that all the signal comes from 12 LAE candidates with EW_rest < 400 \AA, and 2 of them contribute about half of the signal. Following-up spectroscopic observations show that the two are a low-redshift emission line galaxy and a Lyman break galaxy at z = 4.4. Excluding these two and combined with ECDF-S data, we derive a 3-sigma upper limit on the average luminosity of $L_{0.5-2 keV}$ $<$ 2.4 $\times 10^{42}$ ergs/s for z ~ 4.5 LAEs. If the average X-ray emission is due to star formation, it corresponds to a star-formation rate (SFR) of < 180--530 M$_\sun$ per yr. We use this SFR_X as an upper limit of the unobscured SFR to constrain the escape fraction of Lyman-alpha photons, and find a lower limit of f_esc > 3-10%. However, our upper limit on the SFR_X is ~7 times larger than the upper limit on SFR_X on z ~ 3.1 LAEs in the same field, and at least 30 times higher than the SFR estimated from Lyman-alpha emission. From the average X-ray to Lyman-alpha line ratio, we estimate that fewer than 3.2% (6.3%) of our LAEs could be high redshift type 1 (type 2) AGNs, and those hidden AGNs likely show low rest frame EWs.

We report the discovery of a galaxy cluster at z=1.62 located in the Spitzer Wide-Area Infrared Extragalactic survey XMM-LSS field. This structure was selected solely as an overdensity of galaxies with red Spitzer/IRAC colors, satisfying [3.6]-[4.5] > -0.1 AB mag. Photometric redshifts derived from Subaru XMM Deep Survey (BViz-bands), UKIRT Infrared Deep Survey-Ultra-Deep Survey (UKIDSS-UDS, JK-bands), and from the Spitzer Public UDS survey (3.6-8.0 micron) show that this cluster corresponds to a surface density of galaxies at z ~ 1.6 that is more than 20 sigma above the mean at this redshift. We obtained optical spectroscopic observations of galaxies in the cluster region using IMACS on the Magellan telescope. We measured redshifts for seven galaxies in the range z=1.62-1.63 within 2.8 arcmin (<1.4 Mpc) of the astrometric center of the cluster. A posteriori analysis of the XMM data in this field reveal a weak (4 sigma) detection in the [0.5-2 keV] band compatible with the expected thermal emission from such a cluster. The color-magnitude diagram of the galaxies in this cluster shows a prominent red-sequence, dominated by a population of red galaxies with (z-J) > 1.7 mag. The photometric redshift probability distributions for the red galaxies are strongly peaked at z=1.62, coincident with the spectroscopically confirmed galaxies. The rest-frame (U-B) color and scatter of galaxies on the red-sequence are consistent with a mean luminosity-weighted age of 1.2 +/- 0.1 Gyr, yielding a formation redshift z_f = 2.35 +/- 0.10, and corresponding to the last significant star-formation period in these galaxies.

As some of the first known objects to exist in the Universe, Lyman alpha emitting galaxies (LAEs) naturally draw a lot of interest. First discovered over a decade ago, they have allowed us to probe the early Universe, as their strong emission line compensates for their faint continuum light. While initially thought to be indicative of the first galaxies forming in the Universe, recent studies have shown them to be increasingly complex, as some fraction appear evolved, and many LAEs appear to be dusty, which one would not expect from primordial galaxies. Presently, much interest resides in discovering not only the highest redshift galaxies to constrain theories of reionization, but also pushing closer to home, as previous ground-based studies have only found LAEs at z > 3 due to observational limitations. In this review talk I will cover everything from the first theoretical predictions of LAEs, to their future prospects for study, including the HETDEX survey here in Texas.

We report the detection of the Paschen-alpha emission line in the z=2.515 galaxy SMM J163554.2+661225 using Spitzer spectroscopy. SMM J163554.2+661225 is a sub-millimeter-selected infrared (IR)-luminous galaxy maintaining a high star-formation rate (SFR), with no evidence of an AGN from optical or infrared spectroscopy, nor X-ray emission. This galaxy is lensed gravitationally by the cluster Abell 2218, making it accessible to Spitzer spectroscopy. Correcting for nebular extinction derived from the H-alpha and Pa-alpha lines, the dust-corrected luminosity is L(Pa-alpha) = (2.57+/-0.43) x 10^43 erg s^-1, which corresponds to an ionization rate, Q = (1.6+/-0.3) x 10^55 photons s^-1. The instantaneous SFR is 171+/-28 solar masses per year, assuming a Salpeter-like initial mass function. The total IR luminosity derived using 70, 450, and 850 micron data is L(IR) = (5-10) x 10^11 solar luminosities, corrected for gravitational lensing. This corresponds to a SFR=90-180 solar masses per year, where the upper range is consistent with that derived from the Paschen-alpha luminosity. While the L(8 micron) / L(Pa-alpha) ratio is consistent with the extrapolated relation observed in local galaxies and star-forming regions, the rest-frame 24 micron luminosity is significantly lower with respect to local galaxies of comparable Paschen-alpha luminosity. Thus, SMM J163554.2+661225 arguably lacks a warmer dust component (T ~ 70 K), which is associated with deeply embedded star formation, and which contrasts with local galaxies with comparable SFRs. Rather, the starburst is consistent with star-forming local galaxies with intrinsic luminosities, L(IR) ~ 10^10 solar luminosities, but "scaled-up" by a factor of 10-100.

With the coming generation of instruments and telescopes capable of spectroscopy of high redshift galaxies, the spectral synthesis technique in the rest-frame UV and Far-UV range will become one of a few number of tools remaining to study their young stellar populations in detail. The rest-frame UV lines and continuum of high redshift galaxies, observed with visible and infrared telescopes on Earth, can be used for accurate line profile fitting such as PV@1118,1128A, CIII@1176A, and CIV@1550A. These lines are very precise diagnostic tools to estimate ages, metallicities, and masses of stellar populations. Here we discuss the potential for spectral synthesis of rest-frame UV spectra obtained at the Keck telescope. As an example, we study the 8 o'clock arc, a lensed galaxy at z=2.7322. We show that the poor spectral type coverage of the actual UV empirical spectral libraries limits the age and metallicity diagnostic. In order to improve our knowledge of high redshift galaxies using spectral synthesis, UV stellar libraries need to be extended to obtain accurate age, metallicity, and mass estimates likely to be occuring in young stellar populations observed in the early universe.

We investigate the fraction of z ~ 0.3 Lyman alpha emitting galaxies (LAEs) which host active galactic nucleus activity, which is typically from 1 - 5% at 2 < z < 6. Using optical spectroscopy of 23 LAEs at 0.2 < z < 0.45 selected with GALEX UV data, we probe for AGN with a variety of methods, including line widths, diagnostic line ratios, high-ionization emission, infrared activity and X-ray luminosity. We found that our sample of low-redshift LAEs has an AGN fraction of 43 ^{+18}{-26}%, significantly higher than at high redshift. While previous results have shown that low-redshift LAEs have a lower space density than their high-redshift analogs, these results show that star-forming LAEs at low-redshift are rarer still. Accounting for differences in available AGN classification methods, we conclude that rest-frame optical spectroscopy is necessary to identify low-luminosity AGNs in LAEs at high redshift, and that limits on the X-ray flux are necessary to determine the extent to which the AGN contaminates the Lyman alpha fluxes.

We present a large sample of Ly-alpha emitting galaxies spectroscopically confirmed at z~4.5, based on IMACS spectroscopic observations of candidate LAEs in the LALA Cetus field. We identify 110 z~4.5 LAEs based on single line detections with no continuum emission blueward of the line. The Ly-alpha confirmation rate varies from <50% to 76% for candidates selected in different narrowband filters at slightly different redshifts. We find a drop in the LAE density at redshift 4.50+-0.03 from redshift 4.39-4.47 by a factor of 66%, which could be a large scale void in the distribution of star-forming galaxies (~18Mpc along the line of sight and ~80Mpc across). The sample includes many objects with equivalent widths >200A. These large EW candidates are spectroscopically confirmed at the same rate as candidates with more modest EWs. A composite spectrum of all 110 confirmed LAEs shows the characteristic asymmetry of the Ly-alpha line. It also places new and stringent upper limits on the CIV 1549/Ly-alpha and HeII 1640/Ly-alpha line ratios, providing a new upper limit on the fraction of active galactic nuclei in Ly-alpha selected galaxy samples, and on the contribution of Pop III populations. Finally, we calculate the Ly-alpha luminosity function for our z~4.5 sample, which is consistent with those at other redshifts, showing that there is no evolution in Ly-alpha luminosity function from z~3.1-6.6.

We present the results from an optical and near-infrared spectroscopic study of the ultraviolet-luminous z = 2.73 galaxy, the 8 o'clock arc. Due to gravitational lensing, this galaxy is magnified by a factor of > 10, allowing in-depth measurements which are usually unfeasible at such redshifts. In the optical spectra, we measured the systemic redshift of the galaxy, z = 2.7322 +/- 0.0012, using stellar photospheric lines. This differs from the redshift of absorption lines in the interstellar medium, z = 2.7302 +/- 0.0006, implying gas outflows on the order of 160 km/s. With H and K-band near-infrared spectra, we have measured nebular emission lines of Halpha, Hbeta, Hgamma, [N II] and [O III], which have a redshift z = 2.7333 +/- 0.0001, consistent with the derived systemic redshift. From the Balmer decrement, we measured the dust extinction to be A_5500 = 1.17 +/- 0.36 mag. Correcting Halpha for dust extinction and the assumed lensing factor, we measure a star-formation rate of ~ 270 Msol/yr, which is higher than ~ 85% of star-forming galaxies at z ~ 2-3. Using combinations of all detected emission lines, we find that the 8 o'clock arc has a gas-phase metallicity of ~ 0.8 Zsol, showing that enrichment at high-redshift is not rare, even in blue, star-forming galaxies. Studying spectra from two of the arc components separately, we find that one component dominates the dust extinction and star-formation rate, although the metallicities between the two components are similar. We derive the mass via stellar population modeling, and find that the arc has a total stellar mass of ~ 4 x 10^11 Msol, which falls on the mass-metallicity relation at z ~ 2. Finally, we estimate the total gas mass, and find it to be only ~ 12% of the stellar mass, implying that the 8 o'clock arc is likely nearing the end of a starburst.

We present the results of a stellar population analysis of 30 Lyman alpha emitting galaxies (LAEs) at z ~ 0.3, previously discovered with the Galaxy Evolution Explorer (GALEX). With a few exceptions, we can accurately fit model spectral energy distributions to these objects, representing the first time this has been done for a large sample of LAEs at z < 3, a gap of ~ 8 Gyr in the history of the Universe. From the 26/30 LAEs which we can fit, we find an age and stellar mass range of 200 Myr - 10 Gyr and 10^9 - 10^11 Msol, respectively. These objects thus appear to be significantly older and more massive than LAEs at high-redshift. We also find that these LAEs show a mild trend towards higher metallicity than those at high redshift, as well as a tighter range of dust attenuation and interstellar medium geometry. These results suggest that low-redshift LAEs have evolved significantly from those at high redshift.

Recent results have shown that a substantial fraction of high-redshift Lyman alpha galaxies contain considerable amounts of dust. This implies that Lyman alpha galaxies are not primordial, as has been thought in the past. However, this dust has not been directly detected in emission; rather it has been inferred based on extinction estimates from rest-frame ultraviolet (UV) and optical observations. This can be tricky, as both dust and old stars redden galactic spectra at the wavelengths used to infer dust. Measuring dust emission directly from these galaxies is thus a more accurate way to estimate the total dust mass, giving us real physical information on the stellar populations and interstellar medium (ISM) enrichment. New generation instruments such as the Atacama Large Millimeter Array (ALMA) and Sub-Millimeter Array (SMA), should be able to detect dust emission from some of these galaxies in the sub-mm. Using measurements of the UV spectral slopes, we derive far-infrared flux predictions for of a sample of 23 z > 4 Lyman alpha galaxies. We find that in only a few hours, we can detect dust emission from 39 +/- 22% of our Lyman alpha galaxies. Comparing these results to those found from a sample of 21 Lyman break galaxies (LBGs), we find that LBGs are on average 60% more likely to be detected than Lyman alpha galaxies, implying that they are more dusty, and thus indicating an evolutionary difference between these objects. These observations will provide better constraints on dust in these galaxies than those derived from their UV and optical fluxes alone. Undeniable proof of dust in these galaxies could explain the larger than expected Lyman alpha equivalent widths seen in many Lyman alpha galaxies today.

We present stellar population modeling results for 10 newly discovered Lyman alpha emitting galaxies (LAEs), as well as four previously known LAEs at z ~ 4.5 in the Chandra Deep Field - South. We fit stellar population models to these objects in order to learn specifically if there exists more than one class of LAE. Past observational and theoretical evidence has shown that while many LAEs appear to be young, they may be much older, with Lyman alpha EWs enhanced due to resonant scattering of Lyman alpha photons in a clumpy interstellar medium (ISM). Our results show a large range of stellar population age (3 - 500 Myr), stellar mass (1.6 x 10^8 - 5.0 x 10^10 Msol) and dust extinction (A_1200 = 0.3 - 4.5 mag), broadly consistent with previous studies. With such a large number of individually analyzed objects, we have looked at the distribution of stellar population ages in LAEs for the first time, and we find a very interesting bimodality, in that our objects are either very young (< 15 Myr) or old (> 450 Myr). This bimodality may be caused by dust, and it could explain the Lyman alpha duty cycle which has been proposed in the literature. We find that eight of the young objects are best fit with a clumpy ISM. We find that dust geometry appears to play a large role in shaping the SEDs that we observe, and that it may be a major factor in the observed Lyman alpha equivalent width distribution in high redshift Lyman alpha galaxies, although other factors (i.e. outflows) may be in play. We conclude that 12 out of our 14 LAEs are dusty star-forming galaxies, with the other two LAEs being evolved galaxies.

Equivalent widths (EWs) observed in high-redshift Lyman alpha galaxies could be stronger than the EW intrinsic to the stellar population if dust is present residing in clumps in the inter-stellar medium (ISM). In this scenario, continuum photons could be extinguished while the Lyman alpha photons would be resonantly scattered by the clumps, eventually escaping the galaxy. We investigate this radiative transfer scenario with a new sample of six Lyman alpha galaxy candidates in the GOODS CDF-S, selected at z = 4.4 with ground-based narrow-band imaging obtained at CTIO. Grism spectra from the HST PEARS survey confirm that three objects are at z = 4.4, and that another object contains an active galactic nuclei (AGN). If we assume the other five (non-AGN) objects are at z = 4.4, they have rest-frame EWs from 47 -- 190 A. We present results of stellar population studies of these objects, constraining their rest-frame UV with HST and their rest-frame optical with Spitzer. Out of the four objects which we analyzed, three objects were best-fit to contain stellar populations with ages on the order of 1 Myr and stellar masses from 3 - 10 x 10^8 solar masses, with dust in the amount of A_1200 = 0.9 - 1.8 residing in a quasi-homogeneous distribution. However, one object (with a rest EW ~ 150 A) was best fit by an 800 Myr, 6.6 x 10^9 solar mass stellar population with a smaller amount of dust (A_1200 = 0.4) attenuating the continuum only. In this object, the EW was enhanced ~ 50% due to this dust. This suggests that large EW Lyman alpha galaxies are a diverse population. Preferential extinction of the continuum in a clumpy ISM deserves further investigation as a possible cause of the overabundance of large-EW objects that have been seen in narrow-band surveys in recent years.

The 174 ks Chandra Advanced CCD Imaging Spectrometer exposure of the Large Area Lyman Alpha Survey (LALA) Cetus field is the second of the two deep Chandra images on LALA fields. In this paper we present the Chandra X-ray sources detected in the Cetus field, along with an analysis of X-ray source counts, stacked X-ray spectrum, and optical identifications. A total of 188 X-ray sources were detected: 174 in the 0.5-7.0 keV band, 154 in the 0.5-2.0 keV band, and 113 in the 2.0-7.0 keV band. The X-ray source counts were derived and compared with LALA Bootes field (172 ks exposure). Interestingly, we find consistent hard band X-ray source density, but 36+-12% higher soft band X-ray source density in Cetus field. The weighted stacked spectrum of the detected X-ray sources can be fitted by a powerlaw with photon index Gamma = 1.55. Based on the weighted stacked spectrum, we find that the resolved fraction of the X-ray background drops from 72+-1% at 0.5-1.0 keV to 63+-4% at 6.0-8.0 keV. The unresolved spectrum can be fitted by a powerlaw over the range 0.5-7 keV, with a photon index Gamma = 1.22. We also present optical counterparts for 154 of the X-ray sources, down to a limiting magnitude of r' = 25.9 (Vega), using a deep r' band image obtained with the MMT.

We examine the stellar populations of a sample of 98 redshift 4.5 Lyman alpha emitting galaxies using their broadband colors derived from deep photometry at the MMT. These galaxies were selected by narrowband excess from the Large Area Lyman Alpha survey. Twenty-two galaxies are detected in two or more of our MMT filters (g', r', i' and z'). By comparing broad and narrowband colors of these galaxies to synthetic colors from stellar population models, we determine their ages and stellar masses. The highest equivalent width objects have an average age of 4 Myr, consistent with ongoing star formation. The lowest EW objects show an age of 40 - 200 Myr, consistent with the expectation that larger numbers of stars are causing low EWs. We found masses ranging from 2e7 solar masses for the youngest objects in the sample to 2e9 solar masses for the oldest. It is possible that dust effects could produce large equivalent widths even in older populations by allowing the Lyman alpha photons to escape, even while the continuum is extinguished, and we present models for this scenario also.

We present new optical emission-line images of the young SNR 1E 0102-7219 (E0102) in the SMC obtained with the HST Advanced Camera for Surveys (ACS). E0102 is a member of the oxygen-rich class of SNRs showing strong oxygen, neon , and other metal-line emissions in its optical and X-ray spectra, and an absence of H and He. The progenitor of E0102 may have been a Wolf-Rayet star that underwent considerable mass loss prior to exploding as a Type Ib/c or IIL/b SN. The ejecta in this SNR are fast-moving (V > 1000 km/s) and emit as they are compressed and heated in the reverse shock. In 2003, we obtained optical [O III], H-alpha, and continuum images with the ACS Wide Field Camera. The [O III] image captures the full velocity range of the ejecta, and shows considerable high-velocity emission projected in the middle of the SNR that was Doppler-shifted out of the narrow F502N bandpass of a previous Wide Field and Planetary Camera 2 image from 1995. Using these two epochs separated by ~8.5 years, we measure the transverse expansion of the ejecta around the outer rim in this SNR for the first time at visible wavelengths. From proper-motion measurements of 12 ejecta filaments, we estimate a mean expansion velocity for the bright ejecta of ~2000 km/s and an inferred kinematic age for the SNR of \~2050 +/- 600 years. The age we derive from HST data is about twice that inferred by Hughes et al.(2000) from X-ray data, though our 1-sigma error bars overlap. Our proper-motion age is consistent with an independent optical kinematic age derived by Eriksen et al.(2003) using spatially resolved [O III] radial-velocity data. We derive an expansion center that lies very close to X-ray and radio hotspots, which could indicate the presence of a compact remnant (neutron star or black hole).

We present near-infrared magnitudes for all white dwarfs (selected from the catalog of McCook & Sion) contained in the 2 Micron All Sky Survey Second Incremental Data Release(2MASS 2IDR). We show that the near-IR color-color diagram is an effective means of identifying candidate binary stars containing a WD and a low mass main sequence star. The loci of single WDs and WD + red dwarf binaries occupy distinct regions of the near-IR color-color diagram. We recovered all known unresolved WD + red dwarf binaries located in the 2IDR sky coverage, and also identified as many new candidate binaries (47 new candidates out of 95 total). Using observational near-IR data for WDs and M-L dwarfs, we have compared a sample of simulated WD + red dwarf binaries with our 2MASS data. The colors of the simulated binaries are dominated by the low mass companion through the late-M to early-L spectral types. As the spectral type of the companion becomes progressively later, however, the colors of unresolved binaries become progressively bluer. Binaries containing the lowest mass companions will be difficult to distinguish from single WDs solely on the basis of their near-IR colors.