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.

Ly$\alpha$ emission is widely used to detect and confirm high-redshift galaxies and characterize the evolution of the intergalactic medium. However, many galaxies do not display Ly$\alpha$ emission in typical spectroscopic observations, and intrinsic Ly$\alpha$-emitters represent a potentially biased set of high-redshift galaxies. In this work, we analyze a set of 703 galaxies at $2\lesssim z\lesssim3$ with both Ly$\alpha$ spectroscopy and measurements of other rest-frame ultraviolet and optical properties in order to develop an empirical model for Ly$\alpha$ emission from galaxies and understand how the probability of Ly$\alpha$ emission depends on other observables. We consider several empirical proxies for the efficiency of Ly$\alpha$ photon production as well as the subsequent escape of these photons through their local interstellar medium. We find that the equivalent width of metal-line absorption and the O3 ratio of rest-frame optical nebular lines are advantageous empirical proxies for Ly$\alpha$ escape and production, respectively. We develop a new quantity, $X_\mathrm{LIS}^\mathrm{O3}$, that combines these two properties into a single predictor of net Ly$\alpha$ emission, which we find describes $\sim$90% of the observed variance in Ly$\alpha$ equivalent width when accounting for our observational uncertainties. We also construct conditional probability distributions demonstrating that galaxy selection based on measurements of galaxy properties yield samples of galaxies with widely varying probabilities of net Ly$\alpha$ emission. The application of the empirical models and probability distributions described here may be used to infer the selection biases of current galaxy surveys and evaluate the significance of high-redshift Ly$\alpha$ (non-)detections in studies of reionization and the intergalactic medium.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

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

We use photoionization models designed to reconcile the joint rest-UV-optical spectra of high-z star-forming galaxies to self-consistently infer the gas chemistry and nebular ionization and excitation conditions for ~150 galaxies from the Keck Baryonic Structure Survey (KBSS), using only observations of their rest-optical nebular spectra. We find that the majority of z~2-3 KBSS galaxies are moderately O-rich, with an interquartile range in 12+log(O/H)=8.29-8.56, and have significantly sub-solar Fe enrichment, with an interquartile range of [Fe/H]=[-0.79,-0.53], contributing additional evidence in favor of super-solar O/Fe in high-z galaxies. Model-inferred ionization parameter and N/O are strongly correlated with common strong-line indices (such as O32 and N2O2), with the latter exhibiting similar behavior to local extragalactic HII regions. In contrast, diagnostics commonly used for measuring gas-phase O/H (such as N2 and O3N2) show relatively large scatter with the overall amount of oxygen present in the gas and behave differently than observed at z~0. We provide a new calibration for using R23 to measure O/H in typical high-z galaxies, although it is most useful for relatively O-rich galaxies; combining O32 and R23 does not yield a more effective calibration. Finally, we consider implications for the intrinsic correlations between physical conditions across the galaxy sample and find that N/O varies with O/H in high-z galaxies in a manner almost identical to local HII regions. However, we do not find a strong anti-correlation between ionization parameter and metallicity (O/H or Fe/H) in high-z galaxies, which is one of the principal bases for using strong-line ratios to infer oxygen abundance.

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 present Keck/OSIRIS laser guide-star assisted adaptive optics (LGSAO) integral field spectroscopy of [O III] 5007 Angstrom nebular emission from twelve galaxies hosting optically faint (R = 20 - 25; nuLnu ~ 10^44 -10^46 erg/s) active galactic nuclei (AGN) at redshift z ~ 2 - 3. In combination with deep Hubble Space Telescope Wide Field Camera 3 rest-frame optical imaging, Keck/MOSFIRE rest-optical spectroscopy, and Keck/KCWI rest-UV integral field spectroscopy we demonstrate that both the continuum and emission-line structure of these sources exhibit a wide range of morphologies from compact isolated point sources to double-AGN merging systems with extensive ~ 50 kpc tidal tails. One of the twelve galaxies previously known to exhibit a proximate damped LyA system coincident in redshift with the galaxy shows evidence for both an extended [O III] narrow-line emission region and spatially offset LyA emission (with morphologically distinct blueshifted and redshifted components) indicative of large scale gas flows photoionized by the central AGN. We do not find widespread evidence of star formation in the host galaxies surrounding these AGN; the [O III] velocity dispersions tend to be high (sigma = 100 - 500 km/s), the continuum morphologies are much more compact than a mass-matched star forming comparison sample, and the diagnostic nebular emission line ratios are dominated by an AGN-like ionizing spectrum. The sample is most consistent with a population of AGN that radiate at approximately their Eddington limit and photoionize extended [O III] nebulae whose characteristic sizes scale approximately as the square root of the AGN luminosity.

Massive galaxy clusters are now found as early as 3 billion years after the Big Bang, containing stars that formed at even earlier epochs. The high-redshift progenitors of these galaxy clusters, termed 'protoclusters', are identified in cosmological simulations with the highest dark matter overdensities. While their observational signatures are less well defined compared to virialized clusters with a substantial hot intra-cluster medium (ICM), protoclusters are expected to contain extremely massive galaxies that can be observed as luminous starbursts. Recent claimed detections of protoclusters hosting such starbursts do not support the kind of rapid cluster core formation expected in simulations because these structures contain only a handful of starbursting galaxies spread throughout a broad structure, with poor evidence for eventual collapse into a protocluster. Here we report that the source SPT2349-56 consists of at least 14 gas-rich galaxies all lying at z = 4.31 based on sensitive observations of carbon monoxide and ionized carbon. We demonstrate that each of these galaxies is forming stars between 50 and 1000 times faster than our own Milky Way, and all are located within a projected region only $\sim$ 130 kiloparsecs in diameter. This galaxy surface density is more than 10 times the average blank field value (integrated over all redshifts) and $>$1000 times the average field volume density. The velocity dispersion ($\sim$ 410 km s$^{-1}$) of these galaxies and enormous gas and star formation densities suggest that this system represents a galaxy cluster core at an advanced stage of formation when the Universe was only 1.4 billion years old. A comparison with other known protoclusters at high redshifts shows that SPT2349-56 is a uniquely massive and dense system that could be building one of the most massive structures in the Universe today.

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.

On 2017 August 17, gravitational waves were detected from a binary neutron star merger, GW170817, along with a coincident short gamma-ray burst, GRB170817A. An optical transient source, Swope Supernova Survey 17a (SSS17a), was subsequently identified as the counterpart of this event. We present ultraviolet, optical and infrared light curves of SSS17a extending from 10.9 hours to 18 days post-merger. We constrain the radioactively-powered transient resulting from the ejection of neutron-rich material. The fast rise of the light curves, subsequent decay, and rapid color evolution are consistent with multiple ejecta components of differing lanthanide abundance. The late-time light curve indicates that SSS17a produced at least ~0.05 solar masses of heavy elements, demonstrating that neutron star mergers play a role in r-process nucleosynthesis in the Universe.

We study HI and metal-line absorption around $z\approx2$ star-forming galaxies by comparing an analysis of data from the Keck Baryonic Structure Survey to mock spectra generated from the EAGLE cosmological, hydrodynamical simulations. We extract sightlines from the simulations and compare the properties of the absorption by HI, CIV and SiIV around simulated and observed galaxies using pixel optical depths. We mimic the resolution, pixel size, and signal-to-noise ratio of the observations, as well as the distributions of impact parameters and galaxy redshift errors. We find that the EAGLE reference model is in excellent agreement with the observations. In particular, the simulation reproduces the high metal-line optical depths found at small galactocentric distances, the optical depth enhancements out to impact parameters of 2 proper Mpc, and the prominent redshift-space distortions which we find are due to peculiar velocities rather than redshift errors. The agreement is best for halo masses $\sim10^{12.0}$ M$_\odot$, for which the observed and simulated stellar masses also agree most closely. We examine the median ion mass-weighted radial gas velocities around the galaxies, and find that most of the gas is infalling, with the infall velocity depending on halo rather than stellar mass. From this we conclude that the observed redshift-space distortions are predominantly caused by infall rather than outflows.

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.

We present rest-frame optical spectra of 60 faint ($R_{AB}\sim 27$; $L\sim0.1 L_*$) Ly$\alpha$-selected galaxies (LAEs) at $z\approx2.56$. The average LAE is consistent with the extreme low-metallicity end of the continuum-selected galaxy distribution at $z\approx2-3$. In particular, the LAEs have extremely high [OIII] $\lambda$5008/H$\beta$ ratios (log([OIII]/H$\beta$) $\sim$ 0.8) and low [NII] $\lambda$6585/H$\alpha$ ratios (log([NII]/H$\alpha$) $<-1.15$). Using the [OIII] $\lambda$4364 auroral line, we find that the star-forming regions in faint LAEs are characterized by high electron temperatures ($T_e\approx1.8\times10^4$K), low oxygen abundances (12 + log(O/H) $\approx$ 8.04, $Z_{neb}\approx0.22Z_\odot$), and high excitations with respect to more luminous galaxies. Our faintest LAEs have line ratios consistent with even lower metallicities, including six with 12 + log(O/H) $\approx$ 6.9$-$7.4 ($Z_{neb}\approx0.02-0.05Z_\odot$). We interpret these observations in light of new models of stellar evolution (including binary interactions). We find that strong, hard ionizing continua are required to reproduce our observed line ratios, suggesting that faint galaxies are efficient producers of ionizing photons and important analogs of reionization-era galaxies. Furthermore, we investigate physical trends accompanying Ly$\alpha$ emission across the largest current sample of combined Ly$\alpha$ and rest-optical galaxy spectroscopy, including 60 faint LAEs and 368 more luminous galaxies at similar redshifts. We find that Ly$\alpha$ emission is strongly correlated with nebular excitation and ionization and weakly correlated with dust attenuation, suggesting that metallicity plays a strong role in determining the observed properties of these galaxies by modulating their stellar spectra, nebular excitation, and dust content.

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.

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 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 have acquired Hubble Space Telescope (HST) and Very Large Telescope near-infrared spectra and images of supernova (SN) Refsdal after its discovery as an Einstein cross in Fall 2014. The HST light curve of SN Refsdal matches the distinctive, slowly rising light curves of SN 1987A-like supernovae (SNe), and we find strong evidence for a broad H-alpha P-Cygni profile in the HST grism spectrum at the redshift (z = 1.49) of the spiral host galaxy. SNe IIn, powered by circumstellar interaction, could provide a good match to the light curve of SN Refsdal, but the spectrum of a SN IIn would not show broad and strong H-alpha absorption. From the grism spectrum, we measure an H-alpha expansion velocity consistent with those of SN 1987A-like SNe at a similar phase. The luminosity, evolution, and Gaussian profile of the H-alpha emission of the WFC3 and X-shooter spectra, separated by ~2.5 months in the rest frame, provide additional evidence that supports the SN 1987A-like classification. In comparison with other examples of SN 1987A-like SNe, SN Refsdal has a blue B-V color and a high luminosity for the assumed range of potential magnifications. If SN Refsdal can be modeled as a scaled version of SN 1987A, we estimate it would have an ejecta mass of 20+-5 solar masses. The evolution of the light curve at late times will provide additional evidence about the potential existence of any substantial circumstellar material (CSM). Using MOSFIRE and X-shooter spectra, we estimate a subsolar host-galaxy metallicity (8.3+-0.1 dex and <8.4 dex, respectively) near the explosion site.

We present a spectroscopic survey of 318 faint $(R\sim 27$, $L\sim0.1L_*)$, Ly{\alpha}-emission-selected galaxies (LAEs) at 2.5<z<3. A sample of 32 LAEs with rest-frame optical spectra from Keck/MOSFIRE are used to interpret the LAE spectra in the context of their systemic redshifts. We find that the Ly{\alpha} emission of LAEs is typically less spectrally extended than among samples of more luminous continuum-selected galaxies (LBGs) at similar redshifts. Using the MOSFIRE subsample, we find that the peak of the Ly{\alpha} line is shifted by +200 km/s with respect to systemic across a diverse set of galaxies including both LAEs and LBGs. We also find a small number of objects with significantly blueshifted Ly{\alpha} emission, a potential indicator of accreting gas. The Ly{\alpha}-to-H{\alpha} line ratios suggest that the LAEs have Ly{\alpha} escape fractions $f_{\rm esc,Ly{\alpha}} \approx 30$%, significantly higher than typical LBG samples. Using redshifts calibrated by our MOSFIRE sample, we construct composite LAE spectra, finding the first evidence for metal-enriched outflows in such intrinsically-faint high-redshift galaxies. These outflows have smaller continuum covering fractions $(f_c \approx 0.3)$ and velocities $(v_{\rm ave} \approx 100-200$ km/s, $v_{\rm max} \approx 500$ km/s$)$ than those associated with typical LBGs, suggesting that gas covering fraction is a likely driver of the high Ly{\alpha} and Ly-continuum escape fractions of LAEs. Our results suggest a similar scaling of outflow velocity with star formation rate as is observed at lower redshifts $(v_{\rm outflow} \sim {\rm SFR}^{0.25})$ and indicate that a substantial fraction of gas is ejected with $v > v_{esc}$.

We use quasar absorption lines to study the physical conditions in the circumgalactic medium of redshift $z\approx 2.3$ star-forming galaxies taken from the Keck Baryonic Structure Survey (KBSS). In Turner et al. 2014 we used the pixel optical depth technique to show that absorption by HI and the metal ions OVI, NV, CIV, CIII and SiIV is strongly enhanced within $|\Delta v|\lesssim170$ km/s and projected distances $|d|\lesssim180$ proper kpc from sightlines to the background quasars. Here we demonstrate that the OVI absorption is also strongly enhanced at fixed HI, CIV, and SiIV optical depths, and that this enhancement extends out to $\sim350$ km/s. At fixed HI the increase in the median OVI optical depth near galaxies is 0.3-0.7 dex and is detected at 2--3-$\sigma$ confidence for all seven HI bins that have $\log_{10}\tau_{\rm HI}\ge-1.5$. We use ionization models to show that the observed strength of OVI as a function of HI is consistent with enriched, photoionized gas for pixels with $\tau_{\rm HI}\gtrsim10$. However, for pixels with $\tau_{\rm HI} \lesssim 1$ this would lead to implausibly high metallicities at low densities if the gas were photoionized by the background radiation. This indicates that the galaxies are surrounded by gas that is sufficiently hot to be collisionally ionized ($T > 10^5\,$K) and that a substantial fraction of the hot gas has a metallicity $\gtrsim 10^{-1}$ of solar. Given the high metallicity and large velocity extent (out to $\sim1.5\times v_{\rm circ}$) of this gas, we conclude that we have detected hot, metal enriched outflows arising from star-forming galaxies.

In the redshift interval of $2<z<3$, the physical conditions of the inter-stellar medium (ISM) in star-forming galaxies are likely to be different from those in the local Universe because of lower gaseous metallicities, higher gas fractions, and higher star formation activities. In fact, observations suggest that higher electron densities, higher ionization parameters, and harder UV radiation fields are common. In this paper, based on the spectra of H$\alpha$-selected star-forming galaxies at $z=2.5$ taken with Multi-Object Spectrometer for InfraRed Exploration (MOSFIRE) on Keck-1 telescope, we measure electron densities ($n_e$) using the oxygen line ratio ( [OII]$\lambda\lambda$3726,3729), and investigate the relationships between the electron density of ionized gas and other physical properties. As a result, we find that the specific star formation rate (sSFR) and the surface density of SFR ($\Sigma_\mathrm{SFR}$) are correlated with the electron density at $z=2.5$ for the first time. The $\Sigma_\mathrm{SFR}-n_e$ relation is likely to be linked to the star formation law in HII regions (where star formation activity is regulated by interstellar pressure). Moreover, we discuss the mode of star formation in those galaxies. The correlation between sSFR and $\Sigma_\mathrm{SFR}$ suggests that highly star-forming galaxies (with high sSFR) tend to be characterized by higher surface densities of star formation ($\Sigma_\mathrm{SFR}$) and thus higher $n_e$ values as well.

(Abridged) We study metal absorption around 854 $z\approx$2.4 star-forming galaxies taken from the Keck Baryonic Structure Survey (KBSS). The galaxies examined in this work lie in the fields of 15 hyper-luminous background QSOs, with galaxy impact parameters ranging from 35 proper kpc (pkpc) to 2 proper Mpc (pMpc). Using the pixel optical depth technique, we present the first galaxy-centred 2-D maps of the median absorption by OVI, NV, CIV, CIII, and SIV, as well as updated results for HI. At small galactocentric radii we detect a strong enhancement of the absorption relative to randomly located regions that extend out to at least 180 pkpc in the transverse direction, and $\pm$240 km/s along the line-of-sight (LOS, ~1 pMpc in the case of pure Hubble flow) for all ions except NV. For CIV (and HI) we detect a significant enhancement of the absorption signal out to 2 pMpc in the transverse direction, corresponding to the maximum impact parameter in our sample. After normalising the median absorption profiles to account for variations in line strengths and detection limits, in the transverse direction we find no evidence for a sharp drop-off in metals distinct from that of HI. We argue instead that non-detection of some metal line species in the extended circumgalactic medium is consistent with differences in the detection sensitivity. We also present measurements of covering fractions and equivalent widths as a function of projected galaxy distance. Limiting the sample to the 340 galaxies with redshifts measured from nebular emission lines does not decrease the extent of the enhancement along the LOS compared to that in the transverse direction. This rules out redshift errors as the source of the observed redshift-space anisotropy and thus implies that we have detected the signature of gas peculiar velocities from infall, outflows, or virial motions for HI, OVI, CIV, CIII, and CIV.

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 present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

We present multi-epoch observations from the Submillimeter Array (SMA) for a sample of 171 bright blazars, 43 of which were detected by Fermi during the first three months of observations. We explore the correlation between their gamma-ray properties and submillimeter observations of their parsec-scale jets, with a special emphasis on spectral index in both bands and the variability of the synchrotron component. Subclass is determined using a combination of Fermi designation and the Candidate Gamma-Ray Blazar Survey (CGRaBS), resulting in 35 BL Lac objects and 136 flat-spectrum radio quasars (FSRQs) in our total sample. We calculate submillimeter energy spectral indices using contemporaneous observations in the 1 mm and 850 micron bands during the months August-October 2008. The submillimeter light curves are modeled as first-order continuous autoregressive processes, from which we derive characteristic timescales. Our blazar sample exhibits no differences in submillimeter variability amplitude or characteristic timescale as a function of subclass or luminosity. All of the the light curves are consistent with being produced by a single process that accounts for both low and high states, and there is additional evidence that objects may be transitioning between blazar class during flaring epochs.

We present new photometric and spectroscopic observations for 2M 1533+3759 (= NSVS 07826147). It has an orbital period of 0.16177042 day, significantly longer than the 2.3--3.0 hour periods of the other known eclipsing sdB+dM systems. Spectroscopic analysis of the hot primary yields Teff = 29230 +/- 125 K, log g = 5.58 +/- 0.03 and log N(He)/N(H) = -2.37 +/- 0.05. The sdB velocity amplitude is K1 = 71.1 +/- 1.0 km/s. The only detectable light contribution from the secondary is due to the surprisingly strong reflection effect. Light curve modeling produced several solutions corresponding to different values of the system mass ratio, q(M2/M1), but only one is consistent with a core helium burning star, q=0.301. The orbital inclination is 86.6 degree. The sdB primary mass is M1 = 0.376 +/- 0.055 Msun and its radius is R1 = 0.166 +/- 0.007 Rsun. 2M1533+3759 joins PG0911+456 (and possibly also HS2333+3927) in having an unusually low mass for an sdB star. SdB stars with masses significantly lower than the canonical value of 0.48 Msun, down to as low as 0.30 Msun, were theoretically predicted by Han et al. (2002, 2003), but observational evidence has only recently begun to confirm the existence of such stars. The existence of core helium burning stars with masses lower than 0.40--0.43 Msun implies that at least some sdB progenitors have initial main sequence masses of 1.8--2.0 Msun or more, i.e. they are at least main sequence A stars. The secondary is a main sequence M5 star.

We detected a correlation between optical and giant radio pulse emission from the Crab pulsar. Optical pulses coincident with the giant radio pulses were on average 3% brighter than those coincident with normal radio pulses. Combined with the lack of any other pulse profile changes, this result indicates that both the giant radio pulses and the increased optical emission are linked to an increase in the electron-positron plasma density.