Search

We present the results of a Keck and NOEMA spectroscopic survey of 507 galaxies, where we confirm the presence of two massive overdensities at $z = 3.090 - 3.110$ and $z = 3.133 - 3.155$ in the neighborhood of the GOODS-N, each with over a dozen spectroscopically confirmed members. We find that both of these have galaxy overdensities of NIR-detected galaxies of $\delta_{\rm gal, obs} = 6 - 9$ within corrected volumes of $(6 - 7) \times 10^3~{\rm cMpc}^3$. We estimate the properties of the $z = 0$ descendants of these overdensities using a spherical collapse model and find that both should virialize by $z \simeq 0.5 - 0.8$, with total masses of $M_{\rm tot} \simeq (6 - 7) \times 10^{14}~{\rm M}_\odot$. The same spherical collapse calculations, as well as a clustering-of-clusters statistical analysis, suggest a >80% likelihood that the two overdensities will collapse into a single cluster with $M_{\rm tot} = (1.0 - 1.5) \times 10^{15}~{\rm M}_\odot$ by $z \sim 0.1-0.4$. The $z = 3.14$ substructure contains a core of four bright dusty star-forming galaxies with $\Sigma {\rm SFR} = 2700 \pm 700~{\rm M}_\odot~{\rm yr}^{-1}$ in a volume of only 280 ${\rm cMpc}^3$., Comment: 18 pages, 8 figures. Accepted to ApJ

The Hawaii Survey Field SSA22 is the fourth deepest Chandra X-ray field. To allow for the fullest exploration of this field, we present new optical spectroscopy from Keck/DEIMOS and Keck/LRIS, which, in combination with the literature, brings the spectroscopic completeness of the 2--8 keV sample to 62%. We also make optical spectral classifications and estimate photometric redshifts for the sources without spectroscopic redshifts. We then determine hard X-ray luminosity functions (XLFs) for the full sample of Active Galactic Nuclei (AGNs), as well as for the broad-line AGNs (BLAGNs) and the non-BLAGNs separately. Our XLF for the full sample is in good agreement with the literature, showing relatively strong evolution over the redshift range $0.25\le z < 4$. The XLFs for the BLAGNs and the non-BLAGNs imply distinct evolution with redshift, with BLAGNs becoming increasingly dominant at higher redshifts and X-ray luminosities., Comment: 11 pages, 7 figures, Accepted to ApJ

Low-metallicity galaxies may provide key insights into the evolutionary history of galaxies. Galaxies with strong emission lines and high equivalent widths (rest-frame EW(H-beta) > 30 A) are ideal candidates for the lowest metallicity galaxies to z ~ 1. Using a Keck/DEIMOS spectral database of about 18,000 galaxies between z = 0.2 and z = 1, we search for such extreme emission-line galaxies with the goal of determining their metallicities. Using the robust direct Te method, we identify 8 new extremely metal-poor galaxies (XMPGs) with 12 + log O/H < 7.65, including one at 6.949 +/- 0.091, making it the lowest metallicity galaxy reported to date at these redshifts. We also improve upon the metallicities for two other XMPGs from previous work. We investigate the evolution of H-beta using both instantaneous and continuous starburst models, finding that XMPGs are best characterized by continuous starburst models. Finally, we study the dependence on age of the build-up of metals and the emission-line strength., Comment: 13 pages, 11 figures, 1 table, accepted for publication in ApJ V2: Matches published version

We present a search for extreme emission line galaxies (EELGs) at z<1 in the COSMOS and North Ecliptic Pole (NEP) fields with imaging from Subaru/Hyper Suprime-Cam (HSC) and a combination of new and existing spectroscopy. We select EELGs on the basis of substantial excess flux in the z broadband, which is sensitive to H$\alpha$ at 0.3

Ultraluminous Lyman alpha (Lya) emitting galaxies (ULLAEs) with log L (Lya)>43.5 erg/s near the epoch of reionization (z>5) make up the bright end of the LAE luminosity function (LF) and may provide insight into the process of reionization, including the formation of ionized bubbles around these extreme systems. We present a spectroscopic LF for ULLAEs at z=5.7. We used data from the HEROES ~45 sq. deg Subaru Hyper Suprime-Cam survey, which is centered on the North Ecliptic Pole and has both broadband (grizY) and narrowband (NB816 and NB921) imaging, to select candidate ULLAEs based on a NB816 excess and a strong broadband Lyman break. We spectroscopically observed 17 ULLAE candidates with DEIMOS on Keck II. We confirmed 12 as LAEs at z=5.7, 9 of which are ULLAEs. The remaining sources are an AGN at z=5.7, an [OIII]5007 emitter at z=0.63, a red star, and two spectroscopic non-detections. Using the 9 confirmed ULLAEs, we construct a ULLAE LF at z=5.7. After applying a comprehensive incompleteness correction, we compare our new z=5.7 LF with our recent z=6.6 LF and with other LFs from the literature to look for evolution at the ultraluminous end. We find the overall ratio of the z=5.7 to z=6.6 ULLAE comoving number densities to be 1.92 (+1.12, -0.71), which corresponds to a LF offset of 0.28 (+0.20, -0.20) dex., Comment: Accepted for publication in ApJ. 13 pages, 8 Figures, 3 Tables. v2 integrates an erratum to correct the coordinates of NEP5.7AGN. This change has been made directly in the article text, tables, and figures, and does not affect the results or conclusions of the work

We report spectroscopic redshift measurements for three bright submillimeter galaxies (SMGs) near the GOODS-N field, each with SCUBA-2 850 micron fluxes > 10 mJy, using the Northern Extended Millimeter Array (NOEMA). Our molecular linescan observations of these sources, which occupy a ~7 arcmin$^2$ area outside of the HST coverage of the field, reveal that two lie at $z \sim$ 3.14. In the remaining object, we detect line emission consistent with CO(7-6), [C I], and H2O at $z$ = 4.42. The far-infrared spectral energy distributions of these galaxies, constrained by SCUBA-2, NOEMA, and Herschel/SPIRE, indicate instantaneous SFRs $\sim4000 ~{\rm M_{\odot}~yr^{-1}}$ in the $z$ = 4.42 galaxy and $\sim 2500~{\rm M_{\odot}~yr^{-1}}$ in the two $z \sim$ 3.14 galaxies. Based on our sources' CO line luminosities, we estimate $M_{{\rm gas}}\sim10^{11} M_{\odot}$ and find gas depletion timescales of $\tau_{{\rm depl}}\sim 50$ Myr, consistent with findings in other high-redshift SMGs. Finally, we show that the two $z \sim$ 3.14 sources, which alone occupy a volume $\sim$10 Mpc$^3$, very likely mark the location of a protocluster of bright SMGs and less dusty optical sources., Comment: Newest version: 12 pages, 6 figures, 2 tables; added references and revised arguments in Section 4, results unchanged. Accepted to ApJ. First version: 8 pages, 3 figures, 2 tables. Submitted to ApJ Letters

We use the wealth of deep archival optical spectroscopy on the GOODS-South field from Keck, the VLT, and other facilities to select candidate high-redshift Lyman continuum (LyC) leakers in the Hubble Deep UV Legacy Survey (HDUV) dataset. We select sources at $2.35 < z < 3.05$, where the HST/WFC3 F275W filter probes only the redshifted LyC. We find five moderately F275W-bright sources (four detected at $\gtrsim3\sigma$ significance) in this redshift range. However, two of these show evidence in their optical spectra for contamination by foreground galaxies along the line-of-sight. We then perform an F275W error-weighted sum of the fluxes of all 129 galaxies at $2.35 < z < 3.05$ in both the GOODS-N and GOODS-S HDUV areas to estimate the total ionizing flux. The result is dominated by just five candidate F275W-bright LyC sources. Lastly, we examine the contributions to the metagalactic ionizing background, finding that, at the sensitivity of the HDUV F275W data and allowing for the effects of LyC transmission in the intergalactic medium (IGM), star-forming galaxies can match the UV flux required to maintain an ionized IGM at $z \sim 2.5$., Comment: 13 pages, 7 figures; accepted to ApJ

We present the luminosity function (LF) for ultraluminous Ly$\alpha$ emitting galaxies (LAEs) at z = 6.6. We define ultraluminous LAEs (ULLAEs) as galaxies with logL(Ly$\alpha$) > 43.5 erg s$^{-1}$. We select our main sample using the g', r', i', z', and NB921 observations of a wide-area (30 deg$^2$) Hyper Suprime-Cam survey of the North Ecliptic Pole (NEP) field. We select candidates with g', r', i' > 26, NB921 $\leq$ 23.5, and NB921 - z' $\leq$ 1.3. Using the DEIMOS spectrograph on Keck II, we confirm 9 of our 14 candidates as ULLAEs at z = 6.6 and the remaining 5 as an AGN at z = 6.6, two [OIII]$\lambda$5007 emitting galaxies at z = 0.84 and z = 0.85, and two non-detections. This emphasizes the need for full spectroscopic follow-up to determine accurate LFs. In constructing the ULLAE LF at z = 6.6, we combine our 9 NEP ULLAEs with two previously discovered and confirmed ULLAEs in the COSMOS field: CR7 and COLA1. We apply rigorous corrections for incompleteness based on simulations. We compare our ULLAE LF at z = 6.6 with LFs at z = 5.7 and z = 6.6 from the literature. Our data reject some previous LF normalizations and power law indices, but they are broadly consistent with others. Indeed, a comparative analysis of the different literature LFs suggests that none is fully consistent with any of the others, making it critical to determine the evolution from z = 5.7 to z = 6.6 using LFs constructed in exactly the same way at both redshifts., Comment: Accepted for publication in ApJ. 12 pages, 8 Figures, 3 Tables

We analyze the >4-sigma sources in the most sensitive 100 arcmin^2 area (rms <0.56 mJy) of a SCUBA-2 850 micron survey of the GOODS-S and present the 75 band 7 ALMA sources (>4.5-sigma) obtained from high-resolution interferometric follow-up observations. The SCUBA-2---and hence ALMA---samples should be complete to 2.25 mJy. Of the 53 SCUBA-2 sources in this complete sample, only five have no ALMA detections, while 13% (68% confidence range 7-19%) have multiple ALMA counterparts. Color-based high-redshift dusty galaxy selection techniques find at most 55% of the total ALMA sample. In addition to using literature spectroscopic and optical/NIR photometric redshifts, we estimate FIR photometric redshifts based on an Arp 220 template. We identify seven z>4 candidates. We see the expected decline with redshift of the 4.5 micron and 24 micron to 850 micron flux ratios, confirming these as good diagnostics of z>4 candidates. We visually classify 52 ALMA sources, finding 44% (68% confidence range 35-53%) to be apparent mergers. We calculate rest-frame 2-8 keV and 8-28 keV luminosities using the 7 Ms Chandra X-ray image. Nearly all of the ALMA sources detected at 0.5-2 keV are consistent with a known X-ray luminosity to 850 micron flux relation for star-forming galaxies, while most of those detected at 2-7 keV are moderate luminosity AGNs that lie just above the 2-7 keV detection threshold. The latter largely have substantial obscurations of log N_H = 23-24 cm^-2, but two of the high-redshift candidates may even be Compton thick., Comment: 44 pages, 35 figures. Submitted to ApJ

We report on a search for ultraluminous Lyman alpha emitting galaxies (LAEs) at z=6.6 using the NB921 filter on Hyper Suprime-Cam on the Subaru telescope. We searched a 30 degree squared area around the North Ecliptic Pole, which we observed in broadband g', r', i', z', and y' and narrowband NB816 and NB921, for sources with NB921 < 23.5 and z' - NB921 > 1.3. This corresponds to a selection of log L(Ly-alpha) > 43.5 erg/s. We followed up seven candidate LAEs (out of thirteen) with the Keck DEIMOS spectrograph and confirmed five z=6.6 LAEs, one z=6.6 AGN with a broad Ly-alpha line and a strong red continuum, and one low-redshift ([OIII]5007) galaxy. The five ultraluminous LAEs have wider line profiles than lower luminosity LAEs, and one source, NEPLA4, has a complex line profile similar to that of COLA1. In combination with previous results, we show that the line profiles of the z=6.6 ultraluminous LAEs are systematically different than those of lower luminosity LAEs at this redshift. This result suggests that ultraluminous LAEs generate highly ionized regions of the intergalactic medium in their vicinity that allow the full Lyman alpha profile of the galaxy---including any blue wings---to be visible. If this interpretation is correct, then ultraluminous LAEs offer a unique opportunity to determine the properties of the ionized zones around them, which will help in understanding the ionization of the z ~ 7 intergalactic medium. A simple calculation gives a very rough estimate of 0.015 for the escape fraction of ionizing photons, but more sophisticated calculations are needed to fully characterize the uncertainties., Comment: Accepted for publication in ApJ. 13 pages, 13 figures

We present a flux-limited sample of $z\sim0.3$ Ly$\alpha$ emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism spectroscopic data. The published GALEX $z\sim0.3$ LAE sample is pre-selected from continuum-bright objects and thus is biased against high equivalent width (EW) LAEs. We remove this continuum pre-selection and compute the EW distribution and the luminosity function of the Ly$\alpha$ emission line directly from our sample. We examine the evolution of these quantities from $z\sim0.3$ to $2.2$ and find that the EW distribution shows little evidence for evolution over this redshift range. As shown by previous studies, the Ly$\alpha$ luminosity density from star-forming galaxies declines rapidly with declining redshift. However, we find that the decline in Ly$\alpha$ luminosity density from $z=2.2$ to $z=0.3$ may simply mirror the decline seen in the H$\alpha$ luminosity density from $z=2.2$ to $z=0.4$, implying little change in the volumetric Ly$\alpha$ escape fraction. Finally, we show that the observed Ly$\alpha$ luminosity density from AGNs is comparable to the observed Ly$\alpha$ luminosity density from star-forming galaxies at $z=0.3$. We suggest that this significant contribution from AGNs to the total observed Ly$\alpha$ luminosity density persists out to $z\sim2.2$., Comment: 27 pages, 19 figures, accepted for publication in ApJ

We use ultradeep 20 cm data from the Karl G. Jansky Very Large Array and 850 micron data from SCUBA-2 and the Submillimeter Array of an 124 arcmin^2 region of the Chandra Deep Field-north to analyze the high radio power (P_20cm>10^31 erg s^-1 Hz^-1) population. We find that 20 (42+/-9%) of the spectroscopically identified z>0.8 sources have consistent star formation rates (SFRs) inferred from both submillimeter and radio observations, while the remaining sources have lower (mostly undetected) submillimeter fluxes, suggesting that active galactic nucleus (AGN) activity dominates the radio power in these sources. We develop a classification scheme based on the ratio of submillimeter flux to radio power versus radio power and find that it agrees with AGN and star-forming galaxy classifications from Very Long Baseline Interferometry. Our results provide support for an extremely rapid drop in the number of high SFR galaxies above about a thousand solar masses per year (Kroupa initial mass function) and for the locally determined relation between X-ray luminosity and radio power for star-forming galaxies applying at high redshifts and high radio powers. We measure far-infrared (FIR) luminosities and find that some AGNs lie on the FIR-radio correlation, while others scatter below. The AGNs that lie on the correlation appear to do so based on their emission from the AGN torus. We measure a median radio size of 1.0+/-0.3 arcsecond for the star-forming galaxies. The radio sizes of the star-forming galaxies are generally larger than those of the AGNs., Comment: Published in ApJ. 10 pages, 8 figures

In this first paper in the SUPER GOODS series on powerfully star-forming galaxies in the two GOODS fields, we present a deep SCUBA-2 survey of the GOODS-N at both 850 and 450 micron (central rms noise of 0.28 mJy and 2.6 mJy, respectively). In the central region the 850 micron observations cover the GOODS-N to near the confusion limit of ~1.65 mJy, while over a wider 450 arcmin^2 region---well complemented by Herschel far-infrared imaging---they have a median 4-sigma limit of 3.5 mJy. We present >4-sigma catalogs of 186 850 micron and 31 450 micron selected sources. We use interferometric observations from the SMA and the VLA to obtain precise positions for 114 SCUBA-2 sources (28 from the SMA, all of which are also VLA sources). We present new spectroscopic redshifts and include all existing spectroscopic or photometric redshifts. We also compare redshifts estimated using the 20 cm to 850 micron and the 250 micron to 850 micron flux ratios. We show that the redshift distribution increases with increasing flux, and we parameterize the dependence. We compute the star formation history and the star formation rate (SFR) density distribution functions in various redshift intervals, finding that they reach a peak at z=2-3 before dropping to higher redshifts. We show that the number density per unit volume of SFR>500 solar mass per year galaxies measured from the SCUBA-2 sample does not change much relative to that of lower SFR galaxies from UV selected samples over z=2-5, suggesting that, apart from changes in the normalization, the shape in the number density as a function of SFR is invariant over this redshift interval., Comment: Accepted for publication in ApJ. 38 pages, 26 figures

We combine the Herschel Space Observatory PACS and SPIRE photometry with archival WISE photometry to construct the spectral energy distributions (SED) for over 300 local ($z < 0.05$), ultra-hard X-ray (14 - 195 keV) selected active galactic nuclei (AGN) from the Swift Burst Alert Telescope (BAT) 58 month catalogue. Using a simple analytical model that combines an exponentially cut-off powerlaw with a single temperature modified blackbody, we decompose the SEDs into a host-galaxy and AGN component. We calculate dust masses, dust temperatures, and star-formation rates (SFR) for our entire sample and compare them to a stellar mass-matched sample of local non-AGN galaxies. We find AGN host galaxies have systematically higher dust masses, dust temperatures, and SFRs due to the higher prevalence of late-type galaxies to host an AGN, in agreement with previous studies of the Swift/BAT AGN. We provide a scaling to convert X-ray luminosities into 8 - 1000 $\mu$m AGN luminosities, as well as determine the best mid-to-far IR colors for identifying AGN dominated galaxies in the IR regime. We find that for nearly 30 per cent of our sample, the 70 $\mu$m emission contains a significant contribution from the AGN ($> 0.5$), especially at higher luminosities ($L_{14-195\,\rm{keV}} > 10^{42.5}$ ergs s$^{-1}$). Finally, we measure the local SFR-AGN luminosity relationship, finding a slope of 0.18, large scatter (0.37 dex), and no evidence for an upturn at high AGN luminosity. We conclude with a discussion on the implications of our results within the context of galaxy evolution with and without AGN feedback., Comment: 26 pages, 17 figures, 3 Tables, accepted for publication in MNRAS

We present far-infrared (FIR) and submillimeter photometry from the Herschel Space Observatory's Spectral and Photometric Imaging Receiver (SPIRE) for 313 nearby z<0.05 active galactic nuclei (AGN). We selected AGN from the 58 month Swift Burst Alert Telescope (BAT) catalog, the result of an all-sky survey in the 14-195 keV energy band, allowing for a reduction in AGN selection effects due to obscuration and host galaxy contamination. We find 46% (143/313) of our sample is detected at all three wavebands and combined with our PACS observations represents the most complete FIR spectral energy distributions of local, moderate luminosity AGN. We find no correlation between the 250, 350, and 500 micron luminosities with 14-195 keV luminosity, indicating the bulk of the FIR emission is not related to the AGN. However, Seyfert 1s do show a very weak correlation with X-ray luminosity compared to Seyfert 2s and we discuss possible explanations. We compare the SPIRE colors (F250/F350 and F350/F500) to a sample of normal star-forming galaxies, finding the two samples are statistically similar, especially after matching in stellar mass. But a color-color plot reveals a fraction of the Herschel-BAT AGN are displaced from the normal star-forming galaxies due to excess 500 micron emission E500). Our analysis shows E500 is strongly correlated with the 14-195 keV luminosity and 3.4/4.6 micron flux ratio, evidence the excess is related to the AGN. We speculate these sources are experiencing millimeter excess emission originating in the corona of the accretion disk., Comment: 22 pages, 11 figures, accepted for publication in MNRAS. We have created a Dataverse that contains all of the FITS images for the Herschel-BAT AGN as well as the photometry at https://dataverse.harvard.edu/dataverse/herschel-bat-agn . Full versions of Tables 1 and 2 and Figure 1 are available at http://www.astro.umd.edu/~tshimizu/#research

We examine the redshift evolution of standard strong emission-line diagnostics for Hbeta-selected star-forming galaxies using the local SDSS sample and a new z = 0.2 - 2.3 sample obtained from HST WFC3 grism and Keck DEIMOS and MOSFIRE data. We use the SDSS galaxies to show that there is a systematic dependence of the strong emission-line properties on Balmer-line luminosity, which we interpret as showing that both the N/O abundance and the ionization parameter increase with increasing line luminosity. Allowing for the luminosity dependence tightens the diagnostic diagrams and the metallicity calibrations. The combined SDSS and high-redshift samples show that there is no redshift evolution in the line properties once the luminosity correction is applied, i.e., all galaxies with a given L(Hbeta) have similar strong emission-line distributions at all the observed redshifts. We argue that the best metal diagnostic for the high-redshift galaxies may be a luminosity-adjusted version of the [NII]6584/Halpha metallicity relation., Comment: Accepted for publication in ApJ. 16 pages, 15 figures

Our recent estimates of galaxy counts and the luminosity density in the near-infrared (Keenan et al. 2010, 2012) indicated that the local universe may be under-dense on radial scales of several hundred megaparsecs (Mpc). Such a large-scale local under-density could introduce significant biases in the measurement and interpretation of cosmological observables, such as the inferred effects of dark energy on the rate of expansion. In Keenan et al. (2013), we measured the K-band luminosity density as a function of distance from us to test for such a local under-density. We made this measurement over the redshift range 0.01 < z < 0.2 (radial distances D ~ 50-800 Mpc). We found that the shape of the K-band luminosity function is relatively constant as a function of distance and environment. We derive a local (z < 0.07, D < 300 Mpc) K-band luminosity density that agrees well with previously published studies. At z > 0.07, we measure an increasing luminosity density that by z~ 0.1 rises to a value of ~1.5 times higher than that measured locally. This implies that the stellar mass density follows a similar trend. Assuming that the underlying dark matter distribution is traced by this luminous matter, this suggests that the local mass density may be lower than the global mass density of the universe at an amplitude and on a scale that is sufficient to introduce significant biases into the measurement of basic cosmological observables. At least one study has shown that an under-density of roughly this amplitude and scale could resolve the apparent tension between direct local measurements of the Hubble constant and those inferred by Planck team. Other theoretical studies have concluded that such an under-density could account for what looks like an accelerating expansion, even when no dark energy is present., Comment: Proceedings of IAU Symposium 308 "The Zeldovich Universe: Genesis and Growth of the Cosmic Web", 23-28 June 2014, Tallinn, Estonia

We construct a flux-limited sample of 135 candidate z~1 Lya emitters (LAEs) from Galaxy Evolution Explorer (GALEX) grism data using a new data cube search method. These LAEs have luminosities comparable to those at high redshifts and lie within a 7 Gyr gap present in existing LAE samples. We use archival and newly obtained optical spectra to verify the UV redshifts of these LAEs. We use the combination of the GALEX UV spectra, optical spectra, and X-ray imaging data to estimate the active galactic nucleus (AGN) fraction and its dependence on Lya luminosity. We remove the AGNs and compute the luminosity function (LF) from 60 z~1 LAE galaxies. We find that the best fit LF implies a luminosity density increase by a factor of ~1.5 from z~0.3 to z~1 and ~20 from z~1 to z~2. We find a z~1 volumetric Lya escape fraction of 0.7+/-0.4%., Comment: 16 pages, accepted for publication in ApJ

We obtained SMA observations of eight faint (intrinsic 850 $\mu$m fluxes $<$ 2 mJy) submillimeter galaxies (SMGs) discovered in SCUBA images of the massive lensing cluster fields A370, A2390, and A1689 and detected five. In total, we obtain 5 SMA detections, all of which have de-lensed fluxes $<$1 mJy with estimated total infrared luminosities 10$^{10}-10^{12}$ $L_\odot$, comparable to luminous infrared galaxies (LIRGs) and normal star-forming galaxies. Based on the latest number counts, these galaxies contribute $\sim$70% of the 850 $\mu$m extragalactic background light and represent the dominant star-forming galaxy population in the dusty universe. However, only 40$^{+30}_{-16}$% of our faint SMGs would be detected in deep optical or near-infrared surveys, which suggests many of these sources are at high redshifts ($z \gtrsim 3$) or extremely dusty, and they are not included in current star formation history estimates., Comment: 12 pages, 11 figures, ApJ in press

We use the SCUBA-2 submillimeter camera mounted on the JCMT to obtain extremely deep number counts at 450 and 850um. We combine data on two cluster lensing fields, A1689 and A370, and three blank fields, CDF-N, CDF-S, and COSMOS, to measure the counts over a wide flux range at each wavelength. We use statistical fits to broken power law representations to determine the number counts. This allows us to probe to the deepest possible level in the data. At both wavelengths our results agree well with the literature in the flux range over which they have been measured, with the exception of the 850um counts in CDF-S, where we do not observe the counts deficit found by previous single-dish observations. At 450um, we detect significant counts down to ~1mJy, an unprecedented depth at this wavelength. By integrating the number counts above this flux limit, we measure 113.9^{+49.7}_{-28.4} Jydeg^{-2} of the 450um extragalactic background light (EBL). The majority of this contribution is from sources with S_450um between 1-10mJy, and these sources are likely to be the ones that are analogous to the local luminous infrared galaxies (LIRGs). At 850um, we measure 37.3^{+21.1}_{-12.9} Jydeg^{-2} of the EBL. Because of the large systematic uncertainties on the COBE measurements, the percentage of the EBL we resolve could range from 48%-153% (44%-178%) at 450 (850)um. Based on high-resolution SMA observations of around half of the 4sigma 850um sample in CDF-N, we find that 12.5^{+12.1}_{-6.8}% of the sources are blends of multiple fainter sources. This is a low multiple fraction, and we find no significant difference between our original SCUBA-2 850um counts and the multiplicity corrected counts., Comment: update fig14 with bugs fixed, 14 pages, 14 figures, ApJ in press

We present the results from a large 850 micron survey of the sigma Orionis cluster using the SCUBA-2 camera on the James Clerk Maxwell Telescope. The 0.5-degree diameter circular region we surveyed contains 297 young stellar objects with an age estimated at about 3Myr. We detect 9 of these objects, 8 of which have infrared excesses from an inner disc. We also serendipitously detect 3 non-stellar sources at > 5sigma that are likely background submillimetre galaxies. The 9 detected stars have inferred disc masses ranging from 5 to about 17MJup, assuming similar dust properties as Taurus discs and an ISM gas-to-dust ratio of 100. There is a net positive signal toward the positions of the individually undetected infrared excess sources indicating a mean disc mass of 0.5 MJup . Stacking the emission toward those stars without infrared excesses constrains their mean disc mass to less than 0.3MJup, or an equivalent Earth mass in dust. The submillimetre luminosity distribution is significantly different from that in the younger Taurus region, indicating disc mass evolution as star forming regions age and the infrared excess fraction decreases. Submillimeter Array observations reveal CO emission toward 4 sources demonstrating that some, but probably not much, molecular gas remains in these relatively evolved discs. These observations provide new constraints on the dust and gas mass of protoplanetary discs during the giant planet building phase and provide a reference level for future studies of disc evolution., Comment: 10 pages, 6 figures, accepted for publication in MNRAS

Galaxy counts and recent measurements of the luminosity density in the near-infrared (NIR) have indicated the possibility that the local universe may be under-dense on scales of several hundred megaparsecs. The presence of a large-scale under-density in the local universe could introduce significant biases into the interpretation of cosmological observables, and, in particular, into the inferred effects of dark energy on the expansion rate. Here we measure the K-band luminosity density as a function of redshift to test for such a local under-density. In this effort, we combine photometry from UKIDSS and 2MASS with redshifts from the SDSS, 2DFGRS, 6DFGRS, 2MR, and GAMA surveys. We find that the overall shape of the z=0 rest-frame K-band luminosity function (M* = -22.15 +/- 0.04 and alpha = -1.02 +/- 0.03) appears to be relatively constant as a function of environment and redshift out to z ~0.2. We find a local (z < 0.07) luminosity density that is in good agreement with previous studies. Beyond z ~ 0.07 we detect a rising luminosity density that reaches a value ~1.5 times higher than that measured locally at z>0.1. This suggests that the stellar mass density as a function of redshift follows a similar trend. Assuming that luminous matter traces the underlying dark matter distribution, this implies that the local mass density of the universe may be lower than the global value on a scale and amplitude sufficient to introduce significant biases into the determination of basic cosmological observables, such as the expansion rate. An under-density of roughly this scale and amplitude would be sufficient to resolve the apparent tension between direct measurements of the Hubble constant and those inferred by Planck., Comment: Accepted by ApJ. Updated to match final version

We present the results of SCUBA2 observations at 450 micron and 850 micron of the field lensed by the massive cluster A370. With a total survey area > 100 arcmin2 and 1 sigma sensitivities of 3.92 and 0.82 mJy/beam at 450 and 850 micron respectively, we find a secure sample of 20 sources at 450 micron and 26 sources at 850 micron with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte Carlo simulations, we derive the number counts at both wavelengths. The 450 micron number counts probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope at similar wavelengths, and we estimate that ~47-61% of the 450 micron extragalactic background light (EBL) resolved into individual sources with 450 micron fluxes greater than 4.5 mJy. The faint 450 micron sources in the 4 sigma sample have positional accuracies of 3 arcseconds, while brighter sources (signal-to-noise > 6 sigma) are good to 1.4 arcseconds. Using the deep radio map (1 sigma ~ 6 uJy) we find that the percentage of submillimeter sources having secure radio counterparts is 85% for 450 micron sources with intrinsic fluxes > 6 mJy and 67% for 850 micron sources with intrinsic fluxes > 4 mJy. We also find that 67% of the > 4 sigma 450 micron sources are detected at 850 micron, while the recovery rate at 450 micron of > 4 sigma 850 micron sources is 54%. Combined with the source redshifts estimated using millimetric flux ratios, the recovered rate is consistent with the scenario where both 450 micron and 20 cm emission preferentially select lower redshift dusty sources, while 850 micron emission traces a higher fraction of dusty sources at higher redshifts. [Abridge], Comment: 22 pages, 14 figures, ApJ in press. Minor changes to the submitted version

Current studies on Submillimeter Galaxies (SMGs) mostly focus on bright sources with 850 micron flux greater than 2 mJy, and the results have shown that they are likely high redshift mergers with z > 2 and could be a dominant population on star formation in the early Universe. However, bright SMGs only contributes 20-30% of the 850 micron extragalactic background light (EBL), meaning the bulk of the cosmic star formation still hidden by dust and our current understanding is biased. We have started a program to study an unbiased sample of highly-amplified and intrinsically faint SCUBA detected SMGs in the field of massive lensing clusters. Here we report the newly obtained SMA observations at 850 micron on one of our sample source, A2390-5, behind the massive lensing cluster A2390. We successfully detect the source with a flux of 3.95 mJy. Surprisingly, it does not have any counterpart in any other wavelengths even though there are tentative candidates, which implies a very dusty and high-z nature. With less than 1" positional accuracy and the adoption of z = 5, we obtain the amplification factor of 12 using current lensing model, which makes A2390-5 a faint SMG with a de-lensed flux of 0.33 mJy. Together with our previous detection on another faint SMG, both of them have no counterpart in other wavelengths and their properties are very different than previously thought from the single-dish data. We emphasize the importance of direct submillimeter high-resolution studies on faint SMGs, which could be the dominant population of the high-z star formation., Comment: 4 pages, 2 figures, Conference Proceedings for Galaxy Mergers in an Evolving Universe (GMiEU) at Hualien, Taiwan, on 10/23-28, 2011; Website: http://events.asiaa.sinica.edu.tw/workshop/20111023/

In order to find the most extreme dust-hidden high-redshift galaxies, we select 196 extremely red objects in the Ks and IRAC bands (KIEROs, [Ks-4.5um](AB)>1.6) in the 0.06 deg^2 GOODS-N region. This selection avoids the Balmer breaks of galactic spectra at z<4 and picks up red galaxies with strong dust extinction. The photometric redshifts of KIEROs are between 1.5 and 5, with ~70% at z~2-4. KIEROs are very massive, with M*~10^10-10^12 Msun. They are optically faint and usually cannot be picked out by the Lyman break selection. On the other hand, the KIERO selection includes approximately half of the known millimeter and submillimeter galaxies in the GOODS-N. Stacking analyses in the radio, millimeter, and submillimeter all show that KIEROs are much more luminous than average 4.5um selected galaxies. Interestingly, the stacked fluxes for ACS-undetected KIEROs in these wavebands are 2.5-5 times larger than those for ACS-detected KIEROs. With the stacked radio fluxes and the local radio-FIR correlation, we derive mean infrared luminosities of 2-7x10^12 Lsun and mean star formation rates of 300-1200 Msun/yr for KIEROs with redshifts. We do not find evidence of a significant subpopulation of passive KIEROs. The large stellar masses and star formation rates imply that KIEROs are $z>2$ massive galaxies in rapid formation. Our results show that a large sample of dusty ultraluminous sources can be selected in this way and that a large fraction of high-redshift star formation is hidden by dust., Comment: Accepted for publication in ApJ. A typo in the abstract is corrected in this version

We report the first space-based very deep near-infrared continuum observations of a uniform sample of z = 6.5 galaxies with log(L(L_alpha)) = 42.5 - 43 erg/s selected from narrow-band line searches and with spectroscopically confirmed Ly alpha emission. The 1.4 micron HST WFC3 observations are deep enough (AB(1 sigma) = 28.75) to measure individual continuum magnitudes at this redshift for all of the objects. We compare the results with continuum-selected samples at the same redshift and find that Ly alpha emission is present in 24% of all galaxies with M_AB(1350 A) < -20 at z = 6.5. The error in this quantity is dominated by systematic uncertainties, which could be as large as multiplicative factors of three. The Ly alpha galaxies are extended but small (size < 1 kpc), and have star formation rates of approximately 10 Msun/yr. We find a mean L(Ly alpha)/(nu L_nu) at 1400A of 0.08, with the seven objects showing a range from 0.026 to 0.26, implying that there is little sign of destruction of the Ly alpha line. All of the properties of the z = 6.5 sample appear to be very similar to those of Ly alpha emitters at lower redshifts., Comment: Accepted for publication in the Astrophysical Journal Letters

We present optical spectroscopy of two samples of GALEX grism selected Ly alpha emitters (LAEs): one at z=0.195-0.44 and the other at z=0.65-1.25. We have also observed a comparison sample of galaxies in the same redshift intervals with the same UV magnitude distributions but with no detected Ly alpha. We use the optical spectroscopy to eliminate active galactic nuclei (AGNs) and to obtain the optical emission-line properties of the samples. We compare the luminosities of the LAEs in the two redshift intervals and show that there is dramatic evolution in the maximum Ly alpha luminosity over z=0-1. Focusing on the z=0.195-0.44 samples alone, we show that there are tightly defined relations between all of the galaxy parameters and the rest-frame equivalent width (EW) of H alpha. The higher EW(H alpha) sources all have lower metallicities, bluer colors, smaller sizes, and less extinction, consistent with their being in the early stages of the galaxy formation process. We find that 75 +- 12% of the LAEs have EW(H alpha)>100 Angstrom, and, conversely, that 31 +/- 3% of galaxies with EW(H alpha)>100 Angstrom are LAEs. We correct the broadband magnitudes for the emission-line contributions and use spectral synthesis fits to estimate the ages of the galaxies. We find a median age of 1.1x10^{8} yr for the LAE sample and 1.4x10^{9} yr for the UV-continuum sample without detected Ly alpha. The median metallicity of the LAE sample is 12+log(O/H)=8.24, or about 0.4 dex lower than the UV-continuum sample., Comment: to be published in the Astrophysical Journal

We report 850 micron Submillimeter Array (SMA) observations of four gravitationally lensed submillimeter galaxies (SMGs), A370-2, A2390-1, A2390-3 and A2390-4, which were originally discovered with the Submillimeter Common-User Bolometer Array (SCUBA). Our SMA detection of A370-2 with a submillimeter flux of 7.95 +/- 0.60 mJy unambiguously identifies the counterparts to this source at optical and radio wavelengths. A2390-1 is an ultraluminous infrared galaxies with a submillimeter flux of 5.55 +/- 0.92 mJy and a redshift of 1.8 +/- 0.2 computed from submillimeter/radio flux ratio analysis. We resolve A2390-3 into two components, A2390-3a and A2390-3b, with fluxes of 3.15 +/- 0.63 mJy and 1.92 +/- 0.60 mJy, respectively. The structure of the system could be consistent with morphological distortion by gravitational lensing. The lack of counterparts in the optical and infrared indicates a heavily dust-enshrouded nature, and a non-detection in the radio implies that these two sources probably lie at z > 4.7, which would make them among the most distant SMGs known to date. Our non-detection of A2390-4 suggests either that there are multiple fainter submillimeter sources within the SCUBA beam or that the SCUBA detection may have been false. Our precise positions allow us to determine accurate amplifications and fluxes for all our detected sources. Our new results give a shallower power-law fit (-1.10) to the faint-end 850 micron cumulative number counts than previous work. We emphasize the need for high-resolution observations of single-dish detected SMGs in order to measure accurately the faint end of the 850 micron counts., Comment: 8 pages, 6 figures, Accepted to ApJ

We carried out extremely sensitive Submillimeter Array (SMA) 340 GHz continuum imaging on two submillimeter galaxies (SMGs): GOODS 850-11 and GOODS 850-13. The observations reach sub-mJy rms sensitivities and, interestingly, resolve both sources into multiple, physically unrelated SMGs. GOODS 850-11 is resolved into two sources at different redshifts. GOODS 850-13 is resolved into three sources, two with different spectroscopic redshifts and one only with a photometric redshift. All the SMA sources have fluxes in the 3-5 mJy range and all are detected at 1.4 GHz. Three of them are detected by Chandra, and one is a previously unknown X-ray SMG. This is the first time that single-dish SMGs are resolved into multiple unrelated sources and also the first time that the SMA has discovered new SMGs. Our results show that identifications of SMGs at any wavelengths other than the submillimeter itself can be misleading, since such identifications usually only pick up one of the real counterparts. Using simulations that mimic our SCUBA and SMA observations, we find that the number of triple systems detected in our SMA survey is much higher than that expected from the current best-determined number counts. We tentatively attribute this to clustering. We also predict that ALMA will find ~1/3 of >5 mJy 850 um SCUBA sources to be multiple systems. Based on our SMA observations and simulations, we suggest that large samples of existing SMGs should be imaged with sensitive interferometric observations, even if the SMGs were previously thought to be securely identified., Comment: 6 pages, 1 table, and three figures. ApJL accepted

We present a current best estimate of the integrated near-infrared (NIR) extragalactic background light (EBL) attributable to resolved galaxies in J, H, and Ks. Our results in units of nW m-2 sr-1 are 11.7+5.6 -2.6 in J, 11.5+4.5 -1.5 in H and 10.0+2.8 -0.8 in Ks. We derive these new limits by combining our deep wide-field NIR photometry from five widely separated fields with other studies from the literature to create a galaxy counts sample that is highly complete and has good counting statistics out to JHKs ~ 27-28. As part of this effort we present new ultradeep Ks-band galaxy counts from 22 hours of observations with the Multi Object Infrared Camera and Spectrograph (MOIRCS) instrument on the Subaru Telescope. We use this MOIRCS Ks-band mosaic to estimate the total missing flux from sources beyond our detection limits. Our new limits to the NIR EBL are in basic agreement with, but 10 - 20% higher than previous estimates, bringing them into better agreement with estimates of the total NIR EBL (resolved + unresolved sources) obtained from TeV gamma-ray opacity measurements and recent direct measurements of the total NIR EBL. We examine field to field variations in our photometry to show that the integrated light from galaxies is isotropic to within uncertainties, consistent with the expected large-scale isotropy of the EBL. Our data also allow for a robust estimate of the NIR light from Galactic stars, which we find to be 14.7 +/- 2.4 in J, 10.1 +/- 1.9 in H and 7.6 +/- 1.8 in Ks in units of nW m-2 sr-1., Comment: Accepted to ApJ

We have measured the redshift evolution of the density of Lyman limit systems (LLS) in the intergalactic medium over the redshift range 0 < z < 6. We have used two new quasar samples to (1) improve coverage at z ~ 1, with GALEX grism spectrograph observations of 50 quasars with 0.8 < z_em < 1.3, and (2) extend coverage to z ~ 6, with Keck ESI spectra of 25 quasars with 4.17 < z_em < 5.99. Using these samples together with published data, we find that the number density of LLS per unit redshift, n(z), can be well fit by a simple evolution of the form n(z) = n_3.5 [(1+z)/4.5]^gamma, with n_3.5 = 2.80 +/- 0.33 and gamma = 1.94^(+0.36)_(-0.32) for the entire range 0 < z < 6. We have also reanalyzed the evolution of damped Lyman alpha systems (DLAs) in the redshift range 4 < z < 5 using our high-redshift quasar sample. We find a total of 17 DLAs and sub-DLAs, which we have analyzed in combination with published data. The DLAs with log (HI column density) > 20.3 show the same redshift evolution as the LLS. When combined with previous results, our DLA sample is also consistent with a constant Omega_DLA= 9 x 10^(-4) from z = 2 to z = 5. We have used the LLS number density evolution to compute the evolution in the mean free path of ionizing photons. We find a smooth evolution to z ~ 6, very similar in shape to that of Madau, Haardt & Rees (1999) but about a factor of two higher. Recent theoretical models roughly match to the z < 6 data but diverge from the measured power law at z > 6 in different ways, cautioning against extrapolating the fit to the mean free path outside the measured redshift range., Comment: To be published in ApJ. 30 pages, emulateapj style. Corrects typos in text and tables

We present an ultradeep Ks-band image that covers 0.5*0.5 deg^2 centered on the Great Observatories Origins Deep Survey-North (GOODS-N). The image reaches a 5 \sigma depth of Ks(AB) = 24.45 in the GOODS-N region, which is as deep as the GOODS-N Spitzer Infrared Array Camera (IRAC) 3.6 \mu m image. We present a new method of constructing IRAC catalogs that uses the higher spatial resolution Ks image and catalog as priors and iteratively subtracts fluxes from the IRAC images to estimate the IRAC fluxes. Our iterative method is different from the \chi^2 approach adopted by other groups. We verified our results using data taken in two different epochs of observations, as well as by comparing our colors with the colors of stars and with the colors derived from model spectral energy distributions (SEDs) of galaxies at various redshifts. We make available to the community our WIRCam Ks-band image and catalog (94951 objects in 0.25 deg^2), the Interactive Data Language (IDL) pipeline used for reducing the WIRCam images, and our IRAC 3.6 to 8.0 \mu m catalog (16950 objects in 0.06 deg^2 at 3.6 \mu m). With this improved Ks and IRAC catalog and a large spectroscopic sample from our previous work, we study the color-magnitude and color-color diagrams of galaxies. We compare the effectiveness of using Ks and IRAC colors to select active galactic nuclei (AGNs) and galaxies at various redshifts. We also study a color selection of z = 0.65--1.2 galaxies using the Ks, 3.6 \mu m, and 4.5 \mu m bands., Comment: Accepted for publication on ApJS. Online data are available

Ultra strong emission-line galaxies (USELs) with extremely high equivalent widths (EW(H beta) > 30A) can be used to pick out galaxies of extremely low metallicity in the z=0-1 redshift range. Large numbers of these objects are easily detected in deep narrow band searches and, since most have detectable [OIII] 4363, their metallicities determined using the direct method. These large samples hold the possibility for determining if there is a metallicity floor for the galaxy population. Here we describe results of an extensive spectroscopic follow-up of the Kakazu et al. (2007) catalog of 542 USELs using the DEIMOS spectrograph on Keck, with high S/N spectra of 348 galaxies. The two lowest metallicity galaxies in our sample have 12+log(O/H)=6.97+/-0.17 and 7.25+/-0.03 -- values comparable to the lowest metallicity galaxies found to date. We determine an empirical metallicity-R23 parameter relation for our sample, and compare this to the relationship for low redshift galaxies. The determined metallicity-luminosity relation is compared with those of magnitude selected samples in the same redshift range. The emission line selected galaxies show a metal-luminosity relation where the metallicity decreases with luminosity and they appear to define the lower bound of the galaxy metallicity distribution at a given continuum luminosity. We also compute the H alpha luminosity function of the USELs as a function of redshift and use this to compute an upper bound on the Ly alpha emitter luminosity function over the z=0-1 redshift range., Comment: 10 pages, 15 PostScript figures, uses emulateapj.sty, submitted to ApJ

GOODS 850-5 is a hyperluminous radio-faint submillimeter source in the GOODS-N. Although it is generally agreed that GOODS 850-5 is at a high redshift z>~4, its exact redshift is unknown. While its stellar SED suggests z~6, its radio/FIR SED suggests a lower redshift of z~4. To better constrain its stellar SED and redshift, we carried out nano-Jansky sensitivity ultradeep NIR observations between 1.2 and 2.1 um with the HST and the 8 m Subaru Telescope. Even with such great depths we did not detect GOODS 850-5, and the results show that it is an extremely curious source. Between the Ks and 3.6 um bands its spectral slope is >3x that of an ERO, and the flux ratio between the two bands is >8x that of Lyman breaks. It is quite challenging to explain this unusually red color without a Lyman break (which would imply z>17). It requires a large amount (M* ~ 10^11.5 Msun) of reddened old stars at z~6, coexisting with an even more extinguished violent ~2400-4400 Msun/yr starburst, which does not have any associated detectable rest-frame UV radiation. We discuss the discrepancy between the NIR and radio/FIR photometric redshifts. We conclude that GOODS 850-5 is at least at z>4 and is more likely at z>~6. We describe the unusual properties of GOODS 850-5, including its SED and formation history, and we discuss the implications of such massive z>6 galaxies., Comment: 12 pages, 9 figures, ApJ accepted

We present a comprehensive multiwavelength analysis of the bright, long duration gamma-ray burst GRB 070125, comprised of observations in $\gamma$-ray, X-ray, optical, millimeter and centimeter wavebands. Simultaneous fits to the optical and X-ray light curves favor a break on day 3.78, which we interpret as the jet break from a collimated outflow. Independent fits to optical and X-ray bands give similar results in the optical bands but shift the jet break to around day 10 in the X-ray light curve. We show that for the physical parameters derived for GRB 070125, inverse Compton scattering effects are important throughout the afterglow evolution. While inverse Compton scattering does not affect radio and optical bands, it may be a promising candidate to delay the jet break in the X-ray band. Radio light curves show rapid flux variations, which are interpreted as due to interstellar scintillation, and are used to derive an upper limit of $2.4 \times 10^{17}$ cm on the radius of the fireball in the lateral expansion phase of the jet. Radio light curves and spectra suggest a high synchrotron self absorption frequency indicative of the afterglow shock wave moving in a dense medium. Our broadband modeling favors a constant density profile for the circumburst medium over a wind-like profile ($R^{-2}$). However, keeping in mind the uncertainty of the parameters, it is difficult to unambiguously distinguish between the two density profiles. Our broadband fits suggest that \event is a burst with high radiative efficiency ($> 60 %$)., Comment: 50 pages, 33 figures, sty file included, Appeared in 20 Aug 2008 edition of Astrophysical Journal

We report an SMA interferometric identification of a bright submillimeter source, GOODS 850-5. This source is one of the brightest 850 um sources in the GOODS-N but is extremely faint at all other wavelengths. It is not detected in the GOODS HST ACS images and only shows a weak 2 sigma signal at 1.4 GHz. It is detected in the Spitzer IRAC bands and the MIPS 24 um band, however, with very low fluxes. We present evidence in the radio, submillimeter, mid-IR, near-IR, and optical that suggest GOODS 850-5 may be a z>4 galaxy., Comment: 4 pages, ApJ Letters accepted

We describe results of a narrow band search for ultra-strong emission line galaxies (USELs) with EW(H beta) > 30 A. 542 candidate galaxies are found in a half square degree survey using two ~100 Angstrom 8150 A and 9140 A filters with Subaru/SuprimeCam. Followup spectroscopy for randomly selected objects in the sample with KeckII/DEIMOS shows they consist of [OIII] 5007, [OII] 3727, and H alpha selected strong-emission line galaxies at intermediate redshifts (z < 1), and Ly alpha emitting galaxies at high-redshift (z >> 5). We determine the H beta luminosity functions and the star formation density of the USELs, which is 5-10% of the value found from ultraviolet continuum objects at z=0-1, suggesting they correspond to a major epoch in galaxy formation at these redshifts. Many USELs show the temperature-sensitive [OIII] 4363 auroral lines and about a dozen have oxygen abundances characteristic of eXtremely Metal Poor Galaxies (XMPGs). These XMPGs are the most distant known today. Our high yield rate of XMPGs suggests this is a powerful method for finding such populations. The lowest metallicity measured in our sample is 12+log(O/H) = 7.06 (6.78-7.44), close to the minimum metallicity found in local galaxies. The luminosities, metallicities and star formation rates of USELs are consistent with the strong emitters being start-up intermediate mass galaxies and suggest that galaxies are still forming in relatively chemically pristine sites at z < 1., Comment: 18 pages, 20 figures (with 2 color plates), aastex, emulateapj,submitted to the Astrophysical Journal

We present a study of three Lyman-alpha emitting galaxies (LAEs), selected via a narrow-band survey in the GOODS northern field, and spectroscopically confirmed to have redshifts of z ~ 5.65. Using HST ACS and Spitzer IRAC data, we constrain the rest-frame UV-to-optical spectral energy distributions (SEDs) of the galaxies. Fitting stellar population synthesis models to the observed SEDs, we find best-fit stellar populations with masses between ~ 10^9 - 10^10 M_sun and ages between ~ 5 - 100 Myr, assuming a simple starburst star formation history. However, stellar populations as old as 700 Myr are admissible if a constant star formation rate model is considered. Very deep near-IR observations may help to narrow the range of allowed models by providing extra constraints on the rest-frame UV spectral slope. Our narrow-band selected objects and other IRAC-detected z ~ 6 i'-dropout galaxies have similar 3.6 um magnitudes and z' - [3.6] colors, suggesting that they posses stellar populations of similar masses and ages. This similarity may be the result of a selection bias, since the IRAC-detected LAEs and i'-dropouts probably only sample the bright end of the luminosity function. On the other hand, our LAEs have blue i' - z' colors compared to the i'-dropouts, and would have been missed by the i'-dropout selection criterion. A better understanding of the overlap between the LAE and the i'-dropout populations is necessary in order to constrain the properties of the overall high-redshift galaxy population, such as the total stellar mass density at z ~ 6., Comment: 10 pages, 8 figures. Accepted for publication in ApJ

The recent development of large, complete samples which identify high-redshift galaxies at z~5.7 and z~6.5 from deep, wide-field surveys provides detailed information on the earliest galaxies, their numbers, spatial and kinematic distributions, and implications for early reionization of the IGM. In this contribution we present results of spectroscopic studies of z~5.7 and z~6.5 galaxies identified from our deep, Lyman alpha narrowband and multicolor surveys conducted with the SuprimeCam mosaic CCD camera on the 8.3-m Subaru telescope and observed with the DEIMOS multi-object spectrograph on Keck. The luminosity function of the z~6.5 galaxies is shown to be similar to the luminosity function of the z~5.7 galaxy samples, suggesting that a substantial star-forming population is already in place at z~6.5. Comparisons of both individual and stacked spectra of galaxies in these two samples show that the Lyman alpha emission profiles, equivalent widths, and continuum break strengths do not substantially change over this redshift interval. The wide-field nature of the surveys also permits mapping the large-scale distribution of the high-redshift galaxies in spatial structures extending across individual SuprimeCam fields (~60 Mpc). Field-to-field variations in the number of objects at z~6.5 may shortly be able to place constraints on the porosity of the reionization boundary., Comment: 6 pages, 6 figures, to appear in IAU 199 Conf. Proc.: "Probing Galaxies through Quasar Absorption Lines," eds. Williams, Shu, Menard

We report the results of an extensive imaging and spectroscopic survey in the GOODS-North field completed using DEIMOS on the Keck II telescope. Observations of 2018 targets in a magnitude-limited sample of 2911 objects to R=24.4 yield secure redshifts for a sample of 1440 galaxies and AGN plus 96 stars. In addition to redshifts and associated quality assessments, our catalog also includes photometric and astrometric measurements for all targets detected in our R-band imaging survey of the GOODS-North region. We investigate various sources of incompleteness and find the redshift catalog to be 53% complete at its limiting magnitude. The median redshift of z=0.65 is lower than in similar deep surveys because we did not select against low-redshift targets. Comparison with other redshift surveys in the same field, including a complementary Hawaii-led DEIMOS survey, establishes that our velocity uncertainties are as low as 40 km/s for red galaxies and that our redshift confidence assessments are accurate. The distributions of rest-frame magnitudes and colors among the sample agree well with model predictions out to and beyond z=1. We will release all survey data, including extracted 1-D and sky-subtracted 2-D spectra, thus providing a sizable and homogeneous database for the GOODS-North field which will enable studies of large scale structure, spectral indices, internal galaxy kinematics, and the predictive capabilities of photometric redshifts., Comment: 17 pages, 18 figures, submitted to AJ; v2 minor changes; see survey database at http://www2.keck.hawaii.edu/realpublic/science/tksurvey/

We discuss the nature of faint 6.7um galaxies detected with the mid-infrared camera ISOCAM on board the Infrared Space Observatory (ISO). The 23 hour integration on the Hawaii Deep Field SSA13 has provided a sample of 65 sources down to 6uJy at 6.7um. For 57 sources, optical or near-infrared counterparts were found with a statistical method. All four Chandra sources, three SCUBA sources, and one VLA/FIRST source in this field were detected at 6.7um with high significance. Using their optical to mid-infrared colors, we divided the 6.7um sample into three categories: low redshift galaxies with past histories of rapid star formation, high redshift ancestors of these, and other star forming galaxies. Rapidly star forming systems at high redshifts dominate the faintest end. Spectroscopically calibrated photometric redshifts were derived from fits to a limited set of template SEDs. They show a high redshift tail in their distribution with faint (<30uJy) galaxies at z>1. The 6.7um galaxies tend to have brighter K magnitudes and redder I-K colors than the blue dwarf population at intermediate redshifts. Stellar masses of the 6.7um galaxies were estimated from their rest-frame near-infrared luminosities. Massive galaxies (M_star~10e11M_sun) were found in the redshift range of z=0.2-3. Epoch dependent stellar mass functions indicate a decline of massive galaxies' comoving space densities with redshift. Even with such a decrease, the contributions of the 6.7um galaxies to the stellar mass density in the universe are found to be comparable to those expected from UV bright galaxies detected in deep optical surveys., Comment: 31 pages, 15 figures, AJ (accepted), a version with color figures at http://www.ioa.s.u-tokyo.ac.jp/~ysato/pub/3/p3c-ysato.ps.gz

We report results of a deep wide-field narrowband survey for redshift z~5.7 Ly alpha emitters carried out with SuprimeCam on Subaru 8.3-m telescope. Deep narrowband imaging of the SSA22 field through a 120 A bandpass filter centered at 8150 A was combined with deep multicolor RIz' SuprimeCam broadband imaging, and BVRZ imaging taken with CFHT's CFH12K camera to select high-redshift galaxy candidates. Spectroscopic observations were made using the new wide-field multi-object DEIMOS spectrograph on Keck for 22 of the 26 candidate objects. Eighteen objects were identified as z~5.7 Lyman alpha emitters, and a further nineteenth candidate was identified based on an LRIS spectrum. At the 3.3 A resolution of the DEIMOS spectra the asymmetric profile for Ly alpha emission with its steep blue fall-off can be clearly seen. We use this to describe the distribution of equivalent widths and the continuum color break properties for z~5.7 Ly alpha galaxies compared with foreground objects. The large majority (>75%) of Ly alpha lines have rest frame equivalent widths less than 240 A and can be understood in terms of young star forming galaxies with a Salpeter initial mass function for the stars. With narrowband selection criteria of I-N > 0.7 and N<25.05 (AB mags) we find a surface density of Ly alpha emitters of 0.03 per square arcminute per (deltaz=0.1) to a limiting flux just under 2e-17 erg/cm2/s. The luminosity function of the Ly alpha emitters is similar to that at lower redshifts to the lowest measurable luminosity of 1e43 ergs/s as is the universal star formation rate based on their continuum properties. We note that the objects are highly structured in both spatial and spectral properties on the angular scale of the fields (~60 Mpc), and that multiple fields will have to be averaged to accurately measure their ensemble properties., Comment: 12 pages, 20 figures, emulateapj, to appear in the Astronomical Journal (February 2004), version with higher quality figures available at http://www.ifa.hawaii.edu/publications/preprints/Hu_03-95.pdf Fixed problem with cited papers appearing as numbered hypertex links instead of text (only affected astro-ph pdf file)

We investigate the evolution of the universal rest-frame ultraviolet luminosity density from z = 1.5 to the present. We analyze an extensive sample of multicolor data (U', B, V = 24.5) plus spectroscopic redshifts from the Hawaii Survey Fields and the Hubble Deep Field. Our multicolor data allow us to select our sample in the rest-frame ultraviolet (2500 angstrom) over the entire redshift range to z = 1.5. We conclude that the evolution in the luminosity density is a function of the form (1+z)^{1.7\pm1.0} for a flat lambda cosmology and (1+z)^{2.4\pm1.0} for an Einstein-de Sitter cosmology., Comment: 12 pages, 7 figs, 5 tables, submitted to AJ

We present measurements of the extended dark halo profiles of bright early type galaxies at redshifts 0.1 to 0.9 obtained via galaxy-galaxy lensing analysis of images taken at the CFHT using the UH8K CCD mosaic camera. Six half degree fields were observed for a total of 2 hours each in I and V, resulting in catalogs containing ~20 000 galaxies per field. We used V-I color and I magnitude to select bright early type galaxies as the lens galaxies, yielding a sample of massive lenses with fairly well determined redshifts and absolute magnitudes M ~ M_* pm 1. We paired these with faint galaxies lying at angular distances 20" to 60", corresponding to physical radii of 26 to 77 kpc (z = 0.1) and 105 to 315 kpc (z = 0.9), and computed the mean tangential shear of the faint galaxies. The shear falls off with radius roughly as expected for flat rotation curve halos. The shear values were weighted in proportion to the square root of the luminosity of the lens galaxy. Our results give a value for the average mean rotation velocity of an L_* galaxy halo at r~50-200 kpc of v_* = 238^{+27}_{-30} km per sec for a flat lambda (Omega_m0 = 0.3, Omega_l0 = 0.7) cosmology (v_* = 269^{+34}_{-39} km per sec for Einstein-de Sitter), and with little evidence for evolution with redshift. We find a mass-to-light ratio of \~121\pm28h(r/100 kpc) and these halos constitute Omega ~0.04 pm 0.01 (r/100 kpc) of closure density., Comment: 4 pages, to appear in Proceedings of the Workshop "The Shapes of Galaxies and their Halos", Yale, May 2001

We report the discovery of a redshift z=6.56 galaxy lying behind the cluster Abell 370. The object HCM 6A was found in a narrowband imaging survey using a 118 Angstrom bandpass filter centered at 9152 Ang in the LRIS camera on the 10 m Keck II Telescope. Candidate Lyman alpha emitters were identified by the equivalent width of the emission and the absence of lower wavelength flux in ultradeep broadband images. HCM6A is the first galaxy to be confirmed at redshift z>6, and has W_lambda(observed)=190 Angstroms, flux = 2.7 x 10^-17 erg cm^-2 s^-1. Spectra obtained with LRIS confirm the emission line and the continuum break across the line, and show an asymmetric line profile with steep fall-off on the blue side. Deep Subaru near-infrared CISCO images in J, H and K' which extend the sampled continuum to longer wavelengths give a consistent estimate of the continuum flux density in these line-free regions of 2.6+/-0.7 x 10^-30 erg cm^-2 s^-1 Hz^-1. The line width and strength, asymmetric profile, and very deep spectral break are only consistent with the interpretation of the line as a redshifted Lyman alpha feature. From the detailed lensing model of this cluster, we estimate a lensing amplification of 4.5 for this galaxy, which is located slightly over an arcminute from the center of the cluster, for an unlensed flux of 6.5 x 10^-18 erg cm^-2 s^-1. The presence of such a galaxy suggests that the reionizing epoch is beyond z=6.6., Comment: 5 pages, 4 PostScript figures, uses emulateapj5, to appear in the April 1, 2002 issue of the Astrophysical Journal Letters

We analyze the Ly alpha forest properties of a sample of 15 high redshift quasars lying between z = 4.42 and z = 5.75, using high signal-to-noise spectra obtained with ESI on the Keck II 10 m telescope. The distribution of transmissions in the Ly alpha region in this redshift range is shown to be consistent with that found in lambda cold dark matter simulations with a smoothly evolving ionization rate as a function of redshift. The extrapolation of the ionizing flux to z = 6.05 lies a factor of two higher than a 2 sigma upper limit placed by Cen & McDonald (2001) at this redshift, based on the Becker et al. (2001) spectra of the z = 6.28 quasar SDSS 1030+0524. However, the data are also consistent with models in which there is substantial variation of the ionization parameter about the mean value, and in this case, dark gaps such as those seen by Becker et al. become much more likely. We conclude that further data are needed before we can be sure that we have seen the epoch of reionization. We also summarize the damped Ly alpha systems seen in these quasar lines of sight and measure their metallicities and kinematic properties. We argue that the mean DLA metallicity has dropped substantially by z = 5 compared with its value at z < 4., Comment: 17 pages, including 13 figures. Accepted for publication in Astronomical Journal (May 2002)

We present measurements of the extended dark halo profiles of bright early type galaxies at redshifts 0.1 to 0.9 obtained via galaxy-galaxy lensing analysis of images taken at the CFHT using the UH8K CCD mosaic camera. Six half degree fields were observed for a total of 2 hours each in I and V, resulting in catalogs containing ~20 000 galaxies per field. We used V-I color and I magnitude to select bright early type galaxies as the lens galaxies, yielding a sample of massive lenses with fairly well determined redshifts and absolute magnitudes M ~ M_* \pm 1. We paired these with faint galaxies lying at angular distances 20" to 60", corresponding to physical radii of 26 to 77 kpc (z = 0.1) and 105 to 315 kpc (z = 0.9), and computed the mean tangential shear of the faint galaxies. The shear falls off with radius roughly as expected for flat rotation curve halos. The shear values were weighted in proportion to the square root of the luminosity of the lens galaxy. Our results give a value for the average mean rotation velocity of an L_* galaxy halo at r~50-200 kpc of v_* = 238^{+27}_{-30} km per sec for a flat lambda (Omega_m0 = 0.3, Omega_l0 = 0.7) cosmology (v_* = 269^{+34}_{-39} km per sec for Einstein-de Sitter), and with little evidence for evolution with redshift. We compare to halo masses measured by other groups/techniques. We find a mass-to-light ratio of ~121\pm28h(r/100 kpc) and these halos constitute Omega ~0.04 \pm 0.01(r/100 kpc) of closure density. (abridged), Comment: Accepted for publication in ApJ (minor modifications) - 32 pages, 11 figs, 5 tables

The rate of mass accumulation due to galaxy merging depends on the mass, density, and velocity distribution of galaxies in the near neighborhood of a host galaxy. The fractional luminosity in kinematic pairs combines all of these effects in a single estimator which is relatively insensitive to population evolution. Here we use a k-corrected and evolution compensated volume-limited sample drawing about 300 redshifts from CFGRS and 3000 from CNOC2 to measure the rate and redshift evolution of merging. We identify kinematic pairs with projected separations less than either 50 or 100 \hkpc and rest-frame velocity differences of less than 1000\kms. The fractional luminosity in pairs is modeled as f_L(Delta v,r_p,M_r^{ke})(1+z)^{m_L} where [f_L,m_L] are [0.14+/-0.07,0+/-1.4] and [0.37+/-0.7,0.1+/-0.5] for r_p<= 50 and 100\hkpc, respectively (Omega_M=0.2, Omega_Lambda=0.8). The value of m_L is about 0.6 larger if Lambda=0. To convert these redshift space statistics to a merger rate we use the data to derive a conversion factor to physical space pair density, a merger probability and a mean in-spiral time. The resulting mass accretion rate per galaxy (M_1,M_2>= 0.2 M*) is 0.02+/-0.01(1+z)^{0.1+/-0.5} M*~Gyr^{-1}. Present day high-luminosity galaxies therefore have accreted approximately 0.15M* of their mass over the approximately 7 Gyr to redshift one. (abridged), Comment: accepted for publication in ApJ Letters

We review the properties of z>5 galaxies studied with HST and with the Keck telescopes, and discuss the detectability of Lyman alpha emission-line galaxies out to z~6.5 based on these data and ongoing narrowband imaging surveys. The brightest sources may show (R-Z) color breaks, although the high sky background at Z (lambda_{eff} ~ 9200 Ang), makes such observations challenging for typical faint sources. Keck LRIS observations of the z=5 SDSS quasar and z>5 galaxies observed with HST in the HDF show that the strength of the Lyman break is evolving more slowly than extrapolations from models at z~3 would predict., Comment: 12 pages, 7 encapsulated PostScript figures, paspconf.sty. Invited review, to appear in "The Hy-Redshift Universe: Galaxy Formation and Evolution at High Redshift" held in Berkeley, June 21-24, 1999, eds. A.J. Bunker and W.J.M. van Breugel, A.S.P. Conf.Ser. in press (1999)

We report new limits on the Gunn-Peterson effect at a redshift near 5 using spectroscopic observations of the z = 5 Sloan Digital Sky Survey quasar, J033829.31+002156.3, made with the LRIS and HIRES spectrographs on the Keck telescopes. Lower resolution spectrophotometrically calibrated observations made with LRIS over the wavelength region 4500-9600 angstroms were used to obtain a continuum shape and to flux calibrate much higher resolution (R = 36,000) observations made with HIRES. The LRIS data show an Oke D_A index of 0.75. Portions of the HIRES spectrum return to near the extrapolated continuum level. Including both statistical and systematic errors we place an upper limit of tau = 0.1 on the regions of minimum opacity. We argue that, even if this opacity arises in underdense regions of the intergalactic gas, we require a high value of the metagalactic ionizing flux at these redshifts ($J_{\nu} \gg 4 \times 10^{-23}~{\rm ergs cm^{-2} s^{-1} Hz^{-1} sr^{-1}}$ at $z\sim4.72$) to produce a solution which is consistent with even minimum nucleosynthesis estimates of the baryon density. We also report the presence of a Mg II absorption system of extremely high equivalent width ($W_{\lambda,rest}(2796)=1.73$ \AA) at z=2.304., Comment: 15 pages including 2 figures. Also available at http://www.ifa.hawaii.edu/~acowie/gp_5.html To be published in Astrophysical Journal Letters