A Deep Hubble and Keck Search for Definitive Identification of Lyman Continuum Emitters at z~3.1

Narrow-band imaging of the rest-frame Lyman continuum (LyC) of galaxies at z~3.1 has produced a large number of candidate LyC-emitting galaxies. These samples are contaminated by galaxies at lower redshift. To better understand LyC escape, we need an uncontaminated sample of galaxies that emit strongly in the LyC. Here we present deep Hubble imaging of five bright galaxies at z~3.1 that had previously been identified as candidate LyC-emitters with ground-based images. The WFC3 F336W images probe the LyC of galaxies at z>3.06 and provide an order-of-magnitude increase in spatial resolution over ground-based imaging. The non-ionizing UV images often show multiple galaxies (or components) within ~1 of the candidate LyC emission seen from the ground. In each case, only one of the components is emitting light in the F336W filter, which would indicate LyC escape if that component is at z>3.06. We use Keck/NIRSPEC near-IR spectroscopy to measure redshifts of these components to distinguish LyC-emitters from foreground contamination. We find that two candidates are low redshift contaminants, one candidate had a previously misidentified redshift, and the other two cannot be confirmed as LyC-emitters. The level of contamination is consistent with previous estimates. For the galaxies with z>3.06, we derive strong 1 sigma limits on the relative escape fraction between 0.07 and 0.09. We still do not have a sample of definitive LyC-emitters, and a much larger study of low luminosity galaxies is required. The combination of high resolution imaging and deep spectroscopy is critical for distinguishing LyC-emitters from foreground contaminants.

Consequences of bursty star formation on galaxy observables at high redshifts

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}[ u L_{ u}(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.

Ultraviolet Emission Lines in Young Low Mass Galaxies at z~2: Physical Properties and Implications for Studies at z>7

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.

Physical Properties of Emission-Line Galaxies at z ~ 2 from Near-Infrared Spectroscopy with Magellan FIRE

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ~ 2 obtained with the FIRE spectrometer on the Magellan Baade telescope. The sample was selected from the WISP survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 to detect emission-line galaxies over 0.3 < z < 2.3. Our FIRE follow-up spectroscopy (R~5000) over 1.0-2.5 micron permits detailed measurements of physical properties of the z~2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ~5-100 M_sun yr-1. We derive a median metallicity for the sample of ~0.45 Z_sun, and the estimated stellar masses range from ~10^8.5 - 10^9.5 M_sun. The average ionization parameters measured for the sample are typically much higher than what is found for local star-forming galaxies. We derive composite spectra from the FIRE sample, from which we infer typical nebular electron densities of ~100-400 cm^-3. Based on the location of the galaxies and composite spectra on BPT diagrams, we do not find evidence for significant AGN activity in the sample. Most of the galaxies as well as the composites are offset in the BPT diagram toward higher [O III]/H-beta at a given [N II]/H-alpha, in agreement with other observations of z > 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O III]/H-beta versus [S II]/H-alpha diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O III]/H-beta versus [N II]/H-alpha diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf-Rayet stars in starbursting galaxies at z~2. (Abridged)

Ultra-faint Ultraviolet Galaxies at z~2 Behind the Lensing Cluster Abell 1689: the Luminosity Function, Dust Extinction and Star Formation Rate Density

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]

A Deep HST Search for Escaping Lyman Continuum Flux at z~1.3: Evidence for an Evolving Ionizing Emissivity

We have obtained deep Hubble Space Telescope far-UV images of 15 starburst galaxies at z~1.3 in the GOODS fields to search for escaping Lyman continuum photons. These are the deepest far-UV images m_{AB}=28.7, 3sigma, 1 diameter) over this large an area (4.83 arcmin^2) and provide the best escape fraction constraints for any galaxy at any redshift. We do not detect any individual galaxies, with 3sigma limits to the Lyman Continuum (~700 AA) flux 50--149 times fainter (in f_nu) than the rest-frame UV (1500 AA) continuum fluxes. Correcting for the mean IGM attenuation (factor ~2), as well as an intrinsic stellar Lyman Break (~3), these limits translate to relative escape fraction limits of f_{esc,rel}<[0.03,0.21]. The stacked limit is f_{esc,rel}(3sigma)<0.02. We use a Monte Carlo simulation to properly account for the expected distribution of IGM opacities. When including constraints from previous surveys at z~1.3 we find that, at the 95% confidence level, no more than 8% of star--forming galaxies at z~1.3 can have relative escape fractions greater than 0.50. Alternatively, if the majority of galaxies have low, but non-zero, escaping Lyman Continuum, the escape fraction can not be more than 0.04. Both the stacked limits, and the limits from the Monte Carlo simulation suggest that the average ionizing emissivity (relative to non-ionizing UV emissivity) at z~1.3 is significantly lower than has been observed in Lyman Break Galaxies (LBGs) at z~3. If the ionizing emissivity of star-forming galaxies is in fact increasing with redshift, it would help to explain the high photoionization rates seen in the IGM at z>4 and reionization of the intergalactic medium at z>6. [Abridged]

Detection of PAH and Far-Infrared Emission from the Cosmic Eye: Probing the Dust and Star Formation of Lyman Break Galaxies

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.

Spitzer Observations of the z=2.73 Lensed Lyman Break Galaxy, MS1512-cB58

We present Spitzer infrared (IR) photometry and spectroscopy of the lensed Lyman break galaxy (LBG), MS1512-cB58 at z=2.73. The large (factor ~30) magnification allows for the most detailed infrared study of an L*_UV(z=3) LBG to date. Broadband photometry with IRAC (3-10 micron), IRS (16 micron), and MIPS (24, 70 & 160 micron) was obtained as well as IRS spectroscopy spanning 5.5-35 microns. A fit of stellar population models to the optical/near-IR/IRAC photometry gives a young age (~9 Myr), forming stars at ~98 M_sun/yr, with a total stellar mass of ~10^9 M_sun formed thus far. The existence of an old stellar population with twice the stellar mass can not be ruled out. IR spectral energy distribution fits to the 24 and 70 micron photometry, as well as previously obtained submm/mm, data give an intrinsic IR luminosity L_IR = 1-2 x10^11 L_sun and a star formation rate, SFR ~20-40 M_sun/yr. The UV derived star formation rate (SFR) is ~3-5 times higher than the SFR determined using L_IR or L_Halpha because the red UV spectral slope is significantly over predicting the level of dust extinction. This suggests that the assumed Calzetti starburst obscuration law may not be valid for young LBGs. We detect strong line emission from Polycyclic Aromatic Hydrocarbons (PAHs) at 6.2, 7.7, and 8.6 microns. The line ratios are consistent with ratios observed in both local and high redshift starbursts. Both the PAH and rest-frame 8 micron luminosities predict the total L_IR based on previously measured relations in starbursts. Finally, we do not detect the 3.3 micron PAH feature. This is marginally inconsistent with some PAH emission models, but still consistent with PAH ratios measured in many local star-forming galaxies.

High-Redshift QSOs in the SWIRE Survey and the z~3 QSO Luminosity Function

We use a simple optical/infrared (IR) photometric selection of high-redshift QSOs that identifies a Lyman Break in the optical photometry and requires a red IR color to distinguish QSOs from common interlopers. The search yields 100 z~3 (U-dropout) QSO candidates with 19<r<22 over 11.7 deg^2 in the ELAIS-N1 (EN1) and ELAIS-N2 (EN2) fields of the Spitzer Wide-area Infrared Extragalactic (SWIRE) Legacy Survey. The z~3 selection is reliable, with spectroscopic follow-up of 10 candidates confirming they are all QSOs at 2.83<z<3.44. We find that our z~4$ (g-dropout) sample suffers from both unreliability and incompleteness but present 7 previously unidentified QSOs at 3.50<z<3.89. Detailed simulations show our z~3 completeness to be ~80-90% from 3.0<z<3.5, significantly better than the ~30-80% completeness of the SDSS at these redshifts. The resulting luminosity function extends two magnitudes fainter than SDSS and has a faint end slope of beta=-1.42 +- 0.15, consistent with values measured at lower redshift. Therefore, we see no evidence for evolution of the faint end slope of the QSO luminosity function. Including the SDSS QSO sample, we have now directly measured the space density of QSOs responsible for ~70% of the QSO UV luminosity density at z~3. We derive a maximum rate of HI photoionization from QSOs at z~3.2, Gamma = 4.8x10^-13 s^-1, about half of the total rate inferred through studies of the Ly-alpha forest. Therefore, star-forming galaxies and QSOs must contribute comparably to the photoionization of HI in the intergalactic medium at z~3.

New Constraints on the Lyman Continuum Escape Fraction at z~1.3

We examine deep far-ultraviolet (1600 Angstrom) imaging of the Hubble Deep Field-North (HDFN) and the Hubble Ultra Deep Field (HUDF) to search for leaking Lyman continuum radiation from starburst galaxies at z~1.3. There are 21 (primarily sub-L*) galaxies with spectroscopic redshifts between 1.1<z<1.5 and none are detected in the far-UV. We fit stellar population templates to the galaxies optical/near-infrared SEDs to determine the starburst age and level of dust attenuation, giving an accurate estimate of the intrinsic Lyman continuum ratio, f_1500/f_700, and allowing a conversion from f_700 limits to relative escape fractions. We show that previous high-redshift studies may have underestimated the amplitude of the Lyman Break, and thus the relative escape fraction, by a factor of ~2. Once the starburst age and intergalactic HI absorption are accounted for, 18 galaxies in our sample have limits to the relative escape fraction, f_esc,rel < 1.0 with some limits as low as f_esc,rel < 0.10 and a stacked limit of f_esc,rel < 0.08. This demonstrates, for the first time, that most sub-L* galaxies at high redshift do not have large escape fractions. When combined with a similar study of more luminous galaxies at the same redshift we show that, if all star-forming galaxies at z~1 have similar relative escape fractions, the value must be less than 0.14 (3 sigma). We also show that less than 20% (3 sigma) of star-forming galaxies at z~1 have relative escape fractions near unity. These limits contrast with the large escape fractions found at z~3 and suggest that the average escape fraction has decreased between z~3 and z~1. (Abridged)