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.

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

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.

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.

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

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.

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.

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

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.

We studied Lyman-$\alpha$ (Ly$\alpha$) escape in a statistical sample of 43 Green Peas with HST/COS Ly$\alpha$ spectra. Green Peas are nearby star-forming galaxies with strong [OIII]$\lambda$5007 emission lines. Our sample is four times larger than the previous sample and covers a much more complete range of Green Pea properties. We found that about 2/3 of Green Peas are strong Ly$\alpha$ line emitters with rest-frame Ly$\alpha$ equivalent width $>20$ \AA. The Ly$\alpha$ profiles of Green Peas are diverse. The Ly$\alpha$ escape fraction, defined as the ratio of observed Ly$\alpha$ flux to intrinsic Ly$\alpha$ flux, shows anti-correlations with a few Ly$\alpha$ kinematic features -- both the blue peak and red peak velocities, the peak separations, and FWHM of the red portion of the Ly$\alpha$ profile. Using properties measured from SDSS optical spectra, we found many correlations -- Ly$\alpha$ escape fraction generally increases at lower dust reddening, lower metallicity, lower stellar mass, and higher [OIII]/[OII] ratio. We fit their Ly$\alpha$ profiles with the HI shell radiative transfer model and found Ly$\alpha$ escape fraction anti-correlates with the best-fit $N_{HI}$. Finally, we fit an empirical linear relation to predict Ly$\alpha$ escape fraction from the dust extinction and Ly$\alpha$ red peak velocity. The standard deviation of this relation is about 0.3 dex. This relation can be used to isolate the effect of IGM scatterings from Ly$\alpha$ escape and to probe the IGM optical depth along the line of sight of each $z>7$ Ly$\alpha$ emission line galaxy in the JWST era.

For the first time, we present the size evolution of a mass-complete (log(M*/Msol)>10) sample of star-forming galaxies over redshifts z=1-7, selected from the FourStar Galaxy Evolution Survey (ZFOURGE). Observed H-band sizes are measured from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) Hubble Space Telescope (HST)/F160W imaging. Distributions of individual galaxy masses and sizes illustrate that a clear mass-size relation exists up to z~7. At z~7, we find that the average galaxy size from the mass-size relation is more compact at a fixed mass of log(M*/Msol)=10.1, with r_1/2,maj=1.02+/-0.29 kpc, than at lower redshifts. This is consistent with our results from stacking the same CANDELS HST/F160W imaging, when we correct for galaxy position angle alignment. We find that the size evolution of star-forming galaxies is well fit by a power law of the form r_e = 7.07(1 + z)^-0.89 kpc, which is consistent with previous works for normal star-formers at 1<z<4. In order to compare our slope with those derived Lyman break galaxy studies, we correct for different IMFs and methodology and find a slope of -0.97+/-0.02, which is shallower than that reported for the evolution of Lyman break galaxies at z>4 (r_e\propto(1 +z)^-1.2+/-0.06). Therefore, we conclude the Lyman break galaxies likely represent a subset of highly star-forming galaxies that exhibit rapid size growth 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 FourStar galaxy evolution survey (ZFOURGE) is a 45 night legacy program with the FourStar near-infrared camera on Magellan and one of the most sensitive surveys to date. ZFOURGE covers a total of $400\ \mathrm{arcmin}^2$ in cosmic fields CDFS, COSMOS and UDS, overlapping CANDELS. We present photometric catalogs comprising $>70,000$ galaxies, selected from ultradeep $K_s$-band detection images ($25.5-26.5$ AB mag, $5\sigma$, total), and $>80\%$ complete to $K_s<25.3-25.9$ AB. We use 5 near-IR medium-bandwidth filters ($J_1,J_2,J_3,H_s,H_l$) as well as broad-band $K_s$ at $1.05\ - 2.16\ \mu m$ to $25-26$ AB at a seeing of $\sim0.5$". Each field has ancillary imaging in $26-40$ filters at $0.3-8\ \mu m$. We derive photometric redshifts and stellar population properties. Comparing with spectroscopic redshifts indicates a photometric redshift uncertainty $\sigma_z={0.010,0.009}$, and 0.011 in CDFS, COSMOS, and UDS. As spectroscopic samples are often biased towards bright and blue sources, we also inspect the photometric redshift differences between close pairs of galaxies, finding $\sigma_{z,pairs}= 0.01-0.02$ at $1<z<2.5$. We quantify how $\sigma_{z,pairs}$ depends on redshift, magnitude, SED type, and the inclusion of FourStar medium bands. $\sigma_{z,pairs}$ is smallest for bright, blue star-forming samples, while red star-forming galaxies have the worst $\sigma_{z,pairs}$. Including FourStar medium bands reduces $\sigma_{z,pairs}$ by 50\% at $1.5<z<2.5$. We calculate SFRs based on ultraviolet and ultradeep far-IR $Spitzer$/MIPS and Herschel/PACS data. We derive rest-frame $U-V$ and $V-J$ colors, and illustrate how these correlate with specific SFR and dust emission to $z=3.5$. We confirm the existence of quiescent galaxies at $z\sim3$, demonstrating their SFRs are suppressed by $>\times15$.

We build a set of composite galaxy SEDs by de-redshifting and scaling multi-wavelength photometry from galaxies in the ZFOURGE survey, covering the CDFS, COSMOS, and UDS fields. From a sample of ~4000 K_s-band selected galaxies, we define 38 composite galaxy SEDs that yield continuous low-resolution spectra (R~45) over the rest-frame range 0.1-4 um. Additionally, we include far infrared photometry from the Spitzer Space Telescope and the Herschel Space Observatory to characterize the infrared properties of our diverse set of composite SEDs. From these composite SEDs we analyze the rest-frame UVJ colors, as well as the ratio of IR to UV light (IRX) and the UV slope ($\beta$) in the IRX$-\beta$ dust relation at 1<z<3. Blue star-forming composite SEDs show IRX and $\beta$ values consistent with local relations; dusty star-forming galaxies have considerable scatter, as found for local IR bright sources, but on average appear bluer than expected for their IR fluxes. We measure a tight linear relation between rest-frame UVJ colors and dust attenuation for star-forming composites, providing a direct method for estimating dust content from either (U-V) or (V-J) rest-frame colors for star-forming galaxies at intermediate redshifts.

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3\sigma$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

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.

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]

We study the rest-frame ultra-violet sizes of massive (~0.8 x 10^11 M_Sun) galaxies at 3.4<z<4.2, selected from the FourStar Galaxy Evolution Survey (ZFOURGE), by fitting single Sersic profiles to HST/WFC3/F160W images from the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS). Massive quiescent galaxies are very compact, with a median circularized half-light radius r_e = 0.63 +/- 0.18 kpc. Removing 5/16 (31%) sources with signs of AGN activity does not change the result. Star-forming galaxies have r_e = 2.0 +/- 0.60 kpc, 3.2 +/- 1.3 x larger than quiescent galaxies. Quiescent galaxies at z~4 are on average 6.0 +\- 0.17 x smaller than at z~0 and 1.9 +/- 0.7 x smaller than at z~2. Star-forming galaxies of the same stellar mass are 2.4 +/- 0.7 x smaller than at z~0. Overall, the size evolution at 0<z<4 is well described by a powerlaw, with r_e = 5.08 +/- 0.28 (1+z)^(-1.44+/-0.08) kpc for quiescent and r_e = 6.02 +/- 0.28 (1+z)^(-0.72+/-0.05) kpc for star-forming galaxies. Compact star-forming galaxies are rare in our sample: we find only 1/14 (7%) with r_e / (M / 10^11 M_Sun)^0.75 < 1.5, whereas 13/16 (81%) of the quiescent galaxies is compact. The number density of compact quiescent galaxies at z~4 is 1.8 +/- 0.8 x 10^-5 Mpc^-3 and increases rapidly, by >5 x, between 2<z<4. The paucity of compact star-forming galaxies at z~4 and their large rest-frame ultra-violet median sizes suggest that the formation phase of compact cores is very short and/or highly dust obscured.

There is ongoing debate regarding the extent that environment affects galaxy size growth beyond z>1. To investigate the differences in star-forming and quiescent galaxy properties as a function of environment at z=2.1, we create a mass-complete sample of 59 cluster galaxies Spitler et al. (2012) and 478 field galaxies with log(M)>9 using photometric redshifts from the ZFOURGE survey. We compare the mass-size relation of field and cluster galaxies using measured galaxy semi-major axis half-light radii ($r_{1/2,maj}$) from CANDELS HST/F160W imaging. We find consistent mass normalized (log(M)=10.7) sizes for quiescent field galaxies ($r_{1/2,maj}=1.81\pm0.29$ kpc) and quiescent cluster galaxies ($r_{1/2,maj}=2.17\pm0.63$ kpc). The mass normalized size of star-forming cluster galaxies ($r_{1/2,maj}=4.00\pm0.26$ kpc ) is 12% larger (KS test $2.1\sigma$) than star-forming field galaxies ($r_{1/2,maj}=3.57\pm0.10$ kpc). From the mass-color relation we find that quiescent field galaxies with 9.7<log(M)<10.4 are slightly redder (KS test $3.6\sigma$) than quiescent cluster galaxies, while cluster and field quiescent galaxies with log(M)>10.4 have consistent colors. We find that star-forming cluster galaxies are on average 20% redder than star-forming field galaxies at all masses. Furthermore, we stack galaxy images to measure average radial color profiles as a function of mass. Negative color gradients are only present for massive star-forming field and cluster galaxies with log(M)>10.4, the remaining galaxy masses and types have flat profiles. Our results suggest given the observed differences in size and color of star-forming field and cluster galaxies, that the environment has begun to influence/accelerate their evolution. However, the lack of differences between field and cluster quiescent galaxies indicates that the environment has not begun to significantly influence their evolution at z~2.

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

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.

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

We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$\sigma$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.

Our understanding of the redshift $z>3$ galaxy population relies largely on samples selected using the popular "dropout" technique, typically consisting of UV-bright galaxies with blue colors and prominent Lyman breaks. As it is currently unknown if these galaxies are representative of the massive galaxy population, we here use the FourStar Galaxy Evolution (ZFOURGE) Survey to create a stellar mass-limited sample at $z=3-4$. Uniquely, ZFOURGE uses deep near-infrared medium-bandwidth filters to derive accurate photometric redshifts and stellar population properties. The mass-complete sample consists of 57 galaxies with log M $>10.6$, reaching below $M^{\star}$ at $z=3-4$. On average, the massive $z=3-4$ galaxies are extremely faint in the observed optical with median $R_{tot}^{AB}=27.48\pm0.41$ (restframe $M_{1700}=-18.05\pm0.37$). They lie far below the UV luminosity-stellar mass relation for Lyman break galaxies and are about $\sim100\times$ fainter at the same mass. The massive galaxies are red ($R-Ks_{AB}=3.9\pm0.2$; restframe UV-slope $\beta=-0.2\pm0.3$) likely from dust or old stellar ages. We classify the galaxy SEDs by their restframe $U-V$ and $V-J$ colors and find a diverse population: $46^{+6+10}_{-6-17}$% of the massive galaxies are quiescent, $54^{+8+17}_{-8-10}$% are dusty star-forming galaxies, and only $14^{+3+10}_{-3-4}$% resemble luminous blue star forming Lyman break galaxies. This study clearly demonstrates an inherent diversity among massive galaxies at higher redshift than previously known. Furthermore,we uncover a reservoir of dusty star-forming galaxies with $4\times$ lower specific star-formation rates compared to submillimeter-selected starbursts at $z>3$. With $5\times$ higher numbers, the dusty galaxies may represent a more typical mode of star formation compared to submillimeter-bright starbursts.

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.

We report the likely identification of a substantial population of massive M~10^11M_Sun galaxies at z~4 with suppressed star formation rates (SFRs), selected on rest-frame optical to near-IR colors from the FourStar Galaxy Evolution Survey. The observed spectral energy distributions show pronounced breaks, sampled by a set of near-IR medium-bandwidth filters, resulting in tightly constrained photometric redshifts. Fitting stellar population models suggests large Balmer/4000AA breaks, relatively old stellar populations, large stellar masses and low SFRs, with a median specific SFR of 2.9+/-1.8 x 10^-11/yr. Ultradeep Herschel/PACS 100micron, 160micron and Spitzer/MIPS 24micron data reveal no dust-obscured SFR activity for 15/19 (79%) galaxies. Two far-IR detected galaxies are obscured QSOs. Stacking the far-IR undetected galaxies yields no detection, consistent with the SED fit, indicating independently that the average specific SFR is at least 10x smaller than of typical star-forming galaxies at z~4. Assuming all far-IR undetected galaxies are indeed quiescent, the volume density is 1.8+/-0.7 x 10^-5Mpc^-3 to a limit of log10M/M_Sun>10.6, which is 10x and 80x lower than at z = 2 and z = 0.1. They comprise a remarkably high fraction (~35%) of z~4 massive galaxies, suggesting that suppression of star formation was efficient even at very high redshift. Given the average stellar age of 0.8Gyr and stellar mass of 0.8x10^11M_Sun, the galaxies likely started forming stars before z =5, with SFRs well in excess of 100M_Sun/yr, far exceeding that of similarly abundant UV-bright galaxies at z>4. This suggests that most of the star-formation in the progenitors of quiescent z~4 galaxies was obscured by dust.

Using observations from the FourStar Galaxy Evolution Survey (ZFOURGE), we obtain the deepest measurements to date of the galaxy stellar mass function at 0.5 < z < 2.5. ZFOURGE provides well-constrained photometric redshifts made possible through deep medium-bandwidth imaging at 1-2um . We combine this with HST imaging from the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS), allowing for the efficient selection of both blue and red galaxies down to stellar masses ~10^9.5 Msol at z ~ 2.5. The total surveyed area is 316 arcmin^2 distributed over three independent fields. We supplement these data with the wider and shallower NEWFIRM Medium-Band Survey (NMBS) to provide stronger constraints at high masses. Several studies at z<=1 have revealed a steepening of the slope at the low-mass end of the stellar mass function (SMF), leading to an upturn at masses <10^10 Msol that is not well-described by a standard single-Schechter function. We find evidence that this feature extends to at least z ~ 2, and that it can be found in both the star-forming and quiescent populations individually. The characteristic mass (M*) and slope at the lowest masses (alpha) of a double-Schechter function fit to the SMF stay roughly constant at Log(M/Msol) ~ 10.65 and ~-1.5 respectively. The SMF of star-forming galaxies has evolved primarily in normalization, while the change in shape is relatively minor. This is not the case for quiescent galaxies: the depth of our imaging allows us to show for the first time significantly more evolution at Log(M/Msol) < 10.5 than at higher masses. We find that the total mass density (down to 10^9 Msol) in star-forming galaxies has increased by a factor of ~2.2 since z ~ 2.5, whereas in quiescent galaxies it has increased by a factor of ~12 .

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.

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

We present spectroscopic measurements of the [OIII] emission line from two subregions of strong Lyman-alpha emission in a radio-quiet Lyman-alpha blob (LAB). The blob under study is LAB1 (Steidel et al. 2000) at z ~ 3.1, and the [OIII] detections are from the two Lyman break galaxies embedded in the blob halo. The [OIII] measurements were made with LUCIFER on the 8.4m Large Binocular Telescope and NIRSPEC on 10m Keck Telescope. Comparing the redshift of the [OIII] measurements to Lyman-alpha redshifts from SAURON (Weijmans et al. 2010) allows us to take a step towards understanding the kinematics of the gas in the blob. Using both LUCIFER and NIRSPEC we find velocity offsets between the [OIII] and Lyman-alpha redshifts that are modestly negative or consistent with 0 km/s in both subregions studied (ranging from -72 +/- 42 -- +6 +/- 33 km/s). A negative offset means Lyman-alpha is blueshifted with respect to [OIII], a positive offset then implies Lyman-alpha is redshifted with respect to [OIII]. These results may imply that outflows are not primarily responsible for Lyman alpha escape in this LAB, since outflows are generally expected to produce a positive velocity offset (McLinden et al. 2011). In addition, we present an [OIII] line flux upper limit on a third region of LAB1, a region that is unassociated with any underlying galaxy. We find that the [OIII] upper limit from the galaxy-unassociated region of the blob is at least 1.4 -- 2.5 times fainter than the [OIII] flux from one of the LBG-associated regions and has an [OIII] to Lyman-alpha ratio measured at least 1.9 -- 3.4 times smaller than the same ratio measured from one of the LBGs.

We report the first results from the Z-FOURGE survey: the discovery of a candidate galaxy cluster at z = 2.2 consisting of two compact overdensities with red galaxies detected at >20{\sigma} above the mean surface density. The discovery was made possible by a new deep (Ks < 24.8 AB 5{\sigma}) Magellan/FourStar near-IR imaging survey with 5 custom medium-bandwidth filters. The filters pinpoint the location of the Balmer/4000\AA\ break in evolved stellar populations at 1.5 < z < 3.5, yielding significantly more accurate photometric redshifts than possible with broadband imaging alone. The overdensities are within 1' of each other in the COSMOS field and appear to be embedded in a larger structure that contains at least one additional overdensity (\sim10{\sigma}). Considering the global properties of the overdensities, the z = 2.2 system appears to be the most distant example of a galaxy cluster with a population of red galaxies. A comparison to a large {\Lambda}CDM simulation suggests that the system may consist of merging subclusters, with properties in between those of z > 2 protoclusters with more diffuse distributions of blue galaxies and the lower-redshift galaxy clusters with prominent red sequences. The structure is completely absent in public optical catalogs in COSMOS and only weakly visible in a shallower near-IR survey. The discovery showcases the potential of deep near-IR surveys with medium-band filters to advance the understanding of environment and galaxy evolution at z > 1.5.

The study of Ly-alpha emission in the high-redshift universe is a useful probe of the epoch of reionization, as the Ly-alpha line should be attenuated by the intergalactic medium (IGM) at low to moderate neutral hydrogen fractions. Here we present the results of a deep and wide imaging search for Ly-alpha emitters in the COSMOS field. We have used two ultra-narrowband filters (filter width of ~8-9 {\deg}A) on the NEWFIRM camera, installed on the Mayall 4m telescope at Kitt Peak National Observatory, in order to isolate Ly-alpha emitters at z = 7.7; such ultra-narrowband imaging searches have proved to be excellent at detecting Ly-alpha emitters. We found 5-sigma detections of four candidate Ly-alpha emitters in a survey volume of 2.8 x 10^4 Mpc^3 (total survey area ~760 arcmin^2). Each candidate has a line flux greater than 8 x 10^-18 erg s^-1 cm^-2. Using these results to construct a luminosity function and comparing to previously established Ly-alpha luminosity functions at z = 5.7 and z = 6.5, we find no conclusive evidence for evolution of the luminosity function between z = 5.7 and z = 7.7. Statistical Monte Carlo simulations suggest that half of these candidates are real z = 7.7 targets, and spectroscopic follow-up will be required to verify the redshift of these candidates. However, our results are consistent with no strong evolution in the neutral hydrogen fraction of the IGM between z = 5.7 and z = 7.7, even if only one or two of the z = 7.7 candidates are spectroscopically confirmed.

Rapid mass assembly, likely from mergers or smooth accretion, has been predicted to play a vital role in star-formation in high-redshift Lyman-alpha (Lya) emitters. Here we predict the major merger, minor merger, and smooth accreting Lya emitter fraction from z~3 to z~7 using a large dark matter simulation, and a simple physical model that is successful in reproducing many observations over this large redshift range. The central tenet of this model, different from many of the earlier models, is that the star-formation in Lya emitters is proportional to the mass accretion rate rather than the total halo mass. We find that at z~3, nearly 35% of the Lya emitters accrete their mass through major (3:1) mergers, and this fraction increases to about 50% at z~7. This imply that the star-formation in a large fraction of high-redshift Lya emitters is driven by mergers. While there is discrepancy between the model predictions and observed merger fractions, some of this difference (~15%) can be attributed to the mass-ratio used to define a merger in the simulation. We predict that future, deeper observations which use a 3:1 definition of major mergers will find >30% major merger fraction of Lya emitters at redshifts >3.

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.

We present a simple physical model for populating dark matter halos with Lyman Alpha Emiiters(LAEs) and predict the physical properties of LAEs at z~3-7. The central tenet of this model is that the Ly-alpha luminosity is proportional to the star formation rate (SFR) which is directly related to the halo mass accretion rate. The only free parameter in our model is then the star-formation efficiency (SFE). An efficiency of 2.5% provides the best-fit to the Ly-alpha luminosity function (LF) at redshift z=3.1, and we use this SFE to construct Ly-alpha LFs at other redshifts. Our model reproduce the Ly-alpha LFs, stellar ages, SFR ~1-10; Msun/yr, stellar masses ~ 10^7-10^8 Msun and the clustering properties of LAEs at z~3-7. We find the spatial correlation lengths ro ~ 3-6 Mpc/h, in agreement with the observations. Finally, we estimate the field-to-field variation ~ 30% for current volume and flux limited surveys, again consistent with observations. Our results suggest that the star formation, and hence Ly-alpha emission in LAEs is powered by the accretion of new material, and that the physical properties of LAEs do not evolve significantly over a wide range of redshifts. Relating the accreted mass, rather than the total mass of halos, to the Ly-alpha luminosity of LAEs naturally gives rise to the duty cycle of LAEs.