This paper describes the acquisition and analysis of imaging data for the identification of galaxies associated with z~4 damped Lya systems. We present deep BRI images of three fields known to contain four z~4 damped systems. We discuss the reduction and calibration of the data, detail the color criteria used to identify z~4 galaxies, and present a photometric redshift analysis to complement the color selection. We have found no galaxy candidates closer to the QSO than 7'' which could be responsible for the damped Lya systems. Assuming that at least one of the galaxies is not directly beneath the QSO, we set an upper limit on this damped Lya system of L < L*/4. Finally, we have established a web site to release these imaging data to the public.

We present a calculation of the systematic component of the error budget in the photometric redshift technique. We make use of it to describe a simple technique that allows for the assignation of confidence limits to redshift measurements obtained through photometric methods. We show that our technique, through the calculation of a redshift probability function, gives complete information on the probable redshift of an object and its associated confidence intervals. This information can and must be used in the calculation of any observable quantity which makes use of the redshift.

We address the effects of cosmological surface brightness dimming on observations of faint galaxies by examining the distribution of "unobscured" star formation rate intensities versus redshift. We use the star formation rate intensity distribution function to assess the ultraviolet luminosity density versus redshift, based on our photometry and photometric redshift measurements of faint galaxies in the HDF and the HDF--S WFPC2 and NICMOS fields. We find that (1) previous measurements have missed a dominant fraction of the ultraviolet luminosity density of the universe at high redshifts by neglecting cosmological surface brightness dimming effects, which are important at redshifts larger than z = 2, (2) the incidence of the highest intensity star forming regions increases monotonically with redshift, and (3) the ultraviolet luminosity density plausibly increases monotonically with redshift through the highest redshifts observed. By measuring the spectrum of the luminosity density versus redshift, we also find that (4) previous measurements of the ultraviolet luminosity density at redshifts z < 2 must be reduced by a factor 2 to allow for the spectrum of the luminosity density between rest-frame wavelengths 1500 and 2800 A. And by comparing with observations of high-redshift damped Lyman-alpha absorption systems detected toward background QSOs, we further find that (5) the distribution of star formation rate intensities matches the distribution of neutral hydrogen column densities at redshifts z = 2 through 5, which establishes a quantitative connection between high-redshift galaxies and high column density gas and suggests that high-redshift damped Lyman-alpha absorption systems trace lower star formation rate intensity regions of the same galaxies detected in star light in the HDF and HDF--S.

We report the first observations of the clustering environment of damped Lyman alpha absorption systems at z=4. Color selection and photometric redshifts were used to select 44 candidate Lyman-break galaxies brighter than I_AB=25.5 from deep BRI images of the 35 sq. arcmin field containing the quasar BR 0951-04. Multislit spectroscopy of 35 candidate galaxies was performed and 8 of these candidates have been confirmed as z>3.5 Lyman-break galaxies. With only BRI photometry, the photometric redshifts are quite accurate for the spectroscopically confirmed galaxies but have a high rate of misclassification due to color degeneracies between Lyman-break galaxies and low-redshift ellipticals. Both of the z>3.5 galaxies found within 15'' of the quasar line-of-sight appear to be causing absorption systems in the quasar spectrum. We use a battery of statistical tests to look for clustering in the redshift histogram of the z>3.5 galaxies but do not find measurable clustering of these Lyman-break galaxies with the damped Lyman alpha absorbers. With a larger sample of galaxies, our method should determine the cross-correlation between these objects, which probes the bias and hence the mass of the damped Lyman alpha absorbers.

We present a comparison between the spectroscopic catalog of the HDF recently published by Cohen et al and the redshifts measured by our group for those objects using photometric techniques. We perform this comparison in order to characterize the errors associated to the photometric redshift technique. The spectroscopic sample includes over 140 objects in the HDF proper, representing the deepest, cleanest, most complete spectroscopic catalog ever compiled. We study each object for which our redshift and the one measured by Cohen et al disagree. In most cases the photometric evidence is strong enough to call for a careful review of the spectroscopic values, as they seem to be in error. It is possible to characterize the systematic errors associated to our technique, which when combined with the photometric errors allow us to obtain complete information on the redshift of each galaxy and its associated confidence interval, regardless of apparent magnitude. One of the main conclusions of this study is that, to date, all the redshifts from our published catalog that have been checked have been shown to be correct (within the stated confidence limits). This implies that our galaxy template set is a fair representation of the galaxy population at all redshifts (0<z<6) and magnitudes (R<24) explored to date. On the other hand, spectroscopy of faint sources is subject to unknown and uncharacterized systematic errors, that will in turn be transmitted to any photometric redshift technique which uses spectroscopic samples as calibration. Our analysis proves that photometric redshift techniques can and must be used to extend the range of applicability of the spectroscopic redshift measurements. (Abridged)

Observations of distant galaxies are important both for understanding how galaxies form and for probing the physical conditions of the universe at the earliest epochs. It is, however, extremely difficult to identify galaxies at redshift z>5, because these galaxies are faint and exhibit few spectral features. In a previous work, we presented observations that supported the identification of a galaxy at redshift z = 6.68 in a deep STIS field. Here we present new ground-based photometry of the galaxy. We find that the galaxy exhibits moderate detections of flux in the optical B and V images that are inconsistent with the expected absence of flux at wavelength shortward of the redshifted Lyman-alpha emission line of a galaxy at redshift z>5. In addition, the new broad-band imaging data not only show flux measurements of this galaxy that are incompatible with the previous STIS measurement, but also suggest a peculiar spectral energy distribution that cannot be fit with any galaxy spectral template at any redshift. We therefore conclude that the redshift identification of this galaxy remains undetermined.

We present a catalog of photometry and photometric redshifts of 335 faint objects in the HDF-S NICMOS field. The analysis is based on (1) infrared images obtained with the Hubble Space Telescope (HST) using the Near Infrared Camera and Multi-Object Spectrograph (NICMOS) with the F110W, F160W, and F222M filters, (2) an optical image obtained with HST using the Space Telescope Imaging Spectrograph (STIS) with no filter, and (3) optical images obtained with the European Southern Observatory (ESO) Very Large Telescope (VLT) with U, B, V, R, and I filters. The primary utility of the catalog of photometric redshifts is as a survey of faint galaxies detected in the NICMOS F160W and F222M images. The sensitivity of the survey varies significantly with position, reaching a limiting depth of AB(16,000) ~ 28.7 and covering 1.01 arcmin^2 to AB(16,000) = 27 and 1.05 arcmin^2 to AB(16,000) = 26.5. The catalog of photometric redshifts identifies 21 galaxies (or 6% of the total) of redshift z > 5, 8 galaxies (or 2% of the total) of redshift z > 10, and 11 galaxies (or 3% of the total) of best-fit spectral type E/S0, of which 5 galaxies (or 1% of the total) are of redshift z > 1.

We report on some aspects of our efforts to establish properties of the extremely faint galaxy population by applying our photometric redshift technique to the HDF and HDF-S WFPC2 and NICMOS fields. We find that cosmological surface brightness dimming effects play a dominant role in setting what is observed at redshifts z > 2, that the comoving number density of high intrinsic surface brightness regions increases monotonically with increasing redshift, and that previous estimates neglect a significant or dominant fraction of the ultraviolet luminosity density of the universe due to surface brightness effects. The ultraviolet luminosity density of the universe plausibly increases monotonically with increasing redshift to redshifts beyond z = 5.

We report on some aspects of the current status of our efforts to establish properties of faint galaxies by applying our photometric redshift technique to faint galaxies in the HDF and HDF-S WFPC2 and NICMOS fields.

We present an electronic catalog of infrared and optical photometry and photometric redshifts of 323 galaxies in the Hubble Deep Field South NICMOS field at http://www.ess.sunysb.edu/astro/hdfs/home.html. The analysis is based on infrared images obtained with the Hubble Space Telescope using the Near Infrared Camera and Multi-Object Spectrograph and the Space Telescope Imaging Spectrograph together with optical images obtained with the Very Large Telescope. The infrared and optical photometry is measured by means of a new quasi-optimal photometric technique that fits model spatial profiles of the galaxies determined by Pixon image reconstruction techniques to the images. In comparison with conventional methods, the new technique provides higher signal-to-noise-ratio measurements and accounts for uncertainty correlations between nearby, overlapping neighbors. The photometric redshifts are measured by means of our redshift likelihood technique, incorporating six spectrophotometric templates which, by comparison with spectroscopic redshifts of galaxies identified in the Hubble Deep Field North, are known to provide redshift measurements accurate to within an RMS relative uncertainty of (Delta z)/(1 + z) < 0.1 at all redshifts z < 6. The analysis reaches a peak H-band sensitivity threshold of AB(16000) = 28.3 and covers 1.02 acrmin^2 to AB(16000) = 27, 1.27 arcmin^2 to AB(16000) = 26, and 1.44 arcmin^2 to AB(16000) = 25. The analysis identifies galaxies at redshifts ranging from z near 0 through z greater than 10, including 17 galaxies of redshift 5 < z < 10 and five candidate galaxies of redshift z > 10.

Here we describe our attempts to establish statistically complete samples of very high redshift galaxies by obtaining photometric redshifts of galaxies in Medium Deep Survey (MDS) fields and photometric and spectroscopic redshifts of galaxies in very deep STIS slitless spectroscopy fields. On the basis of this analysis, we have identified galaxies of redshift z = 4.92 in an MDS field and of redshift z = 6.68 in a very deep STIS field.

We present new statistical study on the origin of Lyman-alpha absorption systems that are detected in QSO fields. The primary objective is to investigate environmental dependence of the relationship between galaxies and Lyman-alpha absorption systems. We find that galaxies exhibit no environmental dependence in giving rise to Lyman-alpha absorbers in their halos. Conversely, Lyman-alpha absorbers do not exhibit preference as to the environment where they reside.

We report new HST and ground-based observations of a damped Lyman-alpha absorption system toward the QSO 0850+4400. The redshift of the absorption system is z = 0.163770 and the neutral hydrogen column density of the absorption system is log N = 19.81 cm**-2. The absorption system is by far the lowest redshift confirmed damped Lyman-alpha absorption system yet identified, which provides an unprecedented opportunity to examine the nature, impact geometry, and kinematics of the absorbing galaxy in great detail. The observations indicate that the absorption system is remarkable in three respects: First, the absorption system is characterized by weak metal absorption lines and a low metal abundance, possibly less than 4% of the solar metal abundance. This cannot be explained as a consequence of dust, because the neutral hydrogen column density of the absorption system is far too low for obscuration by dust to introduce any significant selection effects. Second, the absorption system is associated with a moderate-luminosity early-type S0 galaxy, although the absorption may actually arise in one of several very faint galaxies detected very close to the QSO line of sight. Third, the absorbing material moves counter to the rotating galaxy disk, which rules out the possibility that the absorption arises in a thin or thick co-rotating gaseous disk. These results run contrary to the expectation that low-redshift damped Lyman-alpha absorption systems generally arise in the gas- and metal-rich inner parts of late-type spiral galaxies. We suggest instead that mounting evidence indicates that low-redshift galaxies of a variety of morphological types may contain significant quantities of low metal abundance gas at large galactocentric distances.