We derive the UV-optical stellar dust attenuation curve of galaxies at z=1.4-2.6 as a function of gas-phase metallicity. We use a sample of 218 star-forming galaxies, excluding those with very young or heavily obscured star formation, from the MOSFIRE Deep Evolution Field (MOSDEF) survey with H$\alpha$, H$\beta$, and [NII]$\lambda 6585$ spectroscopic measurements. We constrain the shape of the attenuation curve by comparing the average flux densities of galaxies sorted into bins of dust obscuration using Balmer decrements, i.e., H$\alpha$-to-H$\beta$ luminosities. The average attenuation curve for the high-metallicity sample (12+log(O/H)>8.5, corresponding to $M_*\gtrsim10^{10.4}\,M_{\odot}$) has a shallow slope, identical to that of the Calzetti local starburst curve, and a significant UV 2175A extinction bump that is $\sim 0.5\times$ the strength of the Milky Way bump. On the other hand, the average attenuation curve of the low-metallicity sample (12+log(O/H) $\sim 8.2-8.5$) has a steeper slope similar to that of the SMC curve, only consistent with the Calzetti slope at the $3\sigma$ level. The UV bump is not detected in the low-metallicity curve, indicating the relative lack of the small dust grains causing the bump at low metallicities. Furthermore, we find that on average the nebular reddening (E(B-V)) is a factor of 2 times larger than that of the stellar continuum for galaxies with low metallicities, while the nebular and stellar reddening are similar for galaxies with higher metallicities. The latter is likely due to a high surface density of dusty clouds embedding the star forming regions but also reddening the continuum in the high-metallicity galaxies.

ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of $M_V \approx -20$ mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of $\sim6.0$ mag (100d)$^{-1}$. Owing to the weakening of HI and the appearance of HeI in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution show significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesised $^{56}$Ni mass of $M_{\rm Ni} \sim0.4\rm M_\odot$ and ejecta with high kinetic energy $E_{\rm kin}= (5-9)\times10^{51} $ erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim0.3\, M_\odot$ and $E_{\rm kin}= 3\times10^{51} $ erg. The high $^{56}$Ni mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high $^{56}$Ni yields. The earliest spectrum shows "flash ionization" features, from which we estimate a mass-loss rate of $ \dot{M}\approx 2\times10^{-4} \, \rm M_\odot~yr^{-1} $. This wind density is too low to power the luminous light curve by ejecta-CSM interaction. We measure expansion velocities as high as $ 17,000 \rm~km~s^{-1}$ for H$\alpha$, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of 1.7-3.1 $\rm M_\odot$ using the [OI] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main sequence mass of 19-24 $\rm M_\odot$.

We present new measurements of the spatial distribution and kinematics of neutral hydrogen in the circumgalactic and intergalactic medium surrounding star-forming galaxies at z ~ 2. Using the spectra of ~ 3000 galaxies with redshifts <z> +/- 0.4 from the Keck Baryonic Structure Survey (KBSS), we assemble a sample of more than 200,000 distinct foreground-background pairs with projected angular separations of 3 - 500 arcsec and spectroscopic redshifts, with <$z_{fg}$> = 2.23 and <$z_{bg}$> = 2.57. The ensemble of sightlines and foreground galaxies is used to construct a 2D map of the mean excess Ly$\alpha$ optical depth relative to the intergalactic mean as a function of projected galactocentric distance (20 < $D_{tran}$/pkpc < 4000) and line-of-sight velocity. We provide information on the line-of-sight kinematics of H I gas as a function of projected distance $D_{tran}$. We compare the map with cosmological zoom-in simulation, finding qualitative agreement between them. A simple two-component (accretion, outflow) analytical model generally reproduces the observed line-of-sight kinematics and projected spatial distribution of H I. The best-fitting model suggests that galaxy-scale outflows with initial velocity $v_{out}$ ~ 600 km/s dominate the kinematics of circumgalactic H I out to $D_{tran}$ ~ 50 kpc, while H I at $D_{tran}$ > 100 kpc is dominated by infall with characteristic $v_{in}$ < $v_c$, where $v_c$ is the circular velocity of the host halo ($M_h$ ~ $10^{12} M_\odot$). Over the impact parameter range 80 < $D_{tran}$/pkpc < 200, the H I line-of-sight velocity range reaches a minimum, with a corresponding flattening in the rest-frame Ly$\alpha$ equivalent width. These observations can be naturally explained as the transition between outflow-dominated and accretion-dominated flows. Beyond $D_{tran}$ ~ 300 kpc, the line of sight kinematics are dominated by Hubble expansion.

We present a joint analysis of rest-UV and rest-optical spectra obtained using Keck/LRIS and Keck/MOSFIRE for a sample of 62 star-forming galaxies at $z\sim2.3$. We divide our sample into 2 bins based on their location in the [OIII]/Hb vs. [NII]/Ha BPT diagram, and perform the first differential study of the rest-UV properties of massive ionizing stars as a function of rest-optical emission-line ratios. Fitting BPASS stellar population synthesis models, including nebular continuum emission, to our rest-UV spectra, we find that high-redshift galaxies offset towards higher [OIII]/Hb and [NII]/Ha have younger ages (log(Age/yr)=$7.20^{+0.57}_{-0.20}$) and lower stellar metallicities ($Z_*=0.0010^{+0.0011}_{-0.0003}$) resulting in a harder ionizing spectrum, compared to the galaxies in our sample that lie on the local BPT star-forming sequence (log(Age/yr)=$8.57^{+0.88}_{-0.84}$, $Z_*=0.0019\pm0.0006$). Additionally, we find that the offset galaxies have an ionization parameter of log(U)=$-3.04^{+0.06}_{-0.11}$ and nebular metallicity of 12+log(O/H)=$8.40^{+0.06}_{-0.07}$, and the non-offset galaxies have log(U)=$-3.11\pm0.08$ and 12+log(O/H)=$8.30^{+0.05}_{-0.06}$. The stellar and nebular metallicities derived for our sample imply that the galaxies offset from the local BPT relation are more alpha-enhanced ($7.28^{+2.52}_{-2.82}$ O/Fe$_\odot$) compared to those consistent with the local sequence ($3.04^{+0.95}_{-0.54}$ O/Fe$_\odot$). However, even galaxies that are entirely consistent with the local nebular excitation sequence appear to be alpha-enhanced -- in contrast with typical local systems. Such differences must be considered when estimating gas-phase O/H at high redshift based on strong emission-line ratios. Specifically, a similarity in the location of high-redshift and local galaxies in the BPT diagram may not be indicative of a similarity in their physical properties.

We present ionized gas kinematics for 681 galaxies at $z\sim 1.4-3.8$ from the MOSFIRE Deep Evolution Field survey, measured using models which account for random galaxy-slit misalignments together with structural parameters derived from CANDELS Hubble Space Telescope (HST) imaging. Kinematics and sizes are used to derive dynamical masses. Baryonic masses are estimated from stellar masses and inferred gas masses from dust-corrected star formation rates (SFRs) and the Kennicutt-Schmidt relation. We measure resolved rotation for 105 galaxies. For the remaining 576 galaxies we use models based on HST imaging structural parameters together with integrated velocity dispersions and baryonic masses to statistically constrain the median ratio of intrinsic ordered to disordered motion, $V/\sigma_{V,0}$. We find that $V/\sigma_{V,0}$ increases with increasing stellar mass and decreasing specific SFR (sSFR). These trends may reflect marginal disk stability, where systems with higher gas fractions have thicker disks. For galaxies with detected rotation we assess trends between their kinematics and mass, sSFR, and baryon surface density ($\Sigma_{\mathrm{bar},e}$). Intrinsic dispersion correlates most with $\Sigma_{\mathrm{bar},e}$ and velocity correlates most with mass. By comparing dynamical and baryonic masses, we find that galaxies at $z\sim 1.4-3.8$ are baryon dominated within their effective radii ($R_E$), with Mdyn/Mbaryon increasing over time. The inferred baryon fractions within $R_E$, $f_{\mathrm{bar}}$, decrease over time, even at fixed mass, size, or surface density. At fixed redshift, $f_{\mathrm{bar}}$ does not appear to vary with stellar mass but increases with decreasing $R_E$ and increasing $\Sigma_{\mathrm{bar},e}$. For galaxies at $z\geq2$, the median inferred baryon fractions generally exceed 100%. We discuss possible explanations and future avenues to resolve this tension.

We perform a comprehensive analysis of the redshift evolution of the rest-UV spectra of star-forming galaxies out to z~5. We combine new z~5 measurements of HI Ly$\alpha$ and low- and high-ionization interstellar metal absorption features with comparable measurements at z~2-4. We measure the equivalent widths of interstellar absorption features using stacked spectra in bins of Ly$\alpha$ equivalent width, performing corrections to Ly$\alpha$ strengths based on a model for the transmission of the intergalactic medium. We find a strong correlation between decreasing low-ionization absorption strength and increasing Ly$\alpha$ emission strength over the redshift range z~2-5, suggesting that both of these quantities are fundamentally linked to neutral gas covering fraction. At the highest Ly$\alpha$ equivalent widths, we observe evolution at $z\sim5$ towards greater Ly$\alpha$ emission strength at fixed low-ionization absorption strength. If we interpret the non-evolving relationship of Ly$\alpha$ emission strength and low-ionization line strength at z~2-4 as primarily reflecting the radiative transfer of Ly$\alpha$ photons, this evolution at z~5 suggests a higher intrinsic production rate of Ly$\alpha$ photons than at lower redshift. Our conclusion is supported by the joint evolution of the relationships among Ly$\alpha$ emission strength, interstellar absorption strength, and dust reddening. We perform additional analysis in bins of stellar mass, star-formation rate, UV luminosity, and age, examining how the relationships between galaxy properties and Ly$\alpha$ emission evolve towards higher redshift. We conclude that increasing intrinsic Ly$\alpha$ photon production and strong detection of nebular CIV emission (signaling lower metallicity) at z~5 indicate an elevated ionized photon production efficiency ($\xi_{\rm ion}$).

We present measurements of [SIII]$\lambda\lambda$9069,9531 for a sample of $z\sim1.5$ star-forming galaxies, the first sample with measurements of these lines at z>0.1. We employ the line ratio S$_{32}$$\equiv$[SIII]$\lambda\lambda$9069,9531/[SII]$\lambda\lambda$6716,6731 as a novel probe of evolving ISM conditions. Since this ratio includes the low-ionization line [SII], it is crucial that the effects of diffuse ionized gas (DIG) on emission-line ratios be accounted for in $z\sim0$ integrated galaxy spectra, or else that comparisons be made to samples of local HII regions in which DIG emission is not present. We find that S$_{32}$ decreases with increasing stellar mass at both $z\sim1.5$ and $z\sim0$, but that the dependence is weak suggesting S$_{32}$ has a very shallow anticorrelation with metallicity, in contrast with O$_{32}$ that displays a strong metallicity dependence. As a result, S$_{32}$ only mildly evolves with redshift at fixed stellar mass. The $z\sim1.5$ sample is systematicallty offset towards lower S$_{32}$ and higher [SII]/H$\alpha$ at fixed [OIII]/H$\beta$ relative to $z=0$ HII regions. By comparing to photoionization model grids, we find that such trends can be explained by a scenario in which the ionizing spectrum is harder at fixed O/H with increasing redshift, but are inconsistent with an increase in ionization parameter at fixed O/H. This analysis demonstrates the advantages of expanding beyond the strongest rest-optical lines for evolutionary studies, and the particular utility of [SIII] for characterizing evolving ISM conditions and stellar compositions. These measurements provide a basis for estimating [SIII] line strengths for high-redshift galaxies, a line that the James Webb Space Telescope will measure out to z~5.5.

We present detections of [OIII]$\lambda$4363 and direct-method metallicities for star-forming galaxies at $z=1.7-3.6$. We combine new measurements from the MOSFIRE Deep Evolution Field (MOSDEF) survey with literature sources to construct a sample of 18 galaxies with direct-method metallicities at $z>1$, spanning $7.5<1$2+log(O/H$)<8.2$ and log(M$_*$/M$_{\odot})=7-10$. We find that strong-line calibrations based on local analogs of high-redshift galaxies reliably reproduce the metallicity of the $z>1$ sample on average. We construct the first mass-metallicity relation at $z>1$ based purely on direct-method O/H, finding a slope that is consistent with strong-line results. Direct-method O/H evolves by $\lesssim0.1$ dex at fixed M$_*$ and SFR from $z\sim0-2.2$. We employ photoionization models to constrain the ionization parameter and ionizing spectrum in the high-redshift sample. Stellar models with super-solar O/Fe and binary evolution of massive stars are required to reproduce the observed strong-line ratios. We find that the $z>1$ sample falls on the $z\sim0$ relation between ionization parameter and O/H, suggesting no evolution of this relation from $z\sim0$ to $z\sim2$. These results suggest that the offset of the strong-line ratios of this sample from local excitation sequences is driven primarily by a harder ionizing spectrum at fixed nebular metallicity compared to what is typical at $z\sim0$, naturally explained by super-solar O/Fe at high redshift caused by rapid formation timescales. Given the extreme nature of our $z>1$ sample, the implications for representative $z\sim2$ galaxy samples at $\sim10^{10}$ M$_{\odot}$ are unclear, but similarities to $z>6$ galaxies suggest that these conclusions can be extended to galaxies in the epoch of reionization.

We present the discovery of PS18kh, a tidal disruption event (TDE) discovered at the center of SDSS J075654.53+341543.6 ($d\simeq322$ Mpc) by the Pan-STARRS Survey for Transients. Our dataset includes pre-discovery survey data from Pan-STARRS, the All-Sky Automated Survey for Supernovae (ASAS-SN), and the Asteroid Terrestrial-impact Last Alert System (ATLAS) as well as high-cadence, multi-wavelength follow-up data from ground-based telescopes and Swift, spanning from 56 days before peak light until 75 days after. The optical/UV emission from PS18kh is well-fit as a blackbody with temperatures ranging from $T\simeq12000$ K to $T\simeq25000$ K and it peaked at a luminosity of $L\simeq8.8\times10^{43}$ ergs s$^{-1}$. PS18kh radiated $E=(3.45\pm0.22)\times10^{50}$ ergs over the period of observation, with $(1.42\pm0.20)\times10^{50}$ ergs being released during the rise to peak. Spectra of PS18kh show a changing, boxy/double-peaked H$\alpha$ emission feature, which becomes more prominent over time. We use models of non-axisymmetric accretion disks to describe the profile of the H$\alpha$ line and its evolution. We find that at early times the high accretion rate leads the disk to emit a wind which modifies the shape of the line profile and makes it bell-shaped. At late times, the wind becomes optically thin, allowing the non-axisymmetric perturbations to show up in the line profile. The line-emitting portion of the disk extends from $r_{\rm in}\sim60r_{\rm g}$ to an outer radius of $r_{\rm out}\sim1400r_{\rm g}$ and the perturbations can be represented either as an eccentricity in the outer rings of the disk or as a spiral arm in the inner disk.

Using data from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we present a census of AGN-driven ionized outflows in a sample of 159 AGNs at $1.4 \le z \le 3.8$. The sample spans AGN bolometric luminosities of $10^{44-47} \mathrm{~erg~s}^{-1}$ and includes both quiescent and star-forming galaxies extending across three orders of magnitude in stellar mass. We identify and characterize outflows from the \hbeta, [OIII], \halpha ~and [NII] emission line spectra. We detect outflows in $17\%$ of the AGNs, seven times more often than in a mass-matched sample of inactive galaxies in MOSDEF. The outflows are fast and galaxy-wide, with velocities of $\sim 400-3500 ~\mathrm{km~s}^{-1}$ and spatial extents of $0.3-11.0$ kpc. The incidence of outflows among AGNs is independent of the stellar mass of the host galaxy, with outflows detected in both star-forming and quiescent galaxies. This suggests that outflows exist across different phases in galaxy evolution. We investigate relations between outflow kinematic, spatial, and energetic properties and both AGN and host galaxy properties. Our results show that AGN-driven outflows are widespread in galaxies along the star-forming main sequence. The mass-loading factors of the outflows are typically $0.1-1$ and increase with AGN luminosity, capable of exceeding unity at $L_\mathrm{AGN} \gtrsim 10^{46.3} \mathrm{~erg~s}^{-1}$. In these more luminous sources the ionized outflow alone is likely sufficient to regulate star formation, and when combined with outflowing neutral and molecular gas may be able to quench star formation in their host galaxies.

Population synthesis models predict that high-mass X-ray binary (HMXB) populations produced in low metallicity environments should be more X-ray luminous, a trend supported by studies of nearby galaxies. This trend may be responsible for the observed increase of the X-ray luminosity ($L_{\mathrm{X}}$) per star formation rate (SFR) with redshift due to the decrease of metallicity ($Z$) at fixed stellar mass as a function of redshift. To test this hypothesis, we use a sample of 79 $z\sim2$ star-forming galaxies with oxygen abundance measurements from the MOSDEF survey, which obtained rest-frame optical spectra for $\sim1500$ galaxies in the CANDELS fields at $1.37<z<3.80$. Using Chandra data from the AEGIS-X Deep, Deep Field North, and Deep Field South surveys, we stack the X-ray data at the galaxy locations in bins of redshift and $Z$ because the galaxies are too faint to be individually detected. In agreement with previous studies, the average $L_{\mathrm{X}}$/SFR of our $z\sim2$ galaxy sample is enhanced by $\approx0.4-0.8$ dex relative to local HMXB $L_{\mathrm{X}}$-SFR scaling relations. Splitting our sample by $Z$, we find that $L_{\mathrm{X}}$/SFR and $Z$ are anti-correlated with 97% confidence. This observed $Z$ dependence for HMXB-dominated galaxies is consistent both with the local $L_{\mathrm{X}}$-SFR-$Z$ relation and a subset of population synthesis models. Although the statistical significance of the observed trends is weak due to the low X-ray statistics, these results constitute the first direct evidence connecting the redshift evolution of $L_{\mathrm{X}}$/SFR and the $Z$ dependence of HMXBs.

We present results on the emission-line properties of 1.3<=z<=2.7 galaxies drawn from the complete MOSFIRE Deep Evolution Field (MOSDEF) survey. Specifically, we use observations of the emission-line diagnostic diagram of [OIII]5007/Hb vs. [SII]6717,6731/Ha, i.e., the "[SII] BPT diagram," to gain insight into the physical properties of high-redshift star-forming regions. High-redshift MOSDEF galaxies are offset towards lower [SII]6717,6731/Ha at fixed [OIII]5007/Hb, relative to local galaxies from the Sloan Digital Sky Survey (SDSS). Furthermore, at fixed [OIII]5007/Hb, local SDSS galaxies follow a trend of decreasing [SII]6717,6731/Ha as the surface density of star formation (Sigma_SFR) increases. We explain this trend in terms of the decreasing fractional contribution from diffuse ionized gas (f_DIG) as Sigma_SFR increases in galaxies, which causes galaxy-integrated line ratios to shift towards the locus of pure HII-region emission. The z~0 relationship between f_DIG and Sigma_SFR implies that high-redshift galaxies have lower f_DIG values than typical local systems, given their significantly higher typical Sigma_SFR. When an appropriate low-redshift benchmark with zero or minimal f_DIG is used, high-redshift MOSDEF galaxies appear offset towards higher [SII]6717,6731/Ha and/or [OIII]5007/Hb. The joint shifts of high-redshift galaxies in the [SII] and [NII] BPT diagrams are best explained in terms of the harder spectra ionizing their star-forming regions at fixed nebular oxygen abundance (expected for chemically-young galaxies), as opposed to large variations in N/O ratios or higher ionization parameters. The evolving mixture of HII regions and DIG is an essential ingredient to our description of the ISM over cosmic time.

We present James Clerk Maxwell Telescope Submillimetre Common-User Bolometer Array 2 (SCUBA-2) 850 & 450 $\mu$m observations ($\sigma_{850}\sim0.5$ mJy, $\sigma_{450}\sim5$ mJy) of the HS1549+19 and HS1700+64 survey fields containing two of the largest known galaxy over-densities at $z=2.85$ and $2.30$, respectively. We detect 56 sub-millimetre galaxies (SMGs) with SNR $> 4$ over $\sim 50$ arcmin$^2$ at 850 $\mu$m with flux densities of 3 - 17 mJy. The number counts indicate over-densities in the $3$-arcmin diameter core region ($\sim 1.5$ Mpc at $z=2.5$) of $6^{+4}_{-2}\times$ (HS1549) and $4^{+6}_{-2}\times$ (HS1700) compared to blank field surveys. Within these core regions, we spectroscopically confirm that approximately one third of the SMGs lie at the protocluster redshifts for both HS1549 and HS1700. We use statistical identifications of other SMGs in the wider fields to constrain an additional four candidate protocluster members in each system. We combine multi wavelength estimates of the star-formation rates (SFRs) from Lyman-break dropout- and narrowband-selected galaxies, and the SCUBA-2 SMGs, to estimate total SFRs of 12,$500\pm2800$ M$_\odot$ yr$^{-1}$ ($4900\pm1200$ M$_\odot$ yr$^{-1}$) in HS1549 (HS1700), and SFR densities (SFRDs) within the central 1.5-Mpc diameter of each protocluster to be $3000\pm900$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ ($1300\pm400$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$) in the HS1549 (HS1700) protocluster, $\sim10^4\times$ larger than the global SFRDs found at their respective epochs, due to the concentration of star-forming galaxies in the small volume of the dense cluster cores. Our results suggest centrally concentrated starbursts within protoclusters may be a relatively common scenario for the build up of mass in rich clusters assembling at $z\gtrsim2$.

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.

We present results from the Keck Baryonic Structure Survey (KBSS) including the first detailed measurements of the column densities, kinematics, and internal energy of metal-bearing gas within the virial radius (35-100 physical kpc) of eight ~$L^*$ galaxies at $z\sim2$. From our full sample of 130 metal-bearing absorbers, we infer that halo gas is kinematically complex when viewed in singly, doubly, and triply ionized species. Broad OVI and CIV absorbers are detected at similar velocities to the lower-ionization gas but with very different kinematic structure indicating that the circumgalactic medium (CGM) is multi-phase. There is a high covering fraction of metal-bearing gas within 100 kpc including highly ionized gas such as OVI; however, observations of a single galaxy probed by a lensed background QSO suggest the size of metal-bearing clouds is small (<400 pc for all but the OVI-bearing gas). The mass in metals found within the halo is substantial, equivalent to $\gtrsim$ 25% of the metal mass within the interstellar medium. The gas kinematics unambiguously show that 70% of galaxies with detected metal absorption have some unbound metal-enriched gas, suggesting galactic winds may commonly eject gas from halos at $z\sim2$. Significant thermal broadening is detected in CGM absorbers which dominates the internal energy of the gas. 40% of the detected gas has temperatures in the range $10^{4.5-5.5}$ K where cooling times are short, suggesting the CGM is dynamic, with constant heating and/or cooling to produce this short-lived thermal phase.

We study the properties of 30 spectroscopically-identified pairs of galaxies observed during the peak epoch of star formation in the universe. These systems are drawn from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 \leq z \leq 3.8$, and are interpreted as early-stage galaxy mergers. Galaxy pairs in our sample are identified as two objects whose spectra were collected on the same Keck/MOSFIRE spectroscopic slit. Accordingly, all pairs in the sample have projected separations $R_{{\rm proj}}\leq 60$ kpc. The velocity separation for pairs was required to be $\Delta v \leq 500 \mbox{ km s}^{-1}$, which is a standard threshold for defining interacting galaxy pairs at low redshift. Stellar mass ratios in our sample range from 1.1 to 550, with 12 ratios closer than or equal to 3:1, the common definition of a "major merger." Studies of merging pairs in the local universe indicate anenhancement in star-formation activity and deficit in gas-phase oxygen abundance relative to isolated galaxies of the same mass. We compare the MOSDEF pairs sample to a control sample of isolated galaxies at the same redshift, finding no measurable SFR enhancement or metallicity deficit at fixed stellar mass for the pairs sample. The lack of significant difference between the average properties of pairs and control samples appears in contrast to results from low-redshift studies, although the small sample size and lower signal-to-noise of the high-redshift data limit definitive conclusions on redshift evolution. These results are consistent with some theoretical works suggesting a reduced differential effect of pre-coalescence mergers on galaxy properties at high redshift -- specifically that pre-coalescence mergers do not drive strong starbursts.

We present a detailed census of galaxies in and around PC217.96+32.3, a spectroscopically confirmed Coma analog at z=3.78. Diverse galaxy types identified in the field include Lya emitters (LAEs), massive star-forming galaxies, and ultra-massive galaxies (log (Mstar/Msun)>= 11) which may have already halted their star formation. The sky distribution of the star-forming galaxies suggests the presence of a significant overdensity (delta_g=8+/-2), which is spatially offset from the previously confirmed members by 3-4 Mpc to the west. Candidate quiescent and post-starburst galaxies are also found in large excess (a factor of ~ 8-15 higher surface density than the field) although their redshifts are less certain. We estimate that the total enclosed mass traced by star-forming galaxy candidates is roughly comparable to that of PC217.96+32.3 traced by the LAEs. We speculate that the true extent of PC217.96+32.3 may be larger than previously known, a half of which is missed by our LAE selection. Alternatively, the newly discovered overdensity may belong to another Coma progenitor not associated with PC217.96+32.3. Expectations from theory suggest that both scenarios are equally unlikely (<1%), particularly in the cosmic volume probed in our survey. If confirmed as a single structure, its total mass will be well in excess of Coma, making this a singularly large cosmic structure rarely seen even in large cosmological simulations. Finally, we find that the protocluster galaxies follow the same SFR-M_star scaling relation as the field galaxies, suggesting that the environmental effect at z~4 is a subtle one at best for normal star-forming galaxies.

We use extensive spectroscopy from the MOSFIRE Deep Evolution Field (MOSDEF) survey to investigate the relationships between rest-frame optical emission line equivalent widths ($W$) and a number of galaxy and ISM characteristics for a sample of $1134$ star-forming galaxies at redshifts $1.4\lesssim z\lesssim 3.8$. We examine how the equivalent widths of [OII]$\lambda\lambda 3727, 3730$, H$\beta$, [OIII]$\lambda\lambda 4960, 5008$, [OIII]$+$H$\beta$, H$\alpha$, and H$\alpha$+[NII]$\lambda\lambda 6550, 6585$, depend on stellar mass, UV slope, age, star-formation rate (SFR) and specific SFR (sSFR), ionization parameter and excitation conditions (O32 and [OIII]/H$\beta$), gas-phase metallicity, and ionizing photon production efficiency ($\xi_{\rm ion}$). The trend of increasing $W$ with decreasing stellar mass is strongest for [OIII] (and [OIII]+H$\beta$). More generally, the equivalent widths of all the lines increase with redshift at a fixed stellar mass or fixed gas-phase metallicity, suggesting that high equivalent width galaxies are common at high redshift. This redshift evolution in equivalent widths can be explained by the increase in SFR and decrease in metallicity with redshift at a fixed stellar mass. Consequently, the dependence of $W$ on sSFR is largely invariant with redshift, particularly when examined for galaxies of a given metallicity. Our results show that high equivalent width galaxies, specifically those with high $W({\rm [OIII]})$, have low stellar masses, blue UV slopes, young ages, high sSFRs, ISM line ratios indicative of high ionization parameters, high $\xi_{\rm ion}$, and low metallicities. As these characteristics are often attributed to galaxies with high ionizing escape fractions, galaxies with high $W$ are likely candidates for the population that dominates cosmic reionization.

We present results of a deep spectroscopic survey designed to quantify the statistics of the escape of ionizing photons from star-forming galaxies at z~3. We measure the ratio of ionizing to non-ionizing UV flux density <f900/f1500>_obs, where f900 is the mean flux density evaluated over the range [880,910] A. We quantify the emergent ratio of ionizing to non-ionizing UV flux density by analyzing high-S/N composite spectra formed from sub-samples with common observed properties and numbers sufficient to reduce the statistical uncertainty in the modeled IGM+CGM correction to obtain precise values of <f900/f1500>_out, including a full-sample average <f900/f1500>_out=$0.057\pm0.006$. We further show that <f900/f1500>_out increases monotonically with Ly$\alpha$ rest equivalent width, inducing an inverse correlation with UV luminosity as a by-product. We fit the composite spectra using stellar spectral synthesis together with models of the ISM in which a fraction f_c of the stellar continuum is covered by gas with column density N(HI). We show that the composite spectra simultaneously constrain the intrinsic properties of the stars (L900/L1500)_int along with f_c, N(HI), E(B-V), and $f_{esc,abs}$, the absolute escape fraction of ionizing photons. We find a sample-averaged $f_{esc,abs} =0.09\pm0.01$, and that subsamples fall along a linear relation $\langle f_{esc,abs}\rangle \sim 0.75[W(Ly\alpha)/110 A]$. We use the FUV luminosity function, the distribution function $n[W(Ly\alpha)]$, and the relationship between $W(Ly\alpha)$ and <f900/f1500>_out to estimate the total ionizing emissivity of $z\sim3$ star-forming galaxies with Muv < -19.5: $\epsilon_{LyC}\sim 6\times10^{24}$ ergs/s/Hz/Mpc$^3$, exceeding the contribution of QSOs by a factor of $\sim 3$, and accounting for $\sim50$% of the total $\epsilon_{LyC}$ at $z\sim3$ estimated using indirect methods.

Using the MOSDEF rest-frame optical spectroscopic survey, we investigate the star-formation histories (SFHs) of different galaxy types, ranging from actively star forming to quiescent at $1.4\leq~z\leq2.6$. SFHs are constrained utilizing stellar continuum spectroscopy, specifically through a combination of Balmer absorption lines, the 4000$\AA$ break, and the equivalent width of the H$\alpha$ emission line. To attain a sufficiently high signal-to-noise ratio (S/N) to conduct these measurements we stack spectra of galaxies with similar spectral types, as determined from their rest-frame $U-V$ and $V-J$ colors. We bin the MOSDEF sample into five spectral types, subdividing the quiescent and star-forming bins to better explore galaxies transitioning between the two. We constrain the average SFHs for each type, finding that quiescent and transitional galaxies in the MOSDEF sample are dominated by an SFH with an average star-formation timescale of $\tau\sim0.1-0.2$~Gyr. These findings contrast with measurements from the low-redshift Universe where, on average, galaxies form their stars over a more extended time period ($\tau>1$~Gyr). Furthermore, our spectral index measurements correlate with mass surface density for all spectral types. Finally, we compare the average properties of the galaxies in our transitional bins to investigate possible paths to quiescence, and speculate on the viability of a dusty post-starburst phase.

We present a detailed analysis of 317 $2.0 \leq z \leq 2.7$ star-forming galaxies from the Keck Baryonic Structure Survey (KBSS). Using complementary spectroscopic observations with Keck/LRIS and Keck/MOSFIRE, as well as spectral energy distribution (SED) fits to broadband photometry, we examine the joint rest-UV and rest-optical properties of the same galaxies, including stellar and nebular dust attenuation, metallicity, and star formation rate (SFR). The inferred parameters of the stellar population (reddening, age, SFR, and stellar mass) are strongly dependent on the details of the assumed stellar population model and the shape of the attenuation curve. Nebular reddening is generally larger than continuum reddening, but with large scatter. Compared to local galaxies, high-redshift galaxies have lower gas-phase metallicities (and/or higher nebular excitation) at fixed nebular reddening, and higher nebular reddening at fixed stellar mass, consistent with gas fractions that increase with redshift. We find that continuum reddening is correlated with $12+\log(\mathrm{O}/\mathrm{H})_{\mathrm{O3N2}}$ at $3.0\sigma$ significance, whereas nebular reddening is correlated with only $1.1\sigma$ significance. This may reflect the dependence of both continuum reddening and O3N2 on the shape of the ionizing radiation field produced by the massive stars. Finally, we show that H$\alpha$-based and SED-based estimates of SFR exhibit significant scatter relative to one another, and agree on average only for particular combinations of spectral synthesis models and attenuation curves. We find that the SMC extinction curve predicts consistent SFRs if we assume the sub-solar ($0.14Z_{\odot}$) binary star models that are favored for high-redshift galaxies.

We report the discovery and spectroscopic confirmation of a very large star-forming Lyman Break galaxy, G6025, at z_spec=3.721+/-0.003. In the rest-frame ~2100A, G6025 subtends ~24 kpc in physical extent when measured from the 1.5-sigma isophote, in agreement with the parametric size measurements which yield the half-light radius of 4.9+/-0.5 kpc and the semi-major axis of 12.5+/-0.1 kpc. G6025 is also very UV-luminous (~5L*(z~4}) and young (~140+/-60 Myr). Despite its unusual size and luminosity, the stellar population parameters and dust reddening (M_star~M*(z~4)$, and E(B-V)=0.18+/-0.05) estimated from the integrated light, are similar to those of smaller galaxies at comparable redshifts. The ground-based morphology and spectroscopy show two dominant components, both located off-center, embedded in more diffuse emission. We speculate that G6025 may be a scaled-up version of chain galaxies seen in deep HST imaging, or alternatively, a nearly equal-mass merger involving two super-L* galaxies in its early stage. G6025 lies close to but not within a known massive protocluster at z=3.78. We find four companions within 6 Mpc from G6025, two of which lie within 1.6 Mpc. While the limited sensitivity of the existing spectroscopy does not allow us to robustly characterize the local environment of G6025, it likely resides in a locally overdense environment. The luminosity, size, and youth of G6025 make it uniquely suited to study the early formation of massive galaxies in the universe.

We investigate the nature of the relation among stellar mass, star-formation rate, and gas-phase metallicity (the M$_*$-SFR-Z relation) at high redshifts using a sample of 260 star-forming galaxies at $z\sim2.3$ from the MOSDEF survey. We present an analysis of the high-redshift M$_*$-SFR-Z relation based on several emission-line ratios for the first time. We show that a M$_*$-SFR-Z relation clearly exists at $z\sim2.3$. The strength of this relation is similar to predictions from cosmological hydrodynamical simulations. By performing a direct comparison of stacks of $z\sim0$ and $z\sim2.3$ galaxies, we find that $z\sim2.3$ galaxies have $\sim0.1$ dex lower metallicity at fixed M$_*$ and SFR. In the context of chemical evolution models, this evolution of the M$_*$-SFR-Z relation suggests an increase with redshift of the mass-loading factor at fixed M$_*$, as well as a decrease in the metallicity of infalling gas that is likely due to a lower importance of gas recycling relative to accretion from the intergalactic medium at high redshifts. Performing this analysis simultaneously with multiple metallicity-sensitive line ratios allows us to rule out the evolution in physical conditions (e.g., N/O ratio, ionization parameter, and hardness of the ionizing spectrum) at fixed metallicity as the source of the observed trends with redshift and with SFR at fixed M$_*$ at $z\sim2.3$. While this study highlights the promise of performing high-order tests of chemical evolution models at high redshifts, detailed quantitative comparisons ultimately await a full understanding of the evolution of metallicity calibrations with redshift.

We present an analysis using the MOSFIRE Deep Evolution Field (MOSDEF) survey on the nature of "MIR-excess" galaxies, which have star formation rates (SFR) inferred from mid-infrared (MIR) data that is substantially elevated relative to that estimated from dust-corrected UV data. We use a sample of $\sim$200 galaxies and AGN at $1.40<z<2.61$ with 24 $\mu$m detections (rest-frame 8$\mu$m) from MIPS/\textit{Spitzer}. We find that the identification of MIR-excess galaxies strongly depends on the methodologies used to estimate IR luminosity ($\rm L_{IR}$) and to correct the UV light for dust attenuation. We find that extrapolations of the SFR from the observed 24 $\mu$m flux, using luminosity-dependent templates based on local galaxies, substantially overestimate $\rm L_{IR}$ in $z\sim2$ galaxies. By including \textit{Herschel} observations and using a stellar mass-dependent, luminosity-independent $\rm L_{IR}$, we obtain more reliable estimates of the SFR and a lower fraction of MIR-excess galaxies. Once stellar mass selection biases are taken into account, we identify $\sim24\%$ of our galaxies as MIR-excess. However, $\rm SFR_{H\alpha}$ is not elevated in MIR-excess galaxies compared to MIR-normal galaxies, indicating that the intrinsic fraction of MIR-excess may be lower. Using X-ray, IR, and optically-selected AGN in MOSDEF, we do not find a higher prevalence for AGN in MIR-excess galaxies relative to MIR-normal galaxies. A stacking analysis of X-ray undetected galaxies does not reveal a harder spectrum in MIR-excess galaxies relative to MIR-normal galaxies. Our analysis indicates that AGN activity does not contribute substantially to the MIR excess and instead implies that it is likely due to the enhanced PAH emission.

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

We present the first comprehensive evolutionary analysis of the rest-frame UV spectroscopic properties of star-forming galaxies at z ~ 2-4. We match samples at different redshifts in UV luminosity and stellar mass, and perform systematic measurements of spectral features and stellar population modeling. By creating composite spectra grouped according to Ly$\alpha$ equivalent width (EW), and various galaxy properties, we study the evolutionary trends among Ly$\alpha$, low- and high-ionization interstellar (LIS and HIS) absorption features, and integrated galaxy properties. We also examine the redshift evolution of Ly$\alpha$ and LIS absorption kinematics, and fine-structure emission EWs. The connections among the strengths of Ly$\alpha$, LIS lines, and dust extinction are redshift-independent, as is the decoupling of Ly$\alpha$ and HIS line strengths, and the bulk outflow kinematics as traced by LIS lines. Stronger Ly$\alpha$ emission is observed at higher redshift at fixed UV luminosity, stellar mass, SFR, and age. Much of this variation in average Ly$\alpha$ strength with redshift, and the variation in Ly$\alpha$ strength at fixed redshift, can be explained in terms of variations in neutral gas covering fraction and/or dust content in the ISM and CGM. However, based on the connection between Ly$\alpha$ and CIII] emission strengths, we additionally find evidence for variations in the intrinsic production rate of Ly$\alpha$ photons at the highest Ly$\alpha$ EWs. The challenge now is to understand the observed evolution in neutral gas covering fraction and dust extinction within a coherent model for galaxy formation, and make robust predictions for the escape of ionizing radiation at z > 6.

We combine spectroscopic measurements of H$\alpha$ and H$\beta$ and UV continuum photometry for a sample of 673 galaxies from the MOSFIRE Deep Evolution Field survey to constrain hydrogen ionizing photon production efficiencies ($\xi_{\rm ion}$, xi_ion) at z=1.4-2.6. We find average log(xi_ion/[Hz erg$^{-1}$])=25.06 (25.34), assuming the Calzetti (SMC) curve for the UV dust correction and a scatter of 0.28 dex in xi_ion distribution. After accounting for observational uncertainties and variations in dust attenuation, we conclude that the remaining scatter in xi_ion is likely dominated by galaxy-to-galaxy variations in stellar populations, including the slope and upper-mass cutoff of the initial mass function, stellar metallicity, star-formation burstiness, and stellar evolution (e.g., single/binary star evolution). Moreover, xi_ion is elevated in galaxies with high ionization states (high [OIII]/[OII]) and low oxygen abundances (low [NII]/H$\alpha$ and high [OIII]/H$\beta$) in the ionized ISM. However, xi_ion does not correlate with the offset from the z~0 star-forming locus in the BPT diagram, suggesting no change in the hardness of ionizing radiation accompanying the offset from the z~0 sequence. We also find that galaxies with blue UV spectral slopes ($\langle\beta\rangle$=-2.1) have elevated xi_ion by a factor of ~2 relative to the average xi_ion of the sample ($\langle\beta\rangle$=-1.4). If these blue galaxies are similar to those at z > 6, our results suggest that a lower Lyman continuum escape fraction is required for galaxies to maintain reionization, compared to the canonical xi_ion predictions from stellar population models. Furthermore, we demonstrate that even with robustly dust-corrected H$\alpha$, the UV dust attenuation can cause on average a ~0.3dex systematic uncertainty in xi_ion calculations.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on broad flux from the nebular emission lines H$\alpha$, [NII], [OIII], H$\beta$, and [SII]. The sample consists of 127 star-forming galaxies at $1.37 < z < 2.61$ and 84 galaxies at $2.95 < z < 3.80$. We decompose the emission lines using narrow ($\text{FWHM} < 275 \ \text{km s}^{-1}$) and broad ($\text{FWHM} > 300 \ \text{km s}^{-1}$) Gaussian components for individual galaxies and stacks. Broad emission is detected at $>3\sigma$ in $<10$% of galaxies and the broad flux accounts for 10-70% of the total flux. We find a slight increase in broad to narrow flux ratio with mass but note that we cannot reliably detect broad emission with $\text{FWHM} < 275 \ \text{km s}^{-1}$, which may be significant at low masses. Notably, there is a correlation between higher signal-to-noise (S/N) spectra and a broad component detection indicating a S/N dependence in our ability to detect broad flux. When placed on the N2-BPT diagram ([OIII]/H$\beta$ vs. [NII]/H$\alpha$) the broad components of the stacks are shifted towards higher [OIII]/H$\beta$ and [NII]/$\alpha$ ratios compared to the narrow component. We compare the location of the broad components to shock models and find that the broad component could be due to shocks, but we do not rule out other possibilities such as the presence of an AGN. We estimate the mass loading factor (mass outflow rate/star formation rate) assuming the broad component is a photoionized outflow and find that the mass loading factor increases as a function of mass which agrees with previous studies. We show that adding emission from shocked gas to $z\sim0$ SDSS spectra shifts galaxies towards the location of $z\sim2$ galaxies on several emission line diagnostic diagrams.

Using observations from the first two years of the MOSFIRE Deep Evolution Field (MOSDEF) survey, we study 13 AGN-driven outflows detected from a sample of 67 X-ray, IR and/or optically-selected AGN at $z \sim 2$. The AGN have bolometric luminosities of $\sim10^{44}-10^{46} ~\mathrm{erg~s^{-1}}$, including both quasars and moderate-luminosity AGN. We detect blueshifted, ionized gas outflows in the H$\beta$ , [OIII], H$\alpha$ ~and/or [NII] emission lines of $19\%$ of the AGN, while only 1.8\% of the MOSDEF galaxies have similarly-detected outflows. The outflow velocities span $\sim$300 to 1000 km s$^{-1}$. Eight of the 13 outflows are spatially extended on similar scales as the host galaxies, with spatial extents of 2.5 to 11.0 kpc. Outflows are detected uniformly across the star-forming main sequence, showing little trend with the host galaxy SFR. Line ratio diagnostics indicate that the outflowing gas is photoionized by the AGN. We do not find evidence for positive AGN feedback, in either our small MOSDEF sample or a much larger SDSS sample, using the BPT diagram. Given that a galaxy with an AGN is ten times more likely to have a detected outflow, the outflowing gas is photoionzed by the AGN, and estimates of the mass and energy outflow rates indicate that stellar feedback is insufficient to drive at least some of these outflows, they are very likely to be AGN-driven. The outflows have mass-loading factors of the order of unity, suggesting that they help regulate star formation in their host galaxies, though they may be insufficient to fully quench it.

We present the first spectroscopic measurement of multiple rest-frame optical emission lines at $z>4$. During the MOSFIRE Deep Evolution Field (MOSDEF) survey, we observed the galaxy GOODSN-17940 with the Keck I/MOSFIRE spectrograph. The K-band spectrum of GOODSN-17940 includes significant detections of the [OII]$\lambda\lambda 3726,3729$, [NeIII]$\lambda3869$, and H$\gamma$ emission lines and a tentative detection of H$\delta$, indicating $z_{\rm{spec}}=4.4121$. GOODSN-17940 is an actively star-forming $z>4$ galaxy based on its K-band spectrum and broadband spectral energy distribution. A significant excess relative to the surrounding continuum is present in the Spitzer/IRAC channel 1 photometry of GOODSN-17940, due primarily to strong H$\alpha$ emission with a rest-frame equivalent width of $\mbox{EW(H}\alpha)=1200$ \AA. Based on the assumption of $0.5 Z_{\odot}$ models and the Calzetti attenuation curve, GOODSN-17940 is characterized by $M_*=5.0^{+4.3}_{-0.2}\times 10^9 M_{\odot}$. The Balmer decrement inferred from H$\alpha$/H$\gamma$ is used to dust correct the H$\alpha$ emission, yielding $\mbox{SFR(H}\alpha)=320^{+190}_{-140} M_{\odot}\mbox{ yr}^{-1}$. These $M_*$ and SFR values place GOODSN-17940 an order of magnitude in SFR above the $z\sim 4$ star-forming "main sequence." Finally, we use the observed ratio of [NeIII]/[OII] to estimate the nebular oxygen abundance in GOODSN-17940, finding $\mbox{O/H}\sim 0.2 \mbox{ (O/H)}_{\odot}$. Combining our new [NeIII]/[OII] measurement with those from stacked spectra at $z\sim 0, 2, \mbox{ and } 3$, we show that GOODSN-17940 represents an extension to $z>4$ of the evolution towards higher [NeIII]/[OII] (i.e., lower $\mbox{O/H}$) at fixed stellar mass. It will be possible to perform the measurements presented here out to $z\sim 10$ using the James Webb Space Telescope.

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

We report the detection of diffuse Ly$\alpha$ emission, or Ly$\alpha$ halos (LAHs), around star-forming galaxies at $z\approx3.78$ and $2.66$ in the NOAO Deep Wide-Field Survey Bo\"otes field. Our samples consist of a total of $\sim$1400 galaxies, within two separate regions containing spectroscopically confirmed galaxy overdensities. They provide a unique opportunity to investigate how the LAH characteristics vary with host galaxy large-scale environment and physical properties. We stack Ly$\alpha$ images of different samples defined by these properties and measure their median LAH sizes by decomposing the stacked Ly$\alpha$ radial profile into a compact galaxy-like and an extended halo-like component. We find that the exponential scale-length of LAHs depends on UV continuum and Ly$\alpha$ luminosities, but not on Ly$\alpha$ equivalent widths or galaxy overdensity parameters. The full samples, which are dominated by low UV-continuum luminosity Ly$\alpha$ emitters ($M_{\rm UV} \gtrsim -21$), exhibit LAH sizes of 5$\,-\,6\,$kpc. However, the most UV- or Ly$\alpha$-luminous galaxies have more extended halos with scale-lengths of 7$\,-\,9\,$kpc. The stacked Ly$\alpha$ radial profiles decline more steeply than recent theoretical predictions that include the contributions from gravitational cooling of infalling gas and from low-level star formation in satellites. On the other hand, the LAH extent matches what one would expect for photons produced in the galaxy and then resonantly scattered by gas in an outflowing envelope. The observed trends of LAH sizes with host galaxy properties suggest that the physical conditions of the circumgalactic medium (covering fraction, HI column density, and outflow velocity) change with halo mass and/or star-formation rates.

We present results on the variation of 7.7 micron Polycyclic Aromatic Hydrocarbon (PAH) emission in galaxies spanning a wide range in metallicity at z ~ 2. For this analysis, we use rest-frame optical spectra of 476 galaxies at 1.37 < z < 2.61 from the MOSFIRE Deep Evolution Field (MOSDEF) survey to infer metallicities and ionization states. Spitzer/MIPS 24 micron and Herschel/PACS 100 and 160 micron observations are used to derive rest-frame 7.7 micron luminosities (L(7.7)) and total IR luminosities (L(IR)), respectively. We find significant trends between the ratio of L(7.7) to L(IR) (and to dust-corrected SFR) and both metallicity and [OIII]/[OII] (O32) emission-line ratio. The latter is an empirical proxy for the ionization parameter. These trends indicate a paucity of PAH emission in low metallicity environments with harder and more intense radiation fields. Additionally, L(7.7)/L(IR) is significantly lower in the youngest quartile of our sample (ages of 500 Myr) compared to older galaxies, which may be a result of the delayed production of PAHs by AGB stars. The relative strength of L(7.7) to L(IR) is also lower by a factor of ~ 2 for galaxies with masses $M_* < 10^{10}M_{\odot}$, compared to the more massive ones. We demonstrate that commonly-used conversions of L(7.7) (or 24 micron flux density; f(24)) to L(IR) underestimate the IR luminosity by more than a factor of 2 at $M_*$ ~ $10^{9.6-10.0} M_{\odot}$. We adopt a mass-dependent conversion of L(7.7) to L(IR) with L(7.7)/L(IR)= 0.09 and 0.22 for $M_* < 10^{10}$ and $> 10^{10} M_{\odot}$, respectively. Based on the new scaling, the SFR-$M_*$ relation has a shallower slope than previously derived. Our results also suggest a higher IR luminosity density at z ~ 2 than previously measured, corresponding to a ~ 30% increase in the SFR density.

We present a detailed study of the rest-optical (3600-7000 Angstrom) nebular spectra of ~380 star-forming galaxies at z~2-3 obtained with Keck/MOSFIRE as part of the Keck Baryonic Structure Survey (KBSS). The KBSS-MOSFIRE sample is representative of star-forming galaxies at these redshifts, with stellar masses M*=10^9-10^11.5 M_sun and star formation rates SFR=3-1000 M_sun/yr. We focus on robust measurements of many strong diagnostic emission lines for individual galaxies: [O II]3727,3729, [Ne III]3869, H-beta, [O III]4960,5008, [N II]6549,6585, H-alpha, and [S II]6718,6732. Comparisons with observations of typical local galaxies from the Sloan Digital Sky Survey (SDSS) and between subsamples of KBSS-MOSFIRE show that high-redshift galaxies exhibit a number of significant differences in addition to the well-known offset in log([O III]/H-beta) and log([N II]/H-alpha). We argue that the primary difference between H II regions in z~2.3 galaxies and those at z~0 is an enhancement in the degree of nebular excitation, as measured by [O III]/H-beta and R23=log[([O III]+[O II])/H-beta]. At the same time, KBSS-MOSFIRE galaxies are ~10 times more massive than z~0 galaxies with similar ionizing spectra and have higher N/O (likely accompanied by higher O/H) at fixed excitation. These results indicate the presence of harder ionizing radiation fields at fixed N/O and O/H relative to typical z~0 galaxies, consistent with Fe-poor stellar population models that include massive binaries, and highlight a population of massive, high-specific star formation rate galaxies at high-redshift with systematically different star formation histories than galaxies of similar stellar mass today.

Unlike spiral galaxies such as the Milky Way, the majority of the stars in massive elliptical galaxies were formed in a short period early in the history of the Universe. The duration of this formation period can be measured using the ratio of magnesium to iron abundance ([Mg/Fe]), which reflects the relative enrichment by core-collapse and type Ia supernovae. For local galaxies, [Mg/Fe] probes the combined formation history of all stars currently in the galaxy, including younger and metal-poor stars that were added during late-time mergers. Therefore, to directly constrain the initial star-formation period, we must study galaxies at earlier epochs. The most distant galaxy for which [Mg/Fe] had previously been measured is at z~1.4, with [Mg/Fe]=0.45(+0.05,-0.19). A slightly earlier epoch (z~1.6) was probed by stacking the spectra of 24 massive quiescent galaxies, yielding an average [Mg/Fe] of 0.31+/-0.12. However, the relatively low S/N of the data and the use of index analysis techniques for both studies resulted in measurement errors that are too large to allow us to form strong conclusions. Deeper spectra at even earlier epochs in combination with analysis techniques based on full spectral fitting are required to precisely measure the abundance pattern shortly after the major star-forming phase (z>2). Here we report a measurement of [Mg/Fe] for a massive quiescent galaxy at z=2.1. With [Mg/Fe]=0.59+/-0.11, this galaxy is the most Mg-enhanced massive galaxy found so far, having twice the Mg enhancement of similar-mass galaxies today. The abundance pattern of the galaxy is consistent with enrichment exclusively by core-collapse supernovae and with a star-formation timescale of 0.1-0.5 Gyr - characteristics that are similar to population II stars in the Milky Way. With an average past SFR of 600-3000 Msol/yr, this galaxy was among the most vigorous star-forming galaxies in the Universe.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on the identification, selection biases, and host galaxy properties of 55 X-ray, IR and optically-selected active galactic nuclei (AGN) at $1.4 < z < 3.8$. We obtain rest-frame optical spectra of galaxies and AGN and use the BPT diagram to identify optical AGN. We examine the uniqueness and overlap of the AGN identified at different wavelengths. There is a strong bias against identifying AGN at any wavelength in low mass galaxies, and an additional bias against identifying IR AGN in the most massive galaxies. AGN hosts span a wide range of star formation rate (SFR), similar to inactive galaxies once stellar mass selection effects are accounted for. However, we find (at $\sim 2-3\sigma$ significance) that IR AGN are in less dusty galaxies with relatively higher SFR and optical AGN in dusty galaxies with relatively lower SFR. X-ray AGN selection does not display a bias with host galaxy SFR. These results are consistent with those from larger studies at lower redshifts. Within star-forming galaxies, once selection biases are accounted for, we find AGN in galaxies with similar physical properties as inactive galaxies, with no evidence for AGN activity in particular types of galaxies. This is consistent with AGN being fueled stochastically in any star-forming host galaxy. We do not detect a significant correlation between SFR and AGN luminosity for individual AGN hosts, which may indicate the timescale difference between the growth of galaxies and their supermassive black holes.

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

We investigate the burstiness of star formation histories (SFHs) of galaxies at $0.4<z<1$ by using the ratio of star formation rates (SFRs) measured from H$\beta$ and FUV (1500 \AA) (H$\beta$--to--FUV ratio). Our sample contains 164 galaxies down to stellar mass (M*) of $10^{8.5} M_\odot$ in the CANDELS GOODS-N region, where Team Keck Redshift Survey DEIMOS spectroscopy and HST/WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the {\it ratio} of H$\beta$- and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H$\beta$--to--FUV ratio of our sample increases with M* and SFR. The median ratio is $\sim$0.7 at M*$\sim10^{8.5} M_\odot$ (or SFR$\sim 0.5 M_\odot/yr$) and increases to $\sim$1 at M*$\sim10^{10} M_\odot$ (or SFR $\sim 10 M_\odot/yr$). At M*$<10^{9.5} M_\odot$, our median H$\beta$--to--FUV ratio is lower than that of local galaxies at the same M*, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation of our results, and the importance of the burstiness increases as M* decreases. Due to sample selection effects, our H$\beta$--to--FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

We present observations of Q1549-C25, an ~L* star-forming galaxy at z=3.15 for which Lyman-continuum (LyC) radiation is significantly detected in deep Keck/LRIS spectroscopy. We find no evidence for contamination from a lower-redshift interloper close to the line of sight in the high signal-to-noise spectrum of Q1549-C25. Furthermore, the morphology of Q1549-C25 in V_606, J_125, and H_160 Hubble Space Telescope (HST) imaging reveals that the object consists of a single, isolated component within 1". In combination, these data indicate Q1549-C25 as a clean spectroscopic detection of LyC radiation, only the second such object discovered to date at z~3. We model the spectral energy distribution (SED) of Q1549-C25, finding evidence for negligible dust extinction, an age (assuming continuous star formation) of ~1 Gyr, and a stellar mass of M_*=7.9x10^9 M_sun. Although it is not possible to derive strong constraints on the absolute escape fraction of LyC emission, f_esc(LyC), from a single object, we use simulations of intergalactic and circumgalactic absorption to infer f_esc(LyC)>=0.51 at 95% confidence. The combination of deep Keck/LRIS spectroscopy and HST imaging is required to assemble a larger sample of objects like Q1549-C25, and obtain robust constraints on the average f_esc(LyC) at z~3 and beyond.

We present measurements of the electron-temperature based oxygen abundance for a highly star-forming galaxy at z=3.08, COSMOS-1908. This is the highest redshift at which [OIII]$\lambda$4363 has been detected, and the first time that this line has been measured at z>2. We estimate an oxygen abundance of 12+log(O/H)$=8.00^{+0.13}_{-0.14}$. This galaxy is a low-mass ($10^{9.3}$ M$_{\odot}$), highly star-forming ($\sim50$ M$_{\odot}$ yr$^{-1}$) system that hosts a young stellar population ($\sim160$ Myr). We investigate the physical conditions of the ionized gas in COSMOS-1908 and find that this galaxy has a high ionization parameter, little nebular reddening ($E(B-V)_{\rm gas}<0.14$), and a high electron density ($n_e\sim500$ cm$^{-3}$). We compare the ratios of strong oxygen, neon, and hydrogen lines to the direct-method oxygen abundance for COSMOS-1908 and additional star-forming galaxies at z=0-1.8 with [OIII]$\lambda$4363 measurements, and show that galaxies at z$\sim$1-3 follow the same strong-line correlations as galaxies in the local universe. This agreement suggests that the relationship between ionization parameter and O/H is similar for z$\sim$0 and high-redshift galaxies. These results imply that metallicity calibrations based on lines of oxygen, neon, and hydrogen do not strongly evolve with redshift and can reliably estimate abundances out to z$\sim$3, paving the way for robust measurements of the evolution of the mass-metallicity relation to high redshift.

We use a large sample of galaxies at z~3 to establish a relationship between reddening, neutral gas covering fraction (fcov(HI)), and the escape of ionizing photons at high redshift. Our sample includes 933 galaxies at z~3, 121 of which have very deep spectroscopic observations (>7 hrs) in the rest-UV (lambda=850-1300 A) with Keck/LRIS. Based on the high covering fraction of outflowing optically-thick HI indicated by the composite spectra of these galaxies, we conclude that photoelectric absorption, rather than dust attenuation, dominates the depletion of ionizing photons. By modeling the composite spectra as the combination of an unattenuated stellar spectrum including nebular continuum emission with one that is absorbed by HI and reddened by a line-of-sight extinction, we derive an empirical relationship between E(B-V) and fcov(HI). Galaxies with redder UV continua have larger covering fractions of HI characterized by higher line-of-sight extinctions. Our results are consistent with the escape of Lya through gas-free lines-of-sight. Covering fractions based on low-ionization interstellar absorption lines systematically underpredict those deduced from the HI lines, suggesting that much of the outflowing gas may be metal-poor. We develop a model which connects the ionizing escape fraction with E(B-V), and which may be used to estimate the escape fraction for an ensemble of high-redshift galaxies. Alternatively, direct measurements of the escape fraction for our data allow us to constrain the intrinsic 900-to-1500 A flux density ratio to be >0.20, a value that favors stellar population models that include weaker stellar winds, a flatter initial mass function, and/or binary evolution. Lastly, we demonstrate how the framework discussed here may be used to assess the pathways by which ionizing radiation escapes from high-redshift galaxies. [Abridged]

Star-forming galaxies form a sequence in the [OIII]/H-beta vs. [NII]/H-alpha diagnostic diagram, with low metallicity, highly ionized galaxies falling in the upper left corner. Drawing from a large sample of UV-selected star-forming galaxies at z~2 with rest-frame optical nebular emission line measurements from Keck-MOSFIRE, we select the extreme ~5% of the galaxies lying in this upper left corner, requiring log([NII]/H-alpha) <= 1.1 and log([OIII]/H-beta) >= 0.75. These cuts identify galaxies with 12 + log(O/H) <~ 8.0, when oxygen abundances are measured via the O3N2 diagnostic. We study the Ly-alpha properties of the resulting sample of 14 galaxies. The mean (median) rest-frame Ly-alpha equivalent width is 39 (36) A, and 11 of the 14 objects (79%) are Ly-alpha-emitters (LAEs) with W_Lya > 20 A. We compare the equivalent width distribution of a sample of 522 UV-selected galaxies at 2.0<z<2.6 identified without regard to their optical line ratios; this sample has mean (median) Ly-alpha equivalent width -1 (-4) A, and only 9% of these galaxies qualify as LAEs. The extreme galaxies typically have lower attenuation at Ly-alpha than those in the comparison sample, and have ~50% lower median oxygen abundances. Both factors are likely to facilitate the escape of Ly-alpha: in less dusty galaxies Ly-alpha photons are less likely to be absorbed during multiple scatterings, while the harder ionizing spectrum and higher ionization parameter associated with strong, low metallicity star formation may reduce the covering fraction or column density of neutral hydrogen, further easing Ly-alpha escape. The use of nebular emission line ratios may prove useful in the identification of galaxies with low opacity to Ly-alpha photons across a range of redshifts.

We present the first measurements of the shape of the far-ultraviolet (far-UV; lambda=950-1500 A) dust attenuation curve at high redshift (z~3). Our analysis employs rest-frame UV spectra of 933 galaxies at z~3, 121 of which have very deep spectroscopic observations (>7 hrs) at lambda=850-1300 A, with the Low Resolution Imaging Spectrograph on the Keck Telescope. By using an iterative approach in which we calculate the ratios of composite spectra in different bins of continuum color excess, E(B-V), we derive a dust curve that implies a lower attenuation in the far-UV for a given E(B-V) than those obtained with standard attenuation curves. We demonstrate that the UV composite spectra of z~3 galaxies can be modeled well by assuming our new attenuation curve, a high covering fraction of HI, and absorption from the Lyman-Werner bands of H2 with a small (<20%) covering fraction. The low covering fraction of H2 relative to that of the HI and dust suggests that most of the dust in the ISM of typical galaxies at z~3 is unrelated to the catalysis of H2, and is associated with other phases of the ISM (i.e., the ionized and neutral gas). The far-UV dust curve implies a factor of ~2 lower dust attenuation of Lyman continuum (ionizing) photons relative to those inferred from the most commonly assumed attenuation curves for L* galaxies at z~3. Our results may be utilized to assess the degree to which ionizing photons are attenuated in HII regions or, more generally, in the ionized or low column density (N(HI)<10^17.2 cm^-2) neutral ISM of high-redshift galaxies.

We present a combined analysis of rest-frame far-UV (1000-2000 A) and rest-frame optical (3600-7000 A) composite spectra formed from very deep observations of a sample of 30 star-forming galaxies with z=2.4+/-0.1, selected to be representative of the full KBSS-MOSFIRE spectroscopic survey. Since the same massive stars are responsible for the observed FUV continuum and the excitation of the observed nebular emission, a self-consistent stellar population synthesis model must simultaneously match the details of the far-UV stellar+nebular continuum and-- when inserted as the excitation source in photoionization models-- account for all observed nebular emission line ratios. We find that only models including massive star binaries, having low stellar metallicity (Z_*/Z_{sun} ~ 0.1) but relatively high ionized gas-phase oxygen abundances (Z_{neb}/Z_{sun} ~ 0.5), can successfully match all of the observational constraints. We argue that this apparent discrepancy is naturally explained by highly super-solar O/Fe [4-5 times (O/Fe)_{sun}], expected for gas whose enrichment is dominated by the products of core-collapse supernovae. Once the correct ionizing spectrum is identified, photoionization models reproduce all of the observed strong emission line ratios, the direct T_e measurement of O/H, and allow accurate measurement of the gas-phase abundance ratios of N/O and C/O -- both of which are significantly sub-solar but, as for O/Fe, are in remarkable agreement with abundance patterns observed in Galactic thick disk, bulge, and halo stars with similar O/H. High nebular excitation is the rule at high-z (and rare at low-z) because of systematically shorter enrichment timescales (<<1 Gyr): low Fe/O environments produce harder (and longer-lived) stellar EUV spectra at a given O/H, enhanced by dramatic effects on the evolution of massive star binaries.

We present new observations of the field containing the z=3.786 protocluster, PC217.96+32.3. We confirm that it is one of the largest and most overdense high-redshift structures known. Such structures are rare even in the largest cosmological simulations. We used the Mayall/MOSAIC1.1 imaging camera to image a 1.2x0.6 deg area (~150x75 comoving Mpc) surrounding the protocluster's core and discovered 165 candidate Lyman Alpha emitting galaxies (LAEs) and 788 candidate Lyman Break galaxies (LBGs). There are at least 2 overdense regions traced by the LAEs, the largest of which shows an areal overdensity in its core (i.e., within a radius of 2.5 comoving Mpc) of 14+/-7 relative to the average LAE spatial density in the imaged field. Further, the average LAE spatial density in the imaged field is twice that derived by other field LAE surveys. Spectroscopy with Keck/DEIMOS yielded redshifts for 164 galaxies (79 LAEs and 85 LBGs); 65 lie at a redshift of 3.785+/-0.010. The velocity dispersion of galaxies near the core is 350+/-40 km/s, a value robust to selection effects. The overdensities are likely to collapse into systems with present-day masses of >10^{15} solar masses and >6x10^{14} solar masses. The low velocity dispersion may suggest a dynamically young protocluster. We find a weak trend between narrow-band (Lyman Alpha) luminosity and environmental density: the Lyman Alpha luminosity is enhanced on average by 1.35X within the protocluster core. There is no evidence that the Lyman Alpha equivalent width depends on environment. These suggest that star-formation and/or AGN activity is enhanced in the higher density regions of the structure. PC217.96+32.3 is a Coma cluster analog, witnessed in the process of formation.

We present the first direct comparison between Balmer line and panchromatic SED-based SFRs for z~2 galaxies. For this comparison we used 17 star-forming galaxies selected from the MOSFIRE Deep Evolution Field (MOSDEF) survey, with $3\sigma$ detections for H$\alpha$ and at least two IR bands (Spitzer/MIPS 24$\mu$m and Herschel/PACS 100 and 160$\mu$m, and in some cases Herschel/SPIRE 250, 350, and 500$\mu$m). The galaxies have total IR (8-1000$\mu$m) luminosities of $\sim10^{11.4}-10^{12.4}\,\textrm{L}_\odot$ and star-formation rates (SFRs) of $\sim30-250\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We fit the UV-to-far-IR SEDs with flexible stellar population synthesis (FSPS) models - which include both stellar and dust emission - and compare the inferred SFRs with the SFR(H$\alpha$,H$\beta$) values corrected for dust attenuation using Balmer decrements. The two SFRs agree with a scatter of 0.17 dex. Our results imply that the Balmer decrement accurately predicts the obscuration of the nebular lines and can be used to robustly calculate SFRs for star-forming galaxies at z~2 with SFRs up to $\sim200\,\textrm{M}_\odot\,\mathrm{yr^{-1}}$. We also use our data to assess SFR indicators based on modeling the UV-to-mid-IR SEDs or by adding SFR(UV) and SFR(IR), for which the latter is based on the mid-IR only or on the full IR SED. All these SFRs show a poorer agreement with SFR(H$\alpha$,H$\beta$) and in some cases large systematic biases are observed. Finally, we show that the SFR and dust attenuation derived from the UV-to-near-IR SED alone are unbiased when assuming a delayed exponentially declining star-formation history.

We use a sample of 149 spectroscopically confirmed UV-selected galaxies at $z\sim 2$ to investigate the relative dust attenuation of the stellar continuum and the nebular emission lines. For each galaxy in the sample, at least one rest-frame optical emission line (H$\alpha$/[NII]$\lambda6583$ or [OIII]$\lambda5007$) measurement has been taken from the litterature, and 41 galaxies have additional Spitzer/MIPS 24$\mu$m observations that are used to infer infrared luminosities. We use a spectral energy distribution (SED) fitting code that predicts nebular line strengths when fitting the stellar populations of galaxies in our sample, and we perform comparisons between the predictions of our models and the observed/derived physical quantities. We find that on average our code is able to reproduce all the physical quantities (e.g., UV $\beta$ slopes, infrared luminosities, emission line fluxes), but we need to apply a higher dust correction to the nebular emission compared to the stellar emission for the largest SFR ($\log\mathrm{(SFR/M}_\odot\mathrm{yr}^{-1})>1.87$, Salpeter IMF). We find a correlation between SFR and the difference in nebular and stellar color excesses, which could resolve the discrepant results regarding nebular dust correction at $z\sim2$ from previous results.

We present H$\alpha$ gas kinematics for 178 star-forming galaxies at z~2 from the MOSFIRE Deep Evolution Field survey. We have developed models to interpret the kinematic measurements from fixed-angle multi-object spectroscopy, using structural parameters derived from CANDELS HST/F160W imaging. For 35 galaxies we measure resolved rotation with a median $(V/\sigma_{V,0})_{R_E}=2.1$. We derive dynamical masses from the kinematics and sizes and compare them to baryonic masses, with gas masses estimated from dust-corrected H$\alpha$ star formation rates (SFRs) and the Kennicutt-Schmidt relation. When assuming that galaxies with and without observed rotation have the same median $(V/\sigma_{V,0})_{R_E}$, we find good agreement between the dynamical and baryonic masses, with a scatter of $\sigma_{RMS}=0.34$dex and a median offset of $\Delta\log_{10}M=0.04$dex. This comparison implies a low dark matter fraction (8% within an effective radius) for a Chabrier initial mass function (IMF), and disfavors a Salpeter IMF. Moreover, the requirement that $M_{dyn}/M_{baryon}$ should be independent of inclination yields a median value of $(V/\sigma_{V,0})_{R_E}=2.1$ for galaxies without observed rotation. If instead we treat the galaxies without detected rotation as early-type galaxies, the masses are also in reasonable agreement ($\Delta\log_{10}M=-0.07$dex, $\sigma_{RMS}=0.37$dex). The inclusion of gas masses is critical in this comparison; if gas masses are excluded there is an increasing trend of $M_{dyn}/M_{*}$ with higher specific SFR (SSFR). Furthermore, we find indications that $V/\sigma$ decreases with increasing H$\alpha$ SSFR for our full sample, which may reflect disk settling. We also study the Tully-Fisher relation and find that at fixed stellar mass $S_{0.5}=(0.5V_{2.2}^2+\sigma_{V,0}^2)^{1/2}$ was higher at earlier times. At fixed baryonic mass, we observe the opposite trend. [abridged]

We present results on the star-formation rate (SFR) versus stellar mass ($M_*$) relation (i.e., the "main sequence") among star-forming galaxies at 1.37$\leq$$z$$\leq$2.61 using the MOSFIRE Deep Evolution Field (MOSDEF) survey. Based on a sample of 261 galaxies with H$\alpha$ and H$\beta$ spectroscopy, we have estimated robust dust-corrected instantaneous SFRs over a large range in $M_*$ ($\sim10^{9.5}-10^{11.5}M_\odot$). We find a correlation between log(SFR(H$\alpha$)) and log($M_*$) with a slope of 0.65$\pm$0.08 (0.58$\pm$0.10) at 1.4<z<2.6 (2.1<z<2.6). We find that different assumptions for the dust correction, such as using the color-excess of the stellar continuum to correct the nebular lines, sample selection biases against red star-forming galaxies, and not accounting for Balmer absorption can yield steeper slopes of the log(SFR)-log($M_*$) relation. Our sample is immune from these biases as it is rest-frame optically selected, H$\alpha$ and H$\beta$ are corrected for Balmer absorption, and the H$\alpha$ luminosity is dust-corrected using the nebular color-excess computed from the Balmer decrement. The scatter of the log(SFR(H$\alpha$))-log($M_*$) relation, after accounting for the measurement uncertainties, is 0.31 dex at 2.1<z<2.6, which is 0.05 dex larger than the scatter in log(SFR(UV))-log($M_*$). Based on comparisons to a simulated SFR-$M_*$ relation with some intrinsic scatter, we argue that in the absence of direct measurements of galaxy-to-galaxy variations in the attenuation/extinction curves and the IMF, one cannot use the difference in the scatter of the SFR(H$\alpha$)- and SFR(UV)-$M_*$ relations to constrain the stochasticity of star formation in high-redshift galaxies.

Using observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey, we investigate the physical conditions of star-forming regions in $z\sim2.3$ galaxies, specifically the electron density and ionization state. From measurements of the [O II]$\lambda\lambda$3726,3729 and [S II]$\lambda\lambda$6716,6731 doublets, we find a median electron density of $\sim250$ cm$^{-3}$ at $z\sim2.3$, an increase of an order of magnitude compared to measurements of galaxies at $z\sim0$. While $z\sim2.3$ galaxies are offset towards significantly higher O$_{32}$ values relative to local galaxies at fixed stellar mass, we find that the high-redshift sample follows a similar distribution to the low-metallicity tail of the local distribution in the O$_{32}$ vs. R$_{23}$ and O3N2 diagrams. Based on these results, we propose that $z\sim2.3$ star-forming galaxies have the same ionization parameter as local galaxies at fixed metallicity. In combination with simple photoionization models, the position of local and $z\sim2.3$ galaxies in excitation diagrams suggests that there is no significant change in the hardness of the ionizing spectrum at fixed metallicity from $z\sim0$ to $z\sim2.3$. We find that $z\sim2.3$ galaxies show no offset compared to low-metallicity local galaxies in emission line ratio diagrams involving only lines of hydrogen, oxygen, and sulfur, but show a systematic offset in diagrams involving [N II]$\lambda$6584. We conclude that the offset of $z\sim2.3$ galaxies from the local star-forming sequence in the [N II] BPT diagram is primarily driven by elevated N/O at fixed O/H compared to local galaxies. These results suggest that the local gas-phase and stellar metallicity sets the ionization state of star-forming regions at $z\sim0$ and $z\sim2$.

We present $U_{336}V_{606}J_{125}H_{160}$ follow-up $HST$ observations of 16 $z\sim3$ candidate LyC emitters in the HS1549+1919 field. With these data, we obtain high spatial-resolution photometric redshifts of all sub-arcsecond components of the LyC candidates in order to eliminate foreground contamination and identify robust candidates for leaking LyC emission. Of the 16 candidates, we find one object with a robust LyC detection that is not due to foreground contamination. This object (MD5) resolves into two components; we refer to the LyC-emitting component as MD5b. MD5b has an observed 1500\AA\ to 900\AA\ flux-density ratio of $(F_{UV}/F_{LyC})_{obs}=4.0\pm2.0$, compatible with predictions from stellar population synthesis models. Assuming minimal IGM absorption, this ratio corresponds to a relative (absolute) escape fraction of $f_{esc,rel}^{MD5b}=75-100$% ($f_{esc,abs}^{MD5b}=14-19$%). The stellar population fit to MD5b indicates an age of $\lesssim50$Myr, which is in the youngest 10% of the $HST$ sample and the youngest third of typical $z\sim3$ Lyman break galaxies, and may be a contributing factor to its LyC detection. We obtain a revised, contamination-free estimate for the comoving specific ionizing emissivity at $z=2.85$, indicating (with large uncertainties) that star-forming galaxies provide roughly the same contribution as QSOs to the ionizing background at this redshift. Our results show that foreground contamination prevents ground-based LyC studies from obtaining a full understanding of LyC emission from $z\sim3$ star-forming galaxies. Future progress in direct LyC searches is contingent upon the elimination of foreground contaminants through high spatial-resolution observations, and upon acquisition of sufficiently deep LyC imaging to probe ionizing radiation in high-redshift galaxies.

We report the redshift of an unlensed, highly obscured submillimetre galaxy (SMG), HS1700.850.1, the brightest SMG (S850um =19.1 mJy) detected in the JCMT/SCUBA-2 Baryonic Structure Survey, based on the detection of its CO line emission. Using the IRAM PdBI-WIDEX with 3.6GHz band width, we serendipitously detect an emission line at 150.6 GHz. From a search over 14.5 GHz in the 3-mm and 2-mm atmospheric windows, we confirm the identification of this line as CO(5-4) at z = 2.816, meaning that it does not reside in the z~2.30 proto-cluster in this field. Measurement of the 870um source size (<0.85") from the Sub-Millimeter Array (SMA) confirms a compact emission in a S870um =14.5mJy, LIR~10^13 Lsun component, suggesting an Eddington-limited starburst. We use the double-peaked CO line profile measurements along with the SMA size constraints to study the gas dynamics of a HyLIRG, estimating the gas and dynamical masses of HS1700.850.1. While HS1700.850.1 is one of the most extreme galaxies known in the Universe, we find that it occupies a relative void in the Lyman-Break Galaxy distribution in this field. Comparison with other extreme objects at similar epochs (HyLIRG Quasars), and cosmological simulations, suggests such an anti-bias of bright SMGs could be relatively common, with the brightest SMGs rarely occupying the most overdense regions at z=2-4.

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

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

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

We present results on the dust attenuation curve of z~2 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. Our sample consists of 224 star-forming galaxies with nebular spectroscopic redshifts in the range z= 1.36-2.59 and high S/N measurements of, or upper limits on, the H-alpha and H-beta emission lines obtained with Keck/MOSFIRE. We construct composite SEDs of galaxies in bins of specific SFR and Balmer optical depth in order to directly constrain the dust attenuation curve from the UV through near-IR for typical star-forming galaxies at high redshift. Our results imply an attenuation curve that is very similar to the SMC extinction curve at wavelengths redward of 2500 Angstroms. At shorter wavelengths, the shape of the curve is identical to that of the Calzetti relation, but with a lower normalization (R_V). Hence, the new attenuation curve results in SFRs that are ~20% lower, and log stellar masses that are 0.16 dex lower, than those obtained with the Calzetti attenuation curve. Moreover, we find that the difference in the reddening---and the total attenuation---of the ionized gas and stellar continuum correlates strongly with SFR, such that for dust-corrected SFRs larger than 20 Msun/yr assuming a Chabrier IMF, the nebular emission lines suffer an increasing degree of obscuration relative to the continuum. A simple model that can account for these trends is one in which the UV through optical stellar continuum is dominated by a population of less reddened stars, while the nebular line and bolometric luminosities become increasingly dominated by dustier stellar populations for galaxies with large SFRs, as a result of the increased dust enrichment that accompanies such galaxies. Consequently, UV- and SED-based SFRs may underestimate the total SFR at even modest levels of ~20 Msun/yr. [Abridged]

We use a sample of 262 spectroscopically confirmed star-forming galaxies at redshifts $2.08\leq z\leq 2.51$ to compare H$\alpha$, UV, and IR star-formation-rate diagnostics and to investigate the dust properties of the galaxies. At these redshifts, the H$\alpha$ line shifts to the $K_{s}$-band. By comparing $K_{s}$-band photometry to underlying stellar population model fits to other UV, optical, and near-infrared data, we infer the H$\alpha$ flux for each galaxy. We obtain the best agreement between H$\alpha$- and UV-based SFRs if we assume that the ionized gas and stellar continuum are reddened by the same value and that the Calzetti attenuation curve is applied to both. Aided with MIPS 24$\mu$m data, we find that an attenuation curve steeper than the Calzetti curve is needed to reproduce the observed IR/UV ratios of galaxies younger than 100 Myr. Furthermore, using the bolometric star-formation rate inferred from the UV and mid-IR data (SFR$_{IR}$+SFR$_{UV}$), we calculated the conversion between the H$\alpha$ luminosity and SFR to be $(7.5\pm1.3) \times 10^{-42}$ for a Salpeter IMF, which is consistent with the Kennicutt (1998) conversion. The derived conversion factor is independent of any assumption of the dust correction and is robust to stellar population model uncertainties.

We report the blind detection of 12CO emission from a Distant Red Galaxy, HS1700.DRG55. We have used the IRAM PdBI-WIDEX, with its 3.6GHz of instantaneous dual-polarization bandwidth, to target 12CO(3--2) from galaxies lying in the proto-cluster at z=2.300 in the field HS1700+64. If indeed this line in DRG55 is 12CO(3--2), it's detection at 104.9GHz indicates a z_CO=2.296. None of the other eight known z~2.30 proto-cluster galaxies lying within the primary beam (PB) are detected in 12CO, although the limits are ~2x worse towards the edge of the PB where several lie. The optical/near-IR magnitudes of DRG55 (R_AB>27, K_AB=22.3) mean that optical spectroscopic redshifts are difficult with 10m-class telescopes, but near-IR redshifts would be feasible. The 24um-implied SFR (210 M_odot yr-1), stellar mass (~10^11 M-odot) and 12CO line luminosity (3.6x10^10 K km s-1 pc^2) are comparable to other normal 12CO-detected star forming galaxies in the literature, although the galaxy is some ~2 mag (~6x) fainter in the rest-frame UV than 12CO-detected galaxies at z>2. The detection of DRG55 in 12CO complements three other 12CO detected UV-bright galaxies in this proto-cluster from previous studies, and suggests that many optically faint galaxies in the proto-cluster may host substantial molecular gas reservoirs, and a full blind census of 12CO in this overdense environment is warranted.

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

We present results on the excitation properties of z~2.3 galaxies using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) Survey. With its coverage of the full suite of strong rest-frame optical emission lines, MOSDEF provides an unprecedented view of the rest-frame optical spectra of a representative sample of distant star-forming galaxies. We investigate the locations of z~2.3 MOSDEF galaxies in multiple emission-line diagnostic diagrams. These include the [OIII]/Hb vs. [NII]/Ha and [OIII]/Hb vs. [SII]/Ha "BPT" diagrams, as well as the O_32 vs. R_23 excitation diagram. We recover the well-known offset in the star-forming sequence of high-redshift galaxies in the [OIII]/Hb vs. [NII]/Ha BPT diagram relative to SDSS star-forming galaxies. However, the shift for our rest-frame optically selected sample is less significant than for rest-frame-UV selected and emission-line selected galaxies at z~2. Furthermore, we find that the offset is mass-dependent, only appearing within the low-mass half of the z~2.3 MOSDEF sample, where galaxies are shifted towards higher [NII]/Ha at fixed [OIII]/Hb. Within the [OIII]/Hb vs. [SII]/Ha and O_32 vs. R_23 diagrams, we find that z~2.3 galaxies are distributed like local ones, and therefore attribute the shift in the [OIII]/Hb vs. [NII]/Ha BPT diagram to elevated N/O abundance ratios among lower-mass (M_*<10^10 M_sun) high-redshift galaxies. The variation in N/O ratios calls into question the use at high redshift of oxygen abundance indicators based on nitrogen lines, but the apparent invariance with redshift of the excitation sequence in the O_32 vs. R_23 diagram paves the way for using the combination of O_32 and R_23 as an unbiased metallicity indicator over a wide range in redshift. This indicator will allow for an accurate characterization of the shape and normalization of the mass-metallicity relationship over more than 10 Gyr.

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on rest-frame optical AGN identification and completeness at z~2.3. With our sample of 50 galaxies and 10 X-ray and IR-selected AGN with measured H-beta, [OIII], H-alpha, and [NII] emission lines, we investigate the location of AGN in the BPT, MEx (mass-excitation), and CEx (color-excitation) diagrams. We find that the BPT diagram works well to identify AGN at z~2.3 and that the z~0 AGN/star-forming galaxy classifications do not need to shift substantially at z~2.3 to robustly separate these populations. However, the MEx diagram fails to identify all of the AGN identified in the BPT diagram, and the CEx diagram is substantially contaminated at high redshift. We further show that AGN samples selected using the BPT diagram have selection biases in terms of both host stellar mass and stellar population, in that AGN in low mass and/or high specific star formation rate galaxies are difficult to identify using the BPT diagram. These selection biases become increasingly severe at high redshift, such that optically-selected AGN samples at high redshift will necessarily be incomplete. We also find that the gas in the narrow-line region appears to be more enriched than gas in the host galaxy for at least some MOSDEF AGN. However, AGN at z~2 are generally less enriched than local AGN with the same host stellar mass.

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 present results on the z~2.3 mass-metallicity relation (MZR) using early observations from the MOSFIRE Deep Evolution Field (MOSDEF) survey. We use an initial sample of 87 star-forming galaxies with spectroscopic coverage of H\beta, [OIII]\lambda 5007, H\alpha, and [NII]\lambda 6584 rest-frame optical emission lines, and estimate the gas-phase oxygen abundance based on the N2 and O3N2 strong-line indicators. We find a positive correlation between stellar mass and metallicity among individual z~2.3 galaxies using both the N2 and O3N2 indicators. We also measure the emission-line ratios and corresponding oxygen abundances for composite spectra in bins of stellar mass. Among composite spectra, we find a monotonic increase in metallicity with increasing stellar mass, offset ~0.15-0.3 dex below the local MZR. When the sample is divided at the median star-formation rate (SFR), we do not observe significant SFR dependence of the z~2.3 MZR among either individual galaxies or composite spectra. We furthermore find that z~2.3 galaxies have metallicities ~0.1 dex lower at a given stellar mass and SFR than is observed locally. This offset suggests that high-redshift galaxies do not fall on the local "fundamental metallicity relation" among stellar mass, metallicity, and SFR, and may provide evidence of a phase of galaxy growth in which the gas reservoir is built up due to inflow rates that are higher than star-formation and outflow rates. However, robust conclusions regarding the gas-phase oxygen abundances of high-redshift galaxies await a systematic reappraisal of the application of locally calibrated metallicity indicators at high redshift.

We study the Lya profiles of 36 spectroscopically-detected Lya-emitters (LAEs) at z2-3, using Keck MOSFIRE to measure systemic redshifts and velocity dispersions from rest-frame optical nebular emission lines. The sample has a median optical magnitude R=26.0, and ranges from R23 to R>27, corresponding to rest-frame UV absolute magnitudes M_UV-22 to M_UV>-18.2. Dynamical masses range from M_dyn<1.3e8 Msun to M_dyn=6.8e9 Msun, with a median value of M_dyn=6.3e8 Msun. Thirty of the 36 Lya emission lines are redshifted with respect to the systemic velocity with at least 1 sigma significance, and the velocity offset with respect to systemic dv_Lya is correlated with R-band magnitude, M_UV, and the velocity dispersion measured from nebular emission lines with >3 sigma significance: brighter galaxies with larger velocity dispersions tend to have larger values of dv_Lya. We also make use of a comparison sample of 122 UV-color-selected R<25.5 galaxies at z2, all with Lya emission and systemic redshifts measured from nebular emission lines. Using the combined LAE and comparison samples for a total of 158 individual galaxies, we find that dv_Lya is anti-correlated with the Lya equivalent width with 7 sigma significance. Our results are consistent with a scenario in which the Lya profile is determined primarily by the properties of the gas near the systemic redshift; in such a scenario, the opacity to Lya photons in lower mass galaxies may be reduced if large gaseous disks have not yet developed and if the gas is ionized by the harder spectrum of young, low metallicity stars.

We present initial results of a deep near-IR spectroscopic survey covering the 15 fields of the Keck Baryonic Structure Survey (KBSS) using MOSFIRE on the Keck 1 telescope, focusing on a sample of 251 galaxies with redshifts 2.0< z < 2.6, star-formation rates 2 < SFR < 200 M_sun/yr, and stellar masses 8.6 < log(M*/M_sun) < 11.4, with high-quality spectra in both H- and K-band atmospheric windows. We show unambiguously that the locus of z~2.3 galaxies in the "BPT" nebular diagnostic diagram exhibits a disjoint, yet similarly tight, relationship between the ratios [NII]6585/Halpha and [OIII]/Hbeta as compared to local galaxies. Using photoionization models, we argue that the offset of the z~2.3 locus relative to z~ 0 is explained by a combination of harder ionizing radiation field, higher ionization parameter, and higher N/O at a given O/H than applies to most local galaxies, and that the position of a galaxy along the z~2.3 star-forming BPT locus is surprisingly insensitive to gas-phase oxygen abundance. The observed nebular emission line ratios are most easily reproduced by models in which the net ionizing radiation field resembles a blackbody with effective temperature T_eff = 50000-60000 K and N/O close to the solar value at all O/H. We critically assess the applicability of commonly-used strong line indices for estimating gas-phase metallicities, and consider the implications of the small intrinsic scatter in the empirical relationship between excitation-sensitive line indices and stellar mass (i.e., the "mass-metallicity" relation), at z~2.3.

We report the discovery of a large-scale structure containing multiple protoclusters at z=3.78 in the Bo\"otes field. The spectroscopic discovery of five galaxies at z=3.783+/-0.002 lying within 1 Mpc of one another led us to undertake a deep narrow- and broad-band imaging survey of the surrounding field. Within a comoving volume of 72x72x25 Mpc^3, we have identified 65 Lyman alpha emitter (LAE) candidates at z=3.795+/-0.015, and four additional galaxies at z_spec=3.730,3.753,3.780,3.835. The galaxy distribution within the field is highly non-uniform, exhibiting three large (~3-5x) overdensities separated by 8-14 Mpc (physical) and possibly connected by filamentary structures traced by LAEs. The observed number of LAEs in the entire field is nearly twice the average expected in field environments, based on estimates of the Lya luminosity function at these redshifts. We estimate that by z=0 the largest overdensity will grow into a cluster of mass 10^15 Msun; the two smaller overdensities will grow into clusters of mass (2-6)x10^14 Msun. The highest concentration of galaxies is located at the southern end of the image, suggesting that the current imaging may not map the true extent of the large scale structure. Finding three large protocluster candidates within a single 0.3 deg^2 field is highly unusual; expectations from theory suggest that such alignments should occur less than 2% of the time. Searching for and characterizing such structures and accurately measuring their volume space density can therefore place constraints on the theory of structure formation. Such regions can also serve as laboratories for the study of galaxy formation in dense environments.

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.

Several works have now shown that nebular emission can have a significant impact on broad-band photometry of high-redshift galaxies ($z > 3$), and how this can affect parameter estimation from SED fitting. So far relatively small spectroscopic samples have been used to measure this effect. Here we focus on a large spectroscopic sample (N$\sim 2300$) at $z\sim2$, with measured H\alpha\ fluxes for $\sim100$ galaxies. Under appropriate assumptions, our SED fitting code is able to reproduce observed H\alpha\ fluxes, and we infer that $\sim20$% of our sample have parameters significantly affected when nebular emission is taken into account. We also determine how nebular emission can affect selection and parameter estimation of evolved galaxies (Balmer Break galaxies) at high-redshift ($z\sim4$).

We present results from a survey for z~2.85 Lyman-Continuum (LyC) emission in the HS1549+1933 field and place constraints on the amount of ionizing radiation escaping from star-forming galaxies. Using a custom narrowband filter (NB3420) tuned to wavelengths just below the Lyman limit at z>=2.82, we probe the LyC spectral region of 49 Lyman break galaxies (LBGs) and 91 Lya-emitters (LAEs) spectroscopically confirmed at z>=2.82. Four LBGs and seven LAEs are detected in NB3420. Using V-band data probing the rest-frame non-ionizing UV, we observe that many NB3420-detected galaxies exhibit spatial offsets between their LyC and non-ionizing UV emission and are characterized by extremely blue NB3420-V colors, corresponding to low ratios of non-ionizing to ionizing radiation (F_UV/F_LyC) that are in tension with current stellar population synthesis models. We measure average values of F_UV/F_LyC for our LBG and LAE samples, correcting for foreground galaxy contamination and HI absorption in the IGM. We find (F_UV/F_LyC)_corr^LBG=82 +/- 45 and (F_UV/F_LyC)_corr^LAE=7.4 +/- 3.6. These flux-density ratios correspond respectively to relative LyC escape fractions of f_esc,rel^LBG=5-8% and f_esc,rel^LAE=18-49%, absolute LyC escape fractions of f_esc^LBG=1-2% and f_esc^LAE=5-15%, and a comoving LyC emissivity from star-forming galaxies of 8.8-15.0 x 10^24 ergs/s/Hz/Mpc^3. In order to study the differential properties of galaxies with and without LyC detections, we analyze narrowband Lya imaging and rest-frame near-infrared imaging, finding that while LAEs with LyC detections have lower Lya equivalent widths on average, there is no substantial difference in the rest-frame near-infrared colors of LBGs or LAEs with and without LyC detections. These preliminary results are consistent with an orientation-dependent model where LyC emission escapes through cleared paths in a patchy ISM.

We present Keck/MOSFIRE observations of the role of environment in the formation of galaxies at z~2. Using K-band spectroscopy of H-alpha and [N II] emission lines, we have analyzed the metallicities of galaxies within and around a z=2.3 protocluster discovered in the HS1700+643 field. Our main sample consists of 23 protocluster and 20 field galaxies with estimates of stellar masses and gas-phase metallicities based on the N2 strong-line metallicity indicator. With these data we have examined the mass-metallicity relation (MZR) with respect to environment at z~2. We find that field galaxies follow the well-established trend between stellar mass and metallicity, such that more massive galaxies have larger metallicities. The protocluster galaxies, however, do not exhibit a dependence of metallicity on mass, with the low-mass protocluster galaxies showing an enhancement in metallicity compared to field galaxies spanning the same mass range. A comparison with galaxy formation models suggests that the mass-dependent environmental trend we observed can be qualitatively explained in the context of the recycling of "momentum-driven" galaxy wind material. Accordingly, winds are recycled on a shorter timescale in denser environments, leading to an enhancement in metallicity at fixed mass for all but the most massive galaxies. Future hydrodynamical simulations of z~2 overdensities matching the one in the HS1700 field will be crucial for understanding the origin of the observed environmental trend in detail.

We present DEIMOS spectroscopic observations of the most UV-luminous star-forming galaxies at 3.2<z<4.6. Our sample contains galaxies with luminosities of L*<L<7L* and is one of the largest samples to date of the most UV-luminous galaxies at these redshifts. Our data confirm 41 star-forming galaxies at 3.2<z<4.6 and validate the clean selection of the photometric candidates. We find that the fraction of Lya emitting galaxies increases with decreasing UV luminosity. We find strong evidence of large-scale outflows, transporting the neutral/ionized gas in the interstellar medium away from the galaxy. Galaxies exhibiting both interstellar absorption and Lya emission lines show a significant velocity offset between the two features (200-1140 km/s). We find tentative evidence that this measure of the outflow velocity increases with UV luminosity and/or stellar mass. The luminosity- and mass-dependent outflow strengths suggest that the efficiency of feedback and enrichment of the surrounding medium depend on these parameters. We present composite spectra of the absorption-line-only and Lya-emitting subsets of the UV luminous galaxies at z~3.7. The composite spectra are similar to those of lower-z and lower-luminosity LBGs samples, but with some subtle differences. Analyses of the composite spectra suggest that the UV luminous LBGs at z~3.7 may have a higher covering fraction of absorbing gas, and may be older than their lower-z and lower-luminosity counterparts. In addition, we have discovered 5 galaxies that belong to a massive overdensity at z=3.78. Finally, two galaxies each show two distinct sets of interstellar absorption features. The latter may be a sign of a final stage of major merger, or clumpy disk formation. Their presence implies that frequency of such sources among our luminous z~3.7 LBGs may be an order of magnitude higher than in lower redshift and lower luminosity samples.

We present observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

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

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

We use stellar population synthesis modeling to analyze the host galaxy properties of a sample of 33 UV-selected, narrow-lined active galactic nuclei (AGNs) at z ~ 2 - 3. In order to quantify the contribution of AGN emission to host galaxy broadband spectral energy distributions (SEDs), we use the subsample of 11 AGNs with photometric coverage spanning from rest-frame UV through near-IR wavelengths. Modeling the SEDs of these objects with a linear combination of stellar population and AGN templates, we infer the effect of the AGN on derived stellar population parameters. We also estimate the typical bias in derived stellar populations for AGNs lacking rest-frame near-IR wavelength coverage, and develop a method for inferring the true host galaxy properties. We compare AGN host galaxy properties to those of a sample of UV-selected, star-forming non-AGNs in the same redshift range, including a subsample carefully matched in stellar mass. Although the AGNs have higher masses and SFRs than the full non-active sample, their stellar population properties are consistent with those of the mass-selected sample, suggesting that the presence of an AGN is not connected with the cessation of star-formation activity in star-forming galaxies at z ~ 2 - 3. We suggest that a correlation between M_BH and galaxy stellar mass is already in place at this epoch. Assuming a roughly constant Eddington ratio for AGNs at all stellar masses, we are unable to detect the AGNs in low-mass galaxies because they are simply too faint.

We analyze rest-frame optical morphologies and gas-phase kinematics as traced by rest-frame far-UV and optical spectra for a sample of 204 star forming galaxies in the redshift range z ~ 2-3 drawn from the Keck Baryonic Structure Survey (KBSS). We find that spectroscopic properties and gas-phase kinematics are closely linked to morphology: compact galaxies with semi-major axis radii r <~ 2 kpc are substantially more likely than their larger counterparts to exhibit LyA in emission. Although LyA emission strength varies widely within galaxies of a given morphological type, all but one of 19 galaxies with LyA equivalent width W_LyA > 20 Angstroms have compact and/or multiple-component morphologies with r <= 2.5 kpc. The velocity structure of absorption lines in the galactic continuum spectra also varies as a function of morphology. Galaxies of all morphological types drive similarly strong outflows (as traced by the blue wing of interstellar absorption line features), but the outflows of larger galaxies are less highly ionized and exhibit larger optical depth at the systemic redshift that may correspond to a decreasing efficiency of feedback in evacuating gas from the galaxy. This v ~ 0 km/s gas is responsible both for shifting the mean absorption line redshift and attenuating W_LyA (via a longer resonant scattering path) in galaxies with larger rest-optical half light radii. In contrast to galaxies at lower redshifts, there is no evidence for a correlation between outflow velocity and inclination, suggesting that outflows from these puffy and irregular systems may be poorly collimated. (Abbrev.)

We report the Herschel SPIRE detection of dust emission arising from UV-luminous (L>~L*) star-forming galaxies at 3.3<z<4.3. Our sample of 1,913 Lyman Break Galaxy (LBG) candidates is selected over an area of 5.3 deg2 in the Bootes Field of the NOAO Deep Wide-Field Survey. This is one of the largest samples of UV-luminous galaxies at this epoch and enables an investigation of the bright end of the galaxy luminosity function. We divide our sample into three luminosity bins and stack the Herschel SPIRE data to measure the average spectral energy distribution (SED) of LBGs at far-infrared (FIR) wavelengths. We find that these galaxies have average IR luminosities of (3-5)x10^{11}L_sun and 60-70% of their star-formation obscured by dust. The FIR SEDs peak at lambda>100um, suggesting dust temperatures (T_d=27-30K) significantly colder than that of local galaxies of comparable IR luminosities. The observed IR-to-UV luminosity ratio (IRX=L_IR/L_UV) is low (~3-4) compared with that observed for z~2 LBGs (IRX(z~2)~7.1+/-1.1). The correlation between the slope of the UV continuum and IRX for galaxies in the two lower luminosity bins suggests dust properties similar to those of local starburst galaxies. However, the galaxies in the highest luminosity bin appear to deviate from the local relation, suggesting that their dust properties may differ from those of their lower-luminosity and low-redshift counterparts. We speculate that the most UV luminous galaxies at this epoch are being observed in a short-lived and young evolutionary phase.

We present observations of the CO[3-2] emission towards two massive and infrared luminous Lyman Break Galaxies at z = 3.21 and z = 2.92, using the IRAM Plateau de Bure Interferometer, placing first constraints on the molecular gas masses (Mgas) of non-lensed LBGs. Their overall properties are consistent with those of typical (Main-Sequence) galaxies at their redshifts, with specific star formation rates ~1.6 and ~2.2 Gyr^(-1), despite their large infrared luminosities L_IR ~2-3 x 10^12 Lsun derived from Herschel. With one plausible CO detection (spurious detection probability of 10^(-3)) and one upper limit, we investigate the evolution of the molecular gas-to-stellar mass ratio (Mgas/M*) with redshift. Our data suggest that the steep evolution of Mgas/M* of normal galaxies up to z~2 is followed by a flattening at higher redshifts, providing supporting evidence for the existence of a plateau in the evolution of the specific star formation rate at z > 2.5.

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 study dust attenuation at UV wavelengths at high redshift, where the UV is redshifted to the observed visible. In particular, we search for a UV bump and related implications for dust attenuation determinations. We use data in the CDFS, obtained in intermediate and broad band filters by the MUSYC project, to sample the UV rest-frame of 751 galaxies with 0.95<z<2.2. When available, Herschel/PACS data (GOODS-Herschel project), and Spitzer/MIPS measurements, are used to estimate the dust emission. The SED of each source is fit using the CIGALE code. The amount of dust attenuation and the dust attenuation curve are obtained as outputs of the SED fitting process, together with other parameters linked to the SFH. The global amount of dust attenuation at UV wavelengths is found to increase with stellar mass and to decrease as UV luminosity increases. A UV bump at 2175A is securely detected in 20% of the galaxies, and the mean amplitude of the bump for the sample is similar to that observed in the LMC supershell region. This amplitude is found to be lower in galaxies with very high SSFRs, and 90% of the galaxies exhibiting a secure bump are at z<1.5. The attenuation curve is confirmed to be steeper than that of local starburst galaxies for 20$% of the galaxies. The large dispersion found for these two parameters describing the attenuation law is likely to reflect a wide diversity of attenuation laws among galaxies. The relations between dust attenuation, IR-to-UV flux ratio, and the slope of the UV continuum are derived for the mean attenuation curve found for our sample. Deviations from the average trends are found to correlate with the age of the young stellar population and the shape of the attenuation curve.(abriged)

Although relatively common in the local Universe, only one grand-design spiral galaxy has been spectroscopically confirmed to lie at z>2 (HDFX 28; z=2.011), and may prove to be a major merger that simply resembles a spiral in projection. The rarity of spirals has been explained as a result of disks being dynamically 'hot' at z>2 which may instead favor the formation of commonly-observed clumpy structures. Alternatively, current instrumentation may simply not be sensitive enough to detect spiral structures comparable to those in the modern Universe. At redshifts <2, the velocity dispersion of disks decreases, and spiral galaxies are more numerous by z~1. Here we report observations of the grand design spiral galaxy Q2343-BX442 at z=2.18. Spectroscopy of ionized gas shows that the disk is dynamically hot, implying an uncertain origin for the spiral structure. The kinematics of the galaxy are consistent with a thick disk undergoing a minor merger, which can drive the formation of short-lived spiral structure. A duty cycle of < 100 Myr for such tidally-induced spiral structure in a hot massive disk is consistent with their rarity.

A large sample of spectroscopically confirmed galaxies at 1.4<z<3.7, with complementary imaging in the near- and mid-IR from the ground and from Hubble and Spitzer, is used to infer the average star formation histories (SFHs) of typical galaxies from z~7 to 2. For a subset of 302 galaxies at 1.5<z<2.6, we perform a comparison of star formation rates (SFRs) determined from SED modeling (SFRs[SED]) and those calculated from deep Keck UV and Spitzer/MIPS 24 micron imaging (SFRs[IR+UV]). Exponentially declining SFHs yield SFRs[SED] that are 5-10x lower on average than SFRs[IR+UV], indicating that declining SFHs may not be accurate for typical galaxies at z>2. The SFRs of z~2-3 galaxies are directly proportional to their stellar masses M*, with unity slope---a result that is confirmed with Spitzer/IRAC stacks of 1179 UV-faint (R>25.5) galaxies---for M*>5e8 Msun and SFRs >2 Msun/yr. We interpret this result in the context of several systematic biases that can affect determinations of the SFR-M* relation. The average specific SFRs at z~2-3 are similar within a factor of two to those measured at z>4, implying an average SFH where SFRs increase with time. A consequence of these rising SFHs is that (a) a substantial fraction of UV-bright z~2-3 galaxies had faint sub-L* progenitors at z>4; and (b) gas masses must increase with time from z=7 to 2, over which time the net cold gas accretion rate---as inferred from the specific SFR and the Kennicutt-Schmidt relation---is ~2-3x larger than the SFR . However, if we evolve to higher redshift the SFHs and masses of the halos that are expected to host L* galaxies at z~2, we find that <10% of the baryons accreted onto typical halos at z>4 actually contribute to star formation at those epochs. These results highlight the relative inefficiency of star formation even at early cosmic times when galaxies were first assembling. [Abridged]

We present results from the Keck Baryonic Structure Survey (KBSS), a unique spectroscopic survey designed to explore the connection between galaxies and intergalactic baryons. The KBSS is optimized for the redshift range z ~ 2-3, combining S/N ~ 100 Keck/HIRES spectra of 15 hyperluminous QSOs with densely sampled galaxy redshift surveys surrounding each QSO sightline. We perform Voigt profile decomposition of all 6000 HI absorbers within the full Lya forest in the QSO spectra. Here we present the distribution, column density, kinematics, and absorber line widths of HI surrounding 886 star-forming galaxies with 2.0 < z < 2.8 and within 3 Mpc of a QSO sightline. We find that N_HI and the multiplicity of HI components increase rapidly near galaxies. The strongest HI absorbers within ~ 100 physical kpc of galaxies have N_HI ~ 3 dex higher than those near random locations in the IGM. The circumgalactic zone of most enhanced HI absorption (CGM) is found within 300 kpc and 300 km/s of galaxies. Nearly half of absorbers with log(N_HI) > 15.5 are found within the CGM of galaxies meeting our photometric selection, while their CGM occupy only 1.5% of the cosmic volume. The spatial covering fraction, multiplicity of absorption components, and characteristic N_HI remain elevated to transverse distances of 2 physical Mpc. Absorbers with log(N_HI) > 14.5 are tightly correlated with the positions of galaxies, while absorbers with lower N_HI are correlated only on Mpc scales. Redshift anisotropies on Mpc scales indicate coherent infall toward galaxies, while on scales of ~100 physical kpc peculiar velocities of 260 km/s are indicated. The median Doppler widths of absorbers within 1-3 virial radii of galaxies are ~50% larger than randomly chosen absorbers of the same N_HI, suggesting higher gas temperatures and/or increased turbulence likely caused by accretion shocks and/or galactic winds.

We present the results of a 42-orbit HST/WFC3 survey of the rest-frame optical morphologies of star forming galaxies with spectroscopic redshifts in the range z=1.5-3.6. The survey consists of 42 orbits of F160W imaging covering ~65 arcmin^2 distributed widely across the sky and reaching a depth of 27.9 AB for a 5 sigma detection within a 0.2 arcsec radius aperture. Focusing on an optically selected sample of 306 star forming galaxies with stellar masses in the range M* = 10^9 - 10^11 Msun, we find that typical circularized effective half-light radii range from ~ 0.7 - 3.0 kpc and describe a stellar mass - radius relation as early as z ~ 3. While these galaxies are best described by an exponential surface brightness profile, their distribution of axis ratios is strongly inconsistent with a population of inclined exponential disks and is better reproduced by triaxial stellar systems with minor/major and intermediate/major axis ratios ~ 0.3 and 0.7 respectively. While rest-UV and rest-optical morphologies are generally similar for a subset of galaxies with HST/ACS imaging data, differences are more pronounced at higher masses M* > 3 x 10^10 Msun. Finally, we discuss galaxy morphology in the context of efforts to constrain the merger fraction, finding that morphologically-identified mergers/non-mergers generally have insignificant differences in terms of physical observables such as stellar mass and star formation rate, although merger-like galaxies selected according to some criteria have statistically smaller effective radii and correspondingly larger SFR surface density.

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

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

Using a sample of 92 UV continuum-selected, spectroscopically identified galaxies with <z> = 2.65, all of which have been imaged in the Ly-a line with extremely deep narrow-band imaging, we examine galaxy Ly-a emission profiles to very faint surface brightness limits. The galaxies are representative of spectroscopic samples of LBGs at similar redshifts in terms of apparent magnitude, UV luminosity, inferred extinction, and star formation rate, and were selected without regard to Ly-a emission properties. We use extremely deep stacks of UV continuum and Ly-a emission line images to show that all sub-samples exhibit diffuse Ly-a emission to radii of at least 10" (80 physical kpc), including galaxies whose spectra exhibit Ly-a in net absorption. The intensity scaling, but not the surface brightness distribution, is strongly correlated with the emission observed in the central ~1". The characteristic scale length for Ly-a line emission exceeds that of the UV continuum light for the same galaxies by factors of at least 5-10, regardless of the spectral morphology of Ly-a. Including the extended Ly-a halos increases the total Ly-a flux [and rest equivalent width W_0(Lya)] by an average factor of 5. We argue that most, if not all, of the observed Ly-a emission in the diffuse halos originates in the galaxy H II regions and is scattered in our direction by H I gas in the galaxy's circum-galactic medium (CGM). We show that whether or not a galaxy is classified as a giant "Lyman-a Blob" (LAB) depends sensitively on the Ly-a surface brightness threshold reached by an observation. Accounting for diffuse Ly-a halos, all LBGs would be LABs if surveys were routinely sensitive to 10 times lower surface brightness thresholds; also, essentially all LBGs would qualify as LAEs (W_0(Lya) > 20 A).

[Abridged] We investigate the average physical properties and star formation histories of the most UV-luminous star-forming galaxies at z~3.7. Our results are derived from analyses of the average spectral energy distributions (SEDs), constructed from stacked optical to infrared photometry, of a sample of the 1,902 most UV-luminous star-forming galaxies found in 5.3 square degrees of the NOAO Deep Wide-Field Survey. We bin the sample according to UV luminosity, and find that the shape of the average SED in the rest-frame optical and infrared is fairly constant with UV luminosity: i.e., more UV luminous galaxies are, on average, also more luminous at longer wavelengths. In the rest-UV, however, the spectral slope (measured at 0.13-0.28 um) rises steeply with the median UV luminosity from -1.8 at L L* to -1.2 in the brightest bin (L~4-5L*). We use population synthesis analyses to derive the average physical properties of these galaxies and find that: (1) L_UV, and thus star formation rates (SFRs), scale closely with stellar mass such that more UV-luminous galaxies are also more massive; (2) The median ages indicate that the stellar populations are relatively young (200-400 Myr) and show little correlation with UV luminosity; and (3) More UV-luminous galaxies are dustier than their less-luminous counterparts, such that L~4-5L* galaxies are extincted up to A(1600)=2 mag while L L* galaxies have A(1600)=0.7-1.5 mag. Based on these observations, we argue that the average star formation histories of UV-luminous galaxies are better described by models in which SFR increases with time in order to simultaneously reproduce the tight correlation between the observed SFR and stellar mass, and the universally young ages of these galaxies. We demonstrate the potential of measurements of the SFR-M* relation at multiple redshifts to discriminate between simple models of star formation histories.

The recently refurbished Hubble Space Telescope reveals a galaxy from a time when the Universe was just 500 million years old, providing insights into the first throes of galaxy formation and the reionization of the Universe.

Increasingly large samples of galaxies are now being discovered at redshifts z~5-6 and higher. Many of these objects are inferred to be young, low in mass, and relatively unreddened, but detailed analysis of their high quality spectra will not be possible until the advent of future facilities. In this paper we shed light on the physical conditions in a plausibly similar low mass galaxy by presenting the analysis of the rest-frame optical and UV spectra of Q2343-BX418, an L* galaxy at z=2.3 with a very low mass-to-light ratio and unusual properties: BX418 is young (<100 Myr), low mass (M_star ~ 10^9 Msun), low in metallicity (Z ~ 1/6 Zsun), and unreddened (E(B-V)~0.02, UV continuum slope beta=-2.1). We infer a metallicity 12+log(O/H)=7.9 +/- 0.2 from the rest-frame optical emission lines. We also determine the metallicity via the direct, electron temperature method, using the ratio O III] 1661, 1666/[O III] 5007 to determine the electron temperature and finding 12+ log(O/H)=7.8 +/- 0.1. These measurements place BX418 among the most metal-poor galaxies observed in emission at high redshift. The rest-frame UV spectrum contains strong emission from Lya (with rest-frame equivalent width 54 A), He II 1640 (both stellar and nebular), C III] 1907, 1909 and O III] 1661, 1666. The C IV/C III] ratio indicates that the source of ionization is unlikely to be an AGN. Analysis of the He II, O III] and C III] line strengths indicates a very high ionization parameter log U ~ -1, while Lya and the interstellar absorption lines indicate that outflowing gas is highly ionized over a wide range of velocities. It remains to be determined how many of BX418's unique spectral features are due to its global properties, such as low metallicity and dust extinction, and how many are indicative of a short-lived phase in the early evolution of an otherwise normal star-forming galaxy.

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

We present new results on the kinematics and spatial distribution of metal-enriched gas within 125 kpc (physical) of Lyman Break galaxies at redshifts z~2-3. In particular, we demonstrate how rest-UV galaxy spectra can be used to obtain key spatial and spectral information more efficiently than possible with QSO sightlines. After recalibrating the measurement of galaxy systemic redshifts from their UV spectra, we investigate the kinematics of galaxy-scale outflows via the strong interstellar (IS) absorption and Lya emission lines (when present), as well as their dependence on other physical properties of the galaxies. We construct a sample of 512 close (1-15 arcsec) angular pairs of z~2-3 LBGs in which the spectra background galaxies probe the circumgalactic gas surrrounding those in the foreground. The close pairs, together with spectra of the foreground galaxies themselves, sample galactocentric impact parameters b=0-125 kpc (physical) at <z>=2.2. The ensemble provides a spatial map of cool gas as a function of galactocentric distance for a well-characterized population of galaxies. We propose a simple model that simultaneously matches the kinematics, depth, and profile shape of IS absorption and Lya emission lines, as well as the observed variation of absorption line strength (of HI, CII, CIV, SiII, SiIV) versus galactocentric impact parameter. We discuss the results of the observations in the context of "cold accretion", in which cool gas accretes via filamentary streams directly onto the central regions of galaxies. At present, we find little observational support for cool infalling material, whereas evidence supporting the large-scale effects of outflows is strong. Reconciling theory and observation on the subject of gas flows into and out of forming galaxies seems necessary.

We present the results of spectroscopic, narrow-band and X-ray observations of a z=2.30 protocluster in the field of the QSO HS 1700+643. Using a sample of BX/MD galaxies, which are selected to be at z~2.2-2.7 by their rest-frame ultraviolet colours, we find that there are 5 protocluster AGN which have been identified by characteristic emission-lines in their optical/near-IR spectra; this represents an enhancement over the field significant at ~98.5 per cent confidence. Using a ~200 ks Chandra/ACIS-I observation of this field we detect a total of 161 X-ray point sources to a Poissonian false-probability limit of 4x10^{-6} and identify 8 of these with BX/MD galaxies. Two of these are spectroscopically confirmed protocluster members and are also classified as emission-line AGN. When compared to a similarly selected field sample the analysis indicates this is also evidence for an enhancement of X-ray selected BX/MD AGN over the field, significant at ~99 per cent confidence. Deep Lya narrow-band imaging reveals that a total of 4/123 Lya emitters (LAEs) are found to be associated with X-ray sources, with two of these confirmed protocluster members and one highly likely member. We do not find a significant enhancement of AGN activity in this LAE sample over that of the field (result significant at only 87 per cent confidence). The X-ray emitting AGN fractions for the BX/MD and LAE samples are found to be 6.9_{-4.4}^{+9.2} and 2.9_{-1.6}^{+2.9} per cent, respectively, for protocluster AGN with L_{2-10 keV}>4.6x10^{43} erg s^{-1} at z=2.30. These findings are similar to results from the z=3.09 protocluster in the SSA 22 field found by Lehmer et al. (2009), in that both suggest AGN activity is favoured in dense environments at z>2.

We use a sample of 90 spectroscopically-confirmed Lyman Break Galaxies with H-alpha and 24 micron observations to constrain the relationship between rest-frame 8 micron luminosity, L(8), and star formation rate (SFR) for L* galaxies at z~2. We find a tight correlation with 0.24 dex scatter between L8 and L(Ha)/SFR for z~2 galaxies with L(IR)~10^10 - 10^12 Lsun. Employing this relationship with a larger sample of 392 galaxies with spectroscopic redshifts, we find that the UV slope can be used to recover the dust attenuation of the vast majority of L* galaxies at z~2 to within 0.4 dex scatter using the local correlation. Separately, young galaxies with ages <100 Myr appear to follow an extinction curve that is steeper than the one found for local starburst galaxies. Therefore, such young galaxies may be significantly less dusty than inferred previously. Our results provide the first direct evidence, independent of the UV slope, for a correlation between UV and bolometric luminosity at high redshift, in the sense that UV-faint galaxies are also on average less infrared and less bolometrically-luminous than their UV-bright counterparts. Further, as the SFR increases, the UV luminosity reaches a maximum value corresponding to L* at z~2, implying that dust obscuration may be largely responsible for modulating the bright-end of the UV luminosity function. L* galaxies at z~2, while at least an order of magnitude more bolometrically-luminous, exhibit ratios of metals-to-dust that are similar to those of local starbursts. This result is expected if high-redshift galaxies are forming their stars in a less metal-rich environment compared to local galaxies of the same luminosity, thus naturally leading to a redshift evolution in both the luminosity-metallicity and luminosity-obscuration relations. [Abridged]

We present the results of a photometric and spectroscopic survey of 321 Lyman break galaxies (LBGs) at z ~ 3 to investigate systematically the relationship between Lya emission and stellar populations. Lya equivalent widths (EW) were calculated from rest-frame UV spectroscopy and optical/near-infrared/Spitzer photometry was used in population synthesis modeling to derive the key properties of age, dust extinction, star formation rate (SFR), and stellar mass. We directly compare the stellar populations of LBGs with and without strong Lya emission, where we designate the former group (EW > 20 AA) as Lya emitters (LAEs) and the latter group (EW < 20 AA) as non-LAEs. This controlled method of comparing objects from the same UV luminosity distribution represents an improvement over previous studies in which the stellar populations of LBGs and narrowband-selected LAEs were contrasted, where the latter were often intrinsically fainter in broadband filters by an order of magnitude simply due to different selection criteria. Using a variety of statistical tests, we find that Lya equivalent width and age, SFR, and dust extinction, respectively, are significantly correlated in the sense that objects with strong Lya emission also tend to be older, lower in star formation rate, and less dusty than objects with weak Lya emission, or the line in absorption. We accordingly conclude that, within the LBG sample, objects with strong Lya emission represent a later stage of galaxy evolution in which supernovae-induced outflows have reduced the dust covering fraction. We also examined the hypothesis that the attenuation of Lya photons is lower than that of the continuum, as proposed by some, but found no evidence to support this picture.

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

We report on the identification of the highest redshift submm-selected source currently known: LESSJ033229.4-275619. This source was detected in the Large Apex BOlometer CAmera (LABOCA) Extended Chandra Deep Field South (ECDFS) Submillimetre Survey (LESS), a sensitive 870-um survey (~1.2-mJy rms) of the full 30'x30' ECDFS with the LABOCA camera on the Atacama Pathfinder EXperiment (APEX) telescope. The submm emission is identified with a radio counterpart for which optical spectroscopy provides a redshift of z=4.76. We show that the bolometric emission is dominated by a starburst with a star formation rate of ~1000 Msun/yr, although we also identify a moderate luminosity Active Galactic Nucleus (AGN) in this galaxy. Thus it has characteristics similar to those of z~2 submm galaxies (SMGs), with a mix of starburst and obscured AGN signatures. This demonstrates that ultraluminous starburst activity is not just restricted to the hosts of the most luminous (and hence rare) QSOs at z~5, but was also occurring in less extreme galaxies at a time when the Universe was less than 10% of its current age. Assuming that we are seeing the major phase of star formation in this galaxy, then we demonstrate that it would be identified as a luminous distant red galaxy at z~3 and that the current estimate of the space density of z>4 SMGs is only sufficient to produce ~10% of the luminous red galaxy population at these early times. However, this leaves open the possibility that some of these galaxies formed through less intense, but more extended star formation events. If the progenitors of all of the luminous red galaxies at z~3 go through an ultraluminous starburst at z>4 then the required volume density of z>4 SMGs will exceed that predicted by current galaxy formation models by more than an order of magnitude.

We use the deep ground-based optical photometry of the Lyman Break Galaxy (LBG) Survey to derive robust measurements of the faint-end slope (alpha) of the UV LF at redshifts 1.9<z<3.4. Our sample includes >2000 spectroscopic redshifts and ~31000 LBGs in 31 spatially-independent fields over a total area of 3261 arcmin^2. These data allow us to select galaxies to 0.07L* and 0.10L* at z~2 and z~3, respectively. A maximum likelihood analysis indicates steep values of alpha(z=2)=-1.73+/-0.07 and alpha(z=3)=-1.73+/-0.13. This result is robust to luminosity dependent systematics in the Ly-alpha equivalent width and reddening distributions, is similar to the steep values advocated at z>4, and implies that ~93% of the unobscured UV luminosity density at z~2-3 arises from sub-L* galaxies. With a realistic luminosity dependent reddening distribution, faint to moderately luminous galaxies account for >70% and >25% of the bolometric luminosity density and present-day stellar mass density, respectively, when integrated over 1.9<z<3.4. We find a factor of 8-9 increase in the star formation rate density between z~6 and z~2, due to both a brightening of L* and an increasing dust correction proceeding to lower redshifts. The previously observed discrepancy between the integral of the star formation history and stellar mass density measurements at z~2 may be reconciled by invoking a luminosity dependent reddening correction to the star formation history combined with an accounting for the stellar mass contributed by UV-faint galaxies. The steep and relatively constant alpha of the UV LF at z>2 contrasts with the shallower value inferred locally, suggesting that the evolution in the faint-end slope may be dictated simply by the availability of low mass halos capable of supporting star formation at z<2. [Abridged]

Many models that seek to explain the origin of the unresolved X-ray background predict that Compton-thick Active Galactic Nuclei (AGNs) are ubiquitious at high redshift. However, few distant Compton-thick AGNs have been reliably identified to date. Here we present Spitzer-IRS spectroscopy and 3.6-70um photometry of a z=2.2 optically identified AGN (HDF-oMD49) that is formally undetected in the 2Ms Chandra Deep Field-North (CDF-N) survey. The Spitzer-IRS spectrum and spectral energy distribution of this object is AGN dominated, and a comparison of the energetics at X-ray wavelengths to those derived from mid-infrared (mid-IR) and optical spectroscopy shows that the AGN is intrinsically luminous (L_X~3x10^44 erg/s) but heavily absorbed by Compton-thick material (N_H>>10^24 cm^{-2}); i.e., this object is a Compton-thick quasar. Adopting the same approach that we applied to HDF-oMD49, we found a further six objects at z~2-2.5 in the literature that are also X-ray weak/undetected but have evidence for AGN activity from optical and/or mid-IR spectroscopy, and show that all of these sources are also Compton-thick quasars with L_X>10^44 erg/s. On the basis of the definition of Daddi etal. (2007), these Compton-thick quasars would be classified as mid-IR excess galaxies, and our study provides the first spectroscopic confirmation of Compton-thick AGN activity in a subsample of these z~2 mid-IR bright galaxies. Using the four objects that lie in the CDF-N field, we estimate the space-density of Compton-thick quasars [Phi~(0.7-2.5)x10^-5 Mpc^-3 for L_X>10^44 erg/s objects at z~2-2.5] and show that Compton-thick accretion is as ubiquitious as unobscured accretion in the distant Universe.

(Abridged) We have discovered a strong DLA coincident in redshift with the faint QSO Q2343-BX415 (R = 20.2, z_em = 2.57393). Follow-up observations at intermediate spectral resolution reveal that the metal lines associated with this 'proximate' DLA consist of two sets of absorption components. One set is moving towards the quasar with velocities of ~ 150-600 km/s; this gas is highly ionized and does not fully cover the continuum source, suggesting that it is physically close to the active nucleus. The other, which accounts for most of the neutral gas, is blueshifted relative to the QSO, with the strongest component at ~ -160 km/s. We consider the possibility that the PDLA arises in the outflowing interstellar medium of the host galaxy of Q2343-BX415, an interpretation supported by strong C IV and N V absorption at nearby velocities, and by the intense radiation field longward of the Lyman limit implied by the high C II*/H I ratio. If Q2343-BX415 is the main source of these UV photons, then the PDLA is located at either ~ 8 or ~ 37 kpc from the active nucleus. Alternatively, the absorber may be a foreground star-forming galaxy unrelated to the quasar and coincidentally at the same redshift, but our deep imaging and follow-up spectroscopy of the field of Q2343-BX415 has not yet produced a likely candidate. We measure the abundances of 14 elements in the PDLA, finding an overall metallicity of ~ 1/5 solar and a normal pattern of relative element abundances for this metallicity. Thus, in this PDLA there is no evidence for the super-solar metallicities that have been claimed for some proximate, high ionization, systems.

We use a sample of rest-frame UV selected and spectroscopically observed galaxies at redshifts 1.9<z<3.4, combined with ground-based spectroscopic H-alpha and Spitzer MIPS 24 micron data, to derive the most robust measurements of the rest-frame UV, H-alpha, and infrared (IR) luminosity functions (LFs) at these redshifts. Our sample is by far the largest of its kind, with over 2000 spectroscopic redshifts in the range 1.9<z<3.4 and ~15000 photometric candidates in 29 independent fields covering a total area of almost a square degree. Our method for computing the LFs takes into account a number of systematic effects, including photometric scatter, Ly-alpha perturbations to optical colors, and contaminants. Taking into account the latter, we find no evidence for an excess of UV-bright galaxies over what was inferred in early z~3 LBG studies. The UV LF appears to undergo little evolution between z~4 and z~2. Corrected for extinction, the UV luminosity density (LD) at z~2 is at least as large as the value at z~3 and a factor of ~9 larger than the value at z~6, primarily reflecting an increase in the number density of bright galaxies between z~6 and z~2. Our analysis yields the first constraints anchored by extensive spectroscopy on the IR and bolometric LFs for faint and moderately luminous (L[bol]<10^12 L_sun) galaxies. Adding the IR to the emergent UV luminosity, incorporating independent measurements of the LD from ULIRGs, and assuming realistic dust attenuation values for UV-faint galaxies, indicates that galaxies with L[bol]<10^12 L_sun account for ~80% of the bolometric LD and SFRD at z~2-3. Our multi-wavelength constraints on the global SFRD indicate that approximately one-third of the present-day stellar mass density was formed in sub-ultraluminous galaxies between redshifts z=1.9-3.4. [Abridged]

We present results from the first robust investigation of galaxy morphology as a function of environment at z>1.5. Our study is motivated by the fact that star-forming galaxies contained within a protocluster at z=2.3 in the HS1700+64 field have significantly older ages and larger stellar masses on average than those at similar redshifts but more typical environmental densities. In the analysis of HST/ACS images, we apply non-parametric statistics to characterize the rest-frame UV morphologies of a sample of 85 UV-selected star-forming galaxies at z=1.7-2.9, 22 of which are contained in the protocluster. The remaining 63 control-sample galaxies are not in the protocluster but have a similar mean redshift of <z>~2.3. We find no environmental dependence for the distributions of morphological properties. Combining the measured morphologies with the results of population synthesis modeling, we find only weak correlations, if any, between morphological and stellar population properties such as stellar mass, age, extinction and star-formation rate. Given the incomplete census of the protocluster galaxy population, and the lack of correlation between rest-frame UV morphology and star-formation history at z~2 within our sample, the absence of environmental trends in the distribution of morphological properties is not surprising. Additionally, using a larger sample of photometric candidates, we compare morphological distributions for 282 UV-selected and 43 near-IR-selected galaxies. While the difference in the degree of nebulosity between the two samples appears to be a byproduct of the fainter average rest-frame UV surface brightness of the near-IR-selected galaxies, we find that, among the lowest surface brightness galaxies, the near-IR-selected objects have significantly smaller angular sizes (abridged).

We present the rest-frame 8 micron luminosity function (LF) at redshifts z=1 and ~2, computed from Spitzer 24 micron-selected galaxies in the GOODS fields over an area of 291 sq. arcmin. Using classification criteria based on X-ray data and IRAC colours, we identify the AGN in our sample. The rest-frame 8 micron LF for star-forming galaxies at redshifts z=1 and ~2 have the same shape as at z~0, but with a strong positive luminosity evolution. The number density of star-forming galaxies with log_{10}(nu L_nu(8 micron))>11 increases by a factor >250 from redshift z~0 to 1, and is basically the same at z=1 and ~2. The resulting rest-frame 8 micron luminosity densities associated with star formation at z=1 and ~2 are more than four and two times larger than at z~0, respectively. We also compute the total rest-frame 8 micron LF for star-forming galaxies and AGN at z~2 and show that AGN dominate its bright end, which is well-described by a power-law. Using a new calibration based on Spitzer star-forming galaxies at 0<z<0.6 and validated at higher redshifts through stacking analysis, we compute the bolometric infrared (IR) LF for star-forming galaxies at z=1 and ~2. We find that the respective bolometric IR luminosity densities are (1.2+/-0.2) x 10^9 and (6.6^{+1.2}_{-1.0}) x 10^8 L_sun Mpc^{-3}, in agreement with previous studies within the error bars. At z~2, around 90% of the IR luminosity density associated with star formation is produced by luminous and ultraluminous IR galaxies (LIRG and ULIRG), with the two populations contributing in roughly similar amounts. Finally, we discuss the consistency of our findings with other existing observational results on galaxy evolution.

Motivated by the irregular and little-understood morphologies of z ~ 2 - 3 galaxies, we use non-parametric coefficents to quantify the morphologies of 216 galaxies which have been spectroscopically confirmed to lie at redshifts z = 1.8 - 3.4 in the GOODS-N field. Using measurements of ultraviolet (UV) and optical spectral lines, multi-band photometric data, and stellar population models we statistically assess possible correlations between galaxy morphology and physical observables such as stellar mass, star formation rate, and the strength of galaxy-scale outflows. We find evidence that dustier galaxies have more nebulous UV morphologies and that larger, more luminous galaxies may drive stronger outflows, but otherwise conclude that UV morphology is either statistically decoupled from the majority of physical observables or determined by too complex a combination of physical processes to provide characterizations with predictive power. Given the absence of strong correlations between UV morphology and physical parameters such as star formation rates, we are therefore unable to support the hypothesis that morphologically irregular galaxies predominantly represent major galaxy mergers. Comparing galaxy samples, we find that IR-selected BzK galaxies and radio-selected submillimeter galaxies (SMGs) have UV morphologies similar to the optically selected sample, while distant red galaxies (DRGs) are more nebulous.

We present the results of a spectroscopic survey with the Keck I telescope of more than 280 star-forming galaxies and AGN at redshifts 1.4<z<3.0 in the GOODS-N field. Candidates are selected by their UGR colors using the ``BM/BX'' criteria to target redshift 1.4<z<2.5 galaxies and the LBG criteria to target redshift z~3 galaxies; combined these samples account for ~25-30% of the R and K-band counts to R=25.5 and K(AB)=24.4, respectively. The sample of 212 BM/BX galaxies and 74 LBGs is presently the largest spectroscopic sample of galaxies at z>1.4 in GOODS-N. Extensive multi-wavelength data allow us to investigate the stellar populations, stellar masses, bolometric luminosities (L_bol), and extinction of z~2 galaxies. Deep Chandra and Spitzer data indicate that the sample includes galaxies with a wide range in L_bol, from 10^10 L_sun to >10^12 L_sun, and covering 4 orders of magnitude in dust obscuration (L_bol/L_UV). The sample includes galaxies with a large dynamic range in evolutionary state, from very young galaxies (ages <50 Myr) with small stellar masses (M* ~ 10^9 M_sun) to evolved galaxies (ages >2 Gyr) with stellar masses comparable to the most massive galaxies at these redshifts (M* > 10^11 M_sun). Spitzer data indicate that the optical sample includes some fraction of the obscured AGN population at high redshifts: at least 3 of 11 AGN in the z>1.4 sample are undetected in the deep X-ray data but exhibit power-law SEDs longward of ~2 micron (rest-frame) indicative of obscured AGN. The results of our survey indicate that rest-frame UV selection and spectroscopy presently constitute the most time-wise efficient method of culling large samples of high redshift galaxies with a wide range in intrinsic properties, and the data presented here will add significantly to the multi-wavelength legacy of the GOODS survey. [Abridged]

Using H-alpha spectra of 114 rest-frame UV-selected galaxies at z~2, we compare inferred star formation rates (SFRs) with those determined from the UV continuum luminosity. After correcting for extinction using standard techniques based on the UV continuum slope, we find excellent agreement between the indicators, with <SFR_Ha> = 31 Msun/yr and <SFR_UV> = 29 Msun/yr. The agreement between the indicators suggests that the UV luminosity is attenuated by an typical factor of ~4.5 (with a range from no attenuation to a factor of ~100 for the most obscured object in the sample), in good agreement with estimates of obscuration from X-ray, radio and mid-IR data. The H-alpha luminosity is attenuated by a factor of ~1.7 on average, and the maximum H-alpha attenuation is a factor of ~5. In agreement with X-ray and mid-IR studies, we find that the SFR increases with increasing stellar mass and at brighter K magnitudes, to <SFR_Ha> ~ 60 Msun/yr for galaxies with K<20; the correlation between K magnitude and SFR is much stronger than the correlation between stellar mass and SFR. All galaxies in the sample have SFRs per unit area Sigma_SFR in the range observed in local starbursts. We compare the instantaneous SFRs and the past average SFRs as inferred from the ages and stellar masses, finding that for most of the sample, the current SFR is an adequate representation of the past average. There is some evidence that the most massive galaxies (M_star >10^11 Msun) have had higher SFRs in the past.

We present analysis of the near-infrared spectra of 114 rest-frame UV-selected star-forming galaxies at z~2. By combining the H-alpha spectra with photometric measurements from observed 0.3--8 micron, we assess the relationships between kinematics, dynamical masses, inferred gas fractions, and stellar masses and ages. The H-alpha line widths give a mean dynamical mass M_dyn=6.9 +/- 0.6 x 10^10 Msun within a typical radius of ~6 kpc, after excluding AGN. The average dynamical mass is ~2 times larger than the average stellar mass, and the two agree to within a factor of several for most objects. However, ~15% of the sample has M_dyn >> M_star. These objects are best fit by young stellar populations and tend to have high H-alpha equivalent widths, suggesting that they are young starbursts with large gas masses. Rest-frame optical luminosity and velocity dispersion are correlated with 4 sigma significance. Using the local empirical correlation between star formation rate per unit area and gas surface density, we estimate the mass of the gas associated with star formation, and find a mean gas fraction of ~50% and a strong decrease in gas fraction with increasing stellar mass. The masses of gas and stars combined are considerably better correlated with the dynamical masses than are the stellar masses alone, and agree to within a factor of three for 85% of the sample. The combination of kinematic measurements, estimates of gas masses, and stellar population properties suggest that the factor of ~500 range in stellar mass across the sample cannot be fully explained by intrinsic differences in the total masses of the galaxies, which vary by a factor of ~40; the remaining variation is due to the evolution of the stellar population and the conversion of gas into stars. [Abridged]

Using very deep Spitzer/MIPS 24 micron observations, we present an analysis of the bolometric luminosities (L[bol]) and UV extinction properties of more than 200 spectroscopically identified, optically selected (UGR) z~2 galaxies in the GOODS-N field. The large spectroscopic sample is supplemented with near-IR selected (BzK/DRG) galaxies and submm sources at similar redshifts in the same field, providing a representative collection of relatively massive (M*>1e10 Msun) galaxies at high redshifts. We focus on the redshift range 1.5-2.6, where MIPS is sensitive to the strength of the mid-IR PAH features in the galaxy spectra (rest-frame 5-8.5 micron). We demonstrate, using stacked X-ray data and a subset of galaxies with H-alpha measurements, that L(5-8.5) provides a reliable estimate of L(IR) for most star forming galaxies at z~2. The range of L(IR) in the samples considered extends from ~1e10 to >1e12 Lsun, with a mean of 2e11 Lsun. Using 24 micron observations to infer dust extinction in high redshift galaxies, we find that, as in the local universe, the obscuration (L[IR]/L[1600]) is strongly dependent on L(bol), and ranges in value from <1 to \~1000. However, the obscuration is ~10 times smaller at a given L(bol) at z~2 than at z=0. We show that the values of L(IR) and obscuration inferred from the UV spectral slope beta generally agree well with the values inferred from L(5-8.5) for L(bol)<1e12 Lsun. For galaxies with L(bol)>1e12 Lsun, it is common for UV-based estimates to underpredict L(IR) by a factor of ~10-100. Using the specific SFR as a proxy for cold gas fraction, we find a wide range in the evolutionary state of galaxies at z~2, from those which have just begun to form stars to those which have already accumulated most of their stellar mass and are about to become, or already are, passively-evolving. [Abridged]

We use a sample of 87 rest-frame UV-selected star-forming galaxies with mean spectroscopic redshift z=2.26 to study the correlation between metallicity and stellar mass at high redshift. Using stellar masses determined from SED fitting to 0.3-8 micron photometry, we divide the sample into six bins in stellar mass, and construct six composite H-alpha+[NII] spectra from all of the objects in each bin. We estimate the mean oxygen abundance in each bin from the [NII]/H-alpha ratio, and find a monotonic increase in metallicity with increasing stellar mass, from 12+log(O/H) < 8.2 for galaxies with <M_star> = 2.7e9 Msun to 12+log(O/H) = 8.6 for galaxies with <M_star> = 1e11 Msun. We use the empirical relation between star formation rate density and gas density to estimate the gas fractions of the galaxies, finding an increase in gas fraction with decreasing stellar mass. These gas fractions combined with the observed metallicities allow the estimation of the effective yield y_eff as a function of stellar mass; in constrast to observations in the local universe which show a decrease in y_eff with decreasing baryonic mass, we find a slight increase. Such a variation of metallicity with gas fraction is best fit by a model with supersolar yield and an outflow rate ~4 times higher than the star formation rate. We conclude that the mass-metallicity relation at high redshift is driven by the increase in metallicity as the gas fraction decreases through star formation, and is likely modulated by metal loss from strong outflows in galaxies of all masses. There is no evidence for preferential loss of metals from low mass galaxies as has been suggested in the local universe. [Abridged]

We present the discovery of the optical afterglow of GRB 050505 and an optical absorption spectrum obtained with the Keck I 10-m telescope. The spectrum exhibits three redshifted absorption systems with the highest, at z=4.2748, arising in the GRB host galaxy. The host absorption system is marked by a damped Ly-alpha (DLA) feature with a neutral hydrogen column density of logN(HI)=22.05+/-0.10, higher than that of any QSO-DLA detected to date, but similar to several other recent measurements from GRB spectra. In addition, we detect absorption lines from both low- and high-ionization species from which we deduce a metallicity, Z~0.06 Z_solar, with a depletion pattern that is roughly similar to that of the Galactic warm halo, warm disk, or disk+halo. More importantly, we detect strong absorption from SiII* indicating a dense environment, n_H>10^2 cm^-3, in the vicinity of the burst, with a size of about 4 pc. In addition, the CIV absorption system spans a velocity range of about 1000 km/s, which is not detected in any other absorption feature. We show that the most likely interpretation for this wide velocity range is absorption in the wind from the progenitor star. In this context, the lack of corresponding SiIV absorption indicates that the progenitor had a mass of <25 M_solar and a metallicity <0.1 Z_solar, and therefore required a binary companion to eject its hydrogen envelope prior to the GRB explosion. Finally, by extending the GRB-DLA sample to z~4.3 we show that these objects appear to follow a similar metallicity-redshift relation as in QSO-DLAs, but with systematically higher metallicities. It remains to be seen whether this trend is simply due to the higher neutral hydrogen columns in GRB-DLAs, or if it is a manifestation of different star formation properties in GRB-DLAs. [abridged]

We have detected Lyman-alpha emission from a damped Lyman-alpha system (DLA) that lies near the bright quasar HS1549+1919. The DLA has the same redshift as HS1549+1919 and was discovered in the spectrum of a faint QSO that lies 49" away (380 proper kpc). The emission line's luminosity, double-peaked profile, and small spatial separation from the DLA suggest that it may be fluorescent Lyman-alpha emission from gas that is absorbing the nearby QSO's radiation. If this is the case, our observations show that the DLA has a size of at least 1.5" and that the QSO's luminosity one million years ago was similar to its luminosity today. A survey for similar systems within 1' of bright QSOs would put interesting limits on the mean quasar lifetime.

Using the extensive multi-wavelength data in the GOODS-North field, we construct and draw comparisons between samples of optical and near-IR selected star-forming and passively evolving galaxies at redshifts 1.4<z<2.6. We find overlap at the 70-80% level in samples of z~2 star-forming galaxies selected by their optical (UGR) and near-IR (BzK) colors when subjected to common K-band limits. Deep Chandra data indicate a 25% AGN fraction among near-IR selected objects, much of which occurs among near-IR bright objects (K<20; Vega). Stacking the X-ray emission for the star-forming galaxies, we find the SFR distributions of UGR, BzK, and J-K>2.3 galaxies (Distant Red Galaxies; DRGs) are very similar as a function of K, with K<20 galaxies having <SFR>~120 Msun/yr, a factor of two to three higher than those with K>20.5. The absence of X-ray emission from the reddest DRGs and BzK galaxies with z-K>3 indicates they must have declining star formation histories to explain their red colors and low SFRs. While the M/L ratio of passively-evolving galaxies may be larger on average, the Spitzer/IRAC data indicate that their inferred stellar masses do not exceed the range spanned by optically selected galaxies, suggesting that the disparity in current SFR may not indicate a fundamental difference between optical and near-IR selected massive galaxies (M* > 10^11 Msun). We consider the contribution of UGR, BzK, DRG, and submillimeter galaxies (SMGs) to the SFRD at z~2, taking into account sample overlap. The total SFRD in the interval 1.4<z<2.6 of UGR and BzK galaxies to K=22, and DRGs to K=21 is ~ 0.10+/-0.02 Msun/yr/Mpc^3. Optically-selected galaxies to R=25.5 and K=22 account for ~70% of this total. Of the radio-selected SMGs to S(850 micron)~4 mJy with redshifts 1.4<z<2.6, >80% could be selected by the UGR, BzK, and/or DRG criteria.

Absorption-line spectroscopy of 23 background QSOs and numerous background galaxies has let us measure the spatial distribution of metals and neutral hydrogen around 1044 UV-selected galaxies at redshifts 1.8<z<3.3. The typical galaxy is surrounded to radii r~40 proper kpc by gas that has a large velocity spread (dv>260 km/s) and produces very strong absorption lines (N_CIV >> 10^14 cm^-2) in the spectra of background objects. Absorption with an average column density of N_CIV ~ 10^14 cm^-2 extends to 80 kpc, a radius large enough to imply that most strong intergalactic CIV absorption is associated with star-forming galaxies like those in our sample. We find that the galaxy-CIV cross-correlation length increases with CIV column density and is similar to the galaxy-galaxy length (r_0 ~ 4 h^-1 Mpc) for N_CIV > 10^12.5 cm^-2. Distortions in the redshift-space galaxy-CIV correlation function on small scales may imply that some of the CIV systems have large peculiar velocities. Four of the five detected OVI absorption systems in our sample lie within 400 proper kpc of a known galaxy. Strong Lyman-a absorption is produced by the intergalactic gas within 1 h^-1 comoving Mpc of most galaxies, but for a significant minority (~1/3) the absorption is weak or absent. We were unable to identify any statistically significant differences between galaxies with weak nearby HI absorption and the rest, although galaxies with weak absorption may have higher star-formation rates. Galaxies near intergalactic CIV systems appear to reside in relatively dense environments and to have distinctive spectral energy distributions that are characterized by blue colors and young ages. (abridged)

We have discovered a highly significant over-density of galaxies at z=2.300+/-0.015 in the course of a redshift survey designed to select star-forming galaxies in the redshift range z=2.3+/-0.4 in the field of the bright z=2.72 QSO HS1700+643. The structure has a redshift-space galaxy over-density of delta_g,z ~= 7 and an estimated matter over-density in real space of delta_m ~= 1.8, indicating that it will virialize by z~0 with a mass scale of ~= 1.4x10^15 M_sun, that of a rich galaxy cluster. Detailed modeling of the spectral energy distribution -- from the rest-far-UV to the rest-near-IR -- of the 72 spectroscopically confirmed galaxies in this field for which we have obtained K_s and Spitzer/IRAC photometry, allows for a first direct comparison of galaxy properties as a function of large-scale environment at high redshift. We find that galaxies in the proto-cluster environment have mean stellar masses and inferred ages that are ~2 times larger (at z=2.30) than identically UV-selected galaxies outside of the structure, and show that this is consistent with simple theoretical expectations for the acceleration of structure formation in a region that is over-dense on large scales by the observed amount. The proto-cluster environment contains a significant number of galaxies that already appear old, with large stellar masses (>10^11 M_sun), by z=2.3.

We present the broad-band UV through mid-infrared spectral energy distributions (SEDs) of a sample of 72 spectroscopically-confirmed star-forming galaxies at z=2.30+/-0.3. Located in a 72 arcmin-squared field centered on the bright background QSO, HS1700+643, these galaxies were pre-selected to lie at z~2 based solely on their rest-frame UV colors and luminosities, and should be representative of UV-selected samples at high redshift. In addition to deep ground-based photometry spanning from 0.35-2.15 microns, we make use of Spitzer/IRAC data, which probes the rest-frame near-IR at z~2. The range of stellar populations present in the sample is investigated with simple, single-component stellar population synthesis models. Emphasizing stellar mass estimates, which are less subject to systematic uncertainties than other parameters, we find <log(M*/Msun)>=10.32+/-0.51 for the sample. Allowing for the possibility of episodic star formation, we find that typical galaxies in our sample could contain up to three times more stellar mass in an old underlying burst than what was inferred from single-component modeling. In contrast, mass estimates for the most massive galaxies in the sample (M*>10^11 Msun) are fairly insensitive to the method used to model the stellar population. Galaxies in this massive tail, which are also the oldest objects in the sample, could plausibly evolve into the passive galaxies discovered at z~1.7 with near-IR selection techniques. In the general framework of hierarchical galaxy formation and mergers, which implies episodic star-formation histories, galaxies at high redshift may pass in and out of UV-selected and near-IR color-selected samples as they evolve from phases of active star formation to quiescence and back again (Abridged).

We present an analysis of the X-ray emission from a large sample of ultraviolet (UV) selected, star forming galaxies with 0.74<z<1.32 in the Hubble Deep Field North (HDF-N) region. By excluding all sources with significant detected X-ray emission in the 2 Ms Chandra observation we are able to examine the properties of galaxies for which the emission in both UV and X-ray is expected to be predominantly due to star formation. Stacking the X-ray flux from 216 galaxies in the soft and hard bands produces significant detections. The derived mean 2-10 keV rest-frame luminosity is 2.97+/-0.26x10^(40) erg/s, corresponding to an X-ray derived star formation rate (SFR) of 6.0+/-0.6 Msolar/yr. Comparing the X-ray value with the mean UV derived SFR, uncorrected for attenuation, we find that the average UV attenuation correction factor is \~3. By binning the galaxy sample according to UV magnitude and colour, correlations between UV and X-ray emission are also examined. We find a strong positive correlation between X-ray emission and rest-frame UV emission. A correlation between the ratio of X-ray-to-UV emission and UV colour is also seen, such that L(X)/L(UV) increases for redder galaxies. Given that X-ray emission offers a view of star formation regions that is relatively unaffected by extinction, results such as these can be used to evaluate the effects of dust on the UV emission from high-z galaxies. For instance we derive a relationship for estimating UV attenuation corrections as a function of colour excess. The observed relation is inconsistent with the Calzetti et al. (2000) reddening law which over predicts the range in UV attenuation corrections by a factor of ~100 for the UV selected z~1 galaxies in this sample (abridged).

We obtained deep 8.5'x8.5' near-infrared images within four high-redshift survey fields, measured the Ks magnitudes of 300 optically selected galaxies with spectroscopic redshift 1.8<z<2.6, and compared the spatial clustering strength of galaxies with Ks<20.5 and Ks>20.5. We found at greater than 95% confidence that the brighter galaxies cluster more strongly. The best-fit correlation lengths for the bright and faint samples are 10+-3 and 4+-0.8 comoving Mpc/h, respectively (1sigma), although the unusual density of bright QSOs in one of our survey fields may imply that the result is not representative of the universe as a whole. Neglecting this possibility, the correlation length for the optically selected sample with Ks<20.5 agrees well with that reported for comparably bright near-IR-selected samples. The differences in correlation length between optically selected and near-IR-selected samples have been presented as evidence that the two techniques find orthogonal populations of high-redshift galaxies. Our results favor a more nuanced view.

We analyzed the spatial distribution of 28500 photometrically selected galaxies with AB magnitude 23.5<R<25.5 and redshift 1.4<z<3.5 in 21 fields with a total area of 0.81 square degrees. The galaxies were divided into three subsamples, with mean redshifts z=1.7, 2.2, 2.9, according to the UGR selection criteria of Adelberger et al. (2004) and Steidel et al. (2003). Combining the galaxies' measured angular clustering with redshift distributions inferred from 1600 spectroscopic redshifts, we find comoving correlation lengths at the three redshifts of r_0 = 4.5+-0.6, 4.2+-0.5, and 4.0+-0.6 h^-1 Mpc, respectively, and infer a roughly constant correlation function slope of gamma=1.6+-0.1. We derive similar numbers from the 1600 object spectroscopic sample itself with a new statistic, K, that is insensitive to many possible systematics. Galaxies that are bright in R (lambda_rest ~ 1500-2500A) cluster more strongly than fainter galaxies at z=2.9 and z=2.2 but not, apparently, at z=1.7. Comparison to a numerical simulation that is consistent with recent WMAP observations suggests that galaxies in our samples are associated with dark matter halos of mass ~ 10^12, and that a small fraction of the halos contain more than one galaxy that satisfies our selection criteria. Adding recent observations of galaxy clustering at z~0 and z~1 to the simulation results, we conclude that the typical object in our samples will evolve into an elliptical galaxy by redshift z=0 and will already have an early-type spectrum by redshift z=1. We comment briefly on the implied relationship between galaxies in our survey and those selected with other techniques. (slightly abridged)

As part of the In-Orbit Checkout activities for the Spitzer Space Telescope, the IRAC team carried out a deep observation (average integration time ~8 hours) of a field surrounding the bright QSO HS 1700+6416. This field contains several hundred z~3 Lyman-break galaxy candidates, and we report here on their mid-infrared properties, including the IRAC detection rate, flux densities and colors, and the results of fitting population synthesis models to the optical, near-infrared, and IRAC magnitudes. The results of the model-fitting show that previous optical/near-infrared studies of LBGs were not missing large, hidden old stellar populations. The LBG candidates' properties are consistent with those of massive, star-forming galaxies at z~3. Other IRAC sources in the same field have similar properties, so IRAC selection may prove a promising method of finding additional high-redshift galaxies.

Few galaxies have been found between the redshift ranges z < ~1 probed by magnitude-limited surveys and z > ~3 probed by Lyman-break surveys. Comparison of galaxy samples at lower and higher redshift suggests that large numbers of stars were born and the Hubble sequence began to take shape at the intermediate redshifts 1<z<3, but observational challenges have prevented us from observing the process in much detail. We present simple and efficient strategies that can be used to find large numbers of galaxies throughout this important but unexplored redshift range. All the strategies are based on selecting galaxies for spectroscopy on the basis of their colors in ground-based images taken through a small number of optical filters: GRi for redshifts 0.85<z<1.15, GRz for 1<z<1.5, and UGR for 1.4<z<2.1 and 1.9<z<2.7. The performance of our strategies is quantified empirically through spectroscopy of more than 2000 galaxies at 1<z<3.5. We estimate that more than half of the UV-luminosity density at 1<z<3 is produced by galaxies that satisfy our color-selection criteria. Our methodology is described in detail, allowing readers to devise analogous selection criteria for other optical filter sets.

We have examined the stacked radio and X-ray emission from UV-selected galaxies spectroscopically confirmed to lie between redshifts 1.5 < z < 3.0 in the GOODS-North field to determine their average extinction and star formation rates (SFRs). The X-ray and radio data are obtained from the Chandra 2 Msec survey and the Very Large Array, respectively. There is a good agreement between the X-ray, radio, and de-reddened UV estimates of the average SFR for our sample of z~2 galaxies of ~50 solar masses per year, indicating that the locally-calibrated SFR relations appear to be statistically valid from redshifts 1.5 < z < 3.0. We find that UV-estimated SFRs (uncorrected for extinction) underestimate the bolometric SFRs as determined from the 2-10 keV X-ray luminosity by a factor of ~4.5 to 5.0 for galaxies over a large range in redshift from 1.0 < z < 3.5.

We present the first results of a large-scale survey, using the UV/blue sensitive LRIS-B spectrograph on the Keck I telescope, of galaxies in the redshift interval 1.4<z<2.5, often called the `redshift desert' because of historical difficulties in spectroscopically identifying galaxies in that range. We introduce two samples of star forming galaxies, `BX' galaxies at <z>=2.20+/-0.32 and `BM' galaxies at <z>=1.70+/-0.34; currently we have spectroscopically confirmed 749 of the former and 114 of the latter. We also present initial results of deep near-IR photometry and spectroscopy, from which we show that z~2 galaxies are significantly redder in their optical/IR colors compared to similarly UV-selected galaxies at z~3, but that the characteristics of their galaxy-scale outflows are quite similar. We illustrate by example the information which can be deduced on the stellar populations, metallicities, and kinematics of `redshift desert' galaxies from easily accessible rest-frame far-UV and rest-frame optical spectra. Far from being hostile to observations, the universe at z~2 is uniquely suited to providing information on the astrophysics of star-forming galaxies and the intergalactic medium, and the relationship between the two [abridged].

We describe results of optical and near-IR observations of a large spectroscopic sample of star-forming galaxies photometrically-selected to lie in the redshift range 1.4 < z < 2.5, often called the ``redshift desert'' because of historical difficulty in obtaining spectroscopic redshifts in this range. We show that the former ``redshift desert'' is now very much open to observation.

We report on deep near-infrared (NIR) observations of submillimeter-selected galaxies (SMGs) with the Near Infrared Camera (NIRC) on the Keck I telescope. We have identified K-band candidate counterparts for 12 out of 15 sources in the SCUBA Cluster Lens Survey. Three SMGs remain non-detections with K-band limits of K>23 mag, corrected for lensing. Compensating for lensing we find a median magnitude of K=22+/-1 mag for the SMG population, but the range of NIR flux densities spans more than a factor of 400. For SMGs with confirmed counterparts based on accurate positions from radio, CO, and/or millimeter continuum interferometric observations, the median NIR color is J-K=2.6+/-0.6 mag. The NIR-bright SMGs (K<19 mag) have colors of J-K =~ 2 mag, while the faint SMGs tend to be extremely red in the NIR (J-K>3 mag). We argue that a color selection criterion of J-K>~3 mag can be used to help identify counterparts of SMGs that are undetected at optical and radio wavelengths. The number density of sources with J-K>3 mag is 5 arcmin^{-2} at K<22.5 mag, greater than that of SMGs with S(850um)>2 mJy. It is not clear if the excess represents less luminous infrared-bright galaxies with S(850um)<~2 mJy, or if the faint extremely red NIR galaxies represent a different population of sources that could be spatially related to the SMGs.

We report new K-band, radio continuum, and CO (1-0) imaging observations and 850 micron photometric observations of PDS 456, the most luminous QSO in the local universe (z<0.3). The 0.6'' resolution K-band image obtained using the Keck telescope shows three compact K~16.5 (M(K)~ -21) sources at a projected distance of ~10 kpc to the southwest, and the host galaxy of PDS 456 may be interacting or merging with one or more companions. The observations using the OVRO millimeter array has revealed a narrow CO (1-0) line (FWHM = 181 km/s) centered at z=0.1849, and 9 x 10^9 solar mass of molecular gas mass is inferred. Radio continuum luminosity is nearly an order of magnitude larger than expected from its FIR luminosity, and the radio source, unresolved by the 2" beam of the VLA, is dominated by the AGN activity. Our 850 micron photometric observations suggest that the cold dust content of the host galaxy is less than one half of the amount in Arp 220. Its SED has both a QSO-like and a ULIRG-like nature, and the observed IR, X-ray, and gas properties suggest that the AGN activity dominates its luminosity. ULIRGs and IR QSOs form a broad continuous track in the ``star formation efficiency'' plot in the manner consistent with the ULIRG-QSO transition scenario, relating the evolution in the dust processed luminosity with the available fuel (gas and dust) supply. However, the location of PDS 456 is clearly offset from the apparent track traced by the ULIRGs and IR QSOs on this plot, and PDS 456 appears to be a rare, exceptional object, and the duration of the physical process governing its present properties must be short compared with the length of the luminous QSO phase.

We have studied the radio and far-infrared (FIR) emission from 114 galaxies in the 7 nearest clusters (<100 Mpc) with prominent X-ray emission to investigate the impact of the cluster environment on the star formation and AGN activity in the member galaxies. The X-ray selection criterion is adopted to focus on the most massive and dynamically relaxed clusters. A large majority of cluster galaxies show an excess in radio emission over that predicted from the radio-FIR correlation, the fraction of sources with radio excess increases toward cluster cores, and the radial gradient in the FIR/radio flux ratio is a result of radio enhancement. Of the radio-excess sources, 70% are early-type galaxies and the same fraction host an AGN. The galaxy density drops by a factor of 10 from the composite cluster center out to 1.5 Mpc, yet galaxies show no change in FIR properties over this region, and show no indication of mass segregation. We have examined in detail the physical mechanisms that might impact the FIR and radio emission of cluster galaxies. While collisional heating of dust may be important for galaxies in cluster centers, it appears to have a negligible effect on the observed FIR emission for our sample galaxies. The correlations between radio and FIR luminosity and radius could be explained by magnetic compression from thermal ICM pressure. We also find that simple delayed harassment cannot fully account for the observed radio, FIR, and mid-IR properties of cluster galaxies.

We report the redshift measurement for the submillimeter selected galaxy SMMJ04431+0210 (N4) using the Near Infrared Spectrograph on the Keck-II telescope. The data show H-alpha, [NII]6583,6548, and [OIII]5007 lines at a redshift of z=2.51. The high nuclear [NII]/H-alpha line ratio is consistent with a LINER or Type-II AGN. The H-alpha emission is spatially resolved, suggesting the presence of significant star-forming activity outside the nucleus. From imaging with the Near Infrared Camera on the Keck-I telescope, we find an extremely red near-infrared color of J-K=3.2 for N4. Follow-up redshifted CO(3-2) observations with the Owens Valley Millimeter Array constrain the mass of molecular gas to be less than 4x10^{10} Msun, after correcting for lensing. The CO to sub-mm flux limit, the spectroscopic line ratios, and the spectral energy distribution for N4 are all within the range of properties found in other high-redshift sub-mm sources and local ultraluminous infrared galaxies. After the correction for lensing, N4 is the weakest intrinsic sub-mm selected source with a known redshift and represents the first redshift for the <2mJy 850um sources which are responsible for the bulk of the emission from the sub-mm population as a whole. We argue that N4 contains either an AGN or LINER nucleus surrounded by an extended region of active star-formation.

The radio counterparts to the IRAS Redshift Survey galaxies are identified in the NRAO VLA Sky Survey (NVSS) catalog. Our new catalog of the infrared flux-limited (S60 >= 2 Jy) complete sample of 1809 galaxies lists accurate radio positions, redshifts, and 1.4 GHz radio and IRAS fluxes. This sample is six times larger in size and five times deeper in redshift coverage (to z ~ 0.15) compared with those used in earlier studies of the radio and far-infrared (FIR) properties of galaxies in the local volume. The well known radio-FIR correlation is obeyed by the overwhelming majority (>98%) of the infrared-selected galaxies, and the radio AGNs identified by their excess radio emission constitute only about 1% of the sample, independent of the IR luminosity. These FIR-selected galaxies can account for the entire population of late-type field galaxies in the local volume, and their radio continuum may be used directly to infer the extinction-free star formation rate in most cases. Both the 1.4 GHz radio and 60 micron infrared luminosity functions are reasonably well described by linear sums of two Schechter functions, one representing normal, late-type field galaxies and the second representing starbursts and other luminous infrared galaxies. The integrated FIR luminosity density for the local volume is 4.8+/-0.5 x 10^7 L(sun) per Mpc^3, less than 10% of which is contributed by the luminous infrared galaxies with L(FIR) > 10^{11} L(sun) . The inferred extinction-free star formation density for the local volume is 0.015 +/- 0.005 solar mass per year per Mpc^3.