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.

We search for high-redshift (z>4.5) X-ray AGNs in the deep central (off-axis angle <5.7') region of the 7 Ms Chandra Deep Field-South X-ray image. We compile an initial candidate sample from direct X-ray detections. We then probe more deeply in the X-ray data by using pre-selected samples with high spatial resolution NIR/MIR (HST 1.6 micron and Spitzer 4.5 micron) and submillimeter (ALMA 850 micron) observations. The combination of the NIR/MIR and submillimeter pre-selections allows us to find X-ray sources with a wide range of dust properties and spectral energy distributions (SEDs). We use the SEDs from the optical to the submillimeter to determine if previous photometric redshifts were plausible. Only five possible z>5 X-ray AGNs are found, all of which might also lie at lower redshifts. If they do lie at high redshifts, then two are Compton-thick AGNs, and three are ALMA 850 micron sources. We find that (i) the number density of X-ray AGNs is dropping rapidly at high redshifts, (ii) the detected AGNs do not contribute significantly to the photoionization at z>5, and (iii) the measured X-ray light density over z=5-10 implies a very low black hole accretion density with very little growth in the black hole mass density in this redshift range.

The CDF-S is the deepest X-ray image available and will remain so for the near future. We provide a spectroscopic (64.5%; 64% with spectral classifications) and photometric redshift catalog for the full 7 Ms sample, but much of our analysis focuses on the central (off-axis angles <5.7') region, which contains a large, faint ALMA sample of 75 >4.5-sigma 850 micron sources. We measure the 850 micron fluxes at the X-ray positions using the ALMA images, where available, or an ultradeep SCUBA-2 map. We find that the full X-ray sample produces ~10% of the 850 micron extragalactic background light. We separate the submillimeter detected X-ray sources into star-forming galaxies and AGNs using a star formation rate (SFR) versus X-ray luminosity calibration for high SFR galaxies. We confirm this separation using the X-ray photon indices. We measure the X-ray fluxes at the accurate positions of the 75 ALMA sources and detect 70% at >3-sigma in either the 0.5-2 or 2-7 keV bands. However, many of these may produce both their X-ray and submillimeter emission by star formation. Indeed, we find that only 20% of the ALMA sources have intermediate X-ray luminosities (rest-frame 8-28 keV luminosities of 10^42.5-10^44 erg/s), and none has a high X-ray luminosity (>10^44 erg/s). Conversely, after combining the CDF-S with the CDF-N, we find extreme star formation (SFR>300 solar masses per yr) in some intermediate X-ray luminosity sources but not in any high X-ray luminosity sources. We argue that the quenching of star formation in the most luminous AGNs may be a consequence of the clearing of gas in these sources.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

Galaxies grow their supermassive black holes in concert with their stars, although the relationship between these major galactic components is poorly understood. Observations of the cosmic growth of stars and black holes in galaxies suffer from disjoint samples and the strong effects of dust attenuation. The thermal infrared holds incredible potential for simultaneously measuring both the star formation and black hole accretion rates in large samples of galaxies covering a wide range of physical conditions. Spitzer demonstrated this potential at low redshift, and by observing some of the most luminous galaxies at z~2. JWST will apply these methods to normal galaxies at these epochs, but will not be able to generate large spectroscopic samples or access the thermal infrared at high-redshift. An order of magnitude gap in our wavelength coverage will persist between JWST and ALMA. A large, cold infrared telescope can fill this gap to determine when (in cosmic time), and where (within the cosmic web), stars and black holes co-evolve, by measuring these processes simultaneously in statistically complete and unbiased samples of galaxies to z>8. A next-generation radio interferometer will have the resolution and sensitivity to measure star-formation and nuclear accretion in even the dustiest galaxies. Together, the thermal infrared and radio can uniquely determine how stars and supermassive blackholes co-evolve in galaxies over cosmic time.

Deep continuum surveys at radio wavelengths are able to cover large areas, yield high angular resolution, and do not suffer from dust extinction, thus providing a robust way to measure the star formation history of the universe. However, with the current sensitivities of existing radio telescopes, it remains challenging to detect galaxies that dominate the cosmic star formation history even with extremely long integrations. With the ngVLA, a new portion of parameter space will be opened up for radio continuum surveys: deep (~200 nJy/bm), large-area (~1 deg^2), sub-arcsecond surveys at high frequencies (~8 GHz), where the observed radio emission from high-z galaxies should be dominated by free-free emission, providing a robust measurement for the star formation history of the universe. By being able to image the star formation activity, unbiased by dust, for a large, homogeneous sample of galaxies with a wide range of luminosities into the epoch of reionization, such surveys with the ngVLA will be highly complementary to those conducted by JWST, which will only be sensitive to un-obscured star formation.

We use deep F275W imaging from the $Hubble$ Deep UV Legacy Survey (HDUV) and G280 grism spectroscopy from $HST$/WFC3, along with new and archival optical spectra from Keck/DEIMOS, to search for candidate ionizing sources in the GOODS-N field at z ~ 2.5 - 3. Spectroscopic identification of our UV-selected sources are 99% complete to F275W = 25.5 in the region of the UV imaging, and we identify 6 potential ionizing galaxies or AGNs at z ~ 3. By far the brightest of these is a z = 2.583 AGN that totally dominates the ionizing flux in the region, with a specific ionizing volume emissivity at $912~ \AA$ of $\epsilon_{912}$ = $8.3^{27}_{1.4} \times 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$. Based on our spectroscopic data, we find four candidates are contaminated by foreground galaxies at z ~ 0.5 - 0.7. At $\epsilon_{912}$ = $2.2^{7.2}_{0.4} \times 10^{23}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$, the remaining candidate galaxy's contribution to the ionizing background lies well below the flux required to ionize the intergalactic medium at z ~ 2.5 - 3, consistent with previous observations that show AGNs provide the bulk of the ionizing background at these redshifts.

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.

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.

We present a catalog of 10718 objects in the COSMOS field observed through multi-slit spectroscopy with the Deep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck II telescope in the wavelength range ~5500-9800A. The catalog contains 6617 objects with high-quality spectra (two or more spectral features), and 1798 objects with a single spectroscopic feature confirmed by the photometric redshift. For 2024 typically faint objects we could not obtain reliable redshifts. The objects have been selected from a variety of input catalogs based on multi-wavelength observations in the field, and thus have a diverse selection function, which enables the study of the diversity in the galaxy population. The magnitude distribution of our objects is peaked at I_AB~23 and K_AB~21, with a secondary peak at K_AB~24. We sample a broad redshift distribution in the range 0<z<6, with one peak at z~1, and another one around z~4. We have identified 13 redshift spikes at z>0.65 with chance probabilities <4xE-4$, some of which are clearly related to protocluster structures of sizes >10 Mpc. An object-to-object comparison with a multitude of other spectroscopic samples in the same field shows that our DEIMOS sample is among the best in terms of fraction of spectroscopic failures and relative redshift accuracy. We have determined the fraction of spectroscopic blends to about 0.8% in our sample. This is likely a lower limit and at any rate well below the most pessimistic expectations. Interestingly, we find evidence for strong lensing of Ly-alpha background emitters within the slits of 12 of our target galaxies, increasing their apparent density by about a factor of 4.

There is substantial and growing observational evidence from the normalized luminosity density in the near-infrared that the local universe is under-dense on scales of several hundred Megaparsecs. We test whether our parameterization of the observational data of such a "void" is compatible with the latest supernovae type Ia data and with constraints from line-of-sight peculiar velocity motions of galaxy clusters with respect to the cosmic microwave background rest frame, known as the linear kinematic Sunyaev-Zel'dovich (kSZ) effect. Our study is based on the large local void (LLV) radial profile observed by Keenan, Barger, and Cowie (KBC) and a theoretical void description based on the Lema\^itre-Tolman-Bondi model with a nonzero cosmological constant ($\Lambda$LTB). We find consistency with the measured luminosity distance-redshift relation on radial scales relevant to the KBC LLV through a comparison with 217 low-redshift supernovae type Ia over the redshift range $0.0233 < z < 0.15$. We assess the implications of the KBC LLV in light of the tension between "local" and "cosmic" measurements of the Hubble constant, $H_{0}$. We find that when the existence of the KBC LLV is fully accounted for, this tension is reduced from $3.4\sigma$ to $2.75\sigma$. We find that previous linear kSZ constraints, as well as new ones from the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT), are fully compatible with the existence of the KBC LLV.

In this third paper of the Hawaii SCUBA-2 Lensing Cluster Survey series, we present Submillimeter Array (SMA) detections of six intrinsically faint 850 $\mu$m sources detected in SCUBA-2 images of the lensing cluster fields, A1689, A2390, A370, MACS J0717.5+3745, and MACS J1423.8+2404. Two of the SCUBA-2 sources split into doublets, yielding a total of eight SMA detections. The intrinsic 870 $\mu$m flux densities of these submillimeter galaxies (SMGs) are $\sim$ 1 mJy. Five of the eight SMGs are not detected in optical or near-infrared (NIR) images. The NIR-to-submillimeter flux ratios of these faint SMGs suggest that most of them are extremely dusty and/or at very high redshifts. Combining these SMGs and several other samples from the literature, we find a bimodal distribution for the faint sources in the space of submillimeter flux versus NIR-to-submillimeter flux ratio. While most of the SMA-detected lensed sources are very obscured, the other SMGs with similar flux densities are mostly bright in the NIR. Future ALMA observations of a large sample of SCUBA-2 sources in cluster fields will allow us to decide whether the bimodality we observe here really exists.

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

In this second paper of the Hawaii SCUBA-2 Lensing Cluster Survey series, we cross-match SCUBA-2 maps with 3 and 6 GHz images from the Janksy-VLA Frontier Fields Legacy Survey for three cluster fields, MACS J0416.1--2403, MACS J0717.5+3745, and MACS J1149.5+2223. Within the HST coverage, 14 out of 44 850 $\mu$m sources have 3 GHz counterparts, five of which are also detected at 6 GHz. The 850 $\mu$m flux densities of these detected sources span from 0.7 to 4.4 mJy after correcting for lensing amplification. The median redshift of the sample is $z = 1.28^{+0.07}_{-0.09}$, much lower than the typical redshifts ($z = 2-3$) of brighter submillimeter galaxies in the literature. In addition, we find that our sources have lower dust temperatures than those of the brighter submillimeter galaxies. This is also confirmed by an analysis of the ratio between infrared star formation rate and 850 $\mu$m flux density. However, these 14 sources may not represent the general submillimeter population at the same flux range, given that the SCUBA-2 sources without radio counterparts are likely at higher redshifts. Detection of these sources would require deeper radio images or submillimeter interferometry.

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.

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.

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.

We present deep number counts at 450 and 850 $\mu$m using the SCUBA-2 camera on the James Clerk Maxwell Telescope. We combine data for six lensing cluster fields and three blank fields to measure the counts over a wide flux range at each wavelength. Thanks to the lensing magnification, our measurements extend to fluxes fainter than 1 mJy and 0.2 mJy at 450 $\mu$m and 850 $\mu$m, respectively. Our combined data highly constrain the faint end of the number counts. Integrating our counts shows that the majority of the extragalactic background light (EBL) at each wavelength is contributed by faint sources with $L_{\rm IR} < 10^{12} L_{\odot }$, corresponding to luminous infrared galaxies (LIRGs) or normal galaxies. By comparing our result with the 500 $\mu$m stacking of $K$-selected sources from the literature, we conclude that the $K$-selected LIRGs and normal galaxies still cannot fully account for the EBL that originates from sources with $L_{\rm IR} < 10^{12} L_{\odot }$. This suggests that many faint submillimeter galaxies may not be included in the UV star formation history. We also explore the submillimeter flux ratio between the two bands for our 450 $\mu$m and 850 $\mu$m selected sources. At 850 $\mu$m, we find a clear relation between the flux ratio and the observed flux. This relation can be explained by a redshift evolution, where galaxies at higher redshifts have higher luminosities and star formation rates. In contrast, at 450 $\mu$m, we do not see a clear relation between the flux ratio and the observed flux.

We report the detection of the most luminous high-redshift Lyman Alpha Emitting galaxy (LAE) yet seen, with log L(Ly alpha) = 43.9 ergs/s. The galaxy -- COSMOS Lyman alpha 1, or COLA1 -- was detected in a search for ultra-luminous LAEs with Hyper Suprime-Cam on the Subaru telescope. It was confirmed to lie at z = 6.593 based on a Lyman alpha line detection obtained from followup spectroscopy with the DEIMOS spectrograph on Keck2. COLA1 is the first very high-redshift LAE to show a multi-component Lyman alpha line profile with a blue wing, which suggests that it could lie in a highly ionized region of the intergalactic medium and could have significant infall. If this interpretation is correct, then ultra-luminous LAEs like COLA1 offer a unique opportunity to determine the properties of the HII regions around these galaxies which will help in understanding the ionization of the z ~ 7 intergalactic medium.

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.

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.

The Next-Generation Very Large Array (ngVLA) will be critical for understanding how galaxies are built and evolve at the earliest epochs. The sensitivity and frequency coverage will allow for the detection of cold gas and dust in `normal' distant galaxies, including the low-J transitions of molecular gas tracers such as CO, HNC, and HCO+; synchrotron and free-free continuum emission; and even the exciting possibility of thermal dust emission at the highest (z~7) redshifts. In particular, by enabling the total molecular gas reservoirs to be traced to unprecedented sensitivities across a huge range of epochs simultaneously -- something no other radio or submillimeter facility will be capable of -- the detection of the crucial low-J transitions of CO in a diverse body of galaxies will be the cornerstone of ngVLA's contribution to high-redshift galaxy evolution science. The ultra-wide bandwidths will allow a complete sampling of radio SEDs, as well as the detection of emission lines necessary for spectroscopic confirmation of elusive dusty starbursts. The ngVLA will also deliver unique contributions to our understanding of cosmic magnetism and to science accessible through microwave polarimetry. Finally, the superb angular resolution will move the field beyond detection experiments and allow detailed studies of the morphology and dynamics of these systems, including dynamical modeling of disks/mergers, determining the properties of outflows, measuring black hole masses from gas disks, and resolving multiple AGN nuclei. We explore the contribution of a ngVLA to these areas and more, as well as synergies with current and upcoming facilities including ALMA, SKA, large single-dish submillimeter observatories, GMT/TMT, and JWST.

We use ultradeep SCUBA-2 850um observations (~0.37 mJy rms) of the 2 Ms Chandra Deep Field-North (CDF-N) and 4 Ms Chandra Deep Field-South (CDF-S) X-ray fields to examine the amount of dusty star formation taking place in the host galaxies of high-redshift X-ray AGNs. Supplementing with COSMOS, we measure the submillimeter fluxes of the 4-8 keV sources at z>1, finding little flux at the highest X-ray luminosities but significant flux at intermediate luminosities. We determine gray body and MIR luminosities by fitting spectral energy distributions to each X-ray source and to each radio source in an ultradeep Karl G. Jansky Very Large Array (VLA) 1.4 GHz (11.5uJy at 5-sigma) image of the CDF-N. We confirm the FIR-radio and MIR-radio correlations to z=4 using the non-X-ray detected radio sources. Both correlations are also obeyed by the X-ray less luminous AGNs but not by the X-ray quasars. We interpret the low FIR luminosities relative to the MIR for the X-ray quasars as being due to a lack of star formation, while the MIR stays high due to the AGN contribution. We find that the FIR luminosity distributions are highly skewed and the means are dominated by a small number of high-luminosity galaxies. Thus, stacking or averaging analyses will overestimate the level of star formation taking place in the bulk of the X-ray sample. We conclude that most of the host galaxies of X-ray quasars are not strong star formers, perhaps because their star formation is suppressed by AGN feedback.

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.

Far-Infrared (FIR) photometry from the the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory is presented for 313 nearby, hard X-ray selected galaxies from the 58-month Swift Burst Alert Telescope (BAT) Active Galactic catalog. The present data do not distinguish between the FIR luminosity distributions at 70 and 160um for Seyfert 1 and Seyfert 2 galaxies. This result suggests that if the FIR emission is from the nuclear obscuring material surrounding the accretion disk, then it emits isotropically, independent of orientation. Alternatively, a significant fraction of the 70 and 160um could be from star formation, independent of AGN type. Using a non-parametric test for partial correlation with censored data, we find a statistically significant correlation between the AGN intrinsic power (in the 14-195 keV band ) and the FIR emission at 70 and 160um for Seyfert 1 galaxies. We find no correlation between the 14-195 keV and FIR luminosities in Seyfert 2 galaxies. The observed correlations suggest two possible scenarios: (i) if we assume that the FIR luminosity is a good tracer of star formation, then there is a connection between star formation and the AGN at sub-kiloparsec scales, or (ii) dust heated by the AGN has a statistically significant contribution to the FIR emission. Using a Spearman rank-order analysis, the 14-195 keV luminosities for the Seyfert 1 and 2 galaxies are weakly statistically correlated with the F70/F160 ratios.

We investigate whether a subset of high-energy events observed by IceCube may be due to neutrinos from Sagittarius A*. We check both spatial and temporal coincidences of IceCube events with other transient activities of Sagittarius A*. Among the seven IceCube shower events nearest to the galactic center, we have found that event 25 has a time very close to (around three hours after) the brightest X-ray flare of Sagittarius A* observed by the Chandra X-ray Observatory with a p-value of 0.9%. Furthermore, two of the seven events occurred within one day of each other (there is a 1.6% probability that this would occur for a random distribution in time). Thus, the determination that some IceCube events occur at similar times as X-ray flares and others occur in a burst could be the smoking gun that Sagittarius A* is a point source of very high energy neutrinos. We point out that if IceCube galactic center neutrino events originate from charged pion decays, then TeV gamma rays should come from neutral pion decays at a similar rate. We show that the CTA, HAWC, H.E.S.S. and VERITAS experiments should be sensitive enough to test this.

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.

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

We use the James Clerk Maxwell Telescope's SCUBA-2 camera to image a 400 arcmin^2 area surrounding the GOODS-N field. The 850 micron rms noise ranges from a value of 0.49 mJy in the central region to 3.5 mJy at the outside edge. From these data, we construct an 850 micron source catalog to 2 mJy containing 49 sources detected above the 4-sigma level. We use an ultradeep (11.5 uJy at 5-sigma) 1.4 GHz image obtained with the Karl G. Jansky Very Large Array together with observations made with the Submillimeter Array to identify counterparts to the submillimeter galaxies. For most cases of multiple radio counterparts, we can identify the correct counterpart from new and existing Submillimeter Array data. We have spectroscopic redshifts for 62% of the radio sources in the 9 arcmin radius highest sensitivity region (556/894) and 67% of the radio sources in the GOODS-N region (367/543). We supplement these with a modest number of additional photometric redshifts in the GOODS-N region (30). We measure millimetric redshifts from the radio to submillimeter flux ratios for the unidentified submillimeter sample, assuming an Arp 220 spectral energy distribution. We find a radio flux dependent K-z relation for the radio sources, which we use to estimate redshifts for the remaining radio sources. We determine the star formation rates (SFRs) of the submillimeter sources based on their radio powers and their submillimeter and find that they agree well. The radio data are deep enough to detect star-forming galaxies with SFRs >2000 solar masses per year to z~6. We find galaxies with SFRs up to ~6,000 solar masses per year over the redshift range z=1.5-6, but we see evidence for a turn-down in the SFR distribution function above 2000 solar masses per year.

A community meeting on the topic of "Radio Astronomy in the LSST Era" was hosted by the National Radio Astronomy Observatory in Charlottesville, VA (2013 May 6--8). The focus of the workshop was on time domain radio astronomy and sky surveys. For the time domain, the extent to which radio and visible wavelength observations are required to understand several classes of transients was stressed, but there are also classes of radio transients for which no visible wavelength counterpart is yet known, providing an opportunity for discovery. From the LSST perspective, the LSST is expected to generate as many as 1 million alerts nightly, which will require even more selective specification and identification of the classes and characteristics of transients that can warrant follow up, at radio or any wavelength. The LSST will also conduct a deep survey of the sky, producing a catalog expected to contain over 38 billion objects in it. Deep radio wavelength sky surveys will also be conducted on a comparable time scale, and radio and visible wavelength observations are part of the multi-wavelength approach needed to classify and understand these objects. Radio wavelengths are valuable because they are unaffected by dust obscuration and, for galaxies, contain contributions both from star formation and from active galactic nuclei. The workshop touched on several other topics, on which there was consensus including the placement of other LSST "Deep Drilling Fields," inter-operability of software tools, and the challenge of filtering and exploiting the LSST data stream. There were also topics for which there was insufficient time for full discussion or for which no consensus was reached, which included the procedures for following up on LSST observations and the nature for future support of researchers desiring to use LSST data products.

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.

We present deep 450um and 850um observations of a large, uniformly covered 394arcmin^2 area in the COSMOS field obtained with the SCUBA-2 instrument on the James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 4.13mJy at 450um and 0.80mJy at 850um. The differential and cumulative number counts are presented and compared to similar previous works. Individual point sources are identified at >3.6sigma significance, a threshold corresponding to a 3-5% sample contamination rate. We identify 78 sources at 450um and 99 at 850um, with flux densities S450=13-37mJy and S850=2-16mJy. Only 62-76% of 450um sources are 850um detected and 61-81% of 850um sources are 450um detected. The positional uncertainties at 450um are small (1-2.5") and therefore allow a precise identification of multiwavelength counterparts without reliance on detection at 24um or radio wavelengths; we find that only 44% of 450um-selected galaxies and 60% of 850um-sources have 24um or radio counterparts. 450um-selected galaxies peak at <z>=1.95+-0.19 and 850um=selected galaxies peak at <z>=2.16+-0.11. The two samples occupy similar parameter space in redshift and luminosity, while their median SED peak wavelengths differ by ~10-50um (translating to deltaTdust =8-12K, where 450um-selected galaxies are warmer). The similarities of the 450um and 850um populations, yet lack of direct overlap between them, suggests that submillimeter surveys conducted at any single far-infrared wavelength will be significantly incomplete (~>30%) at censusing infrared-luminous star formation at high-z.

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.

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.

We carried out extremely sensitive Submillimeter Array (SMA) 340 GHz (860 micron) continuum imaging of a complete sample of SCUBA 850 micron sources (>4 sigma) with fluxes >3 mJy in the GOODS-N. Using these data and new SCUBA-2 data, we do not detect 4 of the 16 SCUBA sources, and we rule out the original SCUBA fluxes at the 4 sigma level. Three more resolve into multiple fainter SMA galaxies, suggesting that our understanding of the most luminous high-redshift dusty galaxies may not be as reliable as we thought. 10 of the 16 independent SMA sources have spectroscopic redshifts (optical/infrared or CO) to z=5.18. Using a new, ultradeep 20 cm image obtained with the Karl G. Jansky Very Large Array (rms of 2.5 microJy), we find that all 16 of the SMA sources are detected at >5 sigma. Using Herschel far-infrared (FIR) data, we show that the five isolated SMA sources with Herschel detections are well described by an Arp 220 spectral energy distribution template in the FIR. They also closely obey the local FIR-radio correlation, a result that does not suffer from a radio bias. We compute the contribution from the 16 SMA sources to the universal star formation rate (SFR) per comoving volume. With individual SFRs in the range 700-5000 solar masses per year, they contribute ~30% of the extinction-corrected ultraviolet-selected SFR density from z=1 to at least z=5. Star formation histories determined from extinction-corrected ultraviolet populations and from submillimeter galaxy populations only partially overlap, due to the extreme ultraviolet faintness of some submillimeter galaxies.

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]

We present 1.4 GHz catalogs for the cluster fields Abell 370 and Abell 2390 observed with the Very Large Array. These are two of the deepest radio images of cluster fields ever taken. The Abell 370 image covers an area of 40'x40' with a synthesized beam of ~1.7" and a noise level of ~5.7 uJy near field center. The Abell 2390 image covers an area of 34'x34' with a synthesized beam of ~1.4" and a noise level of ~5.6 uJy near field center. We catalog 200 redshifts for the Abell 370 field. We construct differential number counts for the central regions (radius < 16') of both clusters. We find that the faint (S_1.4GHz < 3 mJy) counts of Abell 370 are roughly consistent with the highest blank field number counts, while the faint number counts of Abell 2390 are roughly consistent with the lowest blank field number counts. Our analyses indicate that the number counts are primarily from field radio galaxies. We suggest that the disagreement of our counts can be largely attributed to cosmic variance.

Recent cosmological modeling efforts have shown that a local underdensity on scales of a few hundred Mpc (out to z ~ 0.1), could produce the apparent acceleration of the expansion of the universe observed via type Ia supernovae. Several studies of galaxy counts in the near-infrared (NIR) have found that the local universe appears under-dense by ~25-50% compared with regions a few hundred Mpc distant. Galaxy counts at low redshifts sample primarily L ~ L* galaxies. Thus, if the local universe is under-dense, then the normalization of the NIR galaxy luminosity function (LF) at z>0.1 should be higher than that measured for z<0.1. Here we present a highly complete (> 90%) spectroscopic sample of 1436 galaxies selected in the H-band to study the normalization of the NIR LF at 0.1<z<0.3 and address the question of whether or not we reside in a large local underdensity. We find that for the combination of our six fields, the product phi* L* at 0.1 < z < 0.3 is ~ 30% higher than that measured at lower redshifts. While our statistical errors in this measurement are on the ~10% level, we find the systematics due to cosmic variance may be larger still. We investigate the effects of cosmic variance on our measurement using the COSMOS cone mock catalogs from the Millennium simulation and recent empirical estimates. We find that our survey is subject to systematic uncertainties due to cosmic variance at the 15% level ($1 sigma), representing an improvement by a factor of ~ 2 over previous studies in this redshift range. We conclude that observations cannot yet rule out the possibility that the local universe is under-dense at z<0.1.

We describe a method for obtaining a flux-limited sample of Ly-alpha emitters from GALEX grism data. We show that the multiple GALEX grism images can be converted into a three-dimensional (two spatial axes and one wavelength axis) data cube. The wavelength slices may then be treated as narrowband images and searched for emission-line galaxies. For the GALEX NUV grism data, the method provides a Ly-alpha flux-limited sample over the redshift range z=0.67-1.16. We test the method on the Chandra Deep Field South field, where we find 28 Ly-alpha emitters with faint continuum magnitudes (NUV>22) that are not present in the GALEX pipeline sample. We measure the completeness by adding artificial emitters and measuring the fraction recovered. We find that we have an 80% completeness above a Ly-alpha flux of 10^-15 erg/cm^2/s. We use the UV spectra and the available X-ray data and optical spectra to estimate the fraction of active galactic nuclei in the selection. We report the first detection of a giant Ly-alpha blob at z<1, though we find that these objects are much less common at z=1 than at z=3. Finally, we compute limits on the z~1 Ly-alpha luminosity function and confirm that there is a dramatic evolution in the luminosity function over the redshift range z=0-1.

We use the new 4 Ms exposure of the CDF-S field obtained with the Chandra X-ray satellite to investigate the properties of the faintest X-ray sources over a wide range of redshifts. We use an optimized averaging procedure to investigate the weighted mean X-ray fluxes of optically selected sources in the CDF-S over the redshift range z=0-8 and down to 0.5-2 keV fluxes as low as 5e-19 erg/cm^2/s. None of the samples of sources at high redshifts (z>5) show any significant flux, and at z=6.5 we place an upper limit on the X-ray luminosity of 4e41 erg/s in the rest-frame 3.75-15 keV band for the sample of Bouwens et al. (2006). This is consistent with any X-ray production in the galaxies being solely due to star formation. At lower redshifts we find significant weighted mean X-ray fluxes in many samples of sources over the redshift range z=0-4. We use these to argue that (1) the relation between star formation and X-ray production remains invariant over this redshift range, (2) X-ray sources below the direct detection threshold in the CDF-S are primarily star-forming, and (3) there is full consistency between UV and X-ray estimations of the star formation history.

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.

The Baldwin, Phillips, and Terlevich emission-line ratio diagnostic ([OIII]/H{\beta} versus [NII]/H{\alpha}, hereafter BPT diagram) efficiently separates galaxies whose signal is dominated by star formation <BPT-SF> from those dominated by AGN activity (BPT-AGN). Yet this BPT diagram is limited to z < 0.5, the redshift at which [NII]{\lambda}6584 leaves the optical spectral window. Using the Sloan Digital Sky Survey (SDSS), we construct a new diagnostic, or TBT diagram, that is based on rest-frame g-z color, [NeIII]{\lambda}3869, and [OII]{\lambda}{\lambda}3726 + 3729 and can be used for galaxies out to z < 1.4. The TBT diagram identifies 98.7% of the SDSS BPT-AGN as TBT-AGN and 97% of the SDSS BPT-SF as TBT-SF. Furthermore, it identifies 97% of the OPTX Chandra X-ray selected AGNs as TBT-AGN. This is in contrast to the BPT diagram, which misidentifies 20% of X-ray selected AGNs as BPT-SF. We use the GOODS-N and Lockman Hole galaxy samples, with their accompanying deep Chandra imaging, to perform X-ray and infrared stacking analyses to further validate our TBT-AGN and TBT-SF selections; that is, we verify the dominance of AGN activity in the former and star formation activity in the latter. Finally, we address the inclusion of the majority of the BPT-comp (sources lying between the BPT-SF and BPT-AGN regimes) in our TBT-AGN regime. We find that the stacked BPT-comp source is X-ray hard (<{\Gamma}eff> = 1.0 +/-0.4) and has a high X-ray luminosity to total infrared luminosity ratio. This suggests that, on average, the X-ray signal in BPT-comp is dominated by obscured or low accretion rate AGN activity rather than by star formation, supporting their inclusion in the TBT-AGN regime.

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.

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.

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 for measurements of {\nu}I{\nu} 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 J H Ks ~ 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 lim- its 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, as well as recent model estimates for the total light from galaxies. We exam- ine 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.

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.

We present an atlas of 88 z~5.7 and 30 z~6.5 Ly alpha emitters obtained from a wide-field narrowband survey. We combined deep narrowband imaging in 120A bandpass filters centered at 8150A and 9140A with deep BVRIz broadband imaging to select high-redshift galaxy candidates over an area of 4180 square arcmin. The goal was to obtain a uniform selection of comparable depth over the 7 targeted fields in the two filters. For the GOODS-N region of the HDF-N field, we also selected candidates using a 120A filter centered at 9210A. We made spectroscopic observations with Keck DEIMOS of nearly all the candidates to obtain the final sample of Ly alpha emitters. At the 3.3A resolution of the DEIMOS observations the asymmetric profile for Ly alpha emission with its steep blue fall-off can be clearly seen in the spectra of nearly all the galaxies. We show that the spectral profiles are surprisingly similar for many of the galaxies and that the composite spectral profiles are nearly identical at z=5.7 and z=6.5. We analyze the distributions of line widths and Ly alpha equivalent widths and find that the lines are marginally narrower at the higher redshift, with median values of 0.77A at z=6.5 and 0.92A at z=5.7. The line widths have a dependence on the Ly alpha luminosity of the form L(L alpha)^(0.3). We compare the surface densities and the luminosity functions at the two redshifts and find that there is a multiplicative factor of 2 decrease in the number density of bright Ly alpha emitters from z=5.7 to z=6.5, while the characteristic luminosity is unchanged.

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.

We compare optical and hard X-ray identifications of AGNs using a uniformly selected (above a flux limit of f_2-8 keV = 3.5e-15 erg/cm2/s) and highly optically spectroscopically complete ( > 80% for f_2-8 keV > 1e-14 erg/cm2/s and > 60% below) 2-8 keV sample observed in three Chandra fields (CLANS, CLASXS, and the CDF-N). We find that empirical emission-line ratio diagnostic diagrams misidentify 20-50% of the X-ray selected AGNs that can be put on these diagrams as star formers, depending on which division is used. We confirm that there is a large (2 orders in magnitude) dispersion in the log ratio of the [OIII]5007A to hard X-ray luminosities for the non-broad line AGNs, even after applying reddening corrections to the [OIII] luminosities. We find that the dispersion is similar for the broad-line AGNs, where there is not expected to be much X-ray absorption from an obscuring torus around the AGN nor much obscuration from the galaxy along the line-of-sight if the AGN is aligned with the galaxy. We postulate that the X-ray selected AGNs that are misidentified by the diagnostic diagrams have low [OIII] luminosities due to the complexity of the structure of the narrow-line region, which causes many ionizing photons from the AGN not to be absorbed. This would mean that the [OIII] luminosity can only be used to predict the X-ray luminosity to within a factor of ~3 (one sigma). Despite selection effects, we show that the shapes and normalizations of the [OIII] and transformed hard X-ray luminosity functions show reasonable agreement, suggesting that the [OIII] samples are not finding substantially more AGNs at low redshifts than hard X-ray samples.

We use a complete sample of Lya emission-line selected AGNs obtained from nine deep blank fields observed with the grism spectrographs on the Galaxy Evolution Explorer (GALEX) satellite to measure the normalization and spectral shape of the AGN contribution to the ionizing background (rest-frame wavelengths 700-900 A) at z~1. Our sample consists of 139 sources selected in the redshift range z=0.65-1.25 in the near-ultraviolet (NUV; 2371 A central wavelength) channel. The area covered is 8.2 square degrees to a NUV magnitude of 20.5 (AB) and 0.92 square degrees at the faintest magnitude limit of 21.8. The GALEX AGN luminosity function agrees well with those obtained using optical and X-ray AGN samples, and the measured redshift evolution of the ionizing volume emissivity is similar to that previously obtained by measuring the GALEX far-ultraviolet (FUV; 1528 A central wavelength) magnitudes of an X-ray selected sample. For the first time we are able to construct the shape of the ionizing background at z~1 in a fully self-consistent way.

We present sensitive phase-referenced VLBI results on the radio continuum emission from the z=1.87 luminous submillimeter galaxy (SMG) GOODS 850-3. The observations were carried out at 1.4 GHz using the High Sensitivity Array (HSA). Our sensitive tapered VLBI image of GOODS 850-3 at 0.47 x 0.34 arcsec (3.9 x 2.9 kpc) resolution shows a marginally resolved continuum structure with a peak flux density of 148 \pm 38 uJy/beam, and a total flux density of 168 \pm 73 uJy, consistent with previous VLA and MERLIN measurements. The derived intrinsic brightness temperature is > 5 \pm 2 x 10^3 K. The radio continuum position of this galaxy coincides with a bright and extended near-infrared source that nearly disappears in the deep HST optical image, indicating a dusty source of nearly 9 kpc in diameter. No continuum emission is detected at the full VLBI resolution (13.2 x 7.2 mas, 111 x 61 pc), with a 4-sigma point source upper limit of 26 uJy/beam, or an upper limit to the intrinsic brightness temperature of 4.7 x 10^5 K. The extent of the observed continuum source at 1.4 GHz and the derived brightness temperature limits are consistent with the radio emission (and thus presumably the far-infrared emission) being powered by a major starburst in GOODS 850-3, with a star formation rate of ~2500 M_sun/yr. Moreover, the absence of any continuum emission at the full resolution of the VLBI observations indicates the lack of a compact radio AGN source in this z=1.87 SMG.

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.

We construct a sample of low-redshift Ly-alpha emission-line selected sources from GALEX grism spectroscopy of nine deep fields to study the role of Ly-alpha emission in galaxy populations with cosmic time. Our final sample consists of 119 (141) sources selected in the redshift interval z=0.195-0.44 (z=0.65-1.25) from the FUV (NUV) channel. We classify the Ly-alpha sources as AGNs if high-ionization emission lines are present in their UV spectra and as possible star-forming galaxies otherwise. We classify additional sources as AGNs using line widths for our Ly-alpha emitter (LAE) analysis. These classifications are broadly supported by comparisons with X-ray and optical spectroscopic observations. Defining the GALEX LAE sample in the same way as high-redshift LAE samples, we show that LAEs constitute only about 5% of NUV-continuum selected galaxies at z~0.3. We also show that they are less common at z~0.3 than they are at z~3. We find that the z~0.3 optically-confirmed Ly-alpha galaxies lie below the metallicity-luminosity relation of the z~0.3 NUV-continuum selected galaxies but have similar H-alpha velocity widths at similar luminosities, suggesting that they also lie below the metallicity-mass relation of the NUV-continuum selected galaxies. We show that, on average, the Ly-alpha galaxies have bluer colors, lower extinctions as measured from the Balmer line ratios, and more compact morphologies than the NUV-continuum selected galaxies. Finally, we confirm that the z~2 Lyman break galaxies (LBGs) have relatively low metallicities for their luminosities, and we find that they lie in the same metallicity range as the z~0.3 Ly-alpha galaxies.

We present a deep, wide-field near-infrared (NIR) survey over five widely separated fields at high Galactic latitude covering a total of ~ 3 deg^2 in J, H, and Ks. The deepest areas of the data (~ 0.25 deg^2) extend to a 5 sigma limiting magnitude of JHKs > 24 in the AB magnitude system. Although depth and area vary from field to field, the overall depth and large area of this dataset make it one of the deepest wide-field NIR imaging surveys to date. This paper discusses the observations, data reduction, and bright galaxy counts in these fields. We compare the slope of the bright galaxy counts with the Two Micron All Sky Survey (2MASS) and other counts from the literature and explore the relationship between slope and supergalactic latitude. The slope near the supergalactic equator is sub- Euclidean on average pointing to the possibility of a decreasing average space density of galaxies by ~ 10-15% over scales of ~ 250-350 Mpc. On the contrary, the slope at high supergalactic latitudes is strongly super-Euclidean on average suggesting an increase in the space density of galaxies as one moves from the voids just above and below the supergalactic plane out to distances of ~ 250-350 Mpc. These results suggest that local large scale structure could be responsible for large discrepancies in the measured slope between different studies in the past. In addition, the local universe away from the supergalactic plane appears to be underdense by ~ 25-100% relative to the space densities of a few hundred megaparsecs distant. Subject headings: cosmology: observations and large scale structure of universe-galaxies: fundamental parameters (counts)-infrared: galaxies

We compare the optical spectral types with the X-ray spectral properties for a uniformly selected (sources with fluxes greater than the 3 sigma level and above a flux limit of f_2-8 keV > 3.5x10^-15 erg/cm2/s), highly spectroscopically complete (>80% for f_2-8 keV > 10^-14 erg/cm2/s and >60% below) 2-8 keV X-ray sample observed in three Chandra fields (CLANS, CLASXS, and the CDF-N) that cover ~1.2 deg^2. For our sample of 645 spectroscopically observed sources, we confirm that there is significant overlap of the X-ray spectral properties, as determined by the effective photon indices, Geff, obtained from the ratios of the 0.5-2 keV to 2-8 keV counts, for the different optical spectral types. For example, of the broad-line AGNs (non-broad-line AGNs), 20% +/- 3% (33% +/- 4%) have Geff<1.2 (Geff > 1.2). Thus, one cannot use the X-ray spectral classifications and the optical spectral classifications equivalently. Since it is not understood how X-ray and optical classifications relate to the obscuration of the central engine, we strongly advise against a mixed classification scheme, as it can only complicate the interpretation of X-ray AGN samples. We confirm the dependence of optical spectral type on X-ray luminosity, and for z<1, we find a similar luminosity dependence of Geff. However, this dependence breaks down at higher redshifts due to the highly redshift-dependent nature of Geff. We therefore also caution that any classification scheme which depends on Geff is likely to suffer from serious redshift bias.

We report a highly significant Submillimeter Array (SMA) detection of the prototypical submillimeter source HDF 850.1, which is the brightest submillimeter source in the Hubble Deep Field-North proper. The detection yields an extremely precise position of RA(2000)=12:36:51.99 and Dec(2000)=+62:12:25.83 with a 1-sigma positional uncertainty of 0.17 arcsec. The position is consistent with the location of a millimeter wavelength interferometric detection and with the locations of weak VLA detections at 1.4 and 8.4 GHz, but it is not consistent with any previous optical/near-infrared identifications. The source appears pointlike at the 2 arcsec resolution of the SMA, and the detected flux of 7.8+/-1.0 mJy is consistent with the measured SCUBA fluxes. We tabulate fluxes and limits on HDF 850.1 at other wavelengths. Our redshift estimate for HDF 850.1 based on the radio through mid-infrared measurements is z=4.1. The faintness of the source at optical/near-infrared wavelengths and the high estimated redshift suggest that HDF 850.1 may be an analog of the brighter submillimeter source GOODS 850-5, which is also thought to be at z>4. The fact that a source like HDF 850.1 should have appeared in one of the very first blank-field SCUBA observations ever made suggests that such high-redshift sources are quite common. Thus, we are led to conclude that high-redshift star formation is dominated by giant dusty star-forming galaxies, just as it is at lower redshifts.

We use the largest, most uniform, and most spectroscopically complete to faint X-ray flux limits Chandra sample to date to construct hard 2-8 keV rest-frame X-ray luminosity functions (HXLFs) of spectroscopically identified active galactic nuclei (AGNs) to z~5. In addition, we use a new 2-8 keV local sample selected by the very hard (14-195 keV) SWIFT 9-month Burst Alert Telescope (BAT) survey to construct the local 2-8 keV HXLF. We do maximum likelihood fits of the combined distant plus local sample (as well as of the distant sample alone) over the redshift intervals 0<z<1.2, 0<z<3, and 0<z<5 using a variety of analytic forms, which we compare with the HXLFs. We recommend using our luminosity dependent density evolution (LDDE) model fits of the combined distant plus local sample over 0<z<3 for all the spectroscopically identified sources and for the broad-line AGNs.

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.

We present a table of redshifts for 2907 galaxies and stars in the 145 square arcmin HST ACS GOODS-North, making this the most spectroscopically complete redshift sample obtained to date in a field of this size. We also include the redshifts, where available, in a table containing just under 7000 galaxies from the ACS area with K_s(AB)<24.5 measured from a deep K_s image obtained with WIRCam on the CFHT, as well as in a table containing 1016 sources with NUV(AB)<25 and 478 sources with FUV(AB)<25.5 (there is considerable overlap) measured from the deep GALEX images in the ACS area. Finally, we include the redshifts, where available, in a table containing the 1199 24 micron sources to 80 uJy measured from the wider-area Spitzer GOODS-North. The redshift identifications are greater than 90% complete to magnitudes of F435W(AB)=24.5, F850LP(AB)=23.3, and K_s(AB)=21.5 and to 24 micron fluxes of 250 uJy. An extensive analysis of these data will appear in a parallel paper, but here we determine how efficient color-selection techniques are at identifying high-redshift galaxies and Active Galactic Nuclei. We also examine the feasibility of doing tomography of the intergalactic medium with a 30 m telescope.

We use a wide-field (0.9 square degree) X-ray sample with optical and GALEX ultraviolet observations to measure the contribution of Active Galactic Nuclei (AGNs) to the ionizing flux as a function of redshift. Our analysis shows that the AGN contribution to the metagalactic ionizing background peaks around z=2. The measured values of the ionizing background from the AGNs are lower than previous estimates and confirm that ionization from AGNs is insufficient to maintain the observed ionization of the intergalactic medium (IGM) at z>3. We show that only sources with broad lines in their optical spectra have detectable ionizing flux and that the ionizing flux seen in an AGN is not correlated with its X-ray color. We also use the GALEX observations of the GOODS-N region to place a 2-sigma upper limit of 0.008 on the average ionization fraction fnu(700 A)/fnu(1500 A) for 626 UV selected galaxies in the redshift range z=0.9-1.4. We then use this limit to estimate an upper bound to the galaxy contribution in the redshift range z=0-5. If the z~1.15 ionization fraction is appropriate for higher redshift galaxies, then contributions from the galaxy population are also too low to account for the IGM ionization at the highest redshifts (z>4).

We present the redshift catalogs for the X-ray sources detected in the Chandra Deep Field North (CDF-N), the Chandra Large Area Synoptic X-ray Survey (CLASXS), and the Chandra Lockman Area North Survey (CLANS). The catalogs for the CDF-N and CLASXS fields include redshifts from previous work, while the redshifts for the CLANS field are all new. For fluxes above 10^-14 ergs cm^-2 s^-1 (2-8 keV) we have redshifts for 76% of the sources. We extend the redshift information for the full sample using photometric redshifts. The goal of the OPTX Project is to use these three surveys, which are among the most spectroscopically complete surveys to date, to analyze the effect of spectral type on the shape and evolution of the X-ray luminosity functions and to compare the optical spectral types with the X-ray spectral properties. We also present the CLANS X-ray catalog. The nine ACIS-I fields cover a solid angle of ~0.6 square degrees and reach fluxes of 7x10^-16 ergs cm^-2 s^-1 (0.5-2 keV) and 3.5x10^-15 ergs cm^-2 s^-1 (2-8 keV). We find a total of 761 X-ray point sources. Additionally, we present the optical and infrared photometric catalog for the CLANS X-ray sources, as well as updated optical and infrared photometric catalogs for the X-ray sources in the CLASXS and CDF-N fields. The CLANS and CLASXS surveys bridge the gap between the ultradeep pencil-beam surveys, such as the CDFs, and the shallower, very large-area surveys. As a result, they probe the X-ray sources that contribute the bulk of the 2-8 keV X-ray background and cover the flux range of the observed break in the logN-logS distribution. We construct differential number counts for each individual field and for the full sample.

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.

We present an integrated study of star formation and galactic stellar mass assembly from z=0.05-1.5 and galactic metallicity evolution from z=0.05-0.9 using a very large and highly spectroscopically complete sample selected by rest-frame NIR bolometric flux in the GOODS-N. We assume a Salpeter IMF and fit Bruzual & Charlot (2003) models to compute the galactic stellar masses and extinctions. We determine the expected formed stellar mass density growth rates produced by star formation and compare them with the growth rates measured from the formed stellar mass functions by mass interval. We show that the growth rates match if the IMF is slightly increased from the Salpeter IMF at intermediate masses (~10 solar masses). We investigate the evolution of galaxy color, spectral type, and morphology with mass and redshift and the evolution of mass with environment. We find that applying extinction corrections is critical when analyzing galaxy colors; e.g., nearly all of the galaxies in the green valley are 24um sources, but after correcting for extinction, the bulk of the 24um sources lie in the blue cloud. We find an evolution of the metallicity-mass relation corresponding to a decrease of 0.21+/-0.03 dex between the local value and the value at z=0.77 in the 1e10-1e11 solar mass range. We use the metallicity evolution to estimate the gas mass of the galaxies, which we compare with the galactic stellar mass assembly and star formation histories. Overall, our measurements are consistent with a galaxy evolution process dominated by episodic bursts of star formation and where star formation in the most massive galaxies (>1e11 solar masses) ceases at z<1.5 because of gas starvation. (Abstract abridged)

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.

We report results of a recently-completed pre-Formulation Phase study of SPIRIT, a candidate NASA Origins Probe mission. SPIRIT is a spatial and spectral interferometer with an operating wavelength range 25 - 400 microns. SPIRIT will provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view to accomplish three primary scientific objectives: (1) Learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. Observations with SPIRIT will be complementary to those of the James Webb Space Telescope and the ground-based Atacama Large Millimeter Array. All three observatories could be operational contemporaneously.

We use highly spectroscopically complete observations of the radio sources from the VLA 1.4 GHz survey of the HDF-N region to study the faint radio galaxy population and its evolution. We spectrally classify the sources into four spectral types: absorbers, star formers, Seyfert galaxies, and broad-line AGNs, and we analyze their properties by type. We supplement the spectroscopic redshifts with photometric redshifts measured from the rest-frame UV to MIR spectral energy distributions. Using deep X-ray observations of the field, we do not confirm the existence of an X-ray-radio correlation for star-forming galaxies. We also do not observe any correlations between 1.4 GHz flux and R magnitude or redshift. We find that the radio powers of the host galaxies rise dramatically with increasing redshift, while the optical properties of the host galaxies show at most small changes. Assuming that the locally determined FIR-radio correlation holds at high redshifts, we estimate total FIR luminosities for the radio sources. We note that the FIR luminosity estimates for any radio-loud AGNs will be overestimates. Considering only the radio sources with quasar-like bolometric luminosities, we find a maximum ratio of candidate highly-obscured AGNs to X-ray-luminous (>10^42 ergs/s) sources of about 1.9. We use source-stacking analyses to measure the X-ray surface brightnesses of various X-ray and radio populations. We find the contributions to the 4-8 keV light from our candidate highly-obscured AGNs to be very small, and hence these sources are unable to account for the light that has been suggested may be missing at these energies.

We present a 1.4-GHz catalog of 810 radio sources (560 sources in the complete sample) with 1.8" resolution found within a 17' radius in the SSA13 field (RA=13h12m,DEC=42d38'). The radio image from the VLA has an rms noise level of 4.82 microJy/beam at the field center, and Subaru optical images in r-band (6300A) and z-band (9200A) have a three-sigma detection magnitude of 26.1 and 24.9, respectively. 88% of the radio sources are identified with an optical counterpart, and there is significantly more reddening for objects fainter than 24-mag. The radio and optical parameters are tabulated, and source morphologies are displayed by radio contours overlaying optical false-color images. The radio structures show a wealth of complexity and these are classified into a small number of categories. About one-third of the radio sources are larger than 1.2" and their orientation is often similar to that of the associated galaxy or binary-galaxy system. The density of sources in the SSA13 field above 75 microJy is 0.40 per square arcmin, with a slope of -2.43 in the differential counts. The radio spectral index may steepen for sources below 75 microJy. We estimate that at most 40% of the microJansky radio sources are dominated by AGN processes.

We use new deep near-infrared (NIR) and mid-infrared (MIR) observations to analyze the 850$~\mu$m image of the GOODS HDF-N region. We show that much of the submillimeter background at this wavelength is picked out by sources with $H(AB)$ or 3.6um (AB)<23.25 (1.8 uJy). These sources contribute an 850um background of 24\pm2 Jy deg^-2. This is a much higher fraction of the measured background (31-45 Jy deg^-2) than is found with current 20cm or 24um samples. Roughly one-half of these NIR-selected sources have spectroscopic identifications, and we can assign robust photometric redshifts to nearly all of the remaining sources using their UV to MIR spectral energy distributions. We use the redshift and spectral type information to show that a large fraction of the 850um background light comes from sources with z=0-1.5 and that the sources responsible have intermediate spectral types. Neither the elliptical galaxies, which have no star formation, nor the bluest galaxies, which have little dust, contribute a significant amount of 850um light, despite the fact that together they comprise approximately half of the galaxies in the sample. The galaxies with intermediate spectral types have a mean flux of 0.40\pm0.03 mJy at 850um and 9.1\pm0.3 uJy at 20cm. The redshift distribution of the NIR-selected 850um light lies well below that of the much smaller amount of light traced by the more luminous, radio- selected submillimeter sources. We therefore require a revised star-formation history with a lower star-formation rate at high redshifts. We use a stacking analysis of the 20cm light in the NIR sample to show that the star-formation history of the total 850um sample is relatively flat down to z~1 and that half of the total star formation occurs at redshifts z<1.4.

We present deep Hubble Space Telescope and Spitzer Space Telescope observations of the host galaxy of GRB 050904 at z=6.295. The host is detected in the H-band and marginally at 3.6 micron. From these detections, and limits in the z'-band and 4.5 micron, we infer an extinction-corrected absolute magnitude, M(UV)=-20.7 mag, or ~L*, a substantial star formation rate of 15 solar masses per year, and a stellar mass of a few 10^9 solar masses. A comparison to the published sample of spectroscopically-confirmed galaxies at z>5.5 reveals that the host of GRB 050904 would evade detection and/or confirmation in any of the current surveys due to the lack of detectable Ly-alpha emission, which is likely the result of dust extinction (A[1200]~1.5 mag). This suggests that not all luminous starburst galaxies at z~6 are currently being accounted for. Most importantly, using the metallicity of Z~0.05 solar inferred from the afterglow absorption spectrum, our observations indicate for the first time that the observed evolution in the mass- and luminosity-metallicity relations from z=0 to z~2 continues on to z>6. The ease of measuring redshifts and metallicities from the afterglow emission suggests that in tandem with the next generation ground- and space-based telescopes, a GRB mission with dedicated near-IR follow-up can provide unique information on the evolution of stars and galaxies through the epoch of re-ionization.

We present the spatial correlation function analysis of non-stellar X-ray point sources in the Chandra Large Area Synoptic X-ray Survey of Lockman Hole Northwest (CLASXS). Our 9 ACIS-I fields cover a contiguous solid angle of 0.4 deg^2 and reach a depth of 3x10^-15 c.g.s in the 2-8 keV band. We supplement our analysis with data from the Chandra Deep Field North (CDFN). The addition of this field allows better probe of the correlation function at small scales. A total of 233 and 252 sources with spectroscopic information are used in the study of the CLASXS and CDFN fields respectively. We calculate both redshift-space and projected correlation functions in comoving coordinates, averaged over the redshift range of 0.1<z<3.0, for both CLASXS and CDFN fields for a standard cosmology with \Omega_{\Lambda} = 0.73, \Omega_{M} = 0.27, and H_0 = 0.71. The correlation function for the CLASXS field over scales of 3 Mpc <s< 200 Mpc can be modeled as a power-law of the form \xi(s) = (s/s_0)^{-\gamma}, with \gamma = 1.6^{+0.4}_{-0.3} and s_0 = 8.0^{+1.4}_{-1.5} Mpc. The redshift-space correlation function for CDFN on scales of 1 Mpc$<s<$~100 Mpc is found to have a similar correlation length $s_0 = 8.55^{+0.75}_{-0.74}$ Mpc, but a shallower slope ($\gamma = 1.3 \pm 0.1$). The real-space correlation functions derived from the projected correlation functions, are found to be $r_0 = 8.1^{+1.2}_{-2.2}$ Mpc, and $\gamma = 2.1 \pm 0.5$ for the CLASXS field, and $r_0 = 5.8^{+1.0}_{-1.5}$ Mpc, $\gamma = 1.38^{+0.12}_{-0.14}$ for the CDFN field. By comparing the real- and redshift-space correlation functions in the combined CLASXS and CDFN samples, we are able to estimate the redshift distortion parameter $\beta = 0.4 \pm 0.2$ at an effective redshift $z = 0.94$.(abridged)

We use the combination of the 2 Ms Chandra X-ray image, new J and H band images, and the Spitzer IRAC and MIPS images of the Chandra Deep Field-North to obtain high spectroscopic and photometric redshift completeness of high and intermediate X-ray luminosity sources in the redshift interval z=2-3. We measure the number densities of z=2-3 active galactic nuclei (AGNs) and broad-line AGNs in the rest-frame 2-8 keV luminosity intervals 10^44-10^45 and 10^43-10^44 ergs/s and compare with previous lower redshift results. We confirm a decline in the number densities of intermediate-luminosity sources at z>1. We also measure the number density of z=2-3 AGNs in the luminosity interval 10^43-10^44.5 ergs/s and compare with previous low and high-redshift results. Again, we find a decline in the number densities at z>1. In both cases, we can rule out the hypothesis that the number densities remain flat to z=2-3 at above the 5-sigma level.

We use highly spectroscopically complete deep and wide-area Chandra surveys to determine the cosmic evolution of hard X-ray-selected AGNs. We determine hard X-ray luminosity functions (HXLFs) for all spectral types and for broad-line AGNs (BLAGNs) alone. At z<1.2, both are well described by pure luminosity evolution. Thus, all AGNs drop in luminosity by almost an order of magnitude over this redshift range. We show that this observed drop is due to AGN downsizing. We directly compare our BLAGN HXLFs with the optical QSO LFs and find that the optical QSO LFs do not probe faint enough to see the downturn in the BLAGN HXLFs. We rule out galaxy dilution as a partial explanation for the observation that BLAGNs dominate the number densities at the higher X-ray luminosities, while optically-narrow AGNs (FWHM<2000 km/s) dominate at the lower X-ray luminosities by measuring the nuclear UV/optical properties of the Chandra sources using the HST ACS GOODS-North data. The UV/optical nuclei of the optically-narrow AGNs are much weaker than expected if they were similar to the BLAGNs. We therefore postulate the need for a luminosity dependent unified model. Alternatively, the BLAGNs and the optically-narrow AGNs could be intrinsically different source populations. We cover both interpretations by constructing composite spectral energy distributions--including long-wavelength data from the MIR to the submillimeter--by spectral type and by X-ray luminosity. We use these to infer the bolometric corrections (from hard X-ray luminosities to bolometric luminosities) needed to map the accretion history. We determine the accreted supermassive black hole mass density for all spectral types and for BLAGNs alone using the observed evolution of the hard X-ray energy density production rate and our inferred bolometric corrections.

We present photometric and spectroscopic observations of the X-ray sources detected in the wide-area, moderately deep Chandra Large Area Synoptic X-ray Survey of the Lockman Hole-Northwest field. We have B, V, R, I, and z' photometry for 521 (99%) of the 525 sources in the X-ray catalog and spectroscopic redshifts for 271 (52%), including 20 stars. We do not find evidence for redshift groupings of the X-ray sources, like those found in the Chandra Deep Field surveys, because of the larger solid angle covered by this survey. We separate the X-ray sources by optical spectral type and examine the colors, apparent and absolute magnitudes, and redshift distributions for the broad-line and non-broad-line active galactic nuclei. Combining our wide-area survey with other Chandra and XMM-Newton hard X-ray surveys, we find a definite lack of luminous, high accretion rate sources at z<1, consistent with previous observations that showed that supermassive black hole growth is dominated at low redshifts by sources with low accretion rates.

We present the X-ray catalog and basic results from our Chandra Large Area Synoptic X-ray Survey (CLASXS) of the Lockman Hole-Northwest field. Our 9 ACIS-I fields cover a contiguous solid angle of ~0.4 sq. deg. and reach fluxes of 5E-16 cgs (0.4-2 keV) and 3E-15 cgs (2-8keV). Our survey bridges the gap between ultradeep pencil-beam surveys, such as the Chandra Deep Fields (CDFs), and shallower, large area surveys, allowing a better probe of the X-ray sources that contribute most of the 2-10 keV cosmic X-ray background (CXB). We find a total of 525 X-ray point sources and 4 extended sources. At ~10E-14 cgs 2-8 keV, our number counts are significantly higher than those of several non-contiguous, large area surveys. On the other hand, the integrated flux from the CLASXS field, combined with ASCA and Chandra ultradeep surveys, is consistent with results from other large area surveys, within the variance of the CXB. Spectral evolution is seen in the hardening of the sources at fluxes below 1E-14 cgs Above 4E1-14 cgs(0.4-8 keV), ~60 of the sources are variable. Four extended sources in CLASXS is consistent with the previously measured LogN-LogS of galaxy clusters. We report the discovery of a gravitational lensing arc associated with one of these sources. (Abridged)

The Hubble Deep Field-North (HDF-N) is one of the best studied extragalactic fields, and ultra-deep optical, radio, X-ray, and mid-infrared wide-field images are available for this area. Here we present an 850 um survey around the HDF-N, covering most of the area imaged by the Advanced Camera for Surveys as a part of the Great Observatories Origins Deep Survey. Our map has 0.4-4 mJy sensitivities (1 sigma) over an area ~110 arcmin^2 and there are 45 sources detected at >3 sigma. After correcting the effects of noise, confusion, incompleteness, and the Eddington bias using Monte Carlo simulations, we find that the detected 850 $\mu$m sources with fluxes greater than 2 mJy have a surface density of 3200^+1900_-1000 deg^-2 and account for about 24% to 34% of the far-infrared extragalactic background light. Using the deep radio interferometric image and the deep X-ray image, we are able to accurately locate ~60% of the bright submillimeter (submm) sources. In addition, by assuming the Arp 220 spectral energy distribution in the submm and radio, we estimate millimetric redshifts for the radio detected submm sources, and redshift lower limits for the ones not detected in the radio. Using the millimetric redshifts of the radio identified sources and spectroscopic and optical photometric redshifts for galaxies around the submm positions, we find a median redshift of 2.0 for 11 possibly identified sources, or a lower limit of 2.4 for the median redshift of our 4 sigma sample.

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.

We use redshift observations of two deep 1.4GHz fields to probe the evolution of the bright end of the radio galaxy luminosity function to z=1.5. We show that the number of galaxies with radio power that would correspond locally to an ultraluminous infrared galaxy (ULIG) evolves rapidly over this redshift range. The optical spectra and X-ray luminosities are used to refine this result by separating the sources with clear active galactic nucleus (AGN) signatures from those that may be dominated by star formation. Both populations show extremely rapid evolution over this redshift range. We find that the number of sources with ULIG radio power and no clear AGN signatures evolves as (1+z)^(7).

We report the results of an extensive spectroscopic survey of galaxies in the roughly 160 square arcminute ACS-GOODS region surrounding the HDF-N. We have identified 787 galaxies or stars with z'<24, R<24.5, or B<25 lying in the region. The spectra were obtained with either the DEIMOS or LRIS spectrographs on the Keck 10m telescopes. The results are compared with photometric redshift estimates and with redshifts from the literature, as well as with the redshifts of a parallel effort led by a group at Keck. Our sample, when combined with the literature data, provides identifications for 1180 sources. We use our results to determine the redshift distributions with magnitude, to analyze the rest-frame color distributions with redshift and spectral type, and to investigate the dependence of the X-ray galaxy properties on the local galaxy density in the redshift interval z=0-1.5. We find the rather surprising result that the galaxy X-ray properties are not strongly dependent on the local galaxy density for galaxies in the same luminosity range.

We have conducted a deep multi-color imaging survey of 0.2 degrees^2 centered on the Hubble Deep Field North (HDF-N). We shall refer to this region as the Hawaii-HDF-N. Deep data were collected in U, B, V, R, I, and z' bands over the central 0.2 degrees^2 and in HK' over a smaller region covering the Chandra Deep Field North (CDF-N). The data were reduced to have accurate relative photometry and astrometry across the entire field to facilitate photometric redshifts and spectroscopic followup. We have compiled a catalog of 48,858 objects in the central 0.2 degrees^2 detected at 5 sigma significance in a 3" aperture in either R or z' band. Number counts and color-magnitude diagrams are presented and shown to be consistent with previous observations. Using color selection we have measured the density of objects at 3<z<7. Our multi-color data indicates that samples selected at z>5.5 using the Lyman break technique suffer from more contamination by low redshift objects than suggested by previous studies.

We investigate how the fraction of broad-line sources in the AGN population changes with X-ray luminosity and redshift. We first construct the rest-frame hard-energy (2-8 keV) X-ray luminosity function (HXLF) at z=0.1-1 using Chandra Lockman Hole-Northwest wide-area data, Chandra Deep Field-North 2 Ms data, other Chandra deep field data, and the ASCA Large Sky Survey data. We find that broad-line AGNs dominate above 3e43 ergs/s and have a mean luminosity of 1.3e44 ergs/s. Type II AGNs can only become an important component of the X-ray population at Seyfert-like X-ray luminosities. We then construct z=0.1-0.5 and z=0.5-1 HXLFs and compare them with both the local HXLF measured from HEAO-1 A2 survey data and the z=1.5-3 HXLF measured from soft-energy (0.5-2 keV) Chandra and ROSAT data. We find that the number density of >1e44 ergs/s sources (quasars) steadily declines with decreasing redshift, while the number density of 1e43-1e44 ergs/s sources peaks at z=0.5-1. Strikingly, however, the number density of broad-line AGNs remains roughly constant with redshift while their average luminosities decline at the lower redshifts, showing another example of cosmic downsizing.

We present an optical and NIR catalog for the X-ray sources in the 2 Ms Chandra observation of the Hubble Deep Field-North region. We have high-quality multicolor images of all 503 X-ray point sources and reliable spectroscopic redshifts for 284. We spectroscopically identify six z>1 type II quasars. Our spectroscopic completeness for the R<24 sources is 87%. The spectroscopic redshift distribution shows two broad redshift spikes that have clearly grown over those originally seen in the 1 Ms exposure. The spectroscopically identified extragalactic sources already comprise 75% of the measured 2-8 keV light. Redshift slices versus 2-8 keV flux show that an impressive 54% of the measured 2-8 keV light arises from sources at z<1 and 68% from sources at z<2. We use seven broadband colors and a Bayesian photometric redshift estimation code to obtain photometric redshifts. The photometric redshifts are within 25% of the spectroscopic redshifts for 94% of the non-broad-line sources with both photometric and spectroscopic measurements. We use our wide wavelength coverage to determine rest-frame colors for the X-ray sources with spectroscopic or photometric redshifts. Many of the X-ray sources have the rest-frame colors of evolved red galaxies and there is very little evolution in these colors with redshift. We also determine absolute magnitudes and find that many of the non-broad-line sources are more luminous than Mstar, even at high redshifts. We therefore infer that deep X-ray observations may provide an effective way of locating Mstar galaxies with colors similar to present-day early-type galaxies to high redshifts. (Abridged)

We present catalogs for the ~2 Ms Chandra Deep Field-North, currently the deepest X-ray observation of the Universe in the 0.5-8.0 keV band. Five hundred and three (503) X-ray sources are detected over an ~448 sq.arcmin area in up to seven bands; 20 of these X-ray sources lie in the Hubble Deep Field-North. Source positions are determined using matched-filter and centroiding techniques; the median positional uncertainty is ~0.3 arcsecs. The X-ray colors of the detected sources indicate a broad variety of source types, although absorbed AGNs (including some possible Compton-thick sources) are clearly the dominant type. We also match lower significance X-ray sources to optical counterparts and provide a list of 79 optically bright R<~23) lower significance Chandra sources. The majority of these sources appear to be starburst and normal galaxies. We investigate the source-free background, determine the maximum photon-limited exposures, and investigate source confusion. These analyses directly show that Chandra can achieve significantly higher sensitivities in an efficient nearly photon-limited manner and be largely free of source confusion. To allow consistent comparisons, we have also produced point-source catalogs for the ~1 Ms Chandra Deep Field-South (CDF-S). Three hundred and twenty-six (326) X-ray sources are included in the main Chandra catalog, and an additional 42 optically bright X-ray sources are included in a lower significance Chandra catalog. We find good agreement with the photometry of the previously published CDF-S catalogs; however, we provide significantly improved positional accuracy (ABRIDGED).

We present the first results from a wide solid angle, moderately deep {\it Chandra} survey of the Lockman Hole North-West region. Our 9 ACIS-I fields cover an effective solid angle of 0.4 deg$^{2}$ and reach a depth of $3 \times 10^{-16}$ \ergpcmsqps in the 0.4--2 keV band and $3 \times 10^{-15}$ \ergpcmsqps in the 2--8 keV band. The best fit logN-logS for the entire field, the largest contiguous {\it Chandra} field yet observed, matches well onto that of the {\it Chandra} Deep Field North. We show that the full range of the `cosmic variance' previously seen in different {\it Chandra} fields is reproduced in this small region of the sky. Counts-in-cells analysis shows that the hard band sources are more strongly correlated than the soft band sources.

The high angular resolution and sensitivity of the Chandra X-ray Observatory has yielded large numbers of faint X-ray sources with measured redshifts in the soft (0.5-2 keV) and hard (2-8 keV) energy bands. Many of these sources show few obvious optical signatures of active galactic nuclei (AGN). We use Chandra observations of the Hubble Deep Field North region, A370, and the Hawaii Survey Fields SSA13 and SSA22, together with the ROSAT Ultra Deep Survey soft sample and the ASCA Large Sky Survey hard sample, to construct rest-frame 2-8 keV luminosity functions versus redshift for all the X-ray sources, regardless of their optical AGN characteristics. At z=0.1-1 most of the 2-8 keV light density arises in sources with luminosities in the 10^42 erg/s to 10^44 erg/s range. We show that the number density of sources in this luminosity range is rising, or is at least constant, with decreasing redshift. Broad-line AGN are the dominant population at higher luminosities, and these sources show the well-known rapid positive evolution with increasing redshift to z~3. We argue that the dominant supermassive black hole formation has occurred at recent times in objects with low accretion mass flow rates rather than at earlier times in more X-ray luminous objects with high accretion mass flow rates.

Deep Chandra X-ray exposures provide an efficient route for locating optically faint active galactic nuclei (AGN) at high redshifts. We use deep multicolor optical data to search for z>5 AGN in the 2 Ms X-ray exposure of the Chandra Deep Field-North. Of the 423 X-ray sources bright enough (z'<25.2) for a color-color analysis, at most one lies at z=5-6 and none at z>6. The z>5 object is spectroscopically confirmed at z=5.19. Only 31 of the 77 sources with z'>25.2 are undetected in the B or V bands at the 2-sigma level and could lie at z>5. There are too few moderate luminosity AGN at z=5-6.5 to ionize the intergalactic medium.

(abridged) We present multiwavelength observations for a large sample of microjansky radio sources detected in ultradeep 1.4GHz maps centered on the Hubble Deep Field-North (HDF-N) and the Hawaii Survey Fields SSA13 and SSA22. Our spectroscopic redshifts for 169 radio sources reveal a flat median redshift distribution, and these sources are hosted by similarly luminous optical L* galaxies, regardless of redshift. This is a serious concern for radio estimates of the local star formation rate density, as a substantial fraction of the ultraviolet luminosity density is generated by sub-L* galaxies at low redshifts. From our submillimeter measurements for 278 radio sources, we find error-weighted mean 850micron fluxes of 1.72$\pm$0.09 mJy for the total sample, 2.37$\pm$0.13 mJy for the optically-faint (I>23.5) subsample, and 1.04$\pm$0.13 mJy for the optically-bright (I<23.5) subsample. We significantly (>3\sigma) detect in the submillimeter 50 of the radio sources, 38 with I>23.5. Spectroscopic redshifts for three of the I<23.5 submillimeter-detected radio sources are in the range z=1.0-3.4, and all show AGN signatures. Using only the submillimeter mapped regions we find that 69\pm9% of the submillimeter-detected radio population are at I>23.5. We also find that 66\pm7% of the S850>5 mJy (>4\sigma) sources are radio-identified. We find that millimetric redshift estimates at low redshifts are best made with a FIR template intermediate between a Milky Way type galaxy and a starburst galaxy, and at high redshifts with an Arp220 template.

We present the first comprehensive search for submillimeter and radio emission from the host galaxies of twenty well-localized gamma-ray bursts (GRBs). With the exception of a single source, all observations were undertaken months to years after the GRB explosions to ensure negligible contamination from the afterglows. We detect the host galaxy of GRB 000418 in both the sub-mm and radio, and the host galaxy of GRB 000210 only in the sub-mm. These observations, in conjunction with the previous detections of the host galaxies of GRB 980703 and GRB 010222, indicate that about 20% of GRB host galaxies are ultra-luminous and have star formation rates of about 500 M_sun/yr. As an ensemble, the non-detected hosts have a star formation rate of about 100 M_sun/yr (5-sigma) based on their radio emission. The detected and ensemble star formation rates exceed the optical values by an order of magnitude, indicating significant dust obscuration. In the same vein, the ratio of bolometric dust luminosity to UV luminosity for the hosts detected in the sub-mm and radio ranges from 20-800, and follows the known trend of increasing obscuration with increasing bolometric luminosity. We also show that, both as a sample and individually, the GRB host galaxies have bluer R-K colors as compared with galaxies selected in the sub-mm in the same redshift range. This possibly indicates that the stellar populations in the GRB hosts are on average younger, supporting the massive stellar progenitor scenario for GRBs, but it is also possible that GRB hosts are on average less dusty.

We use deep near-infrared and submillimeter observations of three massive lensing cluster fields, A370, A851, and A2390, to determine the average submillimeter properties of a K'-selected sample. The 38 Extremely Red Objects (EROs; I-K'>4) with K'<21.25 have a very significant error-weighted mean 850-micron flux of 1.58 +/- 0.13 mJy. The ERO contribution to the 850-micron background is 1.88 +/- 0.16 x 10^4 mJy per square degree, or about half the background light. The 17 Very Red Objects (VROs; 3.5<I-K'<4) are also significantly detected (1.32 +/- 0.19 mJy), bringing the combined VRO and ERO contribution to 2.59 +/- 0.19 x 10^4 mJy per square degree. There is a substantial systematic uncertainty (about a factor of two) in this value due to field-to-field variation, but even with this uncertainty it is clear that a large fraction of the 850-micron background arises from red objects. An analysis of the VRO and ERO number counts shows that half of the population's 850-micron light arises in objects with demagnified magnitudes K'<20 and half in fainter objects. On the basis of the I-J versus J-K' colors of the galaxies, the bulk of the submillimeter signal appears to arise from the dusty starburst galaxies in the red object population rather than from the high-redshift elliptical galaxies.

We present the optical, near-infrared, submillimeter, and radio follow-up catalog of the X-ray selected sources from the 1 Ms Chandra observation of the Hubble Deep Field North region. We have B, V, R, I, and z' magnitudes for the 370 X-ray point sources, HK' magnitudes for 276, and spectroscopic redshifts for 182. We present high-quality spectra for 175 of these. The redshift distribution shows indications of structures at z=0.843 and z=1.0175 (also detected in optical surveys) which could account for a part of the field-to-field variation seen in the X-ray number counts. The flux contributions separated into unit bins of redshift show that the z<1 spectroscopically identified sources already contribute about one-third of the total flux in both the hard and soft bands. We find from ratios of the X-ray counts that the X-ray spectra are well-described by absorption of an intrinsic Gamma=1.8 power-law, with log NH values ranging from 21 to 23.7. We estimate that the Chandra sources that produce 87% of the HEAO-A X-ray background (XRB) at 3 keV produce 57% at 20 keV, provided that at high energies the spectral shape of the sources continues to be well-described by a Gamma=1.8 power-law. However, when the Chandra contributions are renormalized to the BeppoSAX XRB at 3 keV, the shape matches fairly well the observed XRB at both energies. Thus, whether a substantial population of as-yet undetected Compton-thick sources is required to completely resolve the XRB above 10 keV depends critically on how the currently discrepant XRB measurements in the 1-10 keV energy range tie together with the higher energy XRB. (Abridged)

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.

We present deep 850 micron maps of three massive lensing clusters, A370, A851, and A2390, with well-constrained mass models. Our cluster exposure times are more than 2 to 5 times longer than any other published cluster field observations. We catalog the sources and determine the submillimeter number counts. The counts are best determined in the 0.3 to 2 mJy range where the areas are large enough to provide a significant sample. At 0.3 mJy the cumulative counts are 3.3 (1.3,6.3) 10^4 per square degree, where the upper and lower bounds in the brackets are the 90% confidence range. The surface density at these faint count limits enters the realm of significant overlap with other galaxy populations.The corresponding percentage of the extragalactic background light (EBL) residing in this flux range is about 45-65%, depending on the EBL measurement used. Given that 20-30% of the EBL is resolved at flux densities between 2 and 10 mJy, most of the submillimeter EBL is arising in sources above 0.3 mJy. We also performed a noise analysis to obtain an independent estimate of the counts. The upper bounds on the counts determined from the noise analysis closely match the upper limits obtained from the direct counts. The differential counts from this and other surveys can reasonably be described by the parameterization n(S)=3 10^4/(0.7 + S^3) per square degree per mJy with S in mJy, which also integrates to match the EBL.

We present the 2-8 keV number counts from the 1Ms Chandra observation of the Chandra Deep Field North (CDF-N). We combine these with the number counts from a 78ks exposure of the Hawaii Survey Field SSA22 and with the number counts obtained in independent analyses of the CDF-S and the Hawaii Survey Field SSA13 to determine the number counts from 2.10^(-16) to 10^(-13) ergs/cm2/s. Over this flux range the contribution to the X-ray background is 1.1 10^(-11) ergs/cm2/s/square degree. When the contributions above 10^(-13) ergs/cm2/s from BeppoSAX or ASCA observations are included, the total rises to 1.4 10^(-11) ergs/cm2/s/square degree. However, there appears to be substantial field-to-field variation in the counts in excess of the statistical uncertainties. When the statistical and flux calibration uncertainties (in both the background and source measurements) are taken into account, as much as 0.5 10^(-11) ergs/cm2/s/square degree could still be present in an unresolved component.

An extremely deep X-ray survey (about 1 Ms) of the Hubble Deep Field North and its environs (about 450 arcmin^2) has been performed with the Advanced CCD Imaging Spectrometer on board the Chandra X-ray Observatory. This is one of the two deepest X-ray surveys ever performed; for point sources near the aim point it reaches 0.5-2.0 keV and 2-8 keV flux limits of 3 x 10^{-17} erg/cm^2/s and 2 x 10^{-16} erg/cm^2/s, respectively. Here we provide source catalogs along with details of the observations, data reduction, and technical analysis. Observing conditions, such as background, were excellent for almost all of the exposure. We have detected 370 distinct point sources: 360 in the 0.5-8.0 keV band, 325 in the 0.5-2.0 keV band, 265 in the 2-8 keV band, and 145 in the 4-8 keV band. Two new Chandra sources in the HDF-N itself are reported and discussed. Source positions are accurate to within 0.6-1.7 arcsec (at 90% confidence) depending mainly on the off-axis angle. We also detect two highly significant extended X-ray sources and several other likely extended X-ray sources. We present basic number count results for sources located near the center of the field. Source densities of 7100^{+1100}_{-940} deg^{-2} (at 4.2 x 10^{-17} erg/cm^2/s) and 4200^{+670}_{-580} deg^{-2} (at 3.8 x 10^{-16} erg/cm^2/s) are observed in the soft and hard bands, respectively.

We present submillimeter observations for 136 of the 370 X-ray sources detected in the 1 Ms exposure of the Chandra Deep Field North. Ten of the X-ray sources are significantly detected in the submillimeter. The average X-ray source in the sample has a significant 850 micron flux of 1.69+/-0.27 mJy. This value shows little dependence on the 2-8 keV flux from 5e-16 erg/cm^2/s to 1e-14 erg/cm^2/s. The ensemble of X-ray sources contribute about 10% of the extragalactic background light at 850 microns. The submillimeter excess is found to be strongest in the optically faint X-ray sources that are also seen at 20 cm, which is consistent with these X-ray sources being obscured and at high redshift (z>1).

We describe the optical, near-infrared, and radio properties of a sample of hard (2-7 keV) X-ray sources detected in a deep Chandra observation of the field surrounding the Abell 370 cluster. We combine these data with similar observations of the Chandra Deep Field-North and the Hawaii Survey Field SSA13 to obtain a sample of 69 hard X-ray sources (45 are spectroscopically identified) with extremely deep 20 cm observations. We find that about 4% of the >Lstar galaxy population is X-ray luminous at any time and hence that black hole accretion has a duration of about half a Gyr. We find that about 30% of the summed 2-7 keV flux from our total sample is from sources at z<1. We estimate the bolometric luminosities of accretion onto supermassive black holes and the mass inflow rates. The time history of the accretion rate density is evaluated, and its integrated value is reasonably consistent with the value inferred from the local black hole mass to bulge mass ratio.

We describe progress in a programme to study a sample of sources typical of those which contribute a large fraction of the hard X-ray background. The sources are selected from the fields of approximately 10 ks Chandra cluster observations with follow-up in the near-infrared and optical. The X-ray data indicate that many of these are powerful, obscured sources. Two hard objects which are lensed by the cluster are good candidates for X-ray type-2 quasars, with luminosities of the order 10^{44}-10^{45} erg s^{-1} and obscuring column densities of the order 10^{23} cm^{-2}. We find that the sources are bright in the infra-red (K typically 17-18). From Keck optical spectra and photometric redshifts, we find that the host galaxies are consistent with being early-type, massive hosts at redshifts ranging from 0.2 to 3, with median redshift approximately 1.5. The two obscured quasars also have mid-infrared detections and 850-micron upper-limits implying that the surrounding dust is at warm-to-hot temperatures (T of the order of 1000 K).

We present optical spectra and near-infrared imaging of a sample of serendipitous X-ray sources detected in the field of Chandra observations of the A2390 cluster of galaxies. The sources have 0.5-7 keV fluxes of 0.6-8 x 10^{-14} erg cm^{-2} s^{-1} and lie around the break in the 2-10 keV source counts. They therefore are typical of sources dominating the X-ray Background in that band. Twelve of the fifteen targets for which we have optical spectra show emission lines, most of which have soft X-ray spectra. Including photometric redshifts and published spectra we have redshifts for seventeen of the sources, ranging from z~0.2 up to z~3 with a peak between z=1-2. Ten of our sources have hard X-ray spectra indicating a spectral slope flatter than that of a typical unabsorbed quasar. Two hard sources that are gravitationally lensed by the foreground cluster are obscured quasars, with intrinsic 2-10 keV luminosities of 0.2-3 x 10^{45} erg/s and absorbing columns of N_H>10^{23} cm^{-2}. Both were detected by ISO, showing that the absorbed flux emerges in the far-infrared.

We combine deep Subaru near-infrared images of the massive lensing clusters A2390 and A370 with Keck optical data to map the spectral energy distributions (SEDs) of Chandra X-ray sources lying behind the clusters. The three sources behind A2390 are found to have extremely red colors with SEDs consistent with evolved galaxies at redshifts z>1.4. One source has extremely anomalous colors, which we interpret as evidence for a type Sa SED at a redshift around 2.5. The photometric redshift of another source has been confirmed at z=1.467 from near-infrared spectroscopy using the CISCO spectrograph on Subaru. Mapping of optically faint hard X-ray sources may prove to be an extremely efficient way to locate luminous evolved galaxies at high redshifts.

We present results from a 221.9 ks Chandra exposure of the HDF-N and its vicinity, concentrating on the 8.6' X 8.7' area covered by the Caltech Faint Field Galaxy Redshift Survey (the `Caltech area'). The minimum detectable fluxes in the 0.5-2 keV and 2-8 keV bands are 1.3e-16 cgs and 6.5e-16 cgs, respectively and a total of 82 sources are detected. More than 80% of the extragalactic X-ray background in the 2-8 keV band is resolved. Redshifts are available for 96% of the sources with R<23; the redshift range is 0.1-3.5 with most sources having z < 1.5. Eight of the X-ray sources are located in the HDF-N itself, including two not previously reported. A population of X-ray faint, optically bright, nearby galaxies emerges at soft-band fluxes of ~< 3e-16 cgs. We set the tightest constraints to date on the X-ray emission properties of microJy radio sources, mid-infrared sources detected by ISO, and very red (R-K_s > 5.0) objects. Where both the infrared and the X-ray coverage are deepest, 75% of the X-ray sources are detected by ISO; the high X-ray to infrared matching rate bodes well for future sensitive infrared observations of faint X-ray sources. Four of the 33 very red objects that have been identified in the Caltech area by Hogg et al. (2000) are detected in X-rays; these four are among our hardest Chandra sources, and we argue that they contain moderately luminous obscured AGN. Overall, however, the small Chandra detection fraction suggests a relatively small AGN content in the optically selected very red object population. (Abridged)

With recent Chandra observations, at least 60% of the 2-10 keV background is now resolved into discrete sources. Here we present deep optical, NIR, submm, and 20 cm (radio) images, as well as high-quality optical spectra, of a complete sample of 20 hard X-ray sources in a deep Chandra observation of the SSA13 field. The thirteen I<23.5 galaxies have redshifts in the range 0.1 to 2.6. Two are quasars, five show AGN signatures, and six are z<1.5 luminous bulge-dominated galaxies whose spectra show no obvious optical AGN signatures. The seven spectroscopically unidentified sources have colors that are consistent with evolved early galaxies at z=1.5-3. Only one hard X-ray source is significantly detected in an ultradeep submm map; its millimetric redshift is in the range z=1.2-2.4. None of the remaining 19 sources is detected in the submm. These results probably reflect the fact that the 850-micron flux limits obtainable with SCUBA are quite close to the expected fluxes from obscured AGN. The hard X-ray sources have an average L(FIR)/L(2-10 keV)~60, similar to that of local obscured AGN. The same ratio for a sample of submm selected sources is in excess of 1100, suggesting that their FIR light is primarily produced by star formation. Our data show that luminous hard X-ray sources are common in bulge-dominated optically luminous galaxies. We use our measured bolometric corrections with the 2-10 keV EBL to infer the growth of supermassive black holes. Even with a high radiative efficiency of accretion (e=0.1), the black hole mass density required to account for the observed light is comparable to the local black hole mass density. (Abridged)

Deep submillimetre surveys have successfully detected distant, star-forming galaxies, enshrouded in vast quantities of dust and which emit most of their energy at far infrared wavelengths. These luminous galaxies are an important constituent of the Universal star-formation history, and any complete model of galaxy evolution must account for their existence. Although these sources have been tentatively identified with very faint and sometimes very red optical counterparts, their poorly constrained redshift distribution has made their interpretation unclear. In particular, it was not understood if these galaxies had been missed in previous surveys, or if they constituted a truly new class of objects, undetectable at other wavelengths. By utilizing a radio selection technique, we have isolated a sample of 20 sub-mm objects representative of the 850 micron population brighter than 5 mJy with z<~3. We show that these galaxies are so heavily dust obscured that they remain essentially 'invisible' to ultraviolet selection. Furthermore, relying on the radio-submillimeter flux density ratio, we estimate their redshift distribution, finding a median of two. While not a substitute for the thorough followup of blank field sub-mm surveys, our results do shed light on a substantial portion of the luminous sub-mm population with z<~3.

The origin of the hard (2-10 keV) X-ray background has remained mysterious for over 35 years. Most of the soft (0.5-2 keV) X-ray background has been resolved into discrete sources, which are primarily quasars; however, these sources do not have the flat spectral shape required to match the X-ray background spectrum. Here we report the results of an X-ray survey 30 times more sensitive than previous studies in the hard band and four times more sensitive in the soft band. The sources detected in our survey account for at least 75 per cent of the hard X-ray background. The mean X-ray spectrum of these sources is in good agreement with that of the background. The X-ray emission from the majority of the detected sources is unambiguously associated with either the nuclei of otherwise normal bright galaxies or optically faint sources, which could either be active nuclei of dust enshrouded galaxies or the first quasars at very high redshifts.

Direct submm imaging has recently revealed the 850-micron background to be mostly composed of a population of distant ultraluminous infrared galaxies, but identifying the optical/NIR counterparts to these sources has proved difficult due to the poor submm spatial resolution. However, the proportionality of both cm and submm data to the star formation rate suggests that high resolution radio continuum maps with subarcsecond positional accuracy can be exploited to locate submm sources. In this paper we present results from a targeted SCUBA survey of micro-Jansky radio sources in the flanking fields of the Hubble Deep Field. Even with relatively shallow 850-micron SCUBA observations (>6 mJy at 3-sigma), we were successful at making submm detections of optical/NIR-faint (I>24 and K~21-22) radio sources, and our counts closely match the bright counts from submm surveys. Redshift estimates can be made from the ratio of the submm flux to the radio flux across the 100 GHz break in the spectral energy distribution. This millimetric redshift estimation places the bright submm population at z=1-3 where it forms the high redshift tail of the faint radio population. The star formation rate density (SFRD) due to ultraluminous infrared galaxies increases by more than two orders of magnitude from z~0 to z~1-3. The SFRD at high redshift inferred from our >6 mJy submm observations is comparable to that observed in the UV/optical. (Abridged)

We present extensive observations of a sample of distant, submm galaxies detected in the field of the massive cluster lens, A1835, using SCUBA. Taken in conjunction with earlier observations of other submm-selected sources (Ivison et al. 1998; Smail et al. 1999; Soucail et al. 1999) we now have detailed, multi-wavelength observations of 7 examples of the submm population, having exploited the combination of achromatic amplification by cluster lenses and lavish archival datasets. These sources, all clearly at z > 1, illustrate the wide range in the radio and optical properties of distant submm-selected galaxies. We include detailed observations of the first candidate `pure' starburst submm galaxy at high redshift, a z=2.56 interacting galaxy which shows no obvious sign of hosting an AGN. The remaining sources have varying degrees of inferred AGN activity (3 from 7 of the most luminous show some evidence of the presence of an AGN) although even when an AGN is obviously present it is still not apparent if reprocessed radiation from this source dominates the submm emission. In contrast with the variation in the spectral properties, we see relatively homogeneous morphologies for the population, with a large fraction of merging or interacting systems. Our study shows that virtually identical SEDs are seen for galaxies which exhibit strikingly different opt/UV spectral-line characteristics. We conclude that standard opt/UV spectral classifications are misleading when applied to distant, highly obscured galaxies and that we must seek other means of determining the various contributions to the overall energy budget of submm galaxies and hence to the far-IR extragalactic background.

Using a combination of radio and optical imaging at 0.1\arcs to 0.2\arcs resolution with the VLA/MERLIN and HST has led to a breakthrough in our understanding of radio emission from distant (0.1<z<3) starburst galaxies. We have recently isolated a number of high redshift, dusty starburst galaxies that remain invisible in ground based images to I_{AB}=25 and I_{AB}=28.5 in the Hubble Deep Field. These galaxies appear as faint radio sources, often accompanied by very red counterparts (I-K > 4-6) and submillimeter sources with S > 2 mJy at 850 microns as measured with SCUBA on the JCMT. The far-infrared luminosities of these galaxies exceeds even the most intense starbursts found in the local universe (e.g., Arp 220), suggesting they are in the process of converting the bulk of their gass mass into stars. These galaxies, completely absent in optical surveys, constitute 50% - 90% of the star-formation density in the distant Universe. Given the poor sub-mm resolution (15\arcs) of the SCUBA/JCMT images, we use the 0.2\arcs radio imaging as a surrogate in order to understand the nature of the dominant emission mechanism driving the FIR luminosity (AGN vs. star-formation). Upcoming developments in radio instrumentation (the Expanded VLA Array and Square Kilometer Array) will soon increase sensitivity and resolution orders of magnitude, providing a natural complement to parallel developments in sub-mm facilities (e.g., ALMA). With dual radio continuum and sub-mm surveys of the distant Universe, a census of galaxy evolution to the earliest cosmic epochs (z = 5-30) will soon be possible.

We report on the identifications of radio sources from our survey of the Hubble Deep Field and the SSA13 fields, both of which comprise the deepest radio surveys to date at 1.4 GHz and 8.5 GHz respectively. About 80% of the microjansky radio sources are associated with moderate redshift starburst galaxies or AGNs within the I magnitude range of 17 to 24 with a median of I = 22 mag. Thirty-one (20%) of the radio sources are: 1) fainter than $I>$25 mag, with two objects in the HDF $I_{AB}>$28.5, 2) often identified with very red objects $I-K>$4, and 3) not significantly different in radio properties than the brighter objects. We suggest that most of these objects are associated with heavily obscured starburst galaxies with redshifts between 1 and 3. However, other mechanisms are discussed and cannot be ruled out with the present observations.

We use deep ground-based imaging in the near-IR to search for counterparts to the luminous submm sources in the catalog of Smail et al (1998). For the majority of the submm sources the near-IR imaging supports the counterparts originally selected from deep optical images. However, in two cases (10% of the sample) we find a relatively bright near-IR source close to the submm position, sources that were unidentified in the deep HST and ground-based R-band images used in Smail et al (1998). We place limits on colours of these sources from deep high-resolution Keck II imaging and find they have 2-sigma limits of (I-K)>6.8 and (I-K)>6.0 respectively. Both sources thus class as EROs. Using the spectral properties of the submm source in the radio and submm we argue that these EROs are probably the source of the submm emission, rather than the bright spiral galaxies previously identified by Smail et al. (1998). From the surface density of these submm-bright EROs we suggest that this class accounts for the majority of the reddest members of the ERO population, in good agreement with the preliminary conclusions of pointed submm observations of individual EROs. We conclude that the most extreme EROs represent a population of dusty, ultraluminous galaxies at high redshifts; further study of these will provide insights into the nature of star formation in obscured galaxies in the early Universe. The identification of similar counterparts in blank field submm surveys will be extremely difficult owing to their faintness (K~20.5, I>26.5). Finally, we discuss the radio and submm properties of the two submm-bright EROs discovered here and suggest that both galaxies lie at z>2.

The strength of the submillimeter background light shows directly that much of the energy radiated by star formation and AGN is moved to far infrared wavelengths. However, it is only as this background at 850 microns has been resolved with direct submillimeter imaging that we have seen that it is created by a population of ultraluminous (or near ultraluminous) infrared galaxies (ULIGs) which appear to lie at relatively high redshifts (z>1). Mapping the redshift evolution of this major portion of the universal star formation has been difficult because of the poor submillimeter spatial resolution, but this difficulty can be overcome by using extremely deep cm continuum radio observations to obtain precise astrometric information since the bulk of the brighter submillimeter sources have detectable radio counterparts. However, with this precise position information available, we find that most of the submillimeter sources are extremely faint in the optical and near infrared (I>>24 and K=21-22) and inaccessible to optical spectroscopy. Rough photometric redshift estimates can be made from combined radio and submillimeter energy distributions. We shall refer to this procedure as millimetric redshift estimation to distinguish it from photometric estimators in the optical and near IR. These estimators place the bulk of the submillimeter population at z=1-3, where it corresponds to the high redshift tail of the faint cm radio population. While still preliminary, the results suggest that the submillimeter population appears to dominate the star formation in this redshift range by almost an order of magnitude over the mostly distinct populations selected in the optical-ultraviolet.

Identifying the optical/near-infrared (NIR) counterparts to the distant submillimeter (submm) source population has proved difficult due to poor submm spatial resolution. However, the proportionality of both centimeter and submm data to the star formation rate suggests that high resolution radio continuum maps with subarcsecond positional accuracy could be exploited to locate submm sources. We targeted with SCUBA a sample of micro-Jy radio sources in the flanking fields of the Hubble Deep Field selected from the uniform (8 micro-Jy at 1-sigma) 1.4 GHz VLA image of Richards (1999). We find that the majority of bright (>6 mJy) submm sources have detectable radio counterparts. With the precise positions from the radio, we also find that these submm sources are extremely faint in the optical and NIR (I>>24 and K=21-22) and are therefore inaccessible to optical spectroscopy. Redshift estimates can, however, be made from the shape of the spectral energy distribution in the radio and submm. This procedure, which we refer to as millimetric redshift estimation, places the bright submm population at z=1-3, where it forms the high redshift tail of the faint radio population.

In hierarchical models, where spheroidals are primarily produced via the merger of disk galaxies, the number of intrinsically red systems at faint limits will be substantially lower than that expected in "traditional" models where the bulk of star formation was completed at high redshifts. We analyse the optical--near-infrared colour distribution of a large flux-limited sample of field spheroidal galaxies identified morphologically from archival Hubble Space Telescope data. The $I_{814}-HK'$ colour distribution for a sample jointly limited at $I_{814}<$23 mag and $HK'<$19.5 mag is used to constrain the star formation history of elliptical and S0 galaxies. We detect a significant deficit of intrinsically red spheroids. From limited spectroscopic and deeper photometric redshift data, we demonstrate that this deficit results from a strong evolutionary effect. The space density of spheroidals of all colours at faint limits is lower than expected for a population of constant comoving density. We conclude that field ellipticals cannot have formed all of their stars at high redshift, as has been suggested for the bulk of their clustered counterparts. We compare alternative evolutionary histories with our data, including an analytical representation of recent models based on the hierarchical assembly of ellipticals, which provides a better, though not ideal fit. The accumulation of spectroscopic redshifts for our sample will allow a quantitative discrimination between models with extended star formation.

We have utilized various magnitude-limited samples drawn from an extremely deep and highly complete spectroscopic redshift survey of galaxies observed in seven colors in the Hawaii Survey Fields and the Hubble Deep Field to investigate the evolution of the universal rest-frame ultraviolet luminosity density from z = 1 to the present. The multi-color data (U', B, V, R, I, J, HK') enable the sample selection to be made in the rest-frame ultraviolet for the entire redshift range. Due to the large sample size and depth (U_{AB} = 24.75, B_{AB} = 24.75, I_{AB} = 23.5), we are able to accurately determine the luminosity density to z = 1. We do not confirm the very steep evolution reported by Lilly et al. (1996) but instead find a shallower slope, approximately (1+z)^{1.5} for q0 = 0.5, which would imply that galaxy formation is continuing smoothly to the present time rather than peaking at z = 1. Much of the present formation is taking place in smaller galaxies. Detailed comparisons with other recent determinations of the evolution are presented.

Recent deep blank field submillimeter surveys have revealed a population of luminous high redshift galaxies that emit most of their energy in the submillimeter. The results suggest that much of the star formation at high redshift may be hidden to optical observations. In this paper we present wide-area 850-micron SCUBA data on the Hawaii Survey Fields SSA13, SSA17, and SSA22. Combining these new data with our previous deep field data, we establish the 850-micron galaxy counts from 2 mJy to 10 mJy with a >3-sigma detection limit. The area coverage is 104 square arcmin to 8 mJy and 7.7 square arcmin to 2.3 mJy. The differential 850-micron counts are well described by the function n(S)=N_0/(a+S^3.2), where S is the flux in mJy, N_0=3.0 x 10^4 per square degree per mJy, and a=0.4-1.0 is chosen to match the 850-micron extragalactic background light. Between 20 to 30 per cent of the 850-micron background resides in sources brighter than 2 mJy. Using an empirical fit to our >2 mJy data constrained by the EBL at lower fluxes, we argue that the bulk of the 850-micron extragalactic background light resides in sources with fluxes near 1 mJy. The submillimeter sources are plausible progenitors of the present-day spheroidal population.

We present a Keck II LRIS spectroscopic follow-up study of the possible optical counterparts to a flux-limited sample of galaxies selected from an 850-micron survey of massive lensing clusters using the SCUBA bolometer array on the JCMT. These sources represent a population of luminous dusty galaxies responsible for the bulk of the 850-micron background detected by COBE and thus for a substantial fraction of the total far-infrared emission in the Universe. We present reliable redshifts for 20 galaxies and redshift limits for a further 4 galaxies selected from the error-boxes of 14 submm sources. Two other submm detections in the sample have no obvious optical counterparts, and the final submm source was only identified from imaging data after the completion of our spectroscopic observations. The optical identifications for 4 of the submm sources have been confirmed through either their detection in CO at mm-wavelengths (two pairs of galaxies at z=2.55 and z=2.80) or from the characteristics of their spectral energy distributions (two of the central cD galaxies in the lensing clusters). Plausible arguments based on the optical spectral properties (starburst or AGN signatures) of the counterparts allow us to identify a further two likely counterparts at z=1.06 and 1.16. For the remaining 8 cases, it is not always clear which, if any, of the optical sources identified are the true counterparts. Possible counterparts for these have redshifts ranging from z=0.18 to z=2.11. Working with the current identifications, we suggest that the majority of the extragalactic background light in the submm is emitted by sources at z<3 and hence that the peak activity in highly-obscured sources (both AGN and starbursts) lies at relatively modest redshifts. (Abridged)

We report the detection of CO(3-2) emission from the submillimeter-selected luminous galaxy SMM J14011+0252. The optical counterpart of the submillimeter source has been identified as a merger system with spectral characteristics consistent with a starburst at z=2.565. The CO emission confirms the optical identification of the submillimeter source and implies a molecular gas mass of 5x10^{10}h_{75}^{-2} M(sun), after correcting for a lensing amplification factor of 2.75. The large molecular gas mass and the radio emission are consistent with the starburst interpretation of the source. These results are similar to those found for SMM J02399-0136, which was the first submillimeter selected CO source found at high redshift. The CO detections of these two high-redshift submillimeter galaxies suggest the presence of massive reservoirs of molecular gas which is consistent with the inferred high rates of star-formation (1000 M(sun)/yr). These two systems appear to be associated with merger events which may evolve into present day luminous elliptical galaxies.

We present a detailed analysis of the spectroscopic catalog of galaxies in 10 distant clusters from Dressler et al. (1999, D99). We investigate the nature of the different spectral classes defined by D99 including star forming, post-starburst and passive galaxy populations, and reproduce their basic properties using our spectral synthesis model. We attempt to identify the evolutionary pathways between the various spectral classes in order to search for the progenitors of the numerous post-starburst galaxies. The comparison of the spectra of the distant galaxy populations with samples drawn from the local Universe leads us to identify a significant population of dust-enshrouded starburst galaxies, showing both strong Balmer absorption and relatively modest [OII] emission, that we believe are the most likely progenitors of the post-starburst population. We present the differences between the field and cluster galaxies at z=0.4-0.5. We then compare the spectral and the morphological properties of the distant cluster galaxies, exploring the connection between the quenching of star formation inferred from the spectra and the strong evolution of the S0 population discussed by Dressler et al. (1997). We conclude that either two different timescales and/or two different physical processes are responsible for the spectral and the morphological transformation.

We present the results of an HK' wide-field survey encompassing the Hubble Deep Field and its Flanking Fields. Our wide-field survey provides uniform coverage of an area 61.8 square arcmin to a depth equivalent to K=20.1 at 5-sigma. We have also imaged the Hubble Deep Field in HK', providing uniform coverage of an area 7.8 square arcmin to a depth equivalent to K=21.2 at 5-sigma. Using these data in combination with new deep UH 8K V and I imaging obtained on the CFHT, we find only a small population of objects with colors redder than an equivalent I-K=4, the color expected for an evolved elliptical at z>1. We infer that only a fraction of the local field elliptical population with M_K<-23.4 could have formed in single bursts at high redshift.

The global star formation rate in high redshift galaxies, based on optical surveys, shows a strong peak at a redshift of z=1.5, which implies that we have already seen most of the formation. High redshift galaxies may, however, emit most of their energy at submillimeter wavelengths if they contain substantial amounts of dust. The dust would absorb the starlight and reradiate it as far-infrared light, which would be redshifted to the submillimeter range. Here we report a deep survey of two blank regions of sky performed at submillimeter wavelengths (450 and 850-micron). If the sources we detect in the 850-micron band are powered by star formation, then each must be converting more than 100 solar masses of gas per year into stars, which is larger than the maximum star formation rates inferred for most optically-selected galaxies. The total amount of high redshift star formation is essentially fixed by the level of background light, but where the peak occurs in redshift for the submillimeter is not yet established. However, the background light contribution from only the sources detected at 850-micron is already comparable to that from the optically-selected sources. Establishing the main epoch of star formation will therefore require a combination of optical and submillimeter studies.

We present a thorough quantitative analysis of the evolution of the colour-magnitude relation for early-type galaxies in 17 distant clusters with redshifts 0.31 < z < 1.27 using the Kodama & Arimoto (1997) evolutionary model for elliptical galaxies. The model is calibrated to reproduce the colour- magnitude relation for Coma ellipticals at z ~ 0 and gives the evolution of the slope and zero-point as a function of redshift. We find no significant differences between the colour-magnitude relations of the clusters in our sample. The slopes can be reproduced by a single model sequence in which all elliptical galaxies are assumed to be equally old (the maximum age difference allowed for the brightest 3 magnitudes is only 1 Gyr) and to have mean stellar metallicities which vary as a function of galaxy luminosity. The zero-points of the colour-magnitude relations constrain the epoch of major star formation in early-type galaxies to z_f>2-4. This study provides two important constraints for any model of the formation of rich clusters: the uniformity of the ages of the stellar populations in the early-type galaxies and the universality of the metallicity sequence of these galaxies as a function of galaxy mass.

We investigate various probes of luminosity evolution in the rich cluster environment, concentrating in particular on the spheroidal (E/S0) galaxies, using a newly-constructed catalog of morphologically-classified Hubble Space Telescope images of galaxies observed in the cores of 13 clusters with redshifts 0.17<z<0.56. An important distinction of this study compared to earlier work is the availability of new near-infrared ground-based photometry for a substantial subset of our HST fields, which we have used to select and study the various populations. We find no significant change in the characteristic luminosity of the spheroidal populations at redshifts of 0.31 and 0.56. As a more sensitive probe of luminosity evolution, deriving effective metric radii and mean effective surface brightnesses. At the standard condition corresponding to R_e = 1 kpc, we find convincing evidence of evolutionary brightening in both rest-frame B and K light, consistent with model predictions based on the passive evolution of stellar populations.

We present new results of a program to study the detailed morphologies of galaxies in intermediate redshift clusters and hence understand the physical origin of the enhanced star formation in these environments at earlier epochs. Deep, high resolution imagery has been obtained of 3 rich clusters, AC103,AC118 & AC114 at z=0.31, using the WFPC2 on HST. For AC103 & AC118, single pointings covering a central 0.5x0.5Mpc have been obtained; for AC114, 4 pointings covering a 1.2x0.7Mpc area have been obtained, allowing the more outer regions of a distant cluster to be studied. Hubble types plus evidence for dynamical interactions and/or structural abnormalities have been determined visually for all galaxies down to R=22.25 in AC103 & AC118 and R=23.0 in AC114.We find the numbers of spirals (Sa-Sdm) in our clusters to be up to 4x higher than that seen in present-day clusters; only in the virialised core of our most massive regular cluster, AC114, do we see morphological fractions approaching those of the present epoch. Dynamical interactions are also widespread throughout the clusters with ~20% of members showing morphological evidence of this phenomenon. The highest incidence is amongst those blue members either undergoing a starburst or seen <0.5Gyr thereafter with the majority being involved in major mergers. These galaxies, however, are of modest luminosity (L~L*+1mag), destined to become dwarfs once they fade. Cluster members with ongoing star formation typical of nearby spirals are largely late Sb-Sdm Hubble types. Galaxies whose star formation ended 1-2Gyr prior to the epoch of observation are all normal early-type (S0-Sb) disk systems. The red H-delta strong galaxies are a mixture of E & S0 galaxies, consistent with them being old dormant systems which have undergone a secondary star formation episode.

Using traditional morphological classifications of galaxies in 10 intermediate-redshift (z~0.5) clusters observed with WFPC-2 on the Hubble Space Telescope, we derive relations between morphology and local galaxy density similar to that found by Dressler for low-redshift clusters. Taken collectively, the `morphology-density' relationship, M-D, for these more distant, presumably younger clusters is qualitatively similar to that found for the local sample, but a detailed comparison shows two substantial differences: (1) For the clusters in our sample, the M-D relation is strong in centrally concentrated ``regular'' clusters, those with a strong correlation of radius and surface density, but nearly absent for clusters that are less concentrated and irregular, in contrast to the situation for low redshift clusters where a strong relation has been found for both. (2) In every cluster the fraction of elliptical galaxies is as large or larger than in low-redshift clusters, but the S0 fraction is 2-3 times smaller, with a proportional increase of the spiral fraction. Straightforward, though probably not unique, interpretations of these observations are (1) morphological segregation proceeds hierarchically, affecting richer, denser groups of galaxies earlier, and (2) the formation of elliptical galaxies predates the formation of rich clusters, and occurs instead in the loose-group phase or even earlier, but S0's are generated in large numbers only after cluster virialization.

We analyse the detailed distribution of star-forming and post-starburst members in three distant (z = 0.31) galaxy clusters in terms of evolutionary sequences that incorporate secondary bursts of star formation on pre-existing stellar populations. Using the number density of spectroscopically-confirmed members on the EW(H\delta) versus B-R plane from existing data, and for a larger K'-limited sample on the U-I versus I-K' plane from newly-acquired infrared images, we demonstrate that the proportion of cluster members undergoing secondary bursts of star formation during the last ~2 Gyr prior to the epoch of observation is probably as high as 30 per cent of the member galaxies. A key observation leading to this conclusion is the high proportion of H\delta strong galaxies in all three clusters. The evolutionary modelling, whilst necessarily approximate, returns the correct proportions of galaxies in various stages of the star formation cycle both in terms of spectral and colour properties. HST images for the three clusters indicate a high proportion of the active members show signs of interaction, whereas the H\delta strong galaxies appear mainly to be regular spheroidals. We examine results from recent merger simulations in the context of the populations in these clusters and confirm that the merging of individual galaxies, triggered perhaps by the hierarchical assembly of rich clusters at this epoch, is consistent with the star formation cycle identified in our data. The implications of such a high fraction of active objects in cluster cores is briefly discussed.