Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$\sigma$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.

Despite significant progress both observationally and theoretically, the origin of high-ionization nebular HeII emission in galaxies dominated by stellar photoionization remains unclear. Accretion-powered radiation from high-mass X-ray binaries (HMXBs) is still one of the leading proposed explanations for the missing $\mathrm{He^+}$-ionizing photons, but this scenario has yet to be conclusively tested. In this paper, we present nebular line predictions from a grid of photoionization models with input SEDs containing the joint contribution of both stellar atmospheres and a multi-color disk model for HMXBs. This grid demonstrates that HMXBs are inefficient producers of the photons necessary to power HeII, and can only boost this line substantially in galaxies with HMXB populations large enough to power X-ray luminosities of $10^{42}$ erg/s per unit star formation rate (SFR). To test this, we assemble a sample of eleven low-redshift star-forming galaxies with high-quality constraints on both X-ray emission from Chandra and HeII emission from deep optical spectra, including new observations with the MMT. These data reveal that the HMXB populations of these nearby systems are insufficient to account for the observed HeII strengths, with typical X-ray luminosities or upper limits thereon of only $10^{40}$-$10^{41}$ erg/s per SFR. This indicates that HMXBs are not the dominant source of $\mathrm{He^+}$ ionization in these metal-poor star-forming galaxies. We suggest that the solution may instead reside in revisions to stellar wind predictions, softer X-ray sources, or very hot products of binary evolution at low metallicity.

Galaxy evolution is thought to be driven in large part by the flow of gas between galaxies and the circumgalactic medium (CGM), a halo of metal-enriched gas extending out to $\gtrsim100$ kpc from each galaxy. Studying the spatial structure of the CGM holds promise for understanding these gas flow mechanisms; however, the common method using background quasar sightlines provides minimal spatial information. Recent works have shown the utility of extended background sources such as giant gravitationally lensed arcs. Using background lensed arcs from the CSWA 38 lens system, we continuously probed, at a resolution element of about 15 kpc$^2$, the spatial and kinematic distribution of MgII absorption in a star-forming galaxy at $z=0.77$ (stellar mass $\approx 10^{9.7}$ M$_\odot$, star formation rate $\approx 10$ M$_\odot$ yr$^{-1}$) at impact parameters $D=5-40$ kpc. Our results present an anisotropic, optically thick medium whose absorption strength decreases with increasing impact parameter, in agreement with the statistics towards quasars and other gravitational arcs. Furthermore, we detect low line-of-sight velocities ($v\approx-25-80$ km s$^{-1}$) and relatively high velocity dispersion ($\sigma\approx50\pm30$ km s$^{-1}$) in the MgII gas. These measures provide evidence of a mainly pressure-supported, metal-enriched diffuse gas recycling through the CGM rather than an active galactic outflow.

Connecting the observed rest-ultraviolet (UV) luminosities of high-$z$ galaxies to their intrinsic luminosities (and thus star formation rates) requires correcting for the presence of dust. We bypass a common dust-correction approach that uses empirical relationships between infrared (IR) emission and UV colours, and instead augment a semi-empirical model for galaxy formation with a simple -- but self-consistent -- dust model and use it to jointly fit high-$z$ rest-UV luminosity functions (LFs) and colour-magnitude relations ($M_{\mathrm{UV}}$-$\beta$). In doing so, we find that UV colours evolve with redshift (at fixed UV magnitude), as suggested by observations, even in cases without underlying evolution in dust production, destruction, absorption, or geometry. The observed evolution in our model arises due to the reduction in the mean stellar age and rise in specific star formation rates with increasing $z$. The UV extinction, $A_{\mathrm{UV}}$, evolves similarly with redshift, though we find a systematically shallower relation between $A_{\mathrm{UV}}$ and $M_{\mathrm{UV}}$ than that predicted by IRX-$\beta$ relationships derived from $z \sim 3$ galaxy samples. Finally, assuming that high $1600 \r{A}$ transmission ($\gtrsim 0.6$) is a reliable LAE indicator, modest scatter in the effective dust surface density of galaxies can explain the evolution both in $M_{\mathrm{UV}}$-$\beta$ and LAE fractions. These predictions are readily testable by deep surveys with the James Webb Space Telescope.

Extensive surveys with the \textit{Hubble Space Telescope} (HST) over the past decade, targeting some of the most massive clusters in the sky, have uncovered dozens of galaxy-cluster strong lenses. The massive cluster strong-lens scale is typically $\theta_{E}\sim10\arcsec$ to $\sim30-35\arcsec$, with only a handful of clusters known with Einstein radii $\theta_{E}\sim40\arcsec$ or above (for $z_{source}=2$, nominally). Here we report another very large cluster lens, RXC J0032.1+1808 ($z=0.3956$), the second richest cluster in the redMapper cluster catalog and the 85th most massive cluster in the Planck Sunyaev-Zel'dovich catalog. With our Light-Traces-Mass and fully parametric (dPIEeNFW) approaches, we construct strong lensing models based on 18 multiple images of 5 background galaxies newly identified in the \textit{Hubble} data mainly from the \textit{Reionization Lensing Cluster Survey} (RELICS), in addition to a known sextuply imaged system in this cluster. Furthermore, we compare these models to Lenstool and GLAFIC models that were produced independently as part of the RELICS program. All models reveal a large effective Einstein radius of $\theta_{E}\simeq40\arcsec$ ($z_{source}=2$), owing to the obvious concentration of substructures near the cluster center. Although RXC J0032.1+1808 has a very large critical area and high lensing strength, only three magnified high-redshift candidates are found within the field targeted by RELICS. Nevertheless, we expect many more high-redshift candidates will be seen in wider and deeper observations with \textit{Hubble} or \emph{JWST}. Finally, the comparison between several algorithms demonstrates that the total error budget is largely dominated by systematic uncertainties.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We quantify the distribution of [OIII]+H$\beta$ line strengths at z$\simeq$7 using a sample of 22 bright (M$_{\mathrm{UV}}$ $\lesssim$ $-$21) galaxies. We select these systems over wide-area fields (2.3 deg$^2$ total) using a new colour-selection which precisely selects galaxies at z$\simeq$6.63$-$6.83, a redshift range where blue Spitzer/IRAC [3.6]$-$[4.5] colours unambiguously indicate strong [OIII]$+$H$\beta$ emission. These 22 galaxies suggest a log-normal [OIII]$+$H$\beta$ EW distribution with median EW = 692$^{+102}_{-103}$ $\mathrm{\mathring{A}}$ and standard deviation = 0.25$^{+0.06}_{-0.05}$ dex. We find no evidence for strong variation in this EW distribution with UV luminosity. The typical [OIII]+H$\beta$ EW at z$\simeq$7 implied by our sample is considerably larger than that in massive star forming galaxies at z$\simeq$2, consistent with a shift toward larger average sSFR (5 Gyr$^{-1}$) and lower metallicities (0.2 Z$_\odot$). We also find evidence for the emergence of a population with yet more extreme nebular emission ([OIII]+H$\beta$ EW$>$1200 $\mathrm{\mathring{A}}$) that is rarely seen at lower redshifts. These objects have extremely large sSFR ($>$50 Gyr$^{-1}$), as would be expected for systems undergoing a burst or upturn in star formation. While this may be a short-lived phase, our results suggest that 20% of the z$\simeq$7 population has such extreme nebular emission, implying that galaxies likely undergo intense star formation episodes regularly at z$>$6. We argue that this population may be among the most effective ionizing agents in the reionization era, both in terms of photon production efficiency and escape fraction. We furthermore suggest that galaxies passing through this large sSFR phase are likely to be very efficient in forming bound star clusters.

We measure the ionizing photon production efficiency ($\xi_{ion}$) of low-mass galaxies ($10^{7.8}$-$10^{9.8}$ $M_{\odot}$) at $1.4<z<2.7$, allowing us to better understand the contribution of dwarf galaxies to the ionizing background and cosmic reionization. We target galaxies that are magnified by the strong lensing galaxy clusters Abell 1689, MACS J0717, and MACS J1149. We utilize Keck/MOSFIRE spectra to measure optical nebular emission line fluxes and HST imaging to measure the rest-UV and rest-optical photometry. We present two methods of stacking. First, we take the average of the log(L$_{H\alpha}$ /L$_{UV}$) of galaxies in our sample to determine the typical log($\xi_{ion}$). Second, we take the logarithm of the total L$_{H\alpha}$ over the total L$_{UV}$. We prefer the latter as it provides the total ionizing UV luminosity density of galaxies when multiplied by the non-ionizing UV luminosity density from the UV luminosity function. log($\xi_{ion}$) calculated from the second method is $\sim$ 0.2 dex higher than the first method. We do not find any strong dependence between log($\xi_{ion}$) and stellar mass, M$_{UV}$ or UV spectral slope ($\beta$). We report a value of log($\xi_{ion}$) $\sim25.47\pm 0.09$ for our UV-complete sample ($-22<M_{UV}<-17.3$) and $\sim25.37\pm0.11$ for our mass-complete sample ($7.8<\log(M_*)<9.8)$. These values are consistent with measurements of more massive, more luminous galaxies in other high-redshift studies that use the same stacking technique. Our log($\xi_{ion}$) is $0.2-0.3$ dex higher than low-redshift galaxies of similar mass, indicating an evolution in the stellar properties, possibly due to metallicity, age, or the prevalence of binary stars. We also find a correlation between log($\xi_{ion}$) and the equivalent widths of H$\alpha$ and [OIII]$\lambda$5007 fluxes, confirming that these equivalent widths can be used to estimate $\xi_{ion}$.

We identify and characterize compact dwarf starburst (CDS) galaxies in the RESOLVE survey, a volume-limited census of galaxies in the local universe, to probe whether this population contains any residual ``blue nuggets,'' a class of intensely star-forming compact galaxies first identified at high redshift $z$. Our 50 low-$z$ CDS galaxies are defined by dwarf masses (stellar mass $M_* < 10^{9.5}$ M$_{\odot}$), compact bulged-disk or spheroid-dominated morphologies (using a quantitative criterion, $\mu_\Delta > 8.6$), and specific star formation rates above the defining threshold for high-$z$ blue nuggets ($\log$ SSFR [Gyr$^{-1}] > -0.5$). Across redshifts, blue nuggets exhibit three defining properties: compactness relative to contemporaneous galaxies, abundant cold gas, and formation via compaction in mergers or colliding streams. Those with halo mass below $M_{\rm halo} \sim 10^{11.5}$ M$_{\odot}$ may in theory evade permanent quenching and cyclically refuel until the present day. Selected only for compactness and starburst activity, our CDS galaxies generally have $M_{\rm halo} \lesssim 10^{11.5}$ M$_{\odot}$ and gas-to-stellar mass ratio $\gtrsim$1. Moreover, analysis of archival DECaLS photometry and new 3D spectroscopic observations for CDS galaxies reveals a high rate of photometric and kinematic disturbances suggestive of dwarf mergers. The SSFRs, surface mass densities, and number counts of CDS galaxies are compatible with theoretical and observational expectations for redshift evolution in blue nuggets. We argue that CDS galaxies represent a maximally-starbursting subset of traditional compact dwarf classes such as blue compact dwarfs and blue E/S0s. We conclude that CDS galaxies represent a low-$z$ tail of the blue nugget phenomenon formed via a moderated compaction channel that leaves open the possibility of disk regrowth and evolution into normal disk galaxies.

Recent observations have revealed the presence of strong CIII] emission (EW$_{\rm{CIII]}}>20$ \r{A}) in $z>6$ galaxies, the origin of which remains unclear. In an effort to understand the nature of these line emitters, we have initiated a survey targeting CIII] emission in gravitationally-lensed reionization era analogs identified in HST imaging of clusters from the RELICS survey. Here we report initial results on four galaxies selected to have low stellar masses (2-8$\times$10$^7$ M$_\odot$) and J$_{125}$-band flux excesses indicative of intense [OIII]+H$\beta$ emission (EW$_{\rm{[OIII]+H\beta}}$=500-2000 \r{A}), similar to what has been observed at $z>6$. We detect CIII] emission in three of the four sources, with the CIII] EW reaching values seen in the reionization era (EW$_{\rm{CIII]}}\simeq 17-22$ \r{A}) in the two sources with the strongest optical line emission (EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}). We have obtained a Magellan/FIRE near-infrared spectrum of the strongest CIII] emitter in our sample, revealing gas that is both metal poor and highly ionized. Using photoionization models, we are able to simultaneously reproduce the intense CIII] and optical line emission for extremely young (2-3 Myr) and metal poor (0.06-0.08 Z$_\odot$) stellar populations, as would be expected after a substantial upturn in the SFR of a low mass galaxy. The sources in this survey are among the first for which CIII] has been used as the primary means of redshift confirmation. We suggest that it should be possible to extend this approach to $z>6$ with current facilities, using CIII] to measure redshifts of objects with IRAC excesses indicating EW$_{\rm{[OIII]+H\beta}}\simeq 2000$ \r{A}, providing a method of spectroscopic confirmation independent of Ly$\alpha$.

Outflows form an integral component in regulating the gas cycling in and out of galaxies, although their impact on the galaxy hosts is still poorly understood. Here we present an analysis of 405 high mass (log M$_{*}$/M$_{\odot}\geqslant10$), star-forming galaxies (excluding AGN) with low inclinations at $z\sim$0, using stacking techniques of the NaD $\lambda\lambda$5889,5895 A neutral gas tracer in IFU observations from the MaNGA DR15 survey. We detect outflows in the central regions of 78/405 galaxies and determine their extent and power through the construction of stacked annuli. We find outflows are most powerful in central regions and extend out to $\sim$1R$_{e}$, with declining mass outflow rates and loading factors as a function of radius. The stacking of spaxels over key galaxy quantities reveals outflow detections in regions of high $\Sigma_{\text{SFR}}$ ($\gtrsim$0.01 M$_{\odot}$yr$^{-1}$kpc$^{-2}$) and $\Sigma_{M_{*}}$ ($\gtrsim$10$^{7}$ M$_{\odot}$kpc$^{-2}$) along the resolved main sequence. Clear correlations with $\Sigma_{\text{SFR}}$ suggest it is the main regulator of outflows, with a critical threshold of $\sim$0.01 M$_{\odot}$yr$^{-1}$kpc$^{-2}$ needed to escape the weight of the disk and launch them. Furthermore, measurements of the H$\delta$ and D$_{n}$4000 indices reveal virtually identical star formation histories between galaxies with outflows and those without. Finally, through stacking of HI 21 cm observations for a subset of our sample, we find outflow galaxies show reduced HI gas fractions at central velocities compared to their non-detection control counterparts, suggestive of some removal of HI gas, likely in the central regions of the galaxies, but not enough to completely quench the host.

Since 2016, the HI-MaNGA survey has been obtaining neutral hydrogen (which emits a radio 21cm line) measurements for nearby galaxies which are part of a survey known as "MaNGA". This has been done using observations at the Green Bank Radio Telescope (GBT). We use calibration data obtained routinely at the start of observing runs during this entire period, as well as an extra set of calibration data obtained in 2019 to check the standard GBT data reduction calibrations at L-band (1.4 GHz, or 20cm). We find an constant offset of 20% in the calibration relative to the standard values given in the GBT data reduction pipeline for data of this type. This correction will be used to improve measurements for HI-MaNGA.

Strong [OIII]$\lambda\lambda$4959,5007+H$\beta$ emission appears to be typical in star-forming galaxies at z > 6.5. As likely contributors to cosmic reionization, these galaxies and the physical conditions within them are of great interest. At z > 6.5, where Ly$\alpha$ is greatly attenuated by the intergalactic medium, rest-UV metal emission lines provide an alternative measure of redshift and also constraints on the physical properties of star-forming regions and massive stars. We present the first statistical sample of rest-UV line measurements in z $\sim$ 2 galaxies selected as analogs of those in the reionization era based on [OIII]$\lambda\lambda$4959,5007 EW or rest-frame U-B color. Our sample is drawn from the 3D-HST Survey and spans the redshift range 1.36 $\leqslant$ z $\leqslant$ 2.49. We find that the median Ly$\alpha$ and CIII]$\lambda\lambda$1907,1909 EWs of our sample are significantly greater than those of z $\sim$ 2 UV-continuum-selected star-forming galaxies. Measurements from both individual and composite spectra indicate a monotonic, positive correlation between CIII] and [OIII], while a lack of trend is observed between Ly$\alpha$ and [OIII] at [OIII] EW < 1000$\unicode{x212B}$. At higher [OIII] EW, extreme Ly$\alpha$ emission starts to emerge. Using stacked spectra, we find that Ly$\alpha$ and CIII] are significantly enhanced in galaxies with lower metallicity. Two objects in our sample appear comparable to z > 6.5 galaxies with exceptionally strong rest-UV metal line emission. These objects have significant CIV$\lambda\lambda$1548,1550, HeII$\lambda$1640, and OIII]$\lambda\lambda$1661,1665 emission in addition to intense Ly$\alpha$ or CIII]. Detailed characterization of these lower-redshift analogs provides unique insights into the physical conditions in z > 6.5 star-forming regions, motivating future observations of reionization-era analogs at lower redshifts.

We discuss the rest-frame optical emission line spectra of a large (~50) sample of z=3.1 Lyman alpha emitters (LAEs) whose physical properties suggest such sources are promising analogs of galaxies in the reionization era. Reliable Lyman continuum escape fractions have now been determined for a large sample of such LAEs from the Lyman Continuum Escape Survey (LACES) undertaken via deep HST imaging in the SSA22 survey area reported in Fletcher et al. (2019). Using new measures of [OII] emission secured from Keck MOSFIRE spectra we re-examine, for a larger sample, earlier claims that Lyman continuum leakages may correlate with the nebular emission line ratio [OIII]/[OII] as expected for density-bound HII regions. We find that a large [OIII]/[OII] line ratio is indeed a necessary condition for Lyman continuum leakage, strengthening earlier claims made using smaller samples at various redshifts. However, not all LAEs with large [OIII]/[OII] line ratios are leakers and leaking radiation appears not to be associated with differences in other spectral diagnostics. This suggests the detection of leaking radiation is modulated by an additional property, possibly the viewing angle for porous HII regions. We discuss our new results in the context of the striking bimodality of LAE leakers and non-leakers found in the LACES program and the implications for the sources of cosmic reionization.

We present Keck/MOSFIRE and Keck/LRIS spectroscopy of A1689-217, a lensed (magnification $\sim7.9$), star-forming (SFR $\sim$ 16 $\rm M_\odot\ yr^{-1}$), dwarf (log($M_\ast/M_\odot$) = $8.07-8.59$), Ly$\alpha$-emitter ($\rm EW_0\sim138$ A) at $z$ = 2.5918. Dwarf galaxies similar to A1689-217 are common at high redshift and likely responsible for reionization, yet few have been studied with detailed spectroscopy. We report a 4.2$\sigma$ detection of the electron-temperature-sensitive [O III]$\lambda$4363 emission line and use this line to directly measure an oxygen abundance of 12+$\log$(O/H) = 8.06 $\pm$ 0.12 ($\sim1/4\ \rm Z_\odot$). A1689-217 is the lowest mass galaxy at $z>$ 2 with an [O III]$\lambda$4363 detection. Using the rest-optical emission lines, we measure A1689-217's other nebular conditions including electron temperature ($T_e$([O III]) $\sim$ 14,000 K), electron density ($n_e\sim$ 220 $\rm cm^{-3}$) and reddening ($E(B-V)\sim$ 0.39). We study relations between strong-line ratios and direct metallicities with A1689-217 and other galaxies with [O III]$\lambda$4363 detections at $z\sim0-3.1$, showing that the locally-calibrated, oxygen-based, strong-line relations are consistent from $z\sim0-3.1$. We also show additional evidence that the $\rm O_{32}$ vs. $\rm R_{23}$ excitation diagram can be utilized as a redshift-invariant, direct-metallicity-based, oxygen abundance diagnostic out to $z\sim$ 3.1. From this excitation diagram and the strong-line ratio $-$ metallicity plots, we observe that the ionization parameter at fixed O/H is consistent with no redshift evolution. Although A1689-217 is metal-rich for its $M_\ast$ and SFR, we find it to be consistent within the large scatter of the low-mass end of the Fundamental Metallicity Relation.

We investigate the environmental dependence of the local gas-phase metallicity in a sample of star-forming galaxies from the MaNGA survey. Satellite galaxies with stellar masses in the range $9<\log(M_{*}/M_{\odot})<10$ are found to be $\sim 0.05 \, \mathrm{dex}$ higher in metallicity than centrals of similar stellar mass. Within the low-mass satellite population, we find that the interstellar medium (ISM) metallicity depends most strongly on the stellar mass of the galaxy that is central to the halo, though there is no obvious difference in the metallicity gradients. At fixed total stellar mass, the satellites of high mass ($M_{*}>10^{10.5} \, \mathrm{M_{\odot}}$) centrals are $\sim 0.1 \, \mathrm{dex}$ more metal rich than satellites of low-mass ($M_{*} < 10^{10} \, \mathrm{M_{\odot}}$) centrals, controlling for local stellar mass surface density and gas fraction. Fitting a gas-regulator model to the spaxel data, we are able to account for variations in the local gas fraction, stellar mass surface density and local escape velocity-dependent outflows. We find that the best explanation for the metallicity differences is the variation in the average metallicity of accreted gas between different environments that depends on the stellar mass of the dominant galaxies in each halo. This is interpreted as evidence for the exchange of enriched gas between galaxies in dense environments that is predicted by recent simulations.

Ultraviolet (UV) observations of local star-forming galaxies have begun to establish an empirical baseline for interpreting the rest-UV spectra of reionization-era galaxies. However, existing high-ionization emission line measurements at $z>6$ ($\mathrm{W_{C IV, 0}} \gtrsim 20$ {\AA}) are uniformly stronger than observed locally ($\mathrm{W_{C IV, 0}} \lesssim 2$ {\AA}), likely due to the relatively high metallicities ($Z/Z_\odot > 0.1$) typically probed by UV surveys of nearby galaxies. We present new HST/COS spectra of six nearby ($z<0.01$) extremely metal-poor galaxies (XMPs, $Z/Z_\odot \lesssim 0.1$) targeted to address this limitation and provide constraints on the highly-uncertain ionizing spectra powered by low-metallicity massive stars. Our data reveal a range of spectral features, including one of the most prominent nebular C IV doublets yet observed in local star-forming systems and strong He II emission. Using all published UV observations of local XMPs to-date, we find that nebular C IV emission is ubiquitous in very high specific star formation rate systems at low metallicity, but still find equivalent widths smaller than those measured in individual lensed systems at $z>6$. Our moderate-resolution HST/COS data allow us to conduct an analysis of the stellar winds in a local nebular C IV emitter, which suggests that some of the tension with $z>6$ data may be due to existing local samples not yet probing sufficiently high $\mathrm{\alpha/Fe}$ abundance ratios. Our results indicate that C IV emission can play a crucial role in the JWST and ELT era by acting as an accessible signpost of very low metallicity ($Z/Z_\odot < 0.1$) massive stars in assembling reionization-era systems.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

We present the HI-MaNGA programme of HI follow-up for the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. MaNGA, which is part of the Fourth phase of the Sloan Digital Sky Surveys (SDSS-IV), is in the process of obtaining integral field unit (IFU) spectroscopy for a sample of ~10,000 nearby galaxies. We give an overview of the HI 21cm radio follow-up observing plans and progress and present data for the first 331 galaxies observed in the 2016 observing season at the Robert C. Bryd Green Bank Telescope (GBT). We also provide a cross match of the current MaNGA (DR15) sample with publicly available HI data from the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey. The addition of HI data to the MaNGA data set will strengthen the survey's ability to address several of its key science goals that relate to the gas content of galaxies, while also increasing the legacy of this survey for all extragalactic science.

Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z$\sim$4$-$10 by utilizing mock galaxy samples produced by an empirical model, the UniverseMachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with $\gtrsim$5$\sigma$ significance at z$\sim$4$-$10. Halo mass precisions resulting from Cycle 1 angular clustering measurements will be $\sim$0.2 dex for faint (-18 $\gtrsim$ $\mathit{M}_{\mathrm{UV}}^{ }$ $\gtrsim$ -19) galaxies at z$\sim$4$-$10 as well as $\sim$0.3 dex for bright ($\mathit{M}_{\mathrm{UV}}^{ }$ $\sim$ -20) galaxies at z$\sim$4$-$7. Dedicated spectroscopic follow-up over $\sim$150 arcmin$^2$ would improve these precisions by $\sim$0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar-halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z$>$4. We also find that $\sim$1000 individual satellites will be identifiable at z$\sim$4$-$8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z$\lesssim$2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z$\sim$4$-$8 with $\sim$0.1$-$0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization.

Galaxies in the reionization-era have been shown to have prominent [OIII]+H$\beta$ emission. Little is known about the gas conditions and radiation field of this population, making it challenging to interpret the spectra emerging at $z\gtrsim6$. Motivated by this shortcoming, we have initiated a large MMT spectroscopic survey identifying rest-frame optical emission lines in 227 intense [OIII] emitting galaxies at $1.3<z<2.4$. This sample complements the MOSDEF and KBSS surveys, extending to much lower stellar masses ($10^7-10^8 M_\odot$) and larger specific star formation rates ($5-300$ Gyr$^{-1}$), providing a window on galaxies directly following a burst or recent upturn in star formation. The hydrogen ionizing production efficiency ($\xi_{\rm{ion}}$) is found to increase with the [OIII] EW, in a manner similar to that found in local galaxies by Chevallard et al. (2018). We describe how this relationship helps explain the anomalous success rate in identifying Ly$\alpha$ emission in $z\gtrsim7$ galaxies with strong [OIII]+H$\beta$ emission. We probe the impact of the intense radiation field on the ISM using O32 and Ne3O2, two ionization-sensitive indices. Both are found to scale with the [OIII] EW, revealing extreme ionization conditions not commonly seen in older and more massive galaxies. In the most intense line emitters, the indices have very large average values (O32 $=9.1$, Ne3O2 $=0.5$) that have been shown to be linked to ionizing photon escape. We discuss implications for the nature of galaxies most likely to have O32 values associated with significant LyC escape. Finally we consider the optimal strategy for JWST spectroscopic investigations of galaxies at $z\gtrsim10$ where the strongest rest-frame optical lines are no longer visible with NIRSpec.

We present results from the LymAn Continuum Escape Survey (LACES), a Hubble Space Telescope (HST) program designed to characterize the ionizing radiation emerging from a sample of Lyman alpha emitting galaxies at redshift $z\simeq 3.1$. As many show intense [O III] emission characteristic of $z>6.5$ star-forming galaxies, they may represent valuable low redshift analogs of galaxies in the reionization era. Using HST Wide Field Camera 3 / UVIS $F336W$ to image Lyman continuum emission, we investigate the escape fraction of ionizing photons in this sample. For 61 sources, of which 77% are spectroscopically confirmed and 53 have measures of [O III] emission, we detect Lyman continuum leakage in 20%, a rate significantly higher than is seen in individual continuum-selected Lyman break galaxies. We estimate there is a 98% probability that $\leq 2$ of our detections could be affected by foreground contamination. Fitting multi-band spectral energy distributions (SEDs) to take account of the varying stellar populations, dust extinctions and metallicities, we derive individual Lyman continuum escape fractions corrected for foreground intergalactic absorption. We find escape fractions of 15 to 60% for individual objects, and infer an average 20% escape fraction by fitting composite SEDs for our detected samples. Surprisingly however, even a deep stack of those sources with no individual $F336W$ detections provides a stringent upper limit on the average escape fraction of less than 0.5%. We examine various correlations with source properties and discuss the implications in the context of the popular picture that cosmic reionization is driven by such compact, low metallicity star-forming galaxies.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

Diagnosing the physical and chemical conditions within star-forming galaxies (SFGs) is of paramount importance to understanding key components of galaxy formation and evolution: star-formation, gas enrichment, outflows, and accretion. Well established optical emission-line diagnostics used to discern such properties (i.e., metal content, density, strength/shape of ionizing radiation) will be observationally inaccessible for the earliest galaxies, emphasizing the need for robust, reliable interstellar medium (ISM) diagnostics at ultraviolet (UV) wavelengths. Calibrating these UV diagnostics requires a comprehensive comparison of the UV and optical emission lines in nearby SFGs. Optical integral field unit (IFU) surveys have revealed the inhomogeneous nature of the ISM in SFGs, which leads to non-systematic biases in the interpretation of unresolved sources. Spatial variations are especially important to consider at UV wavelengths, where the strongest emission features originate from only the highest excitation regions of the nebula and are challenging to distinguish from competing high-ionization sources (e.g., shocks, AGN, etc.). Since surveys collecting large-scale optical integral field unit (IFU) spectroscopy are already underway, this white paper calls for an IFU or multi-object far-UV (FUV) spectroscopic instrument with high sensitivity, high spatial resolution, and large field of view (FoV). Given the impact of large-scale optical IFU surveys over the past decade, this white paper emphasizes the scientific need for a comparable foundation of spatially-resolved far-UV spectroscopy survey of nearby galaxies that will lay the foundation of diagnostics critical to the interpretation of the distant universe.

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

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

Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Strong gravitational lensing (SL) is a powerful means to map the distribution of dark matter. In this work, we perform a SL analysis of the prominent X-ray cluster RXJ0152.7-1357 (z=0.83, also known as CL 0152.7-1357) in \textit{Hubble Space Telescope} images, taken in the framework of the Reionization Lensing Cluster Survey (RELICS). On top of a previously known $z=3.93$ galaxy multiply imaged by RXJ0152.7-1357, for which we identify an additional multiple image, guided by a light-traces-mass approach we identify seven new sets of multiply imaged background sources lensed by this cluster, spanning the redshift range [1.79-3.93]. A total of 25 multiple images are seen over a small area of ~0.4 $arcmin^2$, allowing us to put relatively high-resolution constraints on the inner matter distribution. Although modestly massive, the high degree of substructure together with its very elongated shape make RXJ0152.7-1357 a very efficient lens for its size. This cluster also comprises the third-largest sample of z~6-7 candidates in the RELICS survey. Finally, we present a comparison of our resulting mass distribution and magnification estimates with those from a Lenstool model. These models are made publicly available through the MAST archive.

We implement a steady, one-dimensional flow model for the X-ray jet of Centaurus A in which entrainment of stellar mass loss is the primary cause of dissipation. Using over 260 ks of new and archival Chandra/ACIS data, we have constrained the temperature, density and pressure distributions of gas in the central regions of the host galaxy of Centaurus A, and so the pressure throughout the length of its jet. The model is constrained by the observed profiles of pressure and jet width, and conserves matter and energy, enabling us to estimate jet velocities, and hence all the other flow properties. Invoking realistic stellar populations within the jet, we find that the increase in its momentum flux exceeds the net pressure force on the jet unless only about one half of the total stellar mass loss is entrained. For self-consistent models, the bulk speed only falls modestly, from ~0.67c to ~0.52c over the range of 0.25-5.94 kpc from the nucleus. The sonic Mach number varies between ~5.3 and 3.6 over this range.

Extremely metal-poor star-forming galaxies (XMPs) represent one of our only laboratories for study of the low-metallicity stars we expect to encounter at early epochs. But as our understanding of the $z>6$ universe has improved, it has become clear that the majority of known XMPs within 100 Mpc host significantly less prominent massive star populations than their reionization-era counterparts, severely limiting their utility as testbeds for interpreting spectral features found at the highest redshifts. Here we present a new photometric selection technique designed to identify nearby XMPs dominated by young stellar populations comparable to those expected in the reionization era. We apply our technique to uncover candidate XMPs in SDSS imaging at magnitudes $16<i'<23$, extending significantly below the completeness limits of the SDSS spectroscopic survey. Spectroscopic observations with the MMT confirm that 32 of the 53 uniformly metal-poor and high specific star formation rate targets we observed have gas-phase oxygen abundances $12+\log\mathrm{O/H}<7.7$ ($Z/Z_\odot<0.1$), including two in the range of the lowest-metallicity galaxies known, $Z/Z_\odot<0.05$. Our observations shed new light onto the longstanding mystery of He II emission in star-forming galaxies: we find that the equivalent width of the He II $\lambda 4686$ high-ionization emission line does not scale with that of H$\beta$ in our sample, suggesting that binary evolution or other processes on $>10$ Myr timescales contribute substantially to the $\mathrm{He^+}$-ionizing photon budget in this metallicity regime. Applying such selection techniques coupled with deep spectroscopy to next-generation photometric surveys like LSST may eventually provide a basis for an empirical understanding of metal-poor massive stars.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

Galactic-scale outflows regulate the stellar mass growth and chemical enrichment of galaxies, yet key outflow properties such as the chemical composition and mass loss rate remain largely unknown. We address these properties with Keck/ESI echellete spectra of nine gravitationally lensed z=2-3 star forming galaxies, probing a range of absorption transitions. Interstellar absorption in our sample is dominated by outflowing material, with typical velocities -150 km/s. Approximately 80% of the total column density is associated with a net outflow. Mass loss rates in the low ionization phase are comparable to or in excess of the star formation rate, with total outflow rates likely higher when accounting for ionized gas. Of order half of the heavy element yield from star formation is ejected in the low ionization phase, confirming that outflows play a critical role in regulating galaxy chemical evolution. Covering fractions vary and are in general non-uniform, with most galaxies having incomplete covering by the low ions across all velocities. Low ion abundance patterns show remarkably little scatter, revealing a distinct "chemical fingerprint" of outflows. Gas phase Si/Fe abundances are significantly super-solar ([Si/Fe]$\gtrsim$0.4) indicating a combination of $\alpha$-enhancement and dust depletion. Derived properties are comparable to the most kinematically broad, metal-rich, and depleted intergalactic absorption systems at similar redshifts, suggesting that these extreme systems are associated with galactic outflows at impact parameters conservatively within a few tens of kpc. We discuss implications of the abundance patterns in z=2-3 galaxies and the role of outflows at this epoch.

We show how the Radon transform (defined as a series of line integrals through an image at different orientations and offsets from the origin) can be used as a simple, non-parametric tool to characterize galaxy velocity fields, specifically their global kinematic position angles (PA_k) and any radial variation or asymmetry in PA_k. This method is fast and easily automated, making it particularly beneficial in an era where IFU and interferometric surveys are yielding samples of thousands of galaxies. We demonstrate the Radon transform by applying it to gas and stellar velocity fields from the first ~2800 galaxies of the SDSS-IV MaNGA IFU survey. We separately classify gas and stellar velocity fields into five categories based on the shape of their radial PA_k profiles. At least half of stellar velocity fields and two-thirds of gas velocity fields are found to show detectable deviations from uniform coplanar circular motion, although most of these variations are symmetric about the center of the galaxy. The behavior of gas and stellar velocity fields is largely independent, even when PA_k profiles for both components are measured over the same radii. We present evidence that one class of symmetric PA_k variations is likely associated with bars and/or oval distortions, while another class is more consistent with warped disks. This analysis sets the stage for more in-depth future studies which explore the origin of diverse kinematic behavior in the galaxy population.

We present a strong-lensing analysis of four massive galaxy clusters imaged with the Hubble Space Telescope in the Reionization Lensing Cluster Survey. We use a Light-Traces-Mass technique to uncover sets of multiply images and constrain the mass distribution of the clusters. These mass models are the first published for Abell S295 and MACS J0159.8-0849, and are improvements over previous models for Abell 697 and MACS J0025.4-1222. Our analysis for MACS J0025.4-1222 and Abell S295 shows a bimodal mass distribution supporting the merger scenarios proposed for these clusters. The updated model for MACS J0025.4-1222 suggests a substantially smaller critical area than previously estimated. For MACS J0159.8-0849 and Abell 697 we find a single peak and relatively regular morphology, revealing fairly relaxed clusters. Despite being less prominent lenses, three of these clusters seem to have lensing strengths, i.e. cumulative area above certain magnification, similar to the Hubble Frontier Fields clusters (e.g., A($\mu>5$) $\sim 1-3$ arcmin$^2$, A($\mu>10$) $\sim 0.5-1.5$ arcmin$^2$), which in part can be attributed to their merging configurations. We make our lens models publicly available through the Mikulski Archive for Space Telescopes. Finally, using Gemini-N/GMOS spectroscopic observations we detect a single emission line from a high-redshift $J_{125}\simeq25.7$ galaxy candidate lensed by Abell 697. While we cannot rule out a lower-redshift solution, we interpret the line as Ly$\alpha$ at $z=5.800\pm 0.001$, in agreement with its photometric redshift and dropout nature. Within this scenario we measure a Ly$\alpha$ rest-frame equivalent width of $52\pm22$ \AA, and an observed Gaussian width of $117\pm 15$ km/s.

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known.

Recent studies have revealed intense UV metal emission lines in a modest sample of z>7 Lyman-alpha emitters, indicating a hard ionizing spectrum is present. If such high ionization features are shown to be common, it may indicate that extreme radiation fields play a role in regulating the visibility of Lyman-alpha in the reionization era. Here we present deep near-infrared spectra of seven galaxies with Lyman-alpha emission at 5.4<z<8.7 (including a newly-confirmed lensed galaxy at z=6.031) and three bright z>7 photometric targets. In nine sources we do not detect UV metal lines. However in the z=8.683 galaxy EGSY8p7, we detect a 4.6 sigma emission line in the narrow spectral window expected for NV 1243. The feature is unresolved (FWHM<90 km/s) and is likely nebular in origin. A deep H-band spectrum of EGSY8p7 reveals non-detections of CIV, He II, and OIII]. The presence of NV requires a substantial flux of photons above 77 eV, pointing to a hard ionizing spectrum powered by an AGN or fast radiative shocks. Regardless of its origin, the intense radiation field of EGSY8p7 may aid the transmission of Lyman-alpha through what is likely a partially neutral IGM. With this new detection, five of thirteen known Lyman-alpha emitters at z>7 have now been shown to have intense UV line emission, suggesting that extreme radiation fields are commonplace among the Lyman-alpha population. Future observations with JWST will eventually clarify the origin of these features and explain their role in the visibility of Lyman-alpha in the reionization era.

Measurements of the galaxy UV luminosity function at z>6 suggest that young stars hosted in low-mass star-forming galaxies produced the bulk of hydrogen-ionizing photons necessary to reionize the intergalactic medium (IGM) by redshift z~6. Whether star-forming galaxies dominated cosmic reionization, however, also depends on their stellar populations and interstellar medium properties, which set, among other things, the production rate of H-ionizing photons, $\xi_\text{ion}^\star$, and the fraction of these escaping into the IGM. Given the difficulty of constraining with existing observatories the physical properties of z>6 galaxies, in this work we focus on a sample of ten nearby objects showing UV spectral features comparable to those observed at z>6. We use the new-generation Beagle tool to model the UV-to-optical photometry and UV/optical emission lines of these Local 'analogues' of high-redshift galaxies, finding that our relatively simple, yet fully self-consistent, physical model can successfully reproduce the different observables considered. Our galaxies span a broad range of metallicities and are characterised by high ionization parameters, low dust attenuation, and very young stellar populations. Through our analysis, we derive a novel diagnostic of the production rate of H-ionizing photons per unit UV luminosity, $\xi_\text{ion}^\star$, based on the equivalent width of the bright $[\text{OIII}] \lambda 4959,5007$ doublet, which does not require measurements of H-recombination lines. This new diagnostic can be used to estimate $\xi_\text{ion}^\star$ from future direct measurements of the $[\text{OIII}] \lambda 4959,5007$ line using JWST/NIRSpec (out to z~9.5), and by exploiting the contamination by $\text{H}\beta + [\text{OIII}] \lambda 4959,5007$ of photometric observations of distant galaxies, for instance from existing Spitzer/IRAC data and from future ones with JWST/NIRCam.

We present a detailed study of a unusually bright, lensed galaxy at $z=5.424$ discovered within the CFHTLS imaging survey. With an observed flux of $i_{\rm AB}=23.0$, J141446.82+544631.9 is one of the brightest galaxies known at $z>5$. It is characterized by strong Ly$\alpha$ emission, reaching a peak in (observed) flux density of $>10^{-16}~{\rm erg}~{\rm s}^{-1}~{\rm cm}^{-2}~{\rm \AA}^{-1}$. A deep optical spectrum from the LBT places strong constraints on NV and CIV emission, disfavouring an AGN source for the emission. However, a detection of the NIV]~$\lambda$1486 emission line indicates a hard ionizing continuum, possibly from hot, massive stars. Resolved imaging from HST deblends the galaxy from a foreground interloper; these observations include narrowband imaging of the Ly$\alpha$ emission, which is marginally resolved on $\sim$few~kpc scales and has EW$_0~\sim$ 260\AA. The Ly$\alpha$ emission extends over ~2000~${\rm km}~{\rm s}^{-1}$ and is broadly consistent with expanding shell models. SED fitting that includes Spitzer/IRAC photometry suggests a complex star formation history that include both a recent burst and an evolved population. J1414+5446 lies 30" from the centre of a known lensing cluster in the CFHTLS; combined with the foreground contribution this leads to a highly uncertain estimate for the lensing magnification in the range $5 \leq \mu \leq 25$. Because of its unusual brightness J1414+5446 affords unique opportunities for detailed study of an individual galaxy near the epoch of reionization and a preview of what can be expected from upcoming wide-area surveys that will yield hundreds of similar objects.

We present the catalog of ~31500 extragalactic HI line sources detected by the completed ALFALFA survey out to z < 0.06 including both high signal-to-noise ratio (> 6.5) detections and ones of lower quality which coincide in both position and recessional velocity with galaxies of known redshift. We review the observing technique, data reduction pipeline, and catalog construction process, focusing on details of particular relevance to understanding the catalog's compiled parameters. We further describe and make available the digital HI line spectra associated with the catalogued sources. In addition to the extragalactic HI line detections, we report nine confirmed OH megamasers and ten OH megamaser candidates at 0.16 < z < 0.22 whose OH line signals are redshifted into the ALFALFA frequency band. Because of complexities in data collection and processing associated with the use of a feed-horn array on a complex single-dish antenna in the terrestrial radio frequency interference environment, we also present a list of suggestions and caveats for consideration by users of the ALFALFA extragalactic catalog for future scientific investigations.

We introduce a new method for determining the influence of galaxies and active galactic nuclei (AGN) on the physical state of the intergalactic medium (IGM) at high redshift and illustrate its potential via a first application to the field of the $z=6.42$ QSO J1148+5251. By correlating the spatial positions of spectroscopically-confirmed Lyman break galaxies (LBGs) with fluctuations in the Lyman alpha forest seen in the high signal-to-noise spectrum of a background QSO, we provide a statistical measure of the typical escape fraction of Lyman continuum photons close to the end of cosmic reionisation. Here we use Keck DEIMOS spectroscopy to locate 7 colour-selected LBGs in the redshift range $5.3\lesssim z\lesssim 6.4$ and confirm a faint $z=5.701$ AGN. We then examine the spatial correlation between this sample and Ly$\alpha$/Ly$\beta$ transmission fluctuations in a Keck ESI spectrum of the QSO. Interpreting the statistical HI proximity effect as arising from faint galaxies clustered around the detected LBGs, we translate the observed mean Ly$\alpha$ transmitted flux around an average detected LBG into a constraint on the mean escape fraction $\langle f_{\rm esc}\rangle\geq0.08$ at $z\simeq6$. We also report evidence of the individual transverse HI proximity effect of a $z=6.177$ luminous LBG via a Ly$\beta$ transmission spike and two broad Ly$\alpha$ transmission spikes around the $z=5.701$ AGN. We discuss the possible origin of such associations which suggest that while faint galaxies are primarily driving reionisation, luminous galaxies and AGN may provide important contributions to the UV background or thermal fluctuations of the IGM at $z\simeq6$. Although a limited sample, our results demonstrate the potential of making progress using this method in resolving one of the most challenging aspects of the contribution of galaxies and AGN to cosmic reionisation.

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

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

Strong gravitational lensing by galaxy clusters has become a powerful tool for probing the high-redshift Universe, magnifying distant and faint background galaxies. Reliable strong lensing (SL) models are crucial for determining the intrinsic properties of distant, magnified sources and for constructing their luminosity function. We present here the first SL analysis of MACS J0308.9+2645 and PLCK G171.9-40.7, two massive galaxy clusters imaged with the Hubble Space Telescope in the framework of the Reionization Lensing Cluster Survey (RELICS). We use the Light-Traces-Mass modeling technique to uncover sets of multiply imaged galaxies and constrain the mass distribution of the clusters. Our SL analysis reveals that both clusters have particularly large Einstein radii ($\theta_E>30"$ for a source redshift of $z_s=2$), providing fairly large areas with high magnifications, useful for high-redshift galaxy searches ($\sim2$ arcmin$^{2}$ with $\mu>5$ to $\sim1$ arcmin$^{2}$ with $\mu>10$, similar to a typical \textit{Hubble Frontier Fields} cluster). We also find that MACS J0308.9+2645 hosts a promising, apparently bright (J$\sim23.2-24.6$ AB), multiply imaged high-redshift candidate at $z\sim6.4$. These images are amongst the brightest high-redshift candidates found in RELICS. Our mass models, including magnification maps, are made publicly available for the community through the Mikulski Archive for Space Telescopes.

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.

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

We present an original phenomenological model to describe the evolution of galaxy number counts, morphologies, and spectral energy distributions across a wide range of redshifts (0.2<z<15) and stellar masses [Log10 M/Msun >6]. Our model follows observed mass and luminosity functions of both star-forming and quiescent galaxies, and reproduces the redshift evolution of colors, sizes, star-formation and chemical properties of the observed galaxy population. Unlike other existing approaches, our model includes a self-consistent treatment of stellar and photoionized gas emission and dust attenuation based on the BEAGLE tool. The mock galaxy catalogs generated with our new model can be used to simulate and optimize extragalactic surveys with future facilities such as the James Webb Space Telescope (JWST), and to enable critical assessments of analysis procedures, interpretation tools, and measurement systematics for both photometric and spectroscopic data. As a first application of this work, we make predictions for the upcoming JWST Advanced Deep Extragalactic Survey (JADES), a joint program of the JWST/NIRCam and NIRSpec Guaranteed Time Observations teams. We show that JADES will detect, with NIRCam imaging, thousands of galaxies at z>6, and tens at z>10 at m_AB<30 (5-sigma) within the 236 arcmin^2 of the survey. The JADES data will enable accurate constraints on the evolution of the UV luminosity function at z>8, and resolve the current debate about the rate of evolution of galaxies at z>8. Ready to use mock catalogs and software to generate new realizations are publicly available as the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) package.

Strong gravitational lensing by galaxy clusters magnifies background galaxies, enhancing our ability to discover statistically significant samples of galaxies at z>6, in order to constrain the high-redshift galaxy luminosity functions. Here, we present the first five lens models out of the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury Program, based on new HST WFC3/IR and ACS imaging of the clusters RXC J0142.9+4438, Abell 2537, Abell 2163, RXC J2211.7-0349, and ACT-CLJ0102-49151. The derived lensing magnification is essential for estimating the intrinsic properties of high-redshift galaxy candidates, and properly accounting for the survey volume. We report on new spectroscopic redshifts of multiply imaged lensed galaxies behind these clusters, which are used as constraints, and detail our strategy to reduce systematic uncertainties due to lack of spectroscopic information. In addition, we quantify the uncertainty on the lensing magnification due to statistical and systematic errors related to the lens modeling process, and find that in all but one cluster, the magnification is constrained to better than 20% in at least 80% of the field of view, including statistical and systematic uncertainties. The five clusters presented in this paper span the range of masses and redshifts of the clusters in the RELICS program. We find that they exhibit similar strong lensing efficiencies to the clusters targeted by the Hubble Frontier Fields within the WFC3/IR field of view. Outputs of the lens models are made available to the community through the Mikulski Archive for Space Telescopes

The study of unobscured active galactic nuclei (AGN) and quasars depends on the reliable decomposition of the light from the AGN point source and the extended host galaxy light. The problem is typically approached using parametric fitting routines using separate models for the host galaxy and the point spread function (PSF). We present a new approach using a Generative Adversarial Network (GAN) trained on galaxy images. We test the method using Sloan Digital Sky Survey (SDSS) r-band images with artificial AGN point sources added which are then removed using the GAN and with parametric methods using GALFIT. When the AGN point source PS is more than twice as bright as the host galaxy, we find that our method, PSFGAN, can recover PS and host galaxy magnitudes with smaller systematic error and a lower average scatter ($49\%$). PSFGAN is more tolerant to poor knowledge of the PSF than parametric methods. Our tests show that PSFGAN is robust against a broadening in the PSF width of $\pm 50\%$ if it is trained on multiple PSF's. We demonstrate that while a matched training set does improve performance, we can still subtract point sources using a PSFGAN trained on non-astronomical images. While initial training is computationally expensive, evaluating PSFGAN on data is more than $40$ times faster than GALFIT fitting two components. Finally, PSFGAN it is more robust and easy to use than parametric methods as it requires no input parameters.

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

We study tidal features (TFs) around galaxies in the REsolved Spectroscopy of a Local VolumE (RESOLVE) survey. Our sample consists of 1048 RESOLVE galaxies that overlap with the DECam Legacy Survey, which reaches an r-band 3$\sigma$ depth of $\sim$27.9 mag arcsec$^{-2}$ for a 100 arcsec$^{2}$ feature. Images were masked, smoothed, and inspected for TFs like streams, shells, or tails/arms. We find TFs in 17$^{\pm 2} \%$ of our galaxies, setting a lower limit on the true frequency. The frequency of TFs in the gas-poor (gas-to-stellar mass ratio $<$ 0.1) subsample is lower than in the gas-rich subsample (13$^{\pm 3} \%$ vs. 19$^{\pm 2} \%$). Within the gas-poor subsample, galaxies with TFs have higher stellar and halo masses, $\sim 3\times$ closer distances to nearest neighbors (in the same group), and possibly fewer group members at fixed halo mass than galaxies without TFs, but similar specific star formation rates. These results suggest TFs in gas-poor galaxies are typically streams/shells from dry mergers or satellite disruption. In contrast, the presence of TFs around gas-rich galaxies does not correlate with stellar or halo mass, suggesting these TFs are often tails/arms from resonant interactions. Similar to TFs in gas-poor galaxies, TFs in gas-rich galaxies imply 1.7x closer nearest neighbors in the same group; however, TFs in gas-rich galaxies are associated with diskier morphologies, higher star formation rates, and higher gas content. In addition to interactions with known neighbors, we suggest that TFs in gas-rich galaxies may arise from accretion of cosmic gas and/or gas-rich satellites below the survey limit.

We present the first spatially-resolved rest-frame UV study of the gravitationally lensed galaxy, the 'Cosmic Horseshoe' (J1148+1930) at z=2.38. Our gravitational lens model shows that the system is made up of four star-forming regions, each ~4-8 kpc^2 in size, from which we extract four spatially exclusive regional spectra. We study the interstellar and wind absorption lines, along with CIII] doublet emission lines, in each region to investigate any variation in emission/absorption line properties. The mapped CIII] emission shows distinct kinematical structure, with velocity offsets of ~+/-50 km/s between regions suggestive of a merging system, and a variation in equivalent width that indicates a change in ionisation parameter and/or metallicity between the regions. Absorption line velocities reveal a range of outflow strengths, with gas outflowing between -200<v(km/s)<-50 relative to the systemic velocity of that region. Interestingly, the strongest gas outflow appears to emanate from the most diffuse star-forming region. The star-formation rates remain relatively constant (~8-16 M_sol/yr), mostly due to large uncertainties in reddening estimates. As such, the outflows appear to be 'global' rather than 'locally' sourced. We measure electron densities with a range of log(Ne)=3.92-4.36 cm^-3, and point out that such high densities may be common when measured using the CIII] doublet due to its large critical density. Overall, our observations demonstrate that while it is possible to trace variations in large scale gas kinematics, detecting inhomogeneities in physical gas properties and their effects on the outflowing gas may be more difficult. This study provides important lessons for the spatially-resolved rest-frame UV studies expected with future observatories, such as JWST.

The most distant galaxies known are at z~10-11, observed 400-500 Myr after the Big Bang. The few z~10-11 candidates discovered to date have been exceptionally small- barely resolved, if at all, by the Hubble Space Telescope. Here we present the discovery of SPT0615-JD, a fortuitous z~10 (z_phot=9.9+/-0.6) galaxy candidate stretched into an arc over ~2.5" by the effects of strong gravitational lensing. Discovered in the Reionization Lensing Cluster Survey (RELICS) Hubble Treasury program and companion S-RELICS Spitzer program, this candidate has a lensed H-band magnitude of 25.7+/-0.1 AB mag. With a magnification of \mu~4-7 estimated from our lens models, the de-lensed intrinsic magnitude is 27.6+/-0.3 AB mag, and the half-light radius is r_e<0.8 kpc, both consistent with other z>9 candidates. The inferred stellar mass (log [M* /M_Sun]=9.7^{+0.7}_{-0.5}) and star formation rate (\log [SFR/M_Sun yr^{-1}]=1.3^{+0.2}_{-0.3}) indicate that this candidate is a typical star-forming galaxy on the z>6 SFR-M* relation. We note that three independent lens models predict two counterimages, at least one of which should be of a similar magnitude to the arc, but these counterimages are not yet detected. Counterimages would not be expected if the arc were at lower redshift. However, the only spectral energy distributions capable of fitting the Hubble and Spitzer photometry well at lower redshifts require unphysical combinations of z~2 galaxy properties. The unprecedented lensed size of this z~10 candidate offers the potential for the James Webb Space Telescope to study the geometric and kinematic properties of a galaxy observed 500 Myr after the Big Bang.

Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.

We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.

The first billion years of the Universe is a pivotal time: stars, black holes (BHs) and galaxies form and assemble, sowing the seeds of galaxies as we know them today. Detecting, identifying and understand- ing the first galaxies and BHs is one of the current observational and theoretical challenges in galaxy formation. In this paper we present a population synthesis model aimed at galaxies, BHs and Active Galactic Nuclei (AGNs) at high redshift. The model builds a population based on empirical relations. Galaxies are characterized by a spectral energy distribution determined by age and metallicity, and AGNs by a spectral energy distribution determined by BH mass and accretion rate. We validate the model against observational constraints, and then predict properties of galaxies and AGN in other wavelength and/or luminosity ranges, estimating the contamination of stellar populations (normal stars and high-mass X-ray binaries) for AGN searches from the infrared to X-rays, and vice-versa for galaxy searches. For high-redshift galaxies, with stellar ages < 1 Gyr, we find that disentangling stellar and AGN emission is challenging at restframe UV/optical wavelengths, while high-mass X-ray binaries become more important sources of confusion in X-rays. We propose a color-color selection in JWST bands to separate AGN vs star-dominated galaxies in photometric observations. We also esti- mate the AGN contribution, with respect to massive, hot, metal-poor stars, at driving high ionization lines, such as C IV and He II. Finally, we test the influence of the minimum BH mass and occupa- tion fraction of BHs in low mass galaxies on the restframe UV/near-IR and X-ray AGN luminosity function.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

We examine the z = 0 group-integrated stellar and cold baryonic (stars + cold atomic gas) mass functions (group SMF and CBMF) and the baryonic collapse efficiency (group cold baryonic to dark matter halo mass ratio) using the RESOLVE and ECO survey galaxy group catalogs and a galform semi-analytic model (SAM) mock catalog. The group SMF and CBMF fall off more steeply at high masses and rise with a shallower low-mass slope than the theoretical halo mass function (HMF). The transition occurs at group-integrated cold baryonic mass M_coldbary ~ 10^11 Msun. The SAM, however, has significantly fewer groups at the transition mass ~ 10^11 Msun and a steeper low-mass slope than the data, suggesting that feedback is too weak in low-mass halos and conversely too strong near the transition mass. Using literature prescriptions to include hot halo gas and potential unobservable galaxy gas produces a group BMF with slope similar to the HMF even below the transition mass. Its normalization is lower by a factor of ~2, in agreement with estimates of warm-hot gas making up the remaining difference. We compute baryonic collapse efficiency with the halo mass calculated two ways, via halo abundance matching (HAM) and via dynamics (extended all the way to three-galaxy groups using stacking). Using HAM, we find that baryonic collapse efficiencies reach a flat maximum for groups across the halo mass range of M_halo ~ 10^11.4-12 Msun, which we label "nascent groups." Using dynamics, however, we find greater scatter in baryonic collapse efficiencies, likely indicating variation in group hot-to-cold baryon ratios. Similarly, we see higher scatter in baryonic collapse efficiencies in the SAM when using its true groups and their group halo masses as opposed to friends-of-friends groups and HAM masses.

Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below $Z_\odot/2$. Such tests are particularly important for interpreting the surprising high-ionization UV line emission detected at $z>6$ in recent years. We present HST/COS ultraviolet spectra of ten nearby metal-poor star-forming galaxies selected to show He II emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance ($7.8<12+\log\mathrm{O/H}<8.5$) and present uniformly large specific star formation rates (sSFR $\sim 10^2$ $\mathrm{Gyr}^{-1}$). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C III] equivalent widths comparable to those seen in systems at $z\sim 6-7$. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below $12+\log\mathrm{O/H}\lesssim 8.0$ ($Z/Z_\odot \lesssim 1/5$) presenting minimal stellar wind features and prominent nebular emission in He II and C IV. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the $\mathrm{He^+}$-ionizing edge relative to H-ionizing flux as metallicity decreases below a fifth solar, well in excess of standard stellar population synthesis predictions. Our results suggest that often neglected sources of energetic radiation such as stripped binary products and very massive O-stars produce a sharper change in the ionizing spectrum with decreasing metallicity than expected. Consequently, nebular emission in C IV and He II powered by these stars may provide useful metallicity constraints in the reionization era.

The ubiquity of Lyman alpha (Ly$\alpha$) emission in a sample of four bright [O III]-strong star-forming galaxies with redshifts above 7 has led to the suggestion that such luminous sources represent a distinct population compared to their fainter, more numerous, counterparts. The presence of Ly$\alpha$ emission within the reionization era could indicate that these sources created early ionized bubbles due to their unusually strong radiation, possibly due to the presence of active galactic nuclei. To test this hypothesis we have secured long integration spectra with XSHOOTER on the VLT for three $z\simeq$7 sources selected to have similar luminosities and prominent excess fluxes in the IRAC 3.6 or 4.5$\mu$m band, usually attributed to strong [O III] emission. We secured additional spectroscopy for one of these galaxies at $z$=7.15 using MOSFIRE at the Keck telescope. For this, the most well-studied source in our sample with the strongest IRAC excess, we detect prominent nebular emission from He II and NV indicative of a non-thermal source. For the other two sources at $z$=6.81 and $z$=6.85, for which no previous optical/near infrared spectroscopy was initially available, Ly$\alpha$ is seen in one and CIII] emission in the other. Although a modest sample, our results further support the hypothesis that the phenomenon of intense [O III] emission is associated preferentially with sources lying in early ionized bubbles. However, even though one of our sources at $z$=7.15 clearly indicates the presence of non-thermal radiation, such ionized bubbles may not uniquely arise in this manner. We discuss the unique advantages of extending such challenging diagnostic studies with JWST.

Recent observations have suggested that the CIII]$\lambda1907/1909$ emission lines could be alternative diagnostic lines for galaxies in the reionization epoch. We use the F128N narrowband filter on the Hubble Space Telescope's ($\it{HST}$) Wide Field Camera 3 (WFC3) to search for CIII] emission in a sample of five galaxies at z = 5.7 in the Subaru Deep Field and the Subaru/XMM-Newton Deep Field. Using the F128N narrowband imaging, together with the broadband imaging, we do not detect CIII] emission for the five galaxies with $J_{\rm{AB}}$ ranging from 24.10 -- 27.00 in our sample. For the brightest galaxy J132416.13+274411.6 in our sample (z = 5.70, $J_{\rm{AB}} = 24.10$), which has a significantly higher signal to noise, we report a CIII] flux of $3.34\pm1.81 \times 10^{-18}$ $\mathrm{erg\ s^{-1}\ cm^{-2}}$, which places a stringent 3-$\rm\sigma$ upper limit of $5.43\times 10^{-18}$ $\mathrm{erg\ s^{-1}\ cm^{-2}}$ on CIII] flux and 6.57 \AA\ on the CIII] equivalent width. Using the stacked image, we put a 3-$\rm\sigma$ upper limit on the mean CIII] flux of $\mathrm{2.55\times10^{-18}\ erg\ s^{-1}\ cm^{-2}}$, and a 3-$\rm\sigma$ upper limit on the mean CIII] equivalent width of 4.20 {\AA} for this sample of galaxies at z = 5.70. Combined with strong CIII] detection reported among high-z galaxies in the literature, our observations suggest that the equivalent widths of CIII] from galaxies at z $>$ 5.70 exhibit a wide range of distribution. Our strong limits on CIII] emission could be used as a guide for future observations in the reionization epoch.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.

We present the Magellan/FIRE detection of highly-ionized CIV 1550 and OIII] 1666 in a deep infrared spectrum of the z=6.11 gravitationally lensed low-mass galaxy RXC J2248.7-4431-ID3, which has previously-known Lyman-alpha. No corresponding emission is detected at the expected location of HeII 1640. The upper limit on HeII paired with detection of OIII] and CIV constrains possible ionization scenarios. Production of CIV and OIII] requires ionizing photons of 2.5-3.5 Ryd, but once in that state their multiplet emission is powered by collisional excitation at lower energies (~0.5 Ryd). As a pure recombination line, HeII emission is powered by 4 Ryd ionizing photons. The data therefore require a spectrum with significant power at 3.5 Ryd but a rapid drop toward 4.0 Ryd. This hard spectrum with a steep drop is characteristic of low-metallicity stellar populations, and less consistent with soft AGN excitation, which features more 4 Ryd photons and hence higher HeII flux. The conclusions based on ratios of metal line detections to Helium non-detection are strengthened if the gas metallicity is low. RXJ2248-ID3 adds to the growing handful of reionization-era galaxies with UV emission line ratios distinct from the general z=2-3 population, in a way that suggests hard ionizing spectra that do not necessarily originate in AGN.

The search for chemically unevolved galaxies remains prevalent in the nearby Universe, mostly because these systems provide excellent proxies for exploring in detail the physics of high-z systems. The most promising candidates are extremely metal-poor galaxies (XMPs), i.e., galaxies with <1/10 solar metallicity. However, due to the bright emission line based search criteria traditionally used to find XMPs, we may not be sampling the full XMP population. In 2014 we reoriented this search using only morphological properties and uncovered a population of ~150 `blue diffuse dwarf (BDD) galaxies', and published a sub-sample of 12 BDD spectra. Here we present optical spectroscopic observations of a larger sample of 51 BDDs, along with their SDSS photometric properties. With our improved statistics, we use direct-method abundances to confirm that BDDs are chemically unevolved (7.43<12+log(O/H)<8.01), with ~20% of our sample classified as being XMP galaxies, and find they are actively forming stars at rates of 1-33x10^-2 M_sol/yr in HII regions randomly embedded in a blue, low-surface brightness continuum. Stellar masses are calculated from population synthesis models and estimated to be in the range log(M_star/M_sol) ~5-9. Unlike other low-metallicity star-forming galaxies, BDDs are in agreement with the mass-metallicity relation at low masses, suggesting they are not accreting large amounts of pristine gas relative to their stellar mass. BDD galaxies appear to be a population of actively star-forming dwarf irregular (dIrr) galaxies who fall within the class of low-surface brightness dIrr galaxies. Their ongoing star-formation and irregular morphology make them excellent analogues for galaxies in the early Universe.

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

We study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

We present the HI mass inventory for the RESOLVE survey, a volume-limited, multi-wavelength census of >1500 z=0 galaxies spanning diverse environments and complete in baryonic mass down to dwarfs of 10^9 Msun. This first 21cm data release provides robust detections or strong upper limits (1.4M_HI < 5 to 10% of stellar mass M_stars) for 94% of RESOLVE. We examine global atomic gas-to-stellar mass ratios (G/S) in relation to galaxy environment using several metrics: group dark matter halo mass M_h , central/satellite designation, relative mass density of the cosmic web, and distance to nearest massive group. We find that at fixed M_stars, satellites have decreasing G/S with increasing M_h starting clearly at M_h = 10^12 Msun, suggesting the presence of starvation and/or stripping mechanisms associated with halo gas heating in intermediate-mass groups. The analogous relationship for centrals is uncertain because halo abundance matching builds in relationships between central G/S, stellar mass, and halo mass, which depend on the integrated group property used as a proxy for halo mass (stellar or baryonic mass). On larger scales G/S trends are less sensitive to the abundance matching method. At fixed M_h < 10^12 Msun, the fraction of gas-poor centrals increases with large-scale structure density. In overdense regions, we identify a rare population of gas-poor centrals in low-mass (M_h < 10^11.4 Msun) halos primarily located within 1.5 times the virial radius of more massive (M_h > 10^12 Msun) halos, suggesting that gas stripping and/or starvation may be induced by interactions with larger halos or the surrounding cosmic web. We find that the detailed relationship between G/S and environment varies when we examine different subvolumes of RESOLVE independently, which we suggest may be a signature of assembly bias.

[Abridged] We present a robust measurement of the rest-frame UV luminosity function (LF) and its evolution during the peak epoch of cosmic star formation at 1<z<3. We use our deep near ultraviolet imaging from WFC3/UVIS on the Hubble Space Telescope (HST) and existing ACS/WFC and WFC3/IR imaging of three lensing galaxy clusters, Abell 2744 and MACSJ0717 from the Hubble Frontier Field survey and Abell 1689. We use photometric redshifts to identify 780 ultra-faint galaxies with $M_{UV}$<-12.5 AB mag at 1<z<3. From these samples, we identified 5 new, faint, multiply imaged systems in A1689. We compute the rest-frame UV LF and find the best-fit faint-end slopes of $\alpha=-1.56\pm0.04$, $\alpha=-1.72\pm0.04$ and $\alpha=-1.94\pm0.06$ at 1.0<z<1.6, 1.6<z<2.2 and 2.2<z<3.0, respectively. Our results demonstrate that the UV LF becomes steeper from z\sim1.3 to z\sim2.6 with no sign of a turnover down to $M_{UV}=-14$ AB mag. We further derive the UV LFs using the Lyman break "dropout" selection and confirm the robustness of our conclusions against different selection methodologies. Because the sample sizes are so large, and extend to such faint luminosities, the statistical uncertainties are quite small, and systematic uncertainties (due to the assumed size distribution, for example), likely dominate. If we restrict our analysis to galaxies and volumes above > 50% completeness in order to minimize these systematics, we still find that the faint-end slope is steep and getting steeper with redshift, though with slightly shallower (less negative) values ($\alpha=-1.55\pm0.06$, $-1.69\pm0.07$ and $-1.79\pm0.08$ for $z\sim1.3$, 1.9 and 2.6, respectively). Finally, we conclude that the faint star-forming galaxies with UV magnitudes of $-18.5<M_{UV}<-12.5$ covered in this study, produce the majority (55%-60%) of the unobscured UV luminosity density at 1<z<3.

The nature of warm, ionized gas outside of galaxies may illuminate several key galaxy evolutionary processes. A serendipitous observation by the MaNGA survey has revealed a large, asymmetric H$\alpha$ complex with no optical counterpart that extends $\approx8"$ ($\approx6.3$ kpc) beyond the effective radius of a dusty, starbursting galaxy. This H$\alpha$ extension is approximately three times the effective radius of the host galaxy and displays a tail-like morphology. We analyze its gas-phase metallicities, gaseous kinematics, and emission-line ratios, and discuss whether this H$\alpha$ extension could be diffuse ionized gas, a gas accretion event, or something else. We find that this warm, ionized gas structure is most consistent with gas accretion through recycled wind material, which could be an important process that regulates the low-mass end of the galaxy stellar mass function.

The fraction of ionizing photons escaping from high-redshift star-forming galaxies remains a key obstacle in evaluating whether galaxies were the primary agents of cosmic reionization. We previously proposed using the covering fraction of low-ionization gas, measured via deep absorption line spectroscopy, as a proxy. We now present a significant update, sampling seven gravitationally-lensed sources at $4<z<5$. We show that the absorbing gas in our sources is spatially inhomogeneous with a median covering fraction of 66\%. Correcting for reddening according to a dust-in-cloud model, this implies an estimated absolute escape fraction of $\simeq19\pm6$\%. With possible biases and uncertainties, collectively we find the average escape fraction could be reduced to no less than 11\%, excluding the effect of spatial variations. For one of our lensed sources, we have sufficient signal/noise to demonstrate the presence of such spatial variations and scatter in its dependence on the Ly$\alpha$ equivalent width, consistent with recent simulations. If this source is typical, our lower limit to the escape fraction could be reduced by a further factor $\simeq$2. Across our sample, we find a modest anti-correlation between the inferred escape fraction and the local star formation rate, consistent with a time delay between a burst and leaking Lyman continuum photons. Our analysis demonstrates considerable variations in the escape fraction consistent with being governed by the small-scale behavior of star-forming regions, whose activities fluctuate over short timescales. This supports the suggestion that the escape fraction may increase toward the reionization era when star formation becomes more energetic and burst-like.

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

In this work, we present galaxy stellar and baryonic (stars plus cold gas) mass functions (SMF and BMF) and their halo mass dependence for two volume-limited data sets. The first, RESOLVE-B, coincides with the Stripe 82 footprint and is extremely complete down to baryonic mass Mbary ~ 10^9.1 Msun, probing the gas-rich dwarf regime below Mbary ~ 10^10 Msun. The second, ECO, covers a ~40 times larger volume (containing RESOLVE-A) and is complete to Mbary ~ 10^9.4 Msun. To construct the SMF and BMF we implement a new "cross-bin sampling" technique with Monte Carlo sampling from the full likelihood distributions of stellar or baryonic mass. Our SMFs exhibit the "plateau" feature starting below Mstar ~ 10^10 Msun that has been described in prior work. However, the BMF fills in this feature and rises as a straight power law below ~10^10 Msun, as gas-dominated galaxies become the majority of the population. Nonetheless, the low-mass slope of the BMF is not as steep as that of the theoretical dark matter halo MF. Moreover, we assign group halo masses by abundance matching, finding that the SMF and BMF separated into four physically motivated halo mass regimes reveal complex structure underlying the simple shape of the overall MFs. In particular, the satellite MFs are depressed below the central galaxy MF "humps" in groups with mass <10^13.5 Msun yet rise steeply in clusters. Our results suggest that satellite destruction and/or stripping are active from the point of nascent group formation. We show that the key role of groups in shaping MFs enables reconstruction of a given survey's SMF or BMF based on its group halo mass distribution.

We discuss spatially resolved emission line spectroscopy secured for a total sample of 15 gravitationally lensed star-forming galaxies at a mean redshift of $z\simeq2$ based on Keck laser-assisted adaptive optics observations undertaken with the recently-improved OSIRIS integral field unit (IFU) spectrograph. By exploiting gravitationally lensed sources drawn primarily from the CASSOWARY survey, we sample these sub-L$^{\ast}$ galaxies with source-plane resolutions of a few hundred parsecs ensuring well-sampled 2-D velocity data and resolved variations in the gas-phase metallicity. Such high spatial resolution data offers a critical check on the structural properties of larger samples derived with coarser sampling using multiple-IFU instruments. We demonstrate how kinematic complexities essential to understanding the maturity of an early star-forming galaxy can often only be revealed with better sampled data. Although we include four sources from our earlier work, the present study provides a more representative sample unbiased with respect to emission line strength. Contrary to earlier suggestions, our data indicates a more diverse range of kinematic and metal gradient behavior inconsistent with a simple picture of well-ordered rotation developing concurrently with established steep metal gradients in all but merging systems. Comparing our observations with the predictions of hydrodynamical simulations suggests that gas and metals have been mixed by outflows or other strong feedback processes, flattening the metal gradients in early star-forming galaxies.

The Smith Cloud is a gaseous high-velocity cloud (HVC) in an advanced state of accretion, only 2.9 kpc below the Galactic plane and due to impact the disk in 27 Myr. It is unique among HVCs in having a known distance (12.4+/-1.3 kpc) and a well-constrained 3D velocity (296 km/s), but its origin has long remained a mystery. Here we present the first absorption-line measurements of its metallicity, using HST/COS UV spectra of three AGN lying behind the Cloud together with Green Bank Telescope 21 cm spectra of the same directions. Using Voigt-profile fitting of the S II 1250, 1253, 1259 triplet together with ionization corrections derived from photoionization modeling, we derive the sulfur abundance in each direction; a weighted average of the three measurements gives [S/H]=-0.28+/-0.14, or 0.53+0.21-0.15 solar metallicity. The finding that the Smith Cloud is metal-enriched lends support to scenarios where it represents recycled Galactic material rather than the remnant of a dwarf galaxy or accreting intergalactic gas. The metallicity and trajectory of the Cloud are both indicative of an origin in the outer disk. However, its large mass and prograde kinematics remain to be fully explained. If the cloud has accreted cooling gas from the corona during its fountain trajectory, as predicted in recent theoretical work, its current mass would be higher than its launch mass, alleviating the mass concern.

Using the magnetofluid unification framework, we show that the accretion disk plasma (embedded in the background geometry of a blackhole) can relax to a class of states known as the Beltrami-Bernoulli (BB) equilibria. Modeling the disk plasma as a Hall MHD system, we find that the space-time curvature can significantly alter the magnetic/velocity decay rate as we move away from the compact object; the velocity profiles in BB states, for example, deviate substantially from the predicted corresponding geodesic velocity profiles. These departures imply a rich interplay of plasma dynamics and general relativity revealed by examining the corresponding Bernoulli condition representing "homogeneity" of total energy. The relaxed states have their origin in the constraints provided by the two helicity invariants of Hall MHD. These helicities conspire to introduce a new oscillatory length scale into the system that is strongly influenced by relativistic and thermal effects.

We study the relationships between galaxy environments and galaxy properties related to disk (re)growth, considering two highly complete samples that are approximately baryonic mass limited into the high-mass dwarf galaxy regime, the Environmental COntext (ECO) catalog (data release herein) and the B-semester region of the REsolved Spectroscopy Of a Local VolumE (RESOLVE) survey. We quantify galaxy environments using both group identification and smoothed galaxy density field methods. We use by-eye and quantitative morphological classifications plus atomic gas content measurements and estimates. We find that blue early-type (E/S0) galaxies, gas-dominated galaxies, and UV-bright disk host galaxies all become distinctly more common below group halo mass ~10^11.5 Msun, implying that this low group halo mass regime may be a preferred regime for significant disk growth activity. We also find that blue early-type and blue late-type galaxies inhabit environments of similar group halo mass at fixed baryonic mass, consistent with a scenario in which blue early types can regrow late-type disks. In fact, we find that the only significant difference in the typical group halo mass inhabited by different galaxy classes is for satellite galaxies with different colors, where at fixed baryonic mass red early and late types have higher typical group halo masses than blue early and late types. More generally, we argue that the traditional morphology-environment relation (i.e., that denser environments tend to have more early types) can be largely attributed to the morphology-galaxy mass relation for centrals and the color-environment relation for satellites.

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

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

We present custom-processed UV, optical, and near-IR photometry for the RESOLVE survey, a volume-limited census of stellar, gas, and dynamical mass within two subvolumes of the nearby universe (RESOLVE-A and -B), complete down to baryonic mass ~10^9.1-9.3 Msun. In contrast to standard pipeline photometry (e.g., SDSS), our photometry uses optimal background subtraction, avoids suppressing color gradients, and includes systematic errors. With these improvements, we measure brighter magnitudes, larger radii, bluer colors, and a real increase in scatter around the red sequence. Combining stellar masses from our photometry with the RESOLVE-A HI mass census, we create volume-limited calibrations of the photometric gas fractions (PGF) technique, which predicts gas-to-stellar mass ratios (G/S) from galaxy colors and optional additional parameters. We analyze G/S-color residuals vs. potential third parameters, finding that axial ratio is the best independent and physically meaningful third parameter. We define a "modified color" from planar fits to G/S as a function of both color and axial ratio. In the complete galaxy population, upper limits on G/S bias linear and planar fits. We therefore model the entire PGF probability density field, enabling iterative statistical modeling of upper limits and prediction of full G/S probability distributions for individual galaxies, with two-component structure for red colors. We use the RESOLVE-B 21cm census to test several PGF calibrations, finding that most systematically under- or overestimate gas masses, but the full probability density method performs well.

We report five Local Volume dwarf galaxies (two of which are presented here for the first time) uncovered during a comprehensive archival search for optical counterparts to ultra-compact high velocity clouds (UCHVCs). The UCHVC population of HI clouds are thought to be candidate gas-rich, low mass halos at the edge of the Local Group and beyond, but no comprehensive search for stellar counterparts to these systems has been presented. Careful visual inspection of all publicly available optical and ultraviolet imaging at the position of the UCHVCs revealed six blue, diffuse counterparts with a morphology consistent with a faint dwarf galaxy beyond the Local Group. Optical spectroscopy of all six candidate dwarf counterparts show that five have an H$\alpha$-derived velocity consistent with the coincident HI cloud, confirming their association, the sixth diffuse counterpart is likely a background object. The size and luminosity of the UCHVC dwarfs is consistent with other known Local Volume dwarf irregular galaxies. The gas fraction ($M_{HI}/M_{star}$) of the five dwarfs are generally consistent with that of dwarf irregular galaxies in the Local Volume, although ALFALFA-Dw1 (associated with ALFALFA UCHVC HVC274.68+74.70$-$123) has a very high $M_{HI}/M_{star}$$\sim$40. Despite the heterogenous nature of our search, we demonstrate that the current dwarf companions to UCHVCs are at the edge of detectability due to their low surface brightness, and that deeper searches are likely to find more stellar systems. If more sensitive searches do not reveal further stellar counterparts to UCHVCs, then the dearth of such systems around the Local Group may be in conflict with $\Lambda$CDM simulations.

The star formation histories (SFHs) of dwarf galaxies are thought to be \emph{bursty}, with large -- order of magnitude -- changes in the star formation rate on timescales similar to O-star lifetimes. As a result, the standard interpretations of many galaxy observables (which assume a slowly varying SFH) are often incorrect. Here, we use the SFHs from hydro-dynamical simulations to investigate the effects of bursty SFHs on sample selection and interpretation of observables and make predictions to confirm such SFHs in future surveys. First, because dwarf galaxies' star formation rates change rapidly, the mass-to-light ratio is also changing rapidly in both the ionizing continuum and, to a lesser extent, the non-ionizing UV continuum. Therefore, flux limited surveys are highly biased toward selecting galaxies in the \emph{burst} phase and very deep observations are required to detect all dwarf galaxies at a given stellar mass. Second, we show that a $\log_{10}[\nu L_{\nu}(1500{\rm \AA})/L_{{\rm H}\alpha}]>2.5$ implies a very recent quenching of star formation and can be used as evidence of stellar feedback regulating star formation. Third, we show that the ionizing continuum can be significantly higher than when assuming a constant SFH, which can affect the interpretation of nebular emission line equivalent widths and direct ionizing continuum detections. Finally, we show that a star formation rate estimate based on continuum measurements only (and not on nebular tracers such as the hydrogen Balmer lines) will not trace the rapid changes in star formation and will give the false impression of a star-forming main sequence with low dispersion.

We report results of a search for CIII] $\lambda \lambda$1907,1909 {\AA} emission using Keck's MOSFIRE spectrograph in a sample of 7 $z_{phot}\sim7-8$ candidates ($H\sim27$) lensed by the Hubble Frontier Field cluster Abell 2744. Earlier work has suggested the promise of using the CIII] doublet for redshift confirmation of galaxies in the reionization era given $Ly\alpha$ ($\lambda$1216 {\AA}) is likely attenuated by the neutral intergalactic medium. The primary challenge of this approach is the feasibility of locating CIII] emission without advanced knowledge of the spectroscopic redshift. With an integration time of 5 hours in the H-band, we reach a $5\sigma$ median flux limit (in between the skylines) of $1.5\times10^{-18}$ ergs cm$^{-2}$ sec$^{-1}$ but no convincing CIII] emission was found. We also incorporate preliminary measurements from two other CLASH/HFF clusters in which, similarly, no line was detected, but these were observed to lesser depth. Using the known distribution of OH emission and the photometric redshift likelihood distribution of each lensed candidate, we present statistical upper limits on the mean total CIII] rest-frame equivalent width for our $z\simeq7-8$ sample. For a signal/noise ratio of 5, we estimate the typical CIII] doublet rest-frame equivalent width is, with 95\% confidence, $<26\pm5$ {\AA}. Although consistent with the strength of earlier detections in brighter objects at $z\simeq6-7$, our study illustrates the necessity of studying more luminous or strongly-lensed examples prior to the launch of the James Webb Space Telescope.

We present Keck/MOSFIRE observations of UV metal lines in four bright gravitationally-lensed z~6-8 galaxies behind the cluster Abell 1703. The spectrum of A1703-zd6, a highly-magnified star forming galaxy with a Lyman-alpha redshift of z=7.045, reveals a confident detection of the nebular CIV emission line (unresolved with FWHM < 125 km/s). UV metal lines are not detected in the three other galaxies. At z~2-3, nebular CIV emission is observed in just 1% of UV-selected galaxies. The presence of strong CIV emission in one of the small sample of galaxies targeted in this paper may indicate hard ionizing spectra are more common at z~7. The total estimated equivalent width of the CIV doublet (38 A) and CIV/Lyman-alpha flux ratio (0.3) are comparable to measurements of narrow-lined AGNs. Photoionization models show that the nebular CIV line can also be reproduced by a young stellar population, with very hot metal poor stars dominating the photon flux responsible for triply ionizing carbon. Regardless of the origin of the CIV, we show that the ionizing spectrum of A1703-zd6 is different from that of typical galaxies at z~2, producing more H ionizing photons per unit 1500A luminosity and a larger flux density at 30-50 eV. If such extreme radiation fields are typical in UV-selected systems at z>7, it would indicate that reionization-era galaxies are more efficient ionizing agents than previously thought. Alternatively, we suggest that the small sample of Lyman-alpha emitters at z>7 may trace a rare population with intense radiation fields capable of ionizing their surrounding hydrogen distribution. Additional constraints on high ionization emission lines in galaxies with and without Lyman-alpha detections will help clarify whether hard ionizing spectra are common in the reionization era.

Extremely metal poor (XMP) galaxies are known to be very rare, despite the large numbers of low-mass galaxies predicted by the local galaxy luminosity function. This paper presents a sub-sample of galaxies that were selected via a morphology-based search on SDSS images with the aim of finding these elusive XMP galaxies. By using the recently discovered extremely metal-poor galaxy, Leo P, as a guide, we obtained a collection of faint, blue systems, each with isolated HII regions embedded in a diffuse continuum, that have remained optically undetected until now. Here we show the first results from optical spectroscopic follow-up observations of 12 of ~100 of these blue, diffuse dwarf (BDD) galaxies yielded by our search algorithm. Oxygen abundances were obtained via the direct method for eight galaxies, and found to be in the range 7.45<12+log(O/H)<8.0, with two galaxies being classified as XMPs. All BDDs were found to currently have a young star-forming population (<10 Myr) and relatively high ionisation parameters of their HII regions. Despite their low luminosities (-11< M_B < -18) and low surface brightnesses (~23-25 mag arcsec^-2), the galaxies were found to be actively star-forming, with current star-formation rates between 0.0003 and 0.078 Msol/yr. From our current subsample, BDD galaxies appear to be a population of non-quiescent dwarf irregular (dIrr) galaxies, or the diffuse counterparts to blue compact galaxies (BCDs) and as such may bridge the gap between these two populations. Our search algorithm demonstrates that morphology-based searches are successful in uncovering more diffuse metal-poor star-forming galaxies, which traditional emission-line based searches overlook.

Motivated by the idea that a subset of HVCs trace dark matter substructure in the Local Group, we search for signs of star formation in the Smith Cloud, a nearby ~2x10^6 Msun HVC currently falling into the Milky Way. Using GALEX NUV and WISE/2MASS NIR photometry, we apply a series of color and apparent magnitude cuts to isolate candidate O and B stars that are plausibly associated with the Smith Cloud. We find an excess of stars along the line of sight to the cloud, but not at a statistically significant level relative to a control region. The number of stars found in projection on the cloud after removing an estimate of the contamination by the Milky Way implies an average star formation rate surface density of 10^(-4.8 +/- 0.3) Msun yr^(-1) kpc^(-2), assuming the cloud has been forming stars at a constant rate since its first passage through the Milky Way ~70 Myr ago. This value is consistent with the star formation rate expected based on the average gas density of the cloud. We also discuss how the newly discovered star forming galaxy Leo P has very similar properties to the Smith Cloud, but its young stellar population would not have been detected at a statistically significant level using our method. Thus, we cannot yet rule out the idea that the Smith Cloud is really a dwarf galaxy.

Strong gravitational lensing provides a powerful means for studying faint galaxies in the distant universe. By magnifying the apparent brightness of background sources, massive clusters enable the detection of galaxies fainter than the usual sensitivity limit for blank fields. However, this gain in effective sensitivity comes at the cost of a reduced survey volume and, in this {\it Letter}, we demonstrate there is an associated increase in the cosmic variance uncertainty. As an example, we show that the cosmic variance uncertainty of the high redshift population viewed through the Hubble Space Telescope Frontier Field cluster Abell 2744 increases from ~35% at redshift z~7 to >~65% at z~10. Previous studies of high redshift galaxies identified in the Frontier Fields have underestimated the cosmic variance uncertainty that will affect the ultimate constraints on both the faint end slope of the high-redshift luminosity function and the cosmic star formation rate density, key goals of the Frontier Field program.

We present the results of a new, ultra-deep, near-infrared imaging survey executed with the Hawk-I imager at the ESO VLT, of which we make all the data public. This survey, named HUGS (Hawk-I UDS and GOODS Survey), provides deep, high-quality imaging in the K and Y bands over the CANDELS UDS and GOODS-South fields. We describe here the survey strategy, the data reduction process, and the data quality. HUGS delivers the deepest and highest quality K-band images ever collected over areas of cosmological interest, and ideally complements the CANDELS data set in terms of image quality and depth. The seeing is exceptional and homogeneous, confined to the range 0.38"-0.43". In the deepest region of the GOODS-S field, (which includes most of the HUDF) the K-band exposure time exceeds 80 hours of integration, yielding a 1-sigma magnitude limit of ~28.0 mag/sqarcsec. In the UDS field the survey matches the shallower depth of the CANDELS images reaching a 1-sigma limit per sq.arcsec of ~27.3mag in the K band and ~28.3mag in the Y-band, We show that the HUGS observations are well matched to the depth of the CANDELS WFC3/IR data, since the majority of even the faintest galaxies detected in the CANDELS H-band images are also detected in HUGS. We present the K-band galaxy number counts produced by combining the HUGS data from the two fields. We show that the slope of the number counts depends sensitively on the assumed distribution of galaxy sizes, with potential impact on the estimated extra-galactic background light (abridged).

Deep spectroscopic observations of z~6.5 galaxies have revealed a marked decline with increasing redshift in the detectability of Lyman-alpha emission. While this may offer valuable insight into the end of the reionisation process, it presents a fundamental challenge to the detailed spectroscopic study of the many hundreds of photometrically-selected distant sources now being found via deep HST imaging, and particularly those bright sources viewed through foreground lensing clusters. In this paper we demonstrate the validity of a new way forward via the convincing detection of an alternative diagnostic line, CIII]1909, seen in spectroscopic exposures of two star forming galaxies at z=6.029 and 7.213. The former detection is based on a 3.5 hour X-shooter spectrum of a bright (J=25.2) gravitationally-lensed galaxy behind the cluster Abell 383. The latter detection is based on a 4.2 hour MOSFIRE spectra of one of the most distant spectroscopically confirmed galaxies, GN-108036, with J=25.2. Both targets were chosen for their continuum brightness and previously-known redshift (based on Lyman-alpha), ensuring that any CIII] emission would be located in a favorable portion of the near-infrared sky spectrum. We compare our CIII] and Lyman-alpha equivalent widths in the context of those found at z~2 from earlier work and discuss the motivation for using lines other than Lyman-alpha to study galaxies in the reionisation era.

We present deep spectroscopy of 17 very low mass (M* ~ 2.0x10^6 Msun to 1.4x10^9 Msun) and low luminosity (M_UV ~ -13.7 to -19.9) gravitationally lensed galaxies in the redshift range z~1.5-3.0. Deep rest-frame ultraviolet spectra reveal large equivalent width emission from numerous lines (NIV], OIII], CIV, Si III], CIII]) which are rarely seen in individual spectra of more massive star forming galaxies. CIII] is detected in 16 of 17 low mass star forming systems with rest-frame equivalent widths as large as 13.5 Angstroms. Nebular CIV emission is present in the most extreme CIII] emitters, requiring an ionizing source capable of producing a substantial component of photons with energies in excess of 47.9 eV. Photoionization models support a picture whereby the large equivalent widths are driven by the increased electron temperature and enhanced ionizing output arising from metal poor gas and stars, young stellar populations, and large ionization parameters. The young ages implied by the emission lines and continuum SEDs indicate that the extreme line emitters in our sample are in the midst of a significant upturn in their star formation activity. The low stellar masses, blue UV colors, and large sSFRs of our sample are similar to those of typical z>6 galaxies. Given the strong attenuation of Ly-alpha in z>6 galaxies we suggest that CIII] is likely to provide our best probe of early star forming galaxies with ground-based spectrographs and one of the most efficient means of confirming z>10 galaxies with the James Webb Space Telescope.

The redshift-dependent fraction of color-selected galaxies revealing Lyman alpha emission has become the most valuable constraint on the evolving neutrality of the early intergalactic medium. However, in addition to resonant scattering by neutral gas, the visibility of Lyman alpha is also dependent on the intrinsic properties of the host galaxy, including its stellar population, dust content and the nature of outflowing gas. Taking advantage of significant progress we have made in determining the line emitting properties of $z \simeq 4-6$ galaxies, we propose an improved method, based on using the measured slopes of the rest-frame ultraviolet continua of galaxies, to interpret the growing body of near-infrared spectra of $z>7$ galaxies in order to take into account these host galaxy dependencies. In a first application of our new method, we demonstrate its potential via a new spectroscopic survey of $7<z<8$ galaxies undertaken with the Keck MOSFIRE spectrograph. Together with earlier published data our data provides improved estimates of the evolving visibility of Lyman alpha, particularly at redshift $z\simeq 8$. As a byproduct, we also present a new line emitting galaxy at a redshift $z=7.62$ which supersedes an earlier redshift record. We discuss the improving constraints on the evolving neutral fraction over $6<z<8$ and the implications for cosmic reionization.

We present a homogeneous, detailed analysis of the spectral energy distribution (SED) of $\sim$ 1700 LBGs from the GOODS-MUSIC catalogue with deep multi-wavelength photometry from $U$ band to 8 $\mu$m to determine stellar mass, age, dust attenuation, and star formation rate. Using our SED fitting tool, which takes into account nebular emission, we explore a wide parameter space. We also explore a set of different star formation histories. Nebular emission is found to significantly affect the determination of the physical parameters for the majority of $z \sim $ 3--6 LBGs. We identify two populations of galaxies by determining the importance of the contribution of emission lines to broadband fluxes. We find that $\sim$ $65\%$ of LBGs show detectable signs of emission lines, whereas $\sim$ $35\%$ show weak or no emission lines. This distribution is found over the entire redshift range. We interpret these groups as actively star forming and more quiescent LBGs, respectively. We find that it is necessary to considerer SED fits with very young ages ($<50$ Myr) to reproduce some colours affected by strong emission lines. Other arguments favouring episodic star formation and relatively short star formation timescales are also discussed. Considering nebular emission generally leads to a younger age, lower stellar mass, higher dust attenuation, higher star formation rate, and a large scatter in the SFR-$M_{\star}$ relation. Our analysis yields a trend of increasing specific star formation rate with redshift, as predicted by recent galaxy evolution models. The physical parameters of approximately two thirds of high redshift galaxies are significantly modified when we account for nebular emission. The SED models which include nebular emission shed new light on the properties of LBGs with numerous important implications.

We present direct constraints on the CO luminosity function at high redshift and the resulting cosmic evolution of the molecular gas density, $\rho_{\rm H2}$(z), based on a blind molecular line scan in the Hubble Deep Field North (HDF-N) using the IRAM Plateau de Bure Interferometer. Our line scan of the entire 3mm window (79-115 GHz) covers a cosmic volume of ~7000 Mpc$^3$, and redshift ranges z<0.45, 1.01<z<1.89 and z>2. We use the rich multiwavelength and spectroscopic database of the HDF-N to derive some of the best constraints on CO luminosities in high redshift galaxies to date. We combine the blind CO detections in our molecular line scan (presented in a companion paper) with stacked CO limits from galaxies with available spectroscopic redshifts (slit or mask spectroscopy from Keck and grism spectroscopy from HST) to give first blind constraints on high-z CO luminosity functions and the cosmic evolution of the H2 mass density $\rho_{\rm H2}$(z) out to redshifts z~3. A comparison to empirical predictions of $\rho_{\rm H2}$(z) shows that the securely detected sources in our molecular line scan already provide significant contributions to the predicted $\rho_{\rm H2}$(z) in the redshift bins <z>~1.5 and <z>~2.7. Accounting for galaxies with CO luminosities that are not probed by our observations results in cosmic molecular gas densities $\rho_{\rm H2}$(z) that are higher than current predictions. We note however that the current uncertainties (in particular the luminosity limits, number of detections, as well as cosmic volume probed) are significant, a situation that is about to change with the emerging ALMA observatory.

We present a molecular line scan in the Hubble Deep Field North (HDF-N) that covers the entire 3mm window (79-115 GHz) using the IRAM Plateau de Bure Interferometer. Our CO redshift coverage spans z<0.45, 1<z<1.9 and all z>2. We reach a CO detection limit that is deep enough to detect essentially all z>1 CO lines reported in the literature so far. We have developed and applied different line searching algorithms, resulting in the discovery of 17 line candidates. We estimate that the rate of false positive line detections is ~2/17. We identify optical/NIR counterparts from the deep ancillary database of the HDF-N for seven of these candidates and investigate their available SEDs. Two secure CO detections in our scan are identified with star-forming galaxies at z=1.784 and at z=2.047. These galaxies have colors consistent with the `BzK' color selection and they show relatively bright CO emission compared with galaxies of similar dust continuum luminosity. We also detect two spectral lines in the submillimeter galaxy HDF850.1 at z=5.183. We consider an additional 9 line candidates as high quality. Our observations also provide a deep 3mm continuum map (1-sigma noise level = 8.6 $\mu$Jy/beam). Via a stacking approach, we find that optical/MIR bright galaxies contribute only to <50% of the SFR density at 1<z<3, unless high dust temperatures are invoked. The present study represents a first, fundamental step towards an unbiased census of molecular gas in `normal' galaxies at high-z, a crucial goal of extragalactic astronomy in the ALMA era.

WFIRST is the highest priority space mission of the Decadal review, however, it is unlikely to begin in this decade primarily due to a anticipated NASA budget that is unlikely to have sufficient resources to fund such a mission. For this reason we present a lower cost mission that accomplishes all of the WFIRST science as described in the Design Reference Mission 1 with a probe class design. This is effort is motivated by a desire to begin WFIRST in a timely manner and within a budget that can fit within the assets available to NASA on a realistic basis. The design utilizes dichroics to form four focal planes all having the same field of view to use the majority of available photons from a 1.2 meter telescope.

We report the discovery of a planetary-mass companion, HD 106906 b, with the new Magellan Adaptive Optics (MagAO) + Clio2 system. The companion is detected with Clio2 in three bands: $J$, $K_S$, and $L^\prime$, and lies at a projected separation of 7.1'' (650 AU). It is confirmed to be comoving with its $13\pm2$ Myr-old F5 host using Hubble Space Telescope/Advanced Camera for Surveys astrometry over a time baseline of 8.3 yr. DUSTY and COND evolutionary models predict the companion's luminosity corresponds to a mass of $11\pm2 M_{Jup}$, making it one of the most widely separated planetary-mass companions known. We classify its Magellan/Folded-Port InfraRed Echellette $J/H/K$ spectrum as L$2.5\pm1$; the triangular $H$-band morphology suggests an intermediate surface gravity. HD 106906 A, a pre-main-sequence Lower Centaurus Crux member, was initially targeted because it hosts a massive debris disk detected via infrared excess emission in unresolved Spitzer imaging and spectroscopy. The disk emission is best fit by a single component at 95 K, corresponding to an inner edge of 15-20 AU and an outer edge of up to 120 AU. If the companion is on an eccentric ($e>0.65$) orbit, it could be interacting with the outer edge of the disk. Close-in, planet-like formation followed by scattering to the current location would likely disrupt the disk and is disfavored. Furthermore, we find no additional companions, though we could detect similar-mass objects at projected separations $>35$ AU. In situ formation in a binary-star-like process is more probable, although the companion-to-primary mass ratio, at $<1\%$, is unusually small.

We identified the z~2 Lyman break galaxies using deep HST ultraviolet (F275W/F336W) imaging of Abell 1689. Because of the imaging depth and the large magnification provided by the cluster, we detect galaxies 100x fainter (-19.5< M_1500 <-13) than previous surveys at this redshift. We are able to calculate the intrinsic sensitivity of the observations as a function of source plane position, allowing determinations of effective volume as a function of luminosity. We fit the faint-end slope of the luminosity function to be alpha = -1.74 +/-0.08, consistent with the values obtained for 2.5 < z < 6. There is no turnover in the luminosity function down to MUV = -13. The trend of increasingly redder UV spectral slopes with luminosity at higher redshifts is observed in our sample, but with redder slopes at all luminosities and average reddening of < E(B - V) >= 0.15. We assume the stars in these galaxies are metal poor (0.2Z_sun) compared to their brighter counterparts (Z_sun), resulting in bluer assumed intrinsic UV slopes and larger derived dust extinction. The total UV luminosity density at z ~ 2 is 4.31x10^26 erg/s/Hz/Mpc^3, more than 70% of which is emitted by galaxies in the luminosity range of our sample. We determine the star formation rate density at z ~ 2 (assuming constant dust extinction correction of 4.2 over all luminosities and a Kroupa IMF) of 0.148 M/yr/Mpc^3, significantly higher than previous determinations because of the additional population of fainter galaxies and the larger dust correction factors.[abridged]

We relate transitions in galaxy structure and gas content to refueling, here defined to include both the external gas accretion and the internal gas processing needed to renew reservoirs for star formation. We analyze two z=0 data sets: a high-quality ~200-galaxy sample (the Nearby Field Galaxy Survey, data release herein) and a volume-limited ~3000-galaxy sample with reprocessed archival data. Both reach down to baryonic masses ~10^9Msun and span void-to-cluster environments. Two mass-dependent transitions are evident: (i) below the "gas-richness threshold" scale (V~125km/s), gas-dominated quasi-bulgeless Sd--Im galaxies become numerically dominant, while (ii) above the "bimodality" scale (V~200km/s), gas-starved E/S0s become the norm. Notwithstanding these transitions, galaxy mass (or V as its proxy) is a poor predictor of gas-to-stellar mass ratio M_gas/M_*. Instead, M_gas/M_* correlates well with the ratio of a galaxy's stellar mass formed in the last Gyr to its preexisting stellar mass, such that the two ratios have numerically similar values. This striking correspondence between past-averaged star formation and current gas richness implies routine refueling of star-forming galaxies on Gyr timescales. We argue that this refueling underlies the tight M_gas/M_* vs. color correlations often used to measure "photometric gas fractions." Furthermore, the threshold and bimodality scale transitions reflect mass-dependent demographic shifts between three refueling regimes --- accretion dominated, processing dominated, and quenched. In this picture, gas-dominated dwarfs are explained not by inefficient star formation but by overwhelming gas accretion, which fuels stellar mass doubling in <~1Gyr. Moreover, moderately gas-rich bulged disks such as the Milky Way are transitional, becoming abundant only in the narrow range between the threshold and bimodality scales.

Earlier work has raised the potential importance of nebular emission in the derivation of the physical characteristics of high redshift Lyman break galaxies. Within certain redshift ranges, and especially at z ~ 6-7, such lines may be strong enough to reduce estimates of the stellar masses and ages of galaxies compared those derived assuming broad-band photometry represents stellar light alone. To test this hypothesis at the highest redshifts where such lines can be probed with ground-based facilities, we examine the near-infrared spectra of a representative sample of 20 3.0 < z < 3.8 Lyman break galaxies using the newly-commissioned MOSFIRE near-infrared spectrograph at the Keck I telescope. We use this data to derive the rest-frame equivalent widths (EW) of [O III] emission and show that these are comparable to estimates derived using the SED fitting technique introduced for sources of known redshift by Stark et al (2013). Although our current sample is modest, its [O III] EW distribution is consistent with that inferred for H\alpha\ based on SED fitting of Stark et al's larger sample of 3.8 < z < 5 galaxies. For a subset of survey galaxies, we use the combination of optical and near-infrared spectroscopy to quantify kinematics of outflows in z ~ 3.5 star-forming galaxies, and discuss the implications for reionization measurements. The trends we uncover underline the dangers of relying purely on broad-band photometry to estimate the physical properties of high redshift galaxies and emphasize the important role of diagnostic spectroscopy.

We present deep spectroscopic observations of a Lyman-alpha emitter (LAE) candidate at z ~ 7.7 using the infrared spectrograph LUCI on the 2 x 8.4m Large Binocular Telescope (LBT). The candidate is the brightest among the four z ~ 7.7 LAE candidates found in a narrow-band imaging survey by Krug et al. 2012. Our spectroscopic data include a total of 7.5 hours of integration with LBT/LUCI and are deep enough to significantly (3.2-4.9 sigma) detect the Lyman-alpha emission line of this candidate, based on its Lyman-alpha flux 1.2 x 10^{-17} erg s^{-1} cm^{-2} estimated from the narrow-band photometry. However, we do not find any convincing signal at the expected position of its Lyman-alpha emission line, suggesting that this source is not an LAE at z ~ 7.7. The non-detection in this work, together with the previous studies of z ~ 7.7 LAEs, puts a strong constraint on the bright-end Lyman-alpha luminosity function (LF) at z ~ 7.7. We find a rapid evolution of the Lyman-alpha LF from z ~ 6.5 to 7.7: the upper limit of the z ~ 7.7 LF is more than 5 times lower than the z ~ 6.5 LF at the bright end (f > 1.0 x 10^{-17} erg s^{-1} cm^{-2}, or L > 6.9 x 10^{42} erg s^{-1}). This is likely caused by an increasing neutral fraction in the IGM that substantially attenuates Lyman-alpha emission at z ~ 7.7.

The fraction of ionizing photons that escape from young star-forming galaxies is one of the largest uncertainties in determining the role of galaxies in cosmic reionization. Yet traditional techniques for measuring this fraction are inapplicable at the redshifts of interest due to foreground screening by the Lyman alpha forest. In an earlier study, we demonstrated a reduction in the equivalent width of low-ionization absorption lines in composite spectra of Lyman break galaxies at z=4 compared to similar measures at z=3. This might imply a lower covering fraction of neutral gas and hence an increase with redshift in the escape fraction of ionizing photons. However, our spectral resolution was inadequate to differentiate between several alternative explanations, including changes with redshift in the outflow kinematics. Here we present higher quality spectra of 3 gravitationally lensed Lyman break galaxies at z=4 with a spectral resolution sufficient to break this degeneracy of interpretation. We present a method for deriving the covering fraction of low-ionization gas as a function of outflow velocity and compare the results with similar quality data taken for galaxies at lower redshift. We find a significant trend of lower covering fractions of low-ionization gas for galaxies with strong \Lya emission. In combination with the demographic trends of \Lya emission with redshift from our earlier work, our results provide new evidence for a reduction in the average H I covering fraction, and hence an increase in the escape fraction of ionizing radiation from Lyman break galaxies, with redshift.

We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope.

To assess how external factors such as local interactions and fresh gas accretion influence the global ISM of galaxies, we analyze the relationship between recent enhancements of central star formation and total molecular-to-atomic (H2/HI) gas ratios, using a broad sample of field galaxies spanning early-to-late type morphologies, stellar masses of 10^(7.2-11.2) Msun, and diverse stages of evolution. We find that galaxies occupy several loci in a "fueling diagram" that plots H2/HI vs. mass-corrected blue-centeredness, a metric tracing the degree to which galaxies have bluer centers than the average galaxy at their stellar mass. Spiral galaxies show a positive correlation between H2/HI and mass-corrected blue-centeredness. When combined with previous results linking mass-corrected blue-centeredness to external perturbations, this correlation suggests a link between local galaxy interactions and molecular gas inflow/replenishment. Intriguingly, E/S0 galaxies show a more complex picture: some follow the same correlation, some are quenched, and a distinct population of blue-sequence E/S0 galaxies (with masses below key transitions in gas richness) defines a separate loop in the fueling diagram. This population appears to be composed of low-mass merger remnants currently in late- or post-starburst states, in which the burst first consumes the H2 while the galaxy center keeps getting bluer, then exhausts the H2, at which point the burst population reddens as it ages. Multiple lines of evidence suggest connected evolutionary sequences in the fueling diagram. In particular, tracking total gas-to-stellar mass ratios within the diagram provides evidence of fresh gas accretion onto low-mass E/S0s emerging from central starbursts. Drawing on a comprehensive literature search, we suggest that virtually all galaxies follow the same evolutionary patterns found in our broad sample.

We present a new determination of the UV galaxy luminosity function (LF) at redshift z ~ 7 and z ~ 8, and a first estimate at z ~ 9. An accurate determination of the form and evolution of the LF at high z is crucial for improving our knowledge of early galaxy evolution and cosmic reionization. Our analysis exploits fully the new, deepest WFC3/IR imaging from our HST UDF12 campaign, and includes a new, consistent analysis of all appropriate, shallower/wider-area HST data. Our new measurement of the evolving LF at z ~ 7-8 is based on a final catalogue of ~600 galaxies, and involves a step-wise maximum likelihood determination based on the redshift probability distribution for each object; this makes full use of the 11-band imaging now available in the HUDF, including the new UDF12 F140W data, and the deep Spitzer IRAC imaging. The final result is a determination of the z ~ 7 LF extending down to M_UV = -16.75, and the z ~ 8 LF down to M_UV = -17.00. Fitting a Schechter function, we find M* = -19.90 (+0.23/-0.28), log phi* = -2.96 (+0.18/-0.23), and a faint-end slope alpha=-1.90 (+0.14/-0.15) at z~7, and M* = -20.12 (+0.37/-0.48), log phi* = -3.35 (+0.28/-0.47), alpha=-2.02 (+0.22/-0.23) at z~8. These results strengthen suggestions that the evolution at z > 7 is more akin to `density evolution' than the apparent `luminosity evolution' seen at z ~ 5-7. We also provide the first meaningful information on the LF at z ~ 9, explore alternative extrapolations to higher z, and consider the implications for the evolution of UV luminosity density. Finally, we provide catalogues (including z_phot, M_UV and all photometry) for the 100 most robust z~6.5-11.9 galaxies in the HUDF used in this analysis. We discuss our results in the context of earlier work and the results of an independent analysis of the UDF12 data based on colour-colour selection (Schenker et al. 2013).

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

Bright gravitationally lensed galaxies provide our most detailed view of galaxies at high redshift. Yet as a result of the small number of ultra-bright z~2 lensed systems with confirmed redshifts, most detailed spectroscopic studies have been limited in their scope. With the goal of increasing the number of bright lensed galaxies available for detailed follow-up, we have undertaken a spectroscopic campaign targeting wide separation (>3 arcsec) galaxy-galaxy lens candidates within the Sloan Digital Sky Survey (SDSS). Building on the earlier efforts of our CASSOWARY survey, we target a large sample of candidate galaxy-galaxy lens systems in SDSS using a well-established search algorithm which identifies blue arc-like structures situated around luminous red galaxies. In this paper, we present a new redshift catalog containing 25 lensed sources in SDSS confirmed through spectroscopic follow-up of candidate galaxy-galaxy lens systems. Included in this new sample are two of the brightest galaxies (r=19.6 and 19.7) galaxies known at z~2, a low metallicity (12 + log (O/H)~8.0) extreme nebular line emitting galaxy at z=1.43, and numerous systems for which detailed follow-up will be possible. The source redshifts span 0.9<z<2.5 (median redshift of 1.9), and their optical magnitudes are in the range 19.6<r<22.3. We present a brief source-by-source discussion of the spectroscopic properties extracted from our confirmatory spectra and discuss some initial science results. With more than 50 gravitationally lensed z>1 galaxies now confirmed within SDSS, it will soon be possible for the first time to develop generalized conclusions from detailed spectroscopic studies of the brightest lensed systems at high redshift.

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of \alpha = -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, \alpha = -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 \HST\ program aimed at extending previous WFC3/IR observations of the UDF by quadrupling the exposure time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized the universe; our observations are designed to provide a robust determination of the star formation density at $z$$\,\gtrsim\,$8, improve measurements of the ultraviolet continuum slope at $z$$\,\sim\,7\,-\,$8, facilitate the construction of new samples of $z$$\,\sim\,9\,-\,$10 candidates, and enable the detection of sources up to $z$$\,\sim\,$12. For this project we committed to combining these and other WFC3/IR imaging observations of the UDF area into a single homogeneous dataset, to provide the deepest near-infrared observations of the sky currently achievable. In this paper we present the observational overview of the project, motivated by its scientific goals, and describe the procedures used in reducing the data as well as the final products that are produced. We have used the most up up-to-date methods for calibrating and combining the images, in particular paying attention to correcting several instrumental effects. We release the full combined mosaics, comprising a single, unified set of mosaics of the UDF, providing the deepest near-infrared blank-field view of the universe obtained to date, reaching magnitudes as deep as AB$\,\sim\,$30 in the near-infrared, and yielding a legacy dataset on this field of lasting scientific value to the community.

We present the results of the deepest search to date for star-forming galaxies beyond a redshift z~8.5 utilizing a new sequence of near-infrared Wide Field Camera 3 images of the Hubble Ultra Deep Field. This `UDF12' campaign completed in September 2012 doubles the earlier exposures with WFC3/IR in this field and quadruples the exposure in the key F105W filter used to locate such distant galaxies. Combined with additional imaging in the F140W filter, the fidelity of high redshift candidates is greatly improved. Using spectral energy distribution fitting techniques on objects selected from a deep multi-band near-infrared stack we find 7 promising z>8.5 candidates. As none of the previously claimed UDF candidates with 8.5<z<10 is confirmed by our deeper multi-band imaging, our campaign has transformed the measured abundance of galaxies in this redshift range. Although we recover the candidate UDFj-39546284 (previously proposed at z=10.3), it is undetected in the newly added F140W image, implying it lies at z=11.9 or is an intense emission line galaxy at z~2.4. Although no physically-plausible model can explain the required line intensity given the lack of Lyman alpha or broad-band UV signal, without an infrared spectrum we cannot rule out an exotic interloper. Regardless, our robust z ~ 8.5 - 10 sample demonstrates a luminosity density that continues the smooth decline observed over 6 < z < 8. Such continuity has important implications for models of cosmic reionization and future searches for z>10 galaxies with JWST.

The physical properties inferred from the SEDs of z>3 galaxies have been influential in shaping our understanding of early galaxy formation and the role galaxies may play in cosmic reionization. Of particular importance is the stellar mass density at early times which represents the integral of earlier star formation. An important puzzle arising from the measurements so far reported is that the specific star formation rates (sSFR) evolve far less rapidly than expected in most theoretical models. Yet the observations underpinning these results remain very uncertain, owing in part to the possible contamination of rest-optical broadband light from strong nebular emission lines. To quantify the contribution of nebular emission to broad-band fluxes, we investigate the SEDs of 92 spectroscopically-confirmed galaxies in the redshift range 3.8<z<5.0 chosen because the H-alpha line lies within the Spitzer/IRAC 3.6 um filter. We demonstrate that the 3.6 um flux is systematically in excess of that expected from stellar continuum, which we derive by fitting the SED with population synthesis models. No such excess is seen in a control sample at 3.1<z<3.6 in which there is no nebular contamination in the IRAC filters. From the distribution of our 3.6 um flux excesses, we derive an H-alpha equivalent width (EW) distribution. The mean rest-frame H-alpha EW we infer at 3.8<z<5.0 (270 A) indicates that nebular emission contributes at least 30% of the 3.6 um flux. Via our empirically-derived EW distribution we correct the available stellar mass densities and show that the sSFR evolves more rapidly at z>4 than previously thought, supporting up to a 5x increase between z~2 and 7. Such a trend is much closer to theoretical expectations. Given our findings, we discuss the prospects for verifying quantitatively the nebular emission line strengths prior to the launch of the James Webb Space Telescope.

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

It has recently been reported that the faintest galaxies at z~6-7 display extremely blue UV continuum slopes, with a UV power-law index beta ~ -3. Such slopes are bluer than previously reported for any other galaxy population, and imply extinction-free, young, very low-metallicity stellar populations with a high ionizing photon escape fraction. Here we undertake a critical study of the evidence for such extreme values of beta, combining three new WFC3/IR-selected samples of galaxies spanning ~2 decades in UV luminosity over the redshift range z~4.5-8. We explore the impact of inclusion/exclusion of less robust high-z candidates, and use the varying depths of the samples to explore the effects of noise and selection bias. Simple data-consistency arguments suggest that artificially blue average values of beta can result when the analysis is extended into the deepest ~ 0.5-mag bin of these WFC3/IR-selected samples, regardless of the actual luminosity or z range probed. By confining attention to robust, well-detected high-z galaxy candidates, we find that the average value of beta is consistent with -2.05 +/- 0.10 for z=5-7, and -22 < M_UV < -18. We create and analyse a set of simulations which demonstrate that a bias towards artifically low/blue average values of beta is indeed expected when the UV slope analysis is extended towards the source detection threshold, and conclude that there is as yet no clear evidence for UV slopes significantly bluer than beta ~ -2, the typical value displayed by the bluest star-forming galaxies at more modest z. A robust measurement of beta for the faintest galaxies at z~7-8 remains a key observational goal, as it provides a fundamental test for high escape fractions from a potentially abundant source of reionizing photons. This goal is achievable with HST, but requires still deeper WFC3/IR imaging in the HUDF.

We present spectroscopic confirmation of ten highly luminous (L >= 2L*) Lyman alpha emitters in the redshift range 6.01<z<6.49 (nine galaxies and one AGN), initially drawn from a sample of fourteen z_phot >= 6 Lyman break galaxies (LBGs) selected from an area of 0.25 square degrees within the UKIDSS Ultra-deep Survey (UDS). Overall, our high rate of spectroscopic confirmation (>= 71%) and low rate of contamination provides a strong vindication of the photometric redshift analysis used to define the original sample. By considering star-formation rate estimates based on the Ly_alpha and UV continuum luminosity we conclude that our sample is consistent with a Ly_alpha escape fraction of ~25%. Moreover, after careful consideration of the potential uncertainties and biases, we find that 40%-50% of our sample of L >= 2L* galaxies at 6.0<z<6.5 display strong Ly_alpha emission (rest-frame equivalent width >= 25 Angs), a fraction which is a factor of ~2 higher than previously reported for L <= L* galaxies at z~6. Our results suggest that, as the epoch of reionization is approached, it is plausible that the Ly_alpha emitter fraction amongst luminous (L >=2 L*) LBGs shows a similarly sharp increase to that observed in their lower-luminosity (L <= L*) counterparts.

We analyse the spectral energy distribution (SED) of U, B, V and i-dropout samples from GOODS-MUSIC and we determine their physical properties, such as stellar age and mass, dust attenuation and star formation rate (SFR). Furthermore, we examine how the strength of Ly{\alpha} emission can be constrained from broad-band SED fits instead of relying in spectroscopy. We use our SED fitting tool including the effects of nebular emission and we explore different star formation histories (SFHs). We find that SEDs are statistically better fitted with nebular emission and exponentially decreasing star formation. Considering this result, stellar mass and star formation rate (SFR) estimations modify the specific SFR (SFR/Mstar) - redshift relation, in comparison to previous studies. Finally, our inferred Ly{\alpha} properties are in good agreement with the available spectroscopic observations.

We present and discuss the mean rest-frame ultraviolet spectrum for a sample of 81 Lyman Break Galaxies (LBGs) selected to be B-band dropouts with a mean redshift of z=3.9 and apparent magnitudes z_AB<26. Most of the individual spectra are drawn from our ongoing survey in the GOODS fields with the Keck DEIMOS spectrograph, and we have augmented our sample with published data taken with FORS2 on the VLT. In general we find similar trends in the spectral diagnostics to those found in the earlier, more extensive survey of LBGs at z=3 undertaken by Shapley et al (2003). Specifically, we find low-ionization absorption lines which trace the presence of neutral outflowing gas are weaker in galaxies with stronger Lyman-alpha emission, bluer UV spectral slopes, lower stellar masses, lower UV luminosities and smaller half-light radii. This is consistent with a physical picture whereby star formation drives outflows of neutral gas which scatters Lyman-alpha and gives rise to strong low-ionization absorption lines, while increasing the stellar mass, size, metallicity, and dust content of galaxies. Typical galaxies are thus expected to have stronger Lyman-alpha emission and weaker low-ionization absorption at earlier times (higher redshifts). Indeed, our mean spectrum at z=4 shows somewhat weaker low-ionization absorption lines than at z=3 and available data at higher redshift indicates a rapid decrease in low-ionization absorption strength with redshift. We argue that the reduced low-ionization absorption is likely caused by a decrease in the covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. Our continuing survey will enable us to extend these diagnostics more reliably to higher redshift and determine the implications for the escape fraction of ionizing photons which governs the role of early galaxies in cosmic reionization. [Abridged]

Determining Lya properties of distant galaxies is of great interest for various astrophysical studies. We examine how the strength of Lya emission can be constrained from broad-band SED fits instead of relying on spectroscopy. We use our SED fitting tool including the effects of nebular emission, considering in particular Lya emission as a free parameter, and we demonstrate our method with simulations of mock galaxies. Using this tool we analyse a large sample of U, B, V, and i dropout galaxies with multi-band photometry. We find significant trends of the fraction of galaxies with Lya emission increasing both with redshift z and towards fainter magnitude (at fixed z), and similar trends for the Lya equivalent width. Our inferred Lya properties are in good agreement with the available spectroscopic observations and other data. These results demonstrate that the strength of Lya emission in distant star-forming galaxies can be inferred quantitatively from broad-band SED fits, at least statistically for sufficiently large samples with a good photometric coverage.

Using deep Keck spectroscopy of Lyman break galaxies selected from infrared imaging data taken with WFC3/IR onboard the Hubble Space Telescope, we present new evidence for a reversal in the redshift-dependent fraction of star forming galaxies with detectable Lyman alpha emission in the redshift range 6.3 < z < 8.8. Our earlier surveys with the DEIMOS spectrograph demonstrated a significant increase with redshift in the fraction of line emitting galaxies over the interval 4 < z < 6, particularly for intrinsically faint systems which dominate the luminosity density. Using the longer wavelength sensitivities of LRIS and NIRSPEC, we have targeted 19 Lyman break galaxies selected using recent WFC3/IR data whose photometric redshifts are in the range 6.3 < z < 8.8 and which span a wide range of intrinsic luminosities. Our spectroscopic exposures typically reach a 5-sigma sensitivity of < 50 A for the rest-frame equivalent width (EW) of Lyman alpha emission. Despite the high fraction of emitters seen only a few hundred million years later, we find only 2 convincing and 1 possible line emitter in our more distant sample. Combining with published data on a further 7 sources obtained using FORS2 on the ESO VLT, and assuming continuity in the trends found at lower redshift, we discuss the significance of this apparent reversal in the redshift-dependent Lyman alpha fraction in the context of our range in continuum luminosity. Assuming all the targeted sources are at their photometric redshift and our assumptions about the Lyman alpha EW distribution are correct, we would expect to find so few emitters in less than 1% of the realizations drawn from our lower redshift samples. Our new results provide further support for the suggestion that, at the redshifts now being probed spectroscopically, we are entering the era where the intergalactic medium is partially neutral.

We present the results of a photometric redshift analysis designed to identify z>6 galaxies from the near-IR HST imaging in three deep fields (HUDF, HUDF09-2 & ERS). By adopting a rigorous set of criteria for rejecting low-z interlopers, and by employing a deconfusion technique to allow the available IRAC imaging to be included in the candidate selection process, we have derived a robust sample of 70 Lyman-break galaxies (LBGs) spanning the redshift range 6.0<z<8.7. Based on our final sample we investigate the distribution of UV spectral slopes (beta), finding a variance-weighted mean value of <beta>=-2.05 +/- 0.09 which, contrary to some previous results, is not significantly bluer than displayed by lower-redshift starburst galaxies. We confirm the correlation between UV luminosity and stellar mass reported elsewhere, but based on fitting galaxy templates featuring a range of star-formation histories, metallicities and reddening we find that, at z>=6, the range in mass-to-light ratio (M*/L_UV) at a given UV luminosity could span a factor of ~50. Focusing on a sub-sample of twenty-one candidates with IRAC detections at 3.6-microns we find that L* LBGs at z~6.5 have a median stellar mass of M* = (2.1 +/- 1.1) x 10^9 Msun and a median specific star-formation rate of 1.9 +/- 0.8 Gyr^-1. Using the same sub-sample we have investigated the influence of nebular continuum and line emission, finding that for the majority of candidates (16 out of 21) the best-fitting stellar-mass estimates are reduced by less than a factor of 2.5. Finally, a detailed comparison of our final sample with the results of previous studies suggests that, at faint magnitudes, several high-redshift galaxy samples in the literature are significantly contaminated by low-redshift interlopers (abridged).

Motivated by recent observational results that focus on high redshift black holes, we explore the effect of scatter and observational biases on the ability to recover the intrinsic properties of the black hole population at high redshift. We find that scatter and selection biases can hide the intrinsic correlations between black holes and their hosts, with 'observable' subsamples of the whole population suggesting, on average, positive evolution even when the underlying population is characterized by no- or negative evolution. We create theoretical mass functions of black holes convolving the mass function of dark matter halos with standard relationships linking black holes with their hosts. Under these assumptions, we find that the local MBH - sigma correlation is unable to fit the z = 6 black hole mass function proposed by Willott et al. (2010), overestimating the number density of all but the most massive black holes. Positive evolution or including scatter in the MBH - sigma correlation makes the discrepancy worse, as it further increases the number density of observable black holes. We notice that if the MBH - sigma correlation at z = 6 is steeper than today, then the mass function becomes shallower. This helps reproducing the mass function of z = 6 black holes proposed by Willott et al. (2010). Alternatively, it is possible that very few halos (of order 1/1000) host an active massive black hole at z = 6, or that most AGN are obscured, hindering their detection in optical surveys. Current measurements of the high redshift black hole mass function might be underestimating the density of low mass black holes if the active fraction or luminosity are a function of host or black hole mass. Finally, we discuss physical scenarios that can possibly lead to a steeper MBH - sigma relation at high redshift.

We report the discovery of a unique $z=6.027$ galaxy, multiply imaged by the cluster Abell 383 and detected in new Hubble Space Telescope ACS and WFC3 imaging, as well as in Warm Spitzer observations. This galaxy was selected as a pair of i-dropouts; its suspected high redshift was confirmed by the measurement of a strong Lyman-alpha line in both images using Keck/DEIMOS. Combining Hubble and Spitzer photometry after correcting for contamination by line emission (estimated to be a small effect), we identify a strong Balmer break of 1.5 magnitudes. Taking into account the magnification factor of 11.4+/-1.9 (2.65+/-0.17 mag) for the brightest image, the unlensed AB magnitude for the source is 27.2+/-0.05 in the H band, corresponding to a 0.4 L* galaxy, and 25.7+/-0.08 at 3.6 um. The UV slope is consistent with beta~2.0, and from the rest-frame UV continuum we measure a current star formation rate of 2.4+/-1.1 Msol/yr. The unlensed half-light radius is measured to be 300 pc, from which we deduce a star-forming surface density of ~10 Msol/yr/kpc2. The Lyman-alpha emission is found to be extended over ~3" along the slit, corresponding to ~5 kpc in the source plane. This can be explained by the presence of a much larger envelope of neutral hydrogen around the star-forming region. Finally, fitting the spectral energy distribution using 7 photometric data points with simple SED models, we derive the following properties: very little reddening, an inferred stellar mass of M*=6e9 Msol, and an inferred age of ~800 Myrs (corresponding to a redshift of formation of ~18). The star-formation rate of this object was likely much stronger in the past than at the time of observation, suggesting that we may be missing a fraction of galaxies at z~6 which have already faded in rest-frame UV wavelengths.

We present and analyse near-infrared spectroscopy for a sample of 28 gravitationally- lensed star-forming galaxies in the redshift range 1.5 < z < 5, observed mostly with the Keck II telescope. With typical magnifications of ~1.5-4 magnitudes, our survey provides a valuable census of star formation rates, gas-phase metallicities and dynamical masses for a representative sample of low luminosity galaxies seen at a formative period in cosmic history. We find less evolution in the mass-metallicity relation compared to earlier work that focused on more luminous systems with z - 2-3, especially in the low mass (- 10^9 Msol) where our sample is - 0.25 dex more metal-rich. We interpret this offset as a result of the lower star formation rates (typically a factor of -10 lower) for a given stellar mass in our sub-luminous systems. Taking this effect into account, we conclude our objects are consistent with a fundamental metallicity relation recently proposed from unlensed observations.

Star formation rate and accummulated stellar mass are two fundamental physical quantities that describe the evolutionary state of a forming galaxy. Two recent attempts to determine the relationship between these quantities, by interpreting a sample of star-forming galaxies at redshift of z~4, have led to opposite conclusions. We use a model galaxy population to investigate possible causes for this discrepancy and conclude that minor errors in the conversion from observables to physical quantities can lead to major misrepresentation when applied without awareness of sample selection. We also investigate, in a general way, the physical origin of the correlation between star formation rate and stellar mass within hierarchical galaxy formation theory.

We consider the relationship between molecular-gas and star-formation surface densities in 19 morphologically defined E/S0s with stellar mass <~ 4x10^10 M_sun, paying particular attention to those found on the blue sequence in color vs. stellar mass parameter space, where spiral galaxies typically reside. While some blue-sequence E/S0s must be young major-merger remnants, many low-mass blue-sequence E/S0s appear much less disturbed, and may be experiencing the milder starbursts associated with inner-disk building as spirals (re)grow. For a sample of eight E/S0s (four blue-, two mid-, and two red-sequence) whose CARMA CO(1-0), Spitzer MIPS 24um, and GALEX FUV emission distributions are spatially resolved on a 750pc scale, we find roughly linear relationships between molecular-gas and star-formation surface densities within all galaxies, with power law indices N = 0.6-1.9 (median 1.2). Adding 11 more blue-sequence E/S0s whose CO(1-0) emission is not as well resolved, we find that most of our E/S0s have global 1-8 kpc aperture-averaged molecular-gas surface densities overlapping the range spanned by the disks and centers of spiral galaxies. While many of our E/S0s fall on the same Schmidt-Kennicutt relation as local spirals, ~80% (predominantly on the blue sequence) are offset towards apparently higher molecular-gas star formation efficiency (i.e., shorter molecular gas depletion time). Possible interpretations of the elevated efficiencies include bursty star formation similar to that in local dwarf galaxies, H2 depletion in advanced starbursts, or simply a failure of the CO(1-0) emission to trace all of the molecular gas.

Star forming galaxies represent a valuable tracer of cosmic history. Recent observational progress with Hubble Space Telescope has led to the discovery and study of the earliest-known galaxies corresponding to a period when the Universe was only ~800 million years old. Intense ultraviolet radiation from these early galaxies probably induced a major event in cosmic history: the reionization of intergalactic hydrogen. New techniques are being developed to understand the properties of these most distant galaxies and determine their influence on the evolution of the universe.

As Lyman-alpha photons are scattered by neutral hydrogen, a change with redshift in the Lyman-alpha equivalent width distribution of distant galaxies offers a promising probe of the degree of ionization in the intergalactic medium and hence when cosmic reionization ended. This simple test is complicated by the fact that Lyman-alpha emission can also be affected by the evolving astrophysical details of the host galaxies. In the first paper in this series, we demonstrated both a luminosity and redshift dependent trend in the fraction of Lyman-alpha emitters seen within color-selected Lyman-break galaxies (LBGs) over the range 3<z<6; lower luminosity galaxies and those at higher redshift show an increased likelihood of strong emission. Here we present the results from much deeper 12.5 hour exposures with the Keck DEIMOS spectrograph focused primarily on LBGs at z~6 which enable us to confirm the redshift dependence of line emission more robustly and to higher redshift than was hitherto possible. We find 54+/-11% of faint z~6 LBGs show strong (W_0>25 A) emission, an increase of 1.6x from a similar sample observed at z~4. With a total sample of 74 z~6 LBGs, we determine the luminosity-dependent Lyman-alpha equivalent width distribution. Assuming continuity in these trends to the new population of z~7 sources located with the Hubble WFC3/IR camera, we predict that unless the neutral fraction rises in the intervening 200 Myr, the success rate for spectroscopic confirmation using Lyman-alpha emission should be high.

The starburst phenomenon can shape the evolution of the host galaxy and the surrounding intergalactic medium. The extent of the evolutionary impact is partly determined by the duration of the starburst, which has a direct correlation with both the amount of stellar feedback and the development of galactic winds, particularly for smaller mass dwarf systems. We measure the duration of starbursts in twenty nearby, ongoing, and "fossil" starbursts in dwarf galaxies based on the recent star formation histories derived from resolved stellar population data obtained with the Hubble Space Telescope. Contrary to the shorter times of 3-10 Myr often cited, the starburst durations we measure range from 450 - 650 Myr in fifteen of the dwarf galaxies and up to 1.3 Gyr in four galaxies; these longer durations are comparable to or longer than the dynamical timescales for each system. The same feedback from massive stars that may quench the flickering SF does not disrupt the overall burst event in our sample of galaxies. While five galaxies present fossil bursts, fifteen galaxies show ongoing bursts and thus the final durations may be longer than we report here for these systems. One galaxy shows a burst that has been ongoing for only 20 Myr; we are likely seeing the beginning of a burst event in this system. Using the duration of the starbursts, we calculate that the bursts deposited 10^(53.9)-10^(57.2) ergs of energy into the interstellar medium through stellar winds and supernovae and produced 3%-26% of the host galaxy's mass.

We use archival HST observations of resolved stellar populations to derive the star formation histories (SFHs) of 18 nearby starburst dwarf galaxies. In this first paper we present the observations, color-magnitude diagrams, and the SFHs of the 18 starburst galaxies, based on a homogeneous approach to the data reduction, differential extinction, and treatment of photometric completeness. We adopt a star formation rate (SFR) threshold normalized to the average SFR of the individual system as a metric for classifying starbursts in SFHs derived from resolved stellar populations. This choice facilitates finding not only currently bursting galaxies but also "fossil" bursts increasing the sample size of starburst galaxies in the nearby (D<8 Mpc) universe. Thirteen of the eighteen galaxies are experiencing ongoing bursts and five galaxies show fossil bursts. From our reconstructed SFHs, it is evident that the elevated SFRs of a burst are sustained for hundreds of Myr with variations on small timescales. A long >100 Myr temporal baseline is thus fundamental to any starburst definition or identification method. The longer lived bursts rule out rapid "self-quenching" of starbursts on global scales. The bursting galaxies' gas consumption timescales are shorter than the Hubble time for all but one galaxy confirming the short-lived nature of starbursts based on fuel limitations. Additionally, we find the strength of the H{\alpha} emission usually correlates with the CMD based SFR during the last 4-10 Myr. However, in four cases, the H{\alpha} emission is significantly less than what is expected for models of starbursts; the discrepancy is due to the SFR changing on timescales of a few Myr. The inherently short timescale of the H{\alpha} emission limits identifying galaxies as starbursts based on the current characteristics which may or may not be representative of the recent SFH of a galaxy.

We study the luminosity gap, dm12, between the first and second ranked galaxies in a sample of 59 massive galaxy clusters, using data from the Hale Telescope, HST, Chandra, and Spitzer. We find that the dm12 distribution, p(dm12), is a declining function of dm12, to which we fitted a straight line: p(dm12) propto -(0.13+/-0.02)dm12. The fraction of clusters with "large" luminosity gaps is p(dm12>=1)=0.37+/-0.08, which represents a 3sigma excess over that obtained from Monte Carlo simulations of a Schechter function that matches the mean cluster galaxy luminosity function. We also identify four clusters with "extreme" luminosity gaps, dm12>=2, giving a fraction of p(dm12>=2)=0.07+0.05-0.03. More generally, large luminosity gap clusters are relatively homogeneous, with elliptical/disky brightest cluster galaxies (BCGs), cuspy gas density profiles (i.e. strong cool cores), high concentrations, and low substructure fractions. In contrast, small luminosity gap clusters are heterogeneous, spanning the full range of boxy/elliptical/disky BCG morphologies, the full range of cool core strengths and dark matter concentrations, and have large substructure fractions. Taken together, these results imply that the amplitude of the luminosity gap is a function of both the formation epoch, and the recent infall history of the cluster. "BCG dominance" is therefore a phase that a cluster may evolve through, and is not an evolutionary "cul-de-sac". We also compare our results with semi-analytic model predictions based on the Millennium Simulation. None of the models are able to reproduce all of the observational results, underlining the inability of current models to match the empirical properties of BCGs. We identify the strength of AGN feedback and the efficiency with which cluster galaxies are replenished after they merge with the BCG in each model as possible causes of these discrepancies. [Abridged]

We present results from a new Keck spectroscopic survey of UV-faint LBGs in the redshift range 3<z<7. Combined with earlier Keck and published ESO VLT data, our sample contains more than 600 dropouts, offering new insight into the nature of sub-L* sources typical of those likely to dominate the cosmic reionisation process. Here we use this sample to characterise the fraction of strong Lya emitters within the continuum-selected dropouts. By quantifying how the "Lya fraction" varies with redshift, we seek to constrain changes in Lya transmission associated with reionisation. In order to distinguish the effects of reionisation from other factors which affect the Lya fraction (e.g. dust, ISM kinematics), we study the luminosity and redshift-dependence of the Lya fraction over 3<z<6, when the IGM is known to be ionised. These results reveal that low luminosity galaxies show strong Lya emission much more frequently than luminous systems, and that at fixed luminosity, the prevalence of strong Lya emission increases moderately with redshift over 3 < z < 6. Based on the correlation between blue UV slopes and strong Lya emitting galaxies in our dataset, we argue that the Lya fraction trends are governed by redshift and luminosity-dependent variations in the dust obscuration, with likely additional contributions from trends in the kinematics and covering fraction of neutral hydrogen. We find a tentative decrease in the Lya fraction at z~7 based on the limited IR spectra for candidate z~7 lensed LBGs, a result which, if confirmed with future surveys, would suggest an increase in the neutral fraction by this epoch. Given the supply of z and Y-drops now available from Hubble WFC3/IR surveys, we show it will soon be possible to significantly improve estimates of the Lya fraction using optical and near-IR spectrographs, thereby extending the study conducted in this paper to 7<z<8.

Massive galaxies in the early Universe have been shown to be forming stars at surprisingly high rates. Prominent examples are dust-obscured galaxies which are luminous when observed at sub-millimeter (sub-mm) wavelengths and which may be forming stars at rates upto 1,000Mo/yr. These intense bursts of star formation are believed to be driven by mergers between gas rich galaxies. However, probing the properties of individual star-forming regions within these galaxies is beyond the spatial resolution and sensitivity of even the largest telescopes at present. Here, we report observations of the sub-mm galaxy SMMJ2135-0102 at redshift z=2.3259 which has been gravitationally magnified by a factor of 32 by a massive foreground galaxy cluster lens. This cosmic magnification, when combined with high-resolution sub-mm imaging, resolves the star-forming regions at a linear scale of just ~100 parsecs. We find that the luminosity densities of these star-forming regions are comparable to the dense cores of giant molecular clouds in the local Universe, but they are ~100x larger and 10^7 times more luminous. Although vigorously star-forming, the underlying physics of the star formation processes at z~2 appears to be similar to that seen in local galaxies even though the energetics are unlike anything found in the present-day Universe.

We have searched for star-forming galaxies at z~7 by applying the Lyman-break technique to newly-released 1.1micron Y-band images from WFC3 on HST. By comparing these images of the Hubble Ultra Deep Field with the ACS z'-band (0.85micron), we identify objects with red colours, (z'-Y)_AB>1.3), consistent with the Ly-alpha forest absorption at z~6.7-8.8. We identify 12 of these z'-drops down to a limiting magnitude Y_AB<28.5 (equivalent to a star formation rate of 1.3M_sun/yr at z=7.1), which are undetected in the other ACS filters. We use the WFC3 J-band image to eliminate contaminant low mass Galactic stars, which typically have redder colours than z~7 galaxies. One of our z'-drops is a probably a T-dwarf star. The z~7 z'-drops have much bluer spectral slopes than Lyman-break galaxies at lower redshift. Our brightest z'-drop is not present in the NICMOS J-band image of the same field taken 5 years before, and is a possible transient object. From the 10 remaining z~7 candidates we determine a lower limit on the star formation rate density of 0.0017M_sun/yr/Mpc^3 for a Salpeter initial mass function, which rises to 0.0025-0.0034M_sun/yr/Mpc^3 after correction for luminosity bias. The star formation rate density is a factor of ~10 less than that at z=3-4, and is about half the value at z~6. While based on a single deep field, our results suggest that this star formation rate density would produce insufficient Lyman continuum photons to reionize the Universe unless the escape fraction of these photons is extremely high (f_esc>0.5), and the clumping factor of the Universe is low. Even then, we need to invoke a large contribution from galaxies below our detection limit. The apparent shortfall in ionizing photons might be alleviated if stellar populations at high redshift are low metallicity or have a top-heavy IMF.

We analyse recently acquired near-infrared Hubble Space Telescope imaging of the GOODS-South field to search for star forming galaxies at z~7.0. By comparing WFC 3 0.98 micron Y-band images with ACS z-band (0.85 micron) images, we identify objects with colours consistent with Lyman break galaxies at z~6.4-7.4. This new data covers an area five times larger than that previously reported in the WFC3 imaging of the Hubble Ultra Deep Field, and affords a valuable constraint on the bright end of the luminosity function. Using additional imaging of the region in the ACS B, V and i-bands from GOODS v2.0 and the WFC3 J-band we attempt to remove any low-redshift interlopers. Our selection criteria yields 6 candidates brighter than Y_AB = 27.0, of which all except one are detected in the ACS z-band imaging and are thus unlikely to be transients. Assuming all 6 candidates are at z~7 this implies a surface density of objects brighter than Y_AB = 27.0 of 0.30\pm0.12 arcmin-2, a value significantly smaller than the prediction from the z ~ 6 luminosity function. This suggests continued evolution of the bright end of the luminosity function between z = 6 to 7, with number densities lower at higher redshift.

We present the results of a systematic search for galaxies in the redshift range z = 6 - 9, within the new, deep, near-infrared imaging of the Hubble Ultra Deep Field provided by the Wide Field Camera 3 (WFC3) on HST. We have performed full SED fitting to the optical+infrared photometry of all high-redshift galaxy candidates detected at greater than 5-sigma in at least one of the WFC3/IR broad-band filters. After rejection of contaminants, the result is a sample of 49 galaxies with primary redshift solutions z > 5.9. Our sample, selected without recourse to specific colour cuts, re-selects all but the faintest one of the 16 z-drops selected by Oesch et al. (2009), recovers all 5 of the Y-drops reported by Bouwens et al. (2009), and adds a further 29 galaxy candidates, of which 12 lie beyond z = 6.3, and 4 lie beyond z = 7. We also present confidence intervals on our photometric redshift estimates, and caution that acceptable low-redshift (z < 2) solutions exist for 28 out of the 37 galaxies at z > 6.3, and for all 8 galaxy candidates at z > 7.5. Nevertheless, the very highest redshift candidates appear to be strongly clustered in the field. We derive new estimates of the ultraviolet galaxy luminosity function at z = 7 and z = 8. Where our results are most robust, at a characteristic luminosity M(1500) ~ -19.5 (AB), we find that the comoving number density of galaxies declines by a factor of ~ 2.5 between z = 6 and z = 7, and by a further factor of ~ 2 by z = 8. These results suggest that it is difficult for the observed population of high-redshift star-forming galaxies to achieve reionisation by z ~ 6 without a significant contribution from galaxies well below the detection limits, plus alterations in the escape fraction of ionising photons and/or continued vigorous star formation at z > 15.

We report the results of a study of the rest-frame UV spectrum of the Cosmic Eye, a luminous Lyman break galaxy at z=3.07331 gravitationally lensed by a factor of 25. The spectrum, recorded with the ESI spectrograph on the Keck II telescope, is rich in absorption features from the gas and massive stars in this galaxy. The interstellar absorption lines are resolved into two components of approximately equal strength and each spanning several hundred km/s in velocity. One component has a net blueshift of -70 km/s relative to the stars and H II regions and presumably arises in a galaxy-scale outflow similar to those seen in most star-forming galaxies at z = 2-3. The other is more unusual in showing a mean redshift of +350 km/s relative to the systemic redshift; possible interpretations include a merging clump, or material ejected by a previous star formation episode and now falling back onto the galaxy, or more simply a chance alignment with a foreground galaxy. In the metal absorption lines, both components only partially cover the OB stars against which they are being viewed. We tentatively associate the redshifted component with the strong damped Lyman alpha line, indicative of a column density N(H I) = (3.0 +/- 0.8) x 10(21) atoms/cm2, and propose that it provides the dust `foreground screen' responsible for the low ratio of far-infrared to UV luminosities of the Cosmic Eye. Compared to other well-studied examples of strongly lensed galaxies, we find that the young stellar population of the Cosmic Eye is essentially indistinguishable from those of the Cosmic Horseshoe and MS 1512-cB58, while the interstellar spectra of all three galaxies are markedly different, attesting to the real complexity of the interplay between starbursts and ambient interstellar matter in young galaxies (abridged).

We present spatially-resolved dynamics for six strongly lensed star-forming galaxies at z=1.7-3.1, each enlarged by a linear magnification factor ~8. Using the Keck laser guide star AO system and the OSIRIS integral field unit spectrograph we resolve kinematic and morphological detail in our sample with an unprecedented fidelity, in some cases achieving spatial resolutions of ~100 pc. With one exception our sources have diameters ranging from 1-7 kpc, star formation rates of 2-40 Msun/yr (uncorrected for extinction) and dynamical masses of 10^(9.7-10.3) Msun. With this exquisite resolution we find that four of the six galaxies display coherent velocity fields consistent with a simple rotating disk model, which can only be recovered with the considerably improved spatial resolution and sampling from the combination of adaptive optics and strong gravitational lensing. Our model fits imply ratios for the systemic to random motion, V sin(i)/sigma, ranging from 0.5-1.3 and Toomre disk parameters Q<1. The large fraction of well-ordered velocity fields in our sample is consistent with data analyzed for larger, more luminous sources at this redshift. Our high resolution data further reveal that all six galaxies contain multiple giant star-forming HII regions whose resolved diameters are in the range 300 pc - 1.0 kpc, consistent with the Jeans length expected in the case of dispersion support. The density of star formation in these regions is ~100 times higher than observed in local spirals; such high values are only seen in the most luminous local starbursts. The global dynamics and demographics of star formation in these HII regions suggest that vigorous star formation is primarily governed by gravitational instability in primitive rotating disks.

We calibrate the Baryonic Tully-Fisher (BTF) Relation using a sample of gas dominated galaxies. These determine the absolute scale of the baryonic mass--rotation speed relation independent of the choice of stellar mass estimator. We find a BTF slope of 3.94 +/- 0.07 (random) +/- 0.08 (systematic) and a zero point of 1.79 +/- 0.26 (random) +/- 0.25 (systematic). We apply this relation to estimate the stellar masses of star dominated galaxies. This procedure reproduces the trend of mass-to-light ratio with color predicted by population synthesis models. The normalization is also correct, consistent with empirical estimates of the IMF used in such models.

The quest to discover the most distant galaxies has developed rapidly in the last decade. We are now exploring redshifts of 6 and beyond, when the Universe was less than a billion years old, an epoch when the previously-neutral intergalactic medium was reionized. The continuing discovery of galaxies at progressively higher and higher redshifts has been driven by the availability of large telescopes on the ground and in space, improvements in detector technology, and new search strategies. Over the past 4 years, the Lyman break technique has been shown to be effective in isolating z~6 star-forming i'-drop galaxies through spectroscopic confirmation with large ground-based telescopes (Keck, Gemini and the ESO VLTs). Narrow-band imaging, notably with the wide field of the Subaru telescope, has also produced samples of Lyman-alpha emitters at these redshifts. A The discovery of this i'-drop galaxy population has been used to infer the global star formation rate density at this epoch (z~6), and we are now beginning to constrain the contribution to reionization of the UV flux from these galaxies. Infrared data from the Spitzer Space Telescope has been used to determine the spectral energy distributions (SEDs) from the rest-frame UV to the optical, and constrain the previous star formation histories, masses and ages. The indications are that much of the stellar mass of these galaxies might have formed in vigorous bursts at z>6. The next big advances would be to test the population synthesis modelling of these z~6 galaxies through spectroscopy of the rest-frame optical (rather than crude broad-band SEDs), and also to push the observational horizon for galaxies further to directly explore star formation during the reionization epoch. JWST is likely to have a profound impact on realising these goals.

We present a detailed study of the spatially resolved kinematics, star-formation and stellar mass in a highly amplified galaxy at z=4.92 behind the lensing cluster MS1358+62. We use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the velocity structure of the galaxy from the nebular [OII] emission. Using a detailed gravitational lens model, we account for lensing amplification factor 12.+/-2.0 and find that this intrinsically L* galaxy has a stellar mass of M*=7+/-2x10^8Mo, a dynamical mass of Mdyn=3+/-1x10^9csc^2(i)Mo (within of 2kpc) and a star-formation rate of 42+/-8Mo/yr. The source-plane UV/optical morphology of this galaxy is dominated by five discrete star-forming regions. Exploiting the dynamical information we derive masses for individual star-forming regions of Mcl~10^(8-9)Mo with sizes of ~200pc. We find that, at a fixed size, the star-formation rate density within these HII regions is approximately two orders of magnitude greater than those observed in local spiral/starburst galaxies, but consistent with the most massive HII regions in the local Universe such as 30Dor. Finally, we compare the spatially resolved nebular emission line velocity with the Ly-alpha and UV-ISM lines and find that this galaxy is surrounded by a galactic scale outflow in which the Ly-alpha appears redshifted by ~150km/s and the UV-ISM lines blue-shifted by -200km/s from the (systemic) nebular emission. The velocity structure of the outflow mirrors that of the nebular emission suggesting the outflow is young (~15Myr), and has yet to burst out of the system. Taken together, these results suggest that this young galaxy is undergoing its first major epoch of mass assembly.

We report the results of a Spitzer infrared study of the Cosmic Eye, a strongly lensed, L*_UV Lyman Break Galaxy (LBG) at z=3.074. We obtained Spitzer IRS spectroscopy as well as MIPS 24 and 70 micron photometry. The Eye is detected with high significance at both 24 and 70 microns and, when including a flux limit at 3.5 mm, we estimate an infrared luminosity of L_IR = 8.3 (+4.7-4.4) x10^11 L_sun assuming a magnification of 28+-3. This L_IR is eight times lower than that predicted from the rest-frame UV properties assuming a Calzetti reddening law. This has also been observed in other young LBGs, and indicates that the dust reddening law may be steeper in these galaxies. The mid-IR spectrum shows strong PAH emission at 6.2 and 7.7 microns, with equivalent widths near the maximum values observed in star-forming galaxies at any redshift. The L_PAH-to-L_IR ratio lies close to the relation measured in local starbursts. Therefore, L_PAH or L_MIR may be used to estimate L_IR and thus, star formation rate, of LBGs, whose fluxes at longer wavelengths are typically below current confusion limits. We also report the highest redshift detection of the 3.3 micron PAH emission feature. The PAH ratio, L_6.2/L_3.3=5.1+- 2.7, and the PAH-to-L_IR ratio, L_3.3/L_IR = 8.5 +- 4.7 x10^-4, are both in agreement with measurements in local starbursts and ULIRGs, suggesting that this line may serve as a good proxy for L_PAH or L_IR at z > 3 with the James Webb Space Telescope.

We present new measurements of the evolution in the Lyman break galaxy (LBG) population between z~4 and z~6. By utilizing the extensive multiwavelength datasets available in the GOODS fields, we identify 2443 B, 506 V, and 137 i'-band dropout galaxies likely to be at z~4, 5, and 6. With the goal of understanding the duration of typical star formation episodes in galaxies at z>4, we examine the distribution of stellar masses and ages as a function of cosmic time. We find that at a fixed rest-UV luminosity, the average stellar masses and ages of galaxies do not increase significantly between z~6 and 4. In order to maintain this near equilibrium in the average properties of high redshift LBGs, we argue that there must be a steady flux of young, newly-luminous objects at each successive redshift. When considered along with the short duty cycles inferred from clustering measurements, these results may suggest that galaxies are undergoing star formation episodes lasting only several hundred million years. In contrast to the unchanging relationship between the average stellar mass and rest-UV luminosity, we find that the number density of massive galaxies increases considerably with time over 4<z<6. Given this rapid increase of UV luminous massive galaxies, we explore the possibility that a significant fraction of massive (M*>1e11 Msun) z~2-3 distant red galaxies (DRGs) were in part assembled in an LBG phase at earlier times. Integrating the growth in the stellar mass function of actively forming LBGs over 4<z<6 down to z~2, we find that z>3 LBGs could have contributed significantly to the quiescent DRG population, indicating that the intense star-forming systems probed by current sub-millimeter observations are not the only route toward the assembly of DRGs at z~2.

The duration of a starburst is a fundamental parameter affecting the evolution of galaxies yet, to date, observational constraints on the durations of starbursts are not well established. Here we study the recent star formation histories (SFHs) of three nearby dwarf galaxies to rigorously quantify the duration of their starburst events using a uniform and consistent approach. We find that the bursts range from ~200 - ~400 Myr in duration resolving the tension between the shorter timescales often derived observationally with the longer timescales derived from dynamical arguments. If these three starbursts are typical of starbursts in dwarf galaxies, then the short timescales (3 - 10 Myr) associated with starbursts in previous studies are best understood as "flickering" events which are simply small components of the larger starburst. In this sample of three nearby dwarfs, the bursts are not localized events. All three systems show bursting levels of star formation in regions of both high and low stellar density. The enhanced star formation moves around the galaxy during the bursts and covers a large fraction of the area of the galaxy. These massive, long duration bursts can significantly affect the structure, dynamics, and chemical evolution of the host galaxy and can be the progenitors of "superwinds" that drive much of the recently chemically enriched material from the galaxy into the intergalactic medium.

Recent studies of galaxies ~2-3 Gyr after the Big Bang have revealed large, turbulent rotating systems. The existence of well-ordered rotation in galaxies during this peak epoch of cosmic star formation may suggest that gas accretion through cold streams is likely to be the dominant mode by which most star-forming galaxies at high redshift since major mergers can completely disrupt the observed velocity fields. However poor spatial resolution and sensitivity have hampered this interpretation, limiting the study to the largest and most luminous galaxies, which may have fundamentally different modes of assembly than more typical star forming galaxies. Here we report observations of a typical star forming galaxy at z=3.07 with a linear resolution of ~100 parsec. This spatial sampling is made possible by the combination of gravitational lensing and laser guide star adaptive optics. We find a well-ordered compact source in which molecular gas is being converted efficiently into stars, likely assembling a spheroidal bulge and disk similar to those seen in spiral galaxies at the present day.

We present the results of a systematic search for gravitationally-lensed continuum Lyman break `drop-outs' beyond a redshift 7 conducted via very deep imaging through six foreground clusters undertaken with the Hubble and Spitzer Space Telescopes. The survey has yielded 10 z-band and 2 J-band drop-out candidates to photometric limits of J_110~=26.2 AB (5sigma). Taking into account the magnifications afforded by our clusters (1-4 magnitudes), we probe the presence of z>7 sources to unlensed limits of J_{110}~=30 AB, fainter than those charted in the Hubble Ultradeep Field. To verify the fidelity of our candidates we conduct a number of tests for instrumental effects which would lead to spurious detections, and carefully evaluate the likelihood of foreground contamination by considering photometric uncertainties in the drop-out signature, the upper limits from stacked IRAC data and the statistics of multiply-imaged sources. Overall, we conclude that we can expect about half of our sample of z-band drop-outs are likely to be at high redshift. An ambitious infrared spectroscopic campaign undertaken with the NIRSPEC spectrograph at the WM Keck Observatory for seven of the most promising candidates failed to detect any Lyman-alpha emission highlighting the challenge of making further progress in this field. While the volume density of high redshift sources will likely remain uncertain until more powerful facilities are available, our data provides the first potentially interesting constraints on the UV luminosity function at z~=7.5 at intrinsically faint limits. We discuss the implications of our results in the context of the hypothesis that the bulk of the reionizing photons in the era 7<z<12 arise in low luminosity galaxies undetected by conventional surveys.

In this article, we apply strong lensing techniques in Abell 1703, a massive X-ray luminous galaxy cluster at z=0.28. Our analysis is based on imaging data both from space and ground in 8 bands, complemented with a spectroscopic survey. Abell 1703 looks rather circular from the general shape of its multiply imaged systems and present a dominant giant elliptical cD galaxy in its centre. This cluster exhibits a remarkable bright 'central ring' formed by 4 bright images at z_{spec}=0.888 located very close to the cD galaxy, providing observational constraints that are potentially very interesting to probe the central mass distribution. The stellar contribution from the cD galaxy (~1.25 10^{12} M_{sun} within 7") is accounted for in our parametric mass modelling, and the underlying smooth dark matter component distribution is described using a generalized NFW profile parametrized with a central logarithmic slope \alpha. We find that within the range where observational constraints are present (from ~5" to ~50"), the slope of the dark matter distribution in Abell1703 is equal to 1.09^{+0.05}_{-0.11} (3\sigma confidence level). The concentration parameter is equal to c_{200} ~ 3.5, and the scale radius is constrained to be larger than the region where observational constraints are available. Within this radius, the 2D mass is equal to M(50")=2.4 10^{14} M_{\sun}. We cannot draw any conclusions on cosmological models at this point since we lack results from realistic numerical simulations containing baryons to make a proper comparison. We advocate the need for a sample of observed and simulated unimodal relaxed galaxy clusters in order to make reliable comparisons, and potentially provide a test of cosmological models.

We report imaging polarimetry of segments B and C of the Jupiter-family Comet 73P/Schwassmann-Wachmann 3 in the I and H bandpasses at solar phase angles of approximately 35 and 85deg. The level of polarization was typical for active comets, but larger than expected for a Jupiter-family comet. The polarimetric color was slightly red (dP/dL = +1.2 +/- 0.4) at a phase angle of ~ 35deg and either neutral or slightly blue at a phase angle of ~ 85deg. Observations during the closest approach from 2006 May 11-13 achieved a resolution of 35 km at the nucleus. Both segments clearly depart from a 1/rho surface brightness for the first 50 - 200 km from the nucleus. Simulations of radiation driven dust dynamics can reproduce some of the observed coma morphology, but only with a wide distribution of initial dust velocities (at least a factor of 10) for a given grain radius. Grain aggregate breakup and fragmentation are able to reproduce the observed profile perpendicular to the Sun-Comet axis, but fit the observations less well along this axis (into the tail). The required fragmentation is significant, with a reduction in the mean grain aggregate size by about a factor of 10. A combination of the two processes could possibly explain the surface brightness profile of the comet.

We present the optical and near-infrared photometry and spectroscopy of four faint T dwarfs newly discovered from the UKIDSS first data release. The sample, drawn from an imaged area of ~136 square degrees to a depth of Y=19.9 (5-sigma, Vega), is located in the SDSS Southern Equatorial Stripe, a region of significant future deep imaging potential. We detail the selection and followup of these objects, three of which are spectroscopically confirmed brown dwarfs ranging from type T2.5 to T7.5, and one is photometrically identified as early T. Their magnitudes range from Y=19.01 to 19.88 with derived distances from 34 to 98 pc, making these among the coldest and faintest brown dwarfs known. The sample brings the total number of T dwarfs found or confirmed by UKIDSS data in this region to nine, and we discuss the projected numbers of dwarfs in the future survey data. We estimate that ~240 early- and late-T dwarfs are discoverable in the UKIDSS LAS data, falling significantly short of published model projections and suggesting that IMFs and/or birthrates may be at the low end of possible models. Thus, deeper optical data has good potential to exploit the UKIDSS survey depth more fully, but may still find the potential Y dwarf sample to be extremely rare.

We have observed the central region of the IR-dark cloud filament associated with IRAS 18507+0121 at millimeter wavelengths in CO(1-0), 13CO(1-0), and C18O(1-0) line emission and with the Spitzer Space Telescope at mid-IR wavelengths. Five massive outflows from two cloud cores were discovered. Three outflows are centered on or near an Ultracompact HII region (G34.4+0.23) while the remaining two outflows originate from the millimeter core G34.4+0.23 MM. Modeling of the SEDs of the mid-IR sources identified 31 young stellar objects in the filament with a combined stellar mass of ~127 +/- 27 Msun. An additional 22 sources were identified as probable cluster members based on the presence of strong 24 micron emission. The total star formation efficiency in the G34.4 cloud filament is estimated to be ~7% while the massive and intermediate mass star formation efficiency in the entire cloud filament is estimated to be roughly 2%. A comparison of the gravitational binding energy with the outflow kinetic energy suggests that the compact core containing G34.4+0.23 MM is being destroyed by its molecular outflows whereas the outflows associated with more massive core surrounding the G34.4 UC HII region are not likely to totally disrupt the cloud. Additionally, a qualitative evaluation of the region appears to suggest that stars in this region may have formed in two stages: first lower mass stars formed and then, a few Myrs later, the more massive stars began to form.

We have carried out an extensive spectroscopic survey with the Keck and VLT telescopes, targeting lensed galaxies in the background of the massive cluster Abell 68. Spectroscopic measurements are obtained for 26 lensed images, including a distant galaxy at z=5.4 . Redshifts have been determined for 5 out of 7 multiply-image systems. Through a careful modeling of the mass distribution in the strongly-lensed regime, we derive a mass estimate of 5.3 x 10^14 Msun within 500 kpc. Our mass model is then used to constrain the redshift distribution of the remaining multiply-imaged and singly-imaged sources. This enables us to examine the physical properties for a subsample of 7 Lyman-alpha emitters at 1.7 < z < 5.5, whose unlensed luminosities of ~ 10^41 ergs/s are fainter than similar objects found in blank fields. Of particular interest is an extended Lyman-alpha emission region surrounding a highly magnified source at z=2.6, detected in VIMOS Integral Field Spectroscopy data. The physical scale of the most distant lensed source at z=5.4 is very small (<300 pc), similar to the lensed z ~ 5.6 emitter reported by Ellis et al. (2001) in Abell 2218. New photometric data available for Abell 2218 allow for a direct comparison between these two unique objects. Our survey illustrates the practicality of using lensing clusters to probe the faint end of the z ~ 2-5 Lyman-alpha luminosity function in a manner that is complementary to blank field narrow-band surveys.

[Abridged] We develop a simple star formation model whose goal is to interpret the emerging body of observational data on star-forming galaxies at z>~6. The efficiency and duty cycle of the star formation activity within dark matter halos are determined by fitting the luminosity functions of Lya emitter and Lyman-break galaxies at redshifts z~5-6. Using our model parameters we predict the likely abundance of star forming galaxies at earlier epochs and compare these to the emerging data in the redshift interval 7<z<10. We find that the abundance of luminous Lyman-break galaxies in the 500 Myr between z~6 and 10 can be naturally explained by the hierarchical assembly of dark matter haloes; there is only marginal evidence for strong physical evolution. In contrast, the first estimates of the abundance of less luminous star forming galaxies at z=9-10 are higher than predicted and, if verified by further data, may suggest a top-heavy stellar mass function at these early epochs. Although these abundances remain uncertain because of the difficulty of spectroscopic confirmation and cosmic variance, even a modest improvement in survey capability with present or upcoming facilities should yield great progress. In this context, we use our model to consider those observational techniques that hold the most promise and make predictions for specific surveys that are, or will soon be, underway. We conclude that narrowband Lya emitter surveys should be efficient on searches at z~7-8; however, such conventional surveys are unlikely to detect sufficient galaxies at z~10 to provide useful constraints. For this reason, gravitational lensing offers the best prospect for probing the z~10 universe prior to JWST.

We discuss new observational constraints on the abundance of faint high redshift Lyman-alpha emitters secured from a deep Keck near-infrared spectroscopic survey which utilizes the strong magnification provided by lensing galaxy clusters at intermediate redshift (z=0.2-0.5). In each of 9 clusters we have undertaken a systematic `blind' search for line emission with NIRSPEC in the J-band within carefully-selected regions which offer very high magnifications (>10-50x) for background sources with redshifts z~10. The high magnification enables the detection of emission at unprecedented limits (10^41 - 10^42 erg/s), much fainter than those of conventional narrow band imaging and other spectroscopic searches. As the comoving volumes probed are small, our survey is designed to address the important question of whether intrinsically feeble star forming galaxies could provide the dominant ionizing flux at z~10. Our survey has yielded six promising (>5 sigma) candidate Lyman-alpha emitters which lie between z=8.7 and z=10.2. While the predicted reionization photon budget depends upon a large number of physical assumptions, our first glimpse at the z~10 Universe suggests that low luminosity star-forming galaxies contribute a significant proportion of the UV photons necessary for cosmic reionization. [Abridged]

We report the discovery using Hubble Space Telescope imaging and Keck spectroscopy of a very bright, highly magnified (~30*) Lyman Break Galaxy (LBG) at z=3.07 in the field of the massive z=0.33 cluster MACSJ2135.2-0102. The system comprises two high-surface brightness arcs with a maximum extent of 3", bracketing a central object which we identify as a massive early-type galaxy at z=0.73. We construct a lens model which reproduces the main features of the system using a combination of a galaxy-scale lens and the foreground cluster. We show that the morphological, spectral and photometric properties of the arcs are consistent with them arising from the lensing of a single ~L* LBG. The most important feature of this system is that the lensing magnification results in an apparent magnitude of r=20.3, making this one of the brightest LBGs known. Such a high magnification provides the opportunity of obtaining very high signal to noise (and potentially spatially resolved) spectroscopy of a high redshift galaxy to study its physical properties. We present initial imaging and spectroscopy demonstrating the basic properties of the system and discuss the opportunities for future observations.

We present a new measurement of the integrated stellar mass per comoving volume at redshift 5 determined via spectral energy fitting drawn from a sample of 214 photometrically-selected galaxies with z'<26.5 in the southern GOODS field. Following procedures introduced by Eyles et al. (2005), we estimate stellar masses for various sub-samples for which reliable and unconfused Spitzer IRAC detections are available. A spectroscopic sample of 14 of the most luminous sources with <z>=4.92 provides a firm lower limit to the stellar mass density of 1e6 Msun/Mpc^3. Several galaxies in this sub-sample have masses of order 10^11 Msun implying significant earlier activity occurred in massive systems. We then consider a larger sample whose photometric redshifts in the publicly-available GOODS-MUSIC catalog lie in the range 4.4 <z 5.6. Before adopting the GOODS-MUSIC photometric redshifts, we check the accuracy of their photometry and explore the possibility of contamination by low-z galaxies and low-mass stars. After excising probable stellar contaminants and using the z'-J color to exclude any remaining foreground red galaxies, we conclude that 196 sources are likely to be at z~5. The implied mass density from the unconfused IRAC fraction of this sample, scaled to the total available, is 6e6 Msun/Mpc^3. We discuss the uncertainties as well as the likelihood that we have underestimated the true mass density. Including fainter and quiescent sources the total integrated density could be as high as 1e7 Msun/Mpc^3. Using the currently available (but highly uncertain) rate of decline in the star formationhistory over 5 <z< 10, a better fit is obtained for the assembled mass at z~5 if we admit significant dust extinction at early times or extend the luminosity function to very faint limits. [abridged]

We measure the ages, stellar masses, and star formation histories of z~6 galaxies, observed within 1Gyr of the Big Bang. We use imaging from HST and Spitzer from the public "Great Observatories Origins Deep Survey", coupled with ground-based near-infrared imaging, to measure the spectral energy distributions from 0.8-5microns, spanning the rest-frame UV and optical. From our sample of ~50 i-drop Lyman-break star-forming galaxies in GOODS-South with z'(AB)<27mag, we focus on ~30 with reliable photometric or spectroscopic redshifts. Half of these are confused with foreground sources at Spitzer resolution, but from the 16 with clean photometry we find that a surprisingly large fraction (40%) have evidence for substantial Balmer-breaks. This indicates the presence of old underlying stellar populations that dominate the stellar masses. For these objects, we find ages of 200-700Myr, implying formation redshifts of 7<z<14, and large stellar masses in the range 1-3x10^10M_sun. Analysis of 7 i-drops that are undetected at 3.6microns indicates that these are younger, considerably less massive systems. We calculate that line contamination should not severely affect our photometry or derived results. Using data out to 8microns, we find little evidence for substantial intrinsic dust reddening. Correcting for incompleteness in our sample, we find a lower limit on the comoving stellar mass density at z~6 to be 2.5x10^6M_sun/Mpc^3. We are able to explore the star formation histories of our selected galaxies, and we infer that the past global star formation rate may have been much higher than that observed at z~6. The associated UV flux we infer at z>7 could have played a major role in reionizing the universe.

We present VLBA observations of 1667 MHz OH maser emission from the massive star formation region G5.89-0.39. The observations were phase referenced allowing the absolute positions of the masers to be obtained. The 1667 MHz masers have radial velocities that span ~50 km/s but show little evidence of tracing the bipolar molecular outflow, as has been claimed in previous studies. We identify 23 Zeeman pairs through comparison of masers in left and right circular polarization. Magnetic field strengths range from -2 mG to +2 mG, and an ordered reversal in magnetic field direction is observed toward the southern region of the UC HII region. We suggest that the velocity and magnetic field structure of the 1667 MHz masers can be explained in the context of a model in which the masers arise in a neutral shell just outside a rapidly exanding ionized shell.

We exploit the gravitational potential of massive cluster lenses to probe the emission line properties of six z=1 galaxies which appear as highly magnified luminous arcs. Using the GMOS integral field spectrograph together with detailed cluster lens models we reconstruct the intrinsic morphologies and two-dimensional velocity fields in these galaxies on scales corresponds to \~0.5kpc (unlensed) at z=1. Four of the galaxies have stable disk-like kinematics, whilst the other two resemble interacting or starburst galaxies. These galaxies lie close to the mean rest-frame I-band Tully-Fisher relation for nearby spirals suggesting a clear preference for hierarchical growth of structure. In the rest-frame B-band, the observations suggest 0.5+/-0.3mag of brightening, consistent with increased star-formation activity at z=1. However, the galaxies with stable disk kinematics have more slowly rising rotation curves than expected from galaxies with similar surface brightness in the local Universe. We suggest that this may arise because the distant galaxies have lower bulge masses than their local counter-parts. Whilst this study is based on only six galaxies, the gain in flux and in spatial resolution achieved via gravitational magnification provides a much more detailed view of the high redshift Universe than possible with conventional surveys.

We present results from a search for high-redshift J--band ``dropout'' galaxies in the portion of the GOODS southern field that is covered by extremely deep imaging from the Hubble Ultradeep Field (HUDF).Using observations at optical, near-infrared and mid-infrared wavelengths from the Hubble and Spitzer Space Telescopes and the ESO-VLT, we search for very massive galaxies at high redshifts and find one particularly remarkable candidate. Its spectral energy distribution is consistent with a galaxy at z ~ 6.5 and a stellar mass of 6x10e11 M(sun) (for a Salpeter IMF). We interpret a prominent photometric break between the near-infrared and Spitzer bandpasses as the 3646A Balmer discontinuity. The best-fitting models have low reddening and ages of several hundred Myr, placing the formation of the bulk of the stars at z > 9. Alternative models of dusty galaxies at z ~ 2.5 are possible but provide significantly poorer fits. The object is detected with Spitzer at 24 micron. This emission originats from an obscured active nucleus or star formation. We present optical and near-infrared spectroscopy which has, thus far, failed to detect any spectral features. This helps limit the solution in which the galaxy is a starburst or active galaxy at z ~ 2.5. If the high-redshift interpretation is correct, this object would be an example of a galaxy that formed by a process strongly resembling traditional models of monolithic collapse, in a way which a very large mass of stars formed within a remarkably short period of time, at very high redshift.

(Abridged) We review recent observations that suggest the star formation rate density is declining over 3$<z<$7 and illustrate the challenges that this poses in explaining the assembled stellar mass in z$\sim$6 galaxies deduced from Spitzer data. A plausible conclusion is a vigorous period of earlier star formation. Prior to JWST and TMT, strong lensing offers a unique probe of the extent of this earlier activity. We discuss the first results of a blind spectroscopic survey of lensing clusters for 8.5$<z<$10 Ly$\alpha$ emitters using NIRSPEC which is achieving a limiting star formation rate of 0.1 $M_{\odot}$ yr$^{-1}$. A companion HST/IRAC survey is targeting lensed $z$ and $J$-band dropouts and probes a $\simeq$1 arcmin$^2$ region 1 magnitude deeper than the UDF/NICMOS observations. Both surveys will constrain the contribution of early, low luminosity, sources to cosmic reionization.

We report the detection of a z~7 galaxy strongly lensed by the massive galaxy cluster Abell 2218 (z=0.175) at 3.6 and 4.5 um using the Spitzer Observatory and at 1.1 um using the Hubble Space Telescope. The new data indicate a refined photometric redshift in the range of 6.6-6.8 depending on the presence of Ly-alpha emission. The spectral energy distribution is consistent with having a significant Balmer break, suggesting that the galaxy is in the poststarburst stage with an age of at least ~50 Myr and quite possibly a few hundred Myr. This suggests the possibility that a mature stellar population is already in place at such a high redshift. Compared with typical Lyman break galaxies at z~3-4, the stellar mass is an order of magnitude smaller (~10^{9} Msun), but the specific star formation rate (star formation rate/M_{star}) is similarly large (> 10^{-9} yr^{-1}), indicating equally vigorous star-forming activity.

We present 2" to 7" resolution 3 mm continuum and CO(J=1-0) line emission and near infrared Ks, H2, and [FeII] images toward the massive star forming region W75 N. The CO emission uncovers a complex morphology of multiple, overlapping outflows. A total flow mass of greater than 255 Msun extends 3 pc from end-to-end and is being driven by at least four late to early-B protostars. More than 10% of the molecular cloud has been accelerated to high velocities by the molecular flows (> 5.2 km/s relative to v{LSR}) and the mechanical energy in the outflowing gas is roughly half the gravitational binding energy of the cloud. The W75 N cluster members represent a range of evolutionary stages, from stars with no apparent circumstellar material to deeply embedded protostars that are actively powering massive outflows. Nine cores of millimeter-wavelength emission highlight the locations of embedded protostars in W75 N. The total mass of gas & dust associated with the millimeter cores ranges from 340 Msun to 11 Msun. The infrared reflection nebula and shocked H2 emission have multiple peaks and extensions which, again, suggests the presence of several outflows. Diffuse H2 emission extends about 0.6 parsecs beyond the outer boundaries of the CO emission while the [FeII] emission is only detected close to the protostars. The infrared line emission morphology suggests that only slow, non-dissociative J-type shocks exist throughout the pc-scale outflows. Fast, dissociative shocks, common in jet-driven low-mass outflows, are absent in W75 N. Thus, the energetics of the outflows from the late to early B protostars in W75 N differ from their low-mass counterparts -- they do not appear to be simply scaled-up versions of low-mass outflows.