Gas surrounding high redshift galaxies has been studied through observations of absorption line systems toward background quasars for decades. However, it has proven difficult to identify and characterize the galaxies associated with these absorbers due to the intrinsic faintness of the galaxies compared to the quasars at optical wavelengths. Utilizing the Atacama Large Millimeter/Submillimeter Array, we report on detections of [CII] 158 micron line and dust continuum emission from two galaxies associated with two such absorbers at a redshift of z~4. Our results indicate that the hosts of these high-metallicity absorbers have physical properties similar to massive star-forming galaxies and are embedded in enriched neutral hydrogen gas reservoirs that extend well beyond the star-forming interstellar medium of these galaxies.

The far-infrared [C II] 158 micrometer fine structure transition is considered to be a dominant coolant in the interstellar medium. For this reason, under the assumption of a thermal steady state, it may be used to infer the heating rate and, in turn, the star formation rate in local, as well as in high redshift systems. In this work, radio and ultraviolet observations of the Galactic interstellar medium are used to understand whether C II is indeed a good tracer of the star formation rate. For a sample of high Galactic latitude sightlines, direct measurements of the temperature indicate the presence of C II in both the cold and the warm phases of the diffuse interstellar gas. The cold gas fraction (~ 10 - 50% of the total neutral gas column density) is not negligible even at high Galactic latitude. It is shown that, to correctly estimate the star formation rate, C II cooling in both the phases should hence be considered. The simple assumption, that the [C II] line originates only from either the cold or the warm phase, significantly underpredicts or overpredicts the star formation rate, respectively. These results are particularly important in the context of the Damped Lyman-alpha systems for which a similar method is often used to estimate the star formation rate. The derived star formation rates in such cases may not be reliable if the temperature of the gas under consideration is not constrained independently.

We present HIRES observations for 30 damped Lyman alpha systems, selected on the basis of their large metal column densities from previous, lower resolution data. The measured metal column densities for Fe, Zn, S, Si, Cr, Mn, and Ni are provided for these 30 systems. Combined with previously observed large metal column density damped Lyman alpha systems, we present a sample of 44 damped Lyman alpha systems observed with high resolution spectrographs (R~30000). These damped Lyman alpha systems probe the most chemically evolved systems at redshifts greater than 1.5. We discuss the context of our sample with the general damped Lyman alpha population, demonstrating that we are probing the top 10% of metal column densities with our sample. In a companion paper, we will present an analysis of the sample's elemental abundances in the context of galactic chemical enrichment.

A new method is used to measure the physical conditions of the gas in damped Lyman-alpha systems (DLAs). Using high resolution absorption spectra of a sample of 80 DLAs, we are able to measure the ratio of the upper to lower fine-structure levels of the ground state of C II and Si II. These ratios are determined solely by the physical conditions of the gas. We explore the allowed physical parameter space using a Monte Carlo Markov Chain method to constrain simultaneously the temperature, neutral hydrogen density, and electron density of each DLA. The results indicate that at least 5 % of all DLAs have the bulk of their gas in a dense, cold phase with typical densities of ~100 cm-3 and temperatures below 500 K. We further find that the typical pressure of DLAs in our sample is log(P/k) = 3.4 [K cm-3], which is comparable to the pressure of the local interstellar medium (ISM), and that the components containing the bulk of the neutral gas can be quite small with absorption sizes as small as a few parsec. We show that the majority of the systems are consistent with having densities significantly higher than expected from a purely canonical WNM, indicating that significant quantities of dense gas (i.e. n_H > 0.1 cm-3) are required to match observations. Finally, we identify 8 systems with positive detections of Si II*. These systems have pressures (P/k) in excess of 20000 K cm-3, which suggest that these systems tag a highly turbulent ISM in young, star-forming galaxies.

Damped Ly-alpha absorbers (DLAs) are a well-studied class of absorption line systems, and yet the properties of their host galaxies remain largely unknown. To investigate the origin of these systems, we have conducted an imaging survey of 32 quasar fields with intervening DLAs between z~1.9-3.8, leveraging a technique that allows us to image galaxies at any small angular separation from the background quasars. In this paper, we present the properties of the targeted DLA sample, new imaging observations of the quasar fields, and the analysis of new and archival spectra of the background quasars.

We present evidence that the cosmological mean metallicity of neutral atomic hydrogen gas shows a sudden decrease at z>4.7 down to <Z>=-2.03(+0.09-0.11), which is 6 sigma deviant from that predicted by a linear fit to the data at lower redshifts. This measurement is made possible by the chemical abundance measurements of 8 new damped Ly-alpha (DLA) systems at z>4.7 observed with the Echellette Spectrograph and Imager on the Keck II telescope, doubling the number of measurements at z>4.7 to 16. Possible explanations for this sudden decrease in metallicity include a change in the physical processes that enrich the neutral gas within disks, or an increase of the covering factor of neutral gas outside disks due to a lower ultra-violet radiation field and higher density at high redshift. The later possibility would result in a new population of presumably lower metallicity DLAs, with an increased contribution to the DLA population at higher redshifts resulting in a reduced mean metallicity. Furthermore, we provide evidence of a possible decrease at z>4.7 in the comoving metal mass density of DLAs, which is flat out to z~4.3. Such a decrease is expected, as otherwise most of the metals from star-forming galaxies would reside in DLAs by z~6. While the metallicity is decreasing at high redshift, the contribution of DLAs to the total metal budget of the universe increases with redshift, with DLAs at z~4.3 accounting for ~20% as many metals as produced by Lyman break galaxies.

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%.

We present the first Keck/OSIRIS infrared IFU observations of a high redshift damped Lyman-alpha (DLA) galaxy detected in the line of sight to a background quasar. By utilizing the Laser Guide Star Adaptive Optics (LGSAO) to reduce the quasar PSF to FWHM~0.15 arcsec, we were able to search for and map the foreground DLA emission free from the quasar contamination. We present maps of the H-alpha and [OIII] $\lambda \lambda$ 5007, 4959 emission of DLA 2222-0946 at a redshift of z ~ 2.35. From the composite spectrum over the H-alpha emission region we measure a star formation rate of 9.5 $\pm$ 1.0 M$_{\odot}$ year$^{-1}$ and a dynamical mass, M$_{dyn}$ = 6.1 x 10$^9$ M$_{\odot}$. The average star formation rate surface density is < \Sigma_{SFR} > = 0.55 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$, with a central peak of 1.7 M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$. Using the standard Kennicutt-Schmidt relation, this corresponds to a gas mass surface density of $\Sigma_{gas}$ = 243 M$_{\odot}$ pc$^{-2}$. Integrating over the size of the galaxy we find a total gas mass of M$_{gas}$ = 4.2 x 10$^9$ M$_{\odot}$. We estimate the gas fraction of DLA 2222-0946 to be $f_{gas}$ ~ 40%. We detect [NII]$\lambda$6583 emission at 2.5 sigma significance with a flux corresponding to a metallicity of 75% solar. Comparing this metallicity with that derived from the low-ion absorption gas ~6 kpc away, ~30% solar, indicates possible evidence for a metallicity gradient or enriched in/outflow of gas. Kinematically, both H-alpha and [OIII] emission show relatively constant velocity fields over the central galactic region. While we detect some red and blueshifted clumps of emission, they do not correspond with rotational signatures that support an edge-on disk interpretation.

We present an overview of a 90-orbit Hubble Space Telescope treasury program to obtain near ultraviolet imaging of the Hubble Ultra Deep Field using the Wide Field Camera 3 UVIS detector with the F225W, F275W, and F336W filters. This survey is designed to: (i) Investigate the episode of peak star formation activity in galaxies at 1<z<2.5; (ii) Probe the evolution of massive galaxies by resolving sub-galactic units (clumps); (iii) Examine the escape fraction of ionizing radiation from galaxies at z~2-3; (iv) Greatly improve the reliability of photometric redshift estimates; and (v) Measure the star formation rate efficiency of neutral atomic-dominated hydrogen gas at z~1-3. In this overview paper, we describe the survey details and data reduction challenges, including both the necessity of specialized calibrations and the effects of charge transfer inefficiency. We provide a stark demonstration of the effects of charge transfer inefficiency on resultant data products, which when uncorrected, result in uncertain photometry, elongation of morphology in the readout direction, and loss of faint sources far from the readout. We agree with the STScI recommendation that future UVIS observations that require very sensitive measurements use the instrument's capability to add background light through a "post-flash". Preliminary results on number counts of UV-selected galaxies and morphology of galaxies at z~1 are presented. We find that the number density of UV dropouts at redshifts 1.7, 2.1, and 2.7 is largely consistent with the number predicted by published luminosity functions. We also confirm that the image mosaics have sufficient sensitivity and resolution to support the analysis of the evolution of star-forming clumps, reaching 28-29th magnitude depth at 5 sigma in a 0.2 arcsecond radius aperture depending on filter and observing epoch.

Using a sample of 100 H I - selected damped Lyman alpha (DLA) systems, observed with the High Resolution Echelle Spectrometer on the Keck I telescope, we present evidence that the scatter in the well-studied correlation between the redshift and metallicity of a DLA is largely due to the existence of a mass-metallicity relationship at each redshift. To describe the fundamental relations that exist between redshift, metallicity and mass, we use a fundamental plane description, which is described by the following equation: [M/H]=(-1.9+-0.5)+(0.74+-0.21)logdv_90-(0.32+-0.06)z. Here, we assert that the velocity width, dv_90, which is defined as the velocity interval containing 90% of the integrated optical depth, traces the mass of the underlying dark matter halo. This description provides two significant improvements over the individual descriptions of the mass-metallicity correlation and metallicity-redshift correlation. Firstly, the fundamental equation reduces the scatter around both relationships by about 20%, providing a more stringent constraint on numerical simulations modeling DLAs. Secondly, it confirms that the dark matter halos that host DLAs satisfy a mass-metallicity relationship at each redshift between redshifts 2 through 5.

New optical spectra obtained with the SMARTS 1.5m telescope and archival IUE far-ultraviolet spectra of the nova-like variable KQ Mon are discussed. The optical spectra reveal Balmer lines in absorption as well as He I absorption superposed on a blue continuum. The 2011 optical spectrum is similar to the KPNO 2.1m IIDS spectrum we obtained 33 years earlier, except that the Balmer and He I absorption is stronger in 2011. Far-ultraviolet IUE spectra reveal deep absorption lines due to C II, Si III, Si IV, C IV, and He II, but no P Cygni profiles indicative of wind outflow. We present the results of the first synthetic spectral analysis of the IUE archival spectra of KQ Mon with realistic optically thick, steady-state, viscous accretion-disk models with vertical structure and high-gravity photosphere models. We find that the photosphere of the white dwarf (WD) contributes very little FUV flux to the spectrum and is overwhelmed by the accretion light of a steady disk. Disk models corresponding to a WD mass of ~0.6 Msun, with an accretion rate of order 10^{-9} Msun/yr and disk inclinations between 60 and 75 deg, yield distances from the normalization in the range of 144 to 165 pc. KQ Mon is discussed with respect to other nova-like variables. Its spectroscopic similarity to the FUV spectra of three definite SW Sex stars suggests that it is likely a member of the SW Sex class and lends support to the possibility that the WD is magnetic.

We use sensitive observations of three high redshift sources; [CII] fine structure and CO(2-1) rotational transitions for the z=6.4 Quasar host galaxy (QSO) J1148+5251, and [CII] and CO(5-4) transitions from the QSO BR1202-0725 and its sub-millimeter companion (SMG) galaxy at z=4.7. We use these observations to place constraints on the quantity Dz = z(CO) - z(CII) for each source where z(CO) and z(CII) are the observed redshifts of the CO rotational transition and [CII] fine structure transition respectively, using a combination of approaches; 1) By modelling the emission line profiles using `shapelets' to compare both the emission redshifts and the line profiles themselves, in order to make inferences about the intrinsic velocity differences between the molecular and atomic gas, and 2) By performing a marginalisation over all model parameters in order to calculate a non-parametric estimate of Dz. We derive 99% confidence intervals for the marginalised posterior of Dz of (-1.9 pm 1.3) x10^-3, (-3 pm 8) x10^-4 and (-2 pm 4) x10^-3 for J1148+5251, and the BR1202-0725 QSO and SMG respectively. We show the [CII] and CO(2-1) line profiles for J1148+5251 are consistent with each other within the limits of the data, whilst the [CII] and CO(5-4) line profiles from the BR1202-0725 QSO and SMG respectively have 65 and >99.9% probabilities of being inconsistent, with the CO(5-4) lines ~ 30% wider than the [CII] lines. Therefore whilst the observed values of Dz can correspond to variations in the quantity Delta F/F with cosmic time, where F=alpha^2/mu, with alpha the fine structure constant, and mu the proton-to-electron mass ratio, of both (-3.3 pm 2.3) x10^-4 for a look back time of 12.9 Gyr and of (-5 pm 15) x10^-5 for a look back time of 12.4 Gyr we propose that they are the result of the two species of gas being spatially separated as indicated by the inconsistencies in their line profiles.

We present chemical abundance measurements for 47 damped Lyman-alpha systems (DLAs), 30 at z>4, observed with the Echellette Spectrograph and Imager and the High Resolution Echelle Spectrometer on the Keck telescopes. HI column densities of the DLAs are measured with Voigt profile fits to the Lyman-alpha profiles, and we find an increased number of false DLA identifications with SDSS at z>4 due to the increased density of the Lyman-alpha forest. Ionic column densities are determined using the apparent optical depth method, and we combine our new metallicity measurements with 195 from previous surveys to determine the evolution of the cosmic metallicity of neutral gas. We find the metallicity of DLAs decreases with increasing redshift, improving the significance of the trend and extending it to higher redshifts, with a linear fit of -0.22+-0.03 dex per unit redshift from z=0.09-5.06. The metallicity 'floor' of ~1/600 solar continues out to z~5, despite our sensitivity for finding DLAs with much lower metallicities. However, this floor is not statistically different from a steep tail to the distribution. We also find that the intrinsic scatter of metallicity among DLAs of ~0.5 dex continues out to z~5. In addition, the metallicity distribution and the alpha/Fe ratios of z>2 DLAs are consistent with being drawn from the same parent population with those of halo stars. It is therefore possible that the halo stars in the Milky Way formed out of gas that commonly exhibits DLA absorption at z>2.

By comparing the widths of absorption lines from OI, SiII and FeII in the redshift z=2.076 single-component damped Lyman alpha absorption system in the spectrum of Q2206-199 we establish that these absorption lines arise in Warm Neutral Medium gas at ~12000 +/- 3000K. This is consistent with thermal equilibrium model estimates of ~ 8000K for the Warm Neutral Medium in galaxies, but not with the presence of a significant cold component. It is also consistent with, but not required by, the absence of CII* fine structure absorption in this system. Some possible implications concerning abundance estimates in narrow-line WNM absorbers are discussed.

Supernova (SN) 2008ax in NGC 4490 was discovered within hours after shock breakout, presenting the rare opportunity to study a core-collapse SN beginning with the initial envelope-cooling phase immediately following shock breakout. We present an extensive sequence of optical and near-infrared spectra, as well as three epochs of optical spectropolarimetry. Our initial spectra, taken two days after shock breakout, are dominated by hydrogen Balmer lines at high velocity. However, by maximum light, He I lines dominated the optical and near-infrared spectra, which closely resembled those of normal Type Ib supernovae (SNe Ib) such as SN 1999ex. This spectroscopic transition defines Type IIb supernovae, but the strong similarity of SN 2008ax to normal SNe Ib beginning near maximum light, including an absorption feature near 6270A due to H-alpha at high velocities, suggests that many objects classified as SNe Ib in the literature may have ejected similar amounts of hydrogen as SN 2008ax, roughly a few x 0.01 M_sun. Early-time spectropolarimetry (6 and 9 days after shock breakout) revealed strong line polarization modulations of 3.4% across H-alpha, indicating the presence of large asphericities in the outer ejecta. The continuum shares a common polarization angle with the hydrogen, helium, and oxygen lines, while the calcium and iron absorptions are oriented at different angles. This is clear evidence of deviations from axisymmetry even in the outer ejecta. Intrinsic continuum polarization of 0.64% only nine days after shock breakout shows that the outer layers of the ejecta were quite aspherical. A single epoch of late-time spectropolarimetry, as well as the shapes of the nebular line profiles, demonstrate that asphericities extended from the outermost layers all the way down to the center of this SN. [Abridged]

We present evidence for spatially extended low surface brightness emission around Lyman break galaxies (LBGs) in the V band image of the Hubble Ultra Deep Field, corresponding to the z~3 rest-frame FUV light, which is a sensitive measure of star formation rates (SFRs). We find that the covering fraction of molecular gas at z~3 is not adequate to explain the emission in the outskirts of LBGs, while the covering fraction of neutral atomic-dominated hydrogen gas at high redshift is sufficient. We develop a theoretical framework to connect this emission around LBGs to the expected emission from neutral H I gas i.e., damped Lyman alpha systems (DLAs), using the Kennicutt-Schmidt (KS) relation. Working under the hypothesis that the observed FUV emission in the outskirts of LBGs is from in situ star formation in atomic-dominated hydrogen gas, the results suggest that the SFR efficiency in such gas at z~3 is between factors of 10 and 50 lower than predictions based on the local KS relation. The total SFR density in atomic-dominated gas at z~3 is constrained to be ~10% of that observed from the inner regions of LBGs. In addition, the metals produced by in situ star formation in the outskirts of LBGs yield metallicities comparable to those of DLAs, which is a possible solution to the "Missing Metals" problem for DLAs. Finally, the atomic-dominated gas in the outskirts of galaxies at both high and low redshifts has similar reduced SFR efficiencies and is consistent with the same power law.

A re-analysis of the Stokes-parameter spectra obtained of the z=0.692 21 cm absorption line toward 3C 286 shows that our original claimed detection of Zeeman splitting by a line-of-sight magnetic field, B_los = 87 microgauss is incorrect. Because of an insidious software error, what we reported as Stokes V is actually Stokes U: the revised Stokes V spectrum indicates a 3-sigma upper limit of B_los < 17 microgauss. The correct analysis reveals an absorption feature in fractional polarization that is offset in velocity from the Stokes I spectrum by -1.9 km/s. The polarization position-angle spectrum shows a dip that is also significantly offset from the Stokes I feature, but at a velocity that differs slightly from the absorption feature in fractional polarization. We model the absorption feature with 3 velocity components against the core-jet structure of 3C 286. Our chisquare minimization fitting results in components with differing (1) ratios of H I column density to spin temperature, (2) velocity centroids, and (3) velocity dispersions. The change in polarization position angle with frequency implies incomplete coverage of the background jet source by the absorber. It also implies a spatial variation of the polarization position angle across the jet source, which is observed at frequencies higher than the 839.4 MHz absorption frequency. The multi-component structure of the gas is best understood in terms of components with spatial scales of ~100 pc comprised of hundreds of low-temperature (T < 200 K) clouds with linear dimensions of about 1 pc.

[ABRIDGED] We examine the velocity structure in the gas associated with \ion{H}{1} in the damped Ly$\alpha$ absorption system at redshift $z=1.7764$ towards the QSO $1331+170$ using 21cm data, optical and STIS spectra. We find at least two, and possibly three, components showing \ion{C}{1} lines. One of these has Doppler parameter $b=0.55${\kms}, corresponding to a kinetic temperature of 220K if the broadening is thermal. We re-examine the H$_2$ analysis undertaken by \citet{Cui05} using the neutral carbon velocity structure, and find a model which is, consistent with a mixture of collisional and background radiation excitation of the observed H$_2$ rotational levels. For singly ionized heavy elements we find eight components covering a velocity range of $\sim 110$ {\kms}. The \ion{H}{1} structure is expected to follow some combination of the singly ionized and neutral gas, but the 21cm absorption profile is considerably different. This may be because of the different extent and brightness distributions of the radio and optical background sources, and so the spin temperature derived by comparing the Ly$\alpha$ and 21cm line strengths has little physical meaning. The neutral and singly ionized heavy element line profiles also show significant differences, and so the dominant components in each appear to be physically distinct. Attempts to use the range of atomic masses to separate thermal and turbulent components of their Doppler widths were not generally successful. The velocity structure in all ionization stages up to $+3$, apart from the neutral heavy elements, is sufficiently complex that it is difficult to separate out the corresponding velocity components for different ionization levels and determine their column densities.

We discuss two mistreatments of damped Lya (DLA) kinematic analysis that were first performed by Haehnelt, Steinmetz, & Rauch (1998; hereafter HSR98) and have recently been repeated by Hong et al. (2010; arXiv:1008.4242v1, arXiv:1008.4242v2; hereafter H10). Each mistreatment led to the improper excising of simulated absorption profiles. Specifically, their analyses are strictly biased against DLA sightlines that have low HI column density log NHI < 20.5, very high NHI values, and (for all NHI) sightlines with low velocity width Dv (<30 km/s for HSR98; <[20-30] km/s for H10). None of these biases exist in the observational analysis. We suspect these mistreatments compromise the results that followed. Hopefully this posting will prevent their repetition in the future.

We probe the physical conditions in high redshift galaxies, specifically, the Damped Lyman-alpha Systems (DLAs) using neutral carbon (CI) fine structure lines and molecular hydrogen (H2). We report five new detections of CI and analyze the CI in an additional 2 DLAs with previously published data. We also present one new detection of H2 in a DLA. We present a new method of analysis that simultaneously constrains \emph{both} the volume density and the temperature of the gas, as opposed to previous studies that a priori assumed a gas temperature. We use only the column density of CI measured in the fine structure states and the assumption of ionization equilibrium in order to constrain the physical conditions in the gas. We present a sample of 11 CI velocity components in 6 DLAs and compare their properties to those derived by the global CII* technique. The resulting median values for this sample are: <n(HI)> = 69 cm^{-3}, <T> = 50 K, and <log(P/k)> = 3.86 cm^{-3} K, with standard deviations, sigma_{n(HI)} = 134 cm^{-3}, sigma_T = 52 K, and sigma_{log(P/k)} = 3.68 cm^{-3} K. This can be compared with the integrated median values for the same DLAs : <n(HI)> = 2.8 cm^{-3}, <T> = 139 K, and <log(P/k)> = 2.57 cm^{-3} K, with standard deviations sigma_{n(HI)} = 3.0 cm^{-3}, sigma_T = 43 K, and sigma_{log(P/k)} = 0.22 cm^{-3} K. Interestingly, the pressures measured in these high redshift CI clouds are similar to those found in the Milky Way. We conclude that the CI gas is tracing a higher-density, higher-pressure region, possibly indicative of post-shock gas or a photodissociation region on the edge of a molecular cloud. We speculate that these clouds may be direct probes of the precursor sites of star formation in normal galaxies at high redshift.

We report two detections of deuterated molecular hydrogen (HD) in QSO absorption-line systems at $z > 2$. Toward J2123-0500, we find $N$(HD) $= 13.84 \pm 0.2$ for a sub-DLA with metallicity $\simeq 0.5Z_{\odot}$ and $N$(H$_2$) = $17.64 \pm 0.15$ at $z = 2.0594$. Toward FJ0812+32, we find $N$(HD) $= 15.38 \pm 0.3$ for a solar-metallicity DLA with $N$(H$_2$) = $19.88 \pm 0.2$ at $z = 2.6265$. These systems have ratios of HD to H$_2$ above that observed in dense clouds within the Milky Way disk and apparently consistent with a simple conversion from the cosmological ratio of D/H. These ratios are not readily explained by any available model of HD chemistry and there are no obvious trends with metallicity or molecular content. Taken together, these two systems and the two published $z > 2$ HD-bearing DLAs indicate that HD is either less effectively dissociated or more efficiently produced in high-redshift interstellar gas, even at low molecular fraction and/or solar metallicity. It is puzzling that such diverse systems should show such consistent HD/H$_2$ ratios. Without clear knowledge of all the aspects of HD chemistry that may help determine the ratio HD/H$_2$, we conclude that these systems are potentially more revealing of gas chemistry than of D/H itself and that it is premature to use such systems to constrain D/H at high-redshift.

We present high resolution echelle spectra of 7 proximate damped Lyman alpha (PDLA) systems whose relative velocity separation from the background quasar is Delta V < 3000 km/s. Combining our sample with a further 9 PDLAs from the literature we compare the chemical properties of the proximate systems with a control sample of intervening DLAs. Taken at face value, the sample of 16 PDLAs exhibits a wide range of metallicities, ranging from Z ~ 1/3 Z_sun down to Z ~ 1/1000 Z_sun, including the DLA with the lowest N(SiII)/N(HI) yet reported in the literature. We find several pieces of evidence that indicate enhanced ionization and the presence of a hard ionizing spectrum in PDLAs which lead to properties that contrast with the intervening DLAs, particularly when the N(HI) is low. The abundances of Zn, Si and S in PDLAs with log N(HI) > 21, where ionization corrections are minimized, are systematically higher than the intervening population by a factor of around 3. We also find possible evidence for a higher fraction of NV absorbers amongst the PDLAs, although the statistics are still modest. 6/7 of our echelle sample show high ionization species (SiIV, CIV, OVI or NV) offset by >100 km/s from the main low ion absorption. We analyse fine-structure transitions of CII* and SiII* to constrain the PDLA distance from the QSO. Lower limits range from tens of kpc up to >160 kpc for the most stringent limit. We conclude that (at least some) PDLAs do exhibit different characteristics relative to the intervening population out to 3000 km/s (and possibly beyond). Nonetheless, the PDLAs appear distinct from lower column density associated systems and the inferred QSO-absorber separations mean they are unlikely to be associated with the QSO host. We speculate that the PDLAs preferentially sample more massive galaxies in more highly clustered regions of the high redshift universe.

The quasar SDSS J153636.22+044127.0, exhibiting peculiar broad emission-line profiles with multiple components, was proposed as a candidate sub-parsec binary supermassive black hole system. More recently, imaging revealed two spatially distinct sources, leading some to suggest the system to be a quasar pair separated by ~5 kpc. We present Palomar and Keck optical spectra of this system from which we identify a third velocity component to the emission lines. We argue that the system is more likely an unusual member of the class of active galactic nuclei (AGNs) known as "double-peaked emitters" than a sub-parsec black hole binary or quasar pair. We find no significant velocity evolution of the two main peaks over the course of 0.95 yr, with a 3-sigma upper limit on any secular change of 70 km/s/yr. We also find that the three velocity components of the emission lines are spatially coincident to within 0.015" along the slit, apparently ruling out the double-quasar hypothesis.

To date, the identification of interactions at z ~ 3 and above has relied on morphological analysis. Here, we present five serendipitous spectroscopic z ~ 3 Lyman break galaxy (LBG) pairs with projected proper separations < 15 h^-1 kpc in our survey of nine separate Keck fields. The data consist of 140 of our highest signal-to-noise ratio LBG spectra and ~500 of our most confident colour-selected LBGs. We show that the pairs are composed of two distinct close and/or interacting LBGs from a detailed analysis of the rest-frame ultraviolet spectra and images. In addition, we show that the pair number and separation distribution is expected from (1) the angular correlation function when applied to our survey and ~2500 colour-selected LBGs from the literature and (2) an analysis of a carefully matched high-resolution hybrid numerical and analytical cosmological simulation. Because the spectroscopic slitlets have random orientations with respect to the close pairs on the sky, the serendipitous pairs provide an unbiased sampling of the underlying close pair fraction. Finally, we discover two Ly-a emitters (LAEs) in our slitlets and find that they reside within 50 projected h^-1 kpc of LBGs. In this work, we uncover a strong relationship between Ly-a emission and pair separation. All confirmed and all candidate LBG pairs with separations of < 15 projected h^-1 kpc exhibit Ly-a in emission and we find an overabundance of Ly-a emission in pairs with < 50 projected h^-1 kpc separations. This relationship suggests a picture in which a measurable fraction of the Ly-a emission of LBGs, and potentially LAEs, is generated via interaction mechanisms. As a result, serendipitous spectroscopic close pairs provide a unique means to help identify and study high-redshift galaxy interactions using conventional ground-based optical data.

We report on an attempt to accurately wavelength calibrate four nights of data taken with the Keck HIRES spectrograph on QSO PHL957, for the purpose of determining whether the fine structure constant was different in the past. Using new software and techniques, we measured the redshifts of various Ni II, Fe II, Si II, etc. lines in a damped Ly-alpha system at z=2.309. Roughly half the data was taken through the Keck iodine cell which contains thousands of well calibrated iodine lines. Using these iodine exposures to calibrate the normal Th-Ar Keck data pipeline output we found absolute wavelength offsets of 500 m/s to 1000 m/s with drifts of more than 500 m/s over a single night, and drifts of nearly 2000 m/s over several nights. These offsets correspond to an absolute redshift of uncertainty of about Delta z=10^{-5} (Delta lambda= 0.02 Ang), with daily drifts of around Delta z=5x10^{-6} (Delta lambda =0.01 Ang), and multiday drifts of nearly Delta z=2x10^{-5} (0.04 Ang). The causes of the wavelength offsets are not known, but since claimed shifts in the fine structure constant would result in velocity shifts of less than 100 m/s, this level of systematic uncertainty makes may make it difficult to use Keck HIRES data to constrain the change in the fine structure constant. Using our calibrated data, we applied both our own fitting software and standard fitting software to measure (Delta alpha)/alpha, but discovered that we could obtain results ranging from significant detection of either sign, to strong null limits, depending upon which sets of lines and which fitting method was used. We thus speculate that the discrepant results on (Delta alpha)/alpha reported in the literature may be due to random fluctuations coming from under-estimated systematic errors in wavelength calibration and fitting procedure.

We present a sample of 407 z~3 Lyman break galaxies (LBGs) to a limiting isophotal u-band magnitude of 27.6 mag in the Hubble Ultra Deep Field (UDF). The LBGs are selected using a combination of photometric redshifts and the u-band drop-out technique enabled by the introduction of an extremely deep u-band image obtained with the Keck I telescope and the blue channel of the LRIS spectrometer. The Keck u-band image, totaling 9 hours of integration time, has a one sigma depth of 30.7 mag arcsec^-2, making it one of the most sensitive u-band images ever obtained. The u-band image also substantially improves the accuracy of photometric redshift measurements of ~50% of the z~3 Lyman break galaxies, significantly reducing the traditional degeneracy of colors between z~3 and z~0.2 galaxies. This sample provides the most sensitive, high-resolution multi-filter imaging of reliably identified z~3 LBGs for morphological studies of galaxy formation and evolution and the star formation efficiency of gas at high redshift.

We present the first direct evidence for cold gas in a high redshift DLA galaxy. We measured several multiplets of weak neutral carbon (CI) transitions in order to perform a curve of growth analysis. A delta chi-squared test constrains the best fit Doppler parameter, b = 0.33_{-0.04}^{+0.05} km/s, and logN(CI) = 13.30 +- 0.2 cm^-2. This Doppler parameter constrains the kinetic temperature of the gas to T <= 78 K (T <= 115 K, 2 sigma). We used the associated CI fine structure lines to constrain the volume density of the gas, n(HI) ~ 40 - 200 cm^-3 (2 sigma), resulting in a lower limit on the cloud size of approximately 0.1 - 1 parsec. While it is difficult to determine the metallicity of the cold component, the absence of Cr II indicates that the cold cloud suffers a high level of dust depletion. Additionally, the large amount of Lyman and Werner-band molecular hydrogen absorption (log N(H2)_{total} = 19.88 cm^-2, f_{H_2} >= 0.06) with an excitation temperature of T_{ex} = 46 K as determined by the rotational J = 0 and J = 1 states, is consistent with the presence of cold gas. We propose that this cloud may be gravitationally confined and may represent a transition gas-phase from primarily neutral atomic gas, to a colder, denser molecular phase that will eventually host star formation.

We continue our recent work to characterize the plasma content of high-redshift damped and sub-damped Lyman-alpha systems (DLAs/sub-DLAs), which represent multi-phase gaseous (proto)galactic disks and halos seen toward a background source. We survey N V absorption in a sample of 91 DLAs and 18 sub-DLAs in the redshift range 1.67<z<4.28 with unblended coverage of the N V doublet, using data from VLT/UVES, Keck/HIRES, and Keck/ESI. In DLAs, we find eight secure N V detections, four marginal detections, and 79 non-detections. The detection rate of N V in DLAs is therefore 13^{+5}_{-4}%. Two sub-DLA N V detections are found among a sample of 18, at a similar detection rate of 11^{+15}_{-7}%. We show that the N V detection rate is a strong function of metallicity, increasing by a factor of ~4 at [N/H]=[NI/HI]>-2.3. The N V and CIV component b-value distributions in DLAs are statistically similar, but the median b(N V) of 18 km/s is lower than the median b(O VI) of 25 km/s. Some ~20% of the N V components have b<10 km/s and thus arise in warm photoionized plasma at log (T/K)<4.92; local sources of ionizing radiation (as opposed to the extragalactic background) are required to keep the cloud sizes physically reasonable. The nature of the remaining ~80% of (broad) N V components is unclear; models of radiatively-cooling collisionally-ionized plasma at log(T/K)=5.2-5.4 are fairly successful in reproducing the observed integrated high-ion column density ratios and the component line widths, but we cannot rule out photoionization by local sources. Finally, we identify several unusual DLAs with extremely low metallicity (<0.01 solar) but strong high-ion absorption [log N(N V)>14 or log N(O VI)>14.2] that present challenges to either galactic inflow or outflow models.

Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.

When, and how, did the first galaxies and supermassive black holes (SMBH) form, and how did they reionization the Universe? First galaxy formation and cosmic reionization are among the last frontiers in studies of cosmic structure formation. We delineate the detailed astrophysical probes of early galaxy and SMBH formation afforded by observations at centimeter through submillimeter wavelengths. These observations include studies of the molecular gas (= the fuel for star formation in galaxies), atomic fine structure lines (= the dominant ISM gas coolant), thermal dust continuum emission (= an ideal star formation rate estimator), and radio continuum emission from star formation and relativistic jets. High resolution spectroscopic imaging can be used to study galaxy dynamics and star formation on sub-kpc scales. These cm and mm observations are the necessary compliment to near-IR observations, which probe the stars and ionized gas, and X-ray observations, which reveal the AGN. Together, a suite of revolutionary observatories planned for the next decade from centimeter to X-ray wavelengths will provide the requisite panchromatic view of the complex processes involved in the formation of the first generation of galaxies and SMBHs, and cosmic reionization.

The magnetic field pervading our Galaxy is a crucial constituent of the interstellar medium: it mediates the dynamics of interstellar clouds, the energy density of cosmic rays, and the formation of stars. The field associated with ionized interstellar gas has been determined through observations of pulsars in our Galaxy. Radio-frequency measurements of pulse dispersion and the rotation of the plane of linear polarization, i.e., Faraday rotation, yield an average value B ~ 3 microGauss. The possible detection of Faraday rotation of linearly polarized photons emitted by high-redshift quasars suggests similar magnetic fields are present in foreground galaxies with redshifts z > 1. As Faraday rotation alone, however, determines neither the magnitude nor the redshift of the magnetic field, the strength of galactic magnetic fields at redshifts z > 0 remains uncertain. Here we report a measurement of a magnetic field of B ~ 84 microGauss in a galaxy at z =0.692, using the same Zeeman-splitting technique that revealed an average value of B = 6 microGauss in the neutral interstellar gas of our Galaxy. This is unexpected, as the leading theory of magnetic field generation, the mean-field dynamo model, predicts large-scale magnetic fields to be weaker in the past rather than stronger.

We present new results on the frequency distribution of projected HI column densities f(N,X), total comoving covering fraction, and integrated mass densities rho_HI of high redshift, HI `disks' from a survey of damped Lya systems (DLAs) in the Sloan Digital Sky Survey, Data Release 5. For the full sample spanning z=2.2 to 5 [738 DLAs], f(N,X) is well fitted by a double power-law with a break column density N_d = 10^(21.55 +/- 0.04) and low/high-end exponents alpha = -2.00 +/- 0.05, -6.4^{+1.1}_{-1.6}. The shape of f(N,X) is invariant during this redshift interval and also follows the projected surface density distribution of present-day HI disks as inferred from 21cm observations. We conclude that HI gas has been distributed in a self-similar fashion for the past 12Gyr. The normalization of f(N,X), in contrast, decreases by a factor of two during the ~2Gyr interval from z=4 to 2.2 giving corresponding decreases in both the total covering fraction and rho_HI. At z~2, these quantities match the present-day values suggesting no evolution during the past ~10Gyr. We argue that the evolution at early times is driven by `violent' processes that removes gas from nearly half the galaxies at z~3 establishing the antecedants of current early-type galaxies. The perceived constancy of rho_HI, meanwhile, implies that HI gas is a necessary but insufficient pre-condition for star formation and that the global star-formation rate is driven by the accretion and condensation of fresh gas from the intergalactic medium.

We investigate the physical processes occuring in the multiphase gas of a damped Ly\alpha system (DLA). We base our analysis on a high quality Keck HIRES spectrum of the QSO J1211+0422 in which a DLA is detected at z=2.377. There is little contamination of the high-ion (OVI, NV, CIV, SiIV) absorption, allowing us to explore the properties of the highly ionized gas and its connection to other gas-phases. The metallicity ([Z/H]=-1.41+/-0.08), HI column density (log N(HI)=20.80+/-0.10), full-width velocity (\Delta(v_ neut)=70 km/s) and relative abundances ([Si/Fe]=+0.23+/-0.05 and [N/Si]=-0.88+/-0.07) of this DLA are not unusual. However, we derive the lowest CII* cooling rate in a DLA, l_c < 10^{-27.8} erg/s per H atom (3\sigma). Using this stringent limit, we show that the neutral gas (confined at |v|<+39 km/s) must be warm and the star formation rate is <7.1x10^{-3} M_odot/yr/kpc^2. Surprisingly, the gas shows strong, complex absorption profiles from highly ionized gas whose kinematics appear connected to each other and the low ions. The total amount of highly and weakly ionized gas is very large with N(HII)/N(HI)>1.5. At |v|>+39 km/s, the gas is fully and highly ionized (H+/H~1, N(CIV)>>N(CII), N(SiIV)>>N(SiII)). Based on ionization models, OVI and NV are generally difficult to produce by hard photons, while SiIV and CIV can be photoionized to a large extent. There is, however, no evidence of OVI-bearing gas at T~10^6 K associated with this DLA. In contrast, there is some evidence for narrow OVI, NV, and CIV components (unexplained by photoionization), implying too low temperatures (T < 10^5 K) for simple collisional ionization models to produce their observed column densities. Stellar feedback is a possible source for producing the high ions, but we cannot rule out accretion of non-pristine material onto the protogalaxy.

The median observed velocity width v_90 of low-ionization species in damped Ly-alpha systems is close to 90 km/s, with approximately 10% of all systems showing v_90 > 210 km/s at z=3. We show that a relative shortage of such high-velocity neutral gas absorbers in state-of-the-art galaxy formation models is a fundamental problem, present both in grid-based and particle-based numerical simulations. Using a series of numerical simulations of varying resolution and box size to cover a wide range of halo masses, we demonstrate that energy from gravitational infall alone is insufficient to produce the velocity dispersion observed in damped Ly-alpha systems, nor does this dispersion arise from an implementation of star formation and feedback in our highest resolution (~ 45 pc) models, if we do not put any galactic winds into our models by hand. We argue that these numerical experiments highlight the need to separate dynamics of different components of the multiphase interstellar medium at z=3.

We report the discovery of a very bright (m_R = 22.2) Lyman break galaxy at z = 3.03 that appears to be a massive system in a late stage of merging. Deep imaging reveals multiple peaks in the brightness profile with angular separations of ~0.''8 (~25 h^-1 kpc comoving). In addition, high signal-to-noise ratio rest-frame UV spectroscopy shows evidence for ~5 components based on stellar photospheric and ISM absorption lines with a velocity dispersion of sigma ~460 km s^-1 for the three strongest components. Both the dynamics and high luminosity, as well as our analysis of a LCDM numerical simulation, suggest a very massive system with halo mass M ~ 10^13 M_solar. The simulation finds that all halos at z = 3 of this mass contain sub-halos in agreement with the properties of these observed components and that such systems typically evolve into M ~ 10^14 M_solar halos in groups and clusters by z = 0. This discovery provides a rare opportunity to study the properties and individual components of z ~ 3 systems that are likely to be the progenitors to brightest cluster galaxies.

We report evidence for a bimodality in damped Ly systems (DLAs). Using [C II] 158 mu cooling rates, lc, we find a distribution with peaks at lc=10^-27.4 and 10^-26.6 ergs s^-1 H^-1 separated by a trough at lc^crit ~= lc < 10^-27.0 ergs s^-1 H^-1. We divide the sample into low cool DLAs with lc < lc^crit and high cool DLAs with lc > lc^crit and find the Kolmogorv-Smirnov probabilities that velocity width, metallicity, dust-to-gas ratio, and Si II equivalent width in the two subsamples are drawn from the same parent population are small. All these quantities are significantly larger in the high cool population, while the H I column densities are indistinguishable in the two populations. We find that heating by X-ray and FUV background radiation is insufficient to balance the cooling rates of either population. Rather, the DLA gas is heated by local radiation fields. The rare appearance of faint, extended objects in the Hubble Ultra Deep Field rules out in situ star formation as the dominant star-formation mode for the high cool population, but is compatible with in situ star formation as the dominant mode for the low cool population. Star formation in the high cool DLAs likely arises in Lyman Break galaxies. We investigate whether these properties of DLAs are analogous to the bimodal properties of nearby galaxies. Using Si II equivalent width as a mass indicator, we construct bivariate distributions of metallicity, lc, and areal SFR versus the mass indicators. Tentative evidence is found for correlations and parallel sequences, which suggest similarities between DLAs and nearby galaxies. We suggest that the transition-mass model provides a plausible scenario for the bimodality we have found. As a result, the bimodality in current galaxies may have originated in DLAs.

We analyzed the spectroscopic and photometric database of the 5th data release of the Sloan Digital Sky Survey (SDSS) to search for evidence of the quasar reddening produced by dust embedded in intervening damped Ly alpha (DLA) systems. From a list of 5164 quasars in the interval of emission redshift 2.25 </= z_e </= 3.5 and SDSS spectra with signal-to-noise ratio SNR >/= 4, we built up an "absorption sample" of 248 QSOs with a single DLA system in the interval of absorption redshift 2.2 < z_a </= 3.5 and a "pool" of 1959 control QSOs without DLA systems or strong metal systems. For each QSO of the absorption sample we extracted from the pool a subset of control QSOs that are closest in redshift and magnitude. The mean color of this subset was used as a zero point to measure the "deviation from the mean color" of individual DLA-QSOs, Delta_i. The colors were measured using "BEST" ugriz SDSS imaging data. The mean color excess of the absorption sample, <E>, was estimated by averaging the individual color deviations Delta_i. We find <E(r-z)> = 27 +/- 9 x 10**(-3) mag and <E(g-z)> = 54 +/- 12 x 10**(-3) mag. These values are representative of the reddening of DLA systems at z_a ~ 2.7 in SDSS QSOs with limiting magnitude r =/~ 20.2. The detection of the mean reddening is confirmed by several statistical tests. Analysis of the results suggests an origin of the reddening in dust embedded in the DLA systems, with an SMC-type extinction curve. By converting the reddening into rest-frame extinction, we derive a mean dust-to-gas ratio <A_V/N(HI)> ~ 2 to 4 x 10**(-23) mag cm^2. This value is ~ -1.25 dex lower than the mean dust-to-gas ratio of the Milky Way, in line with the lower level of metallicity in the present DLA sample.

We analyze the gas kinematics of damped Lya systems (DLAs) hosting high z gamma-ray bursts (GRBs) and those toward quasars (QSO-DLAs) focusing on threestatistics: (1) dv, the velocity interval encompassing 90% of the totaloptical depth, (2,3) Wsi and Wciv, the rest equivalent widths of the SiII1526 and CIV1548 transitions. The dv distributions of the GRB-DLAs and QSO-DLAs are similar, each has median dv~80km/s and a significant tail to several hundred km/s. This suggests comparable galaxy masses for the parent populations of GRB-DLAs and QSO-DLAs and we infer the average dark matter halo mass of GRB galaxies is <~10^{12} Msol. The unique configuration of GRB-DLA sightlines and the presence (and absence) of fine-structure absorption together give special insight into the nature of high z, protogalactic velocity fields. The data support a scenario where the dv statistic reflects dynamics in the interstellar medium (ISM) and Wsi traces motions outside the ISM (e.g. halo gas, galactic-scale winds). The Wsi statistic and gas metallicity [M/H] are tightly correlated, especially for the QSO-DLAs: [M/H]=a + b log(Wsi/1A) with a=-0.92+/-0.05 and b=-1.41+/-0.10. We argue that the Wsi statistic primarily tracks dynamical motions in the halos of high z galaxies and interpret this correlation as a mass-metallicity relation with very similar slope to the trend observed in local, low-metallicity galaxies. Finally, the GRB-DLAs exhibit systematically larger Wsi values (>0.5A) than the QSO-DLAs (<Wsi>~ 0.5A) which may suggest galactic-scale outflows contribute to the largest observed velocity fields.

We publish the Keck/HIRES and Keck/ESI spectra that we have obtained during the first 10 years of Keck observatory operations. Our full sample includes 42 HIRES spectra and 39 ESI spectra along 65 unique sightlines providing abundance measurements on ~85 damped Lya systems. The normalized data can be downloaded from the journal or from our supporting website: http://www.ucolick.org/~xavier/DLA/. The database includes all of the sightlines that have been included in our papers on the chemical abundances, kinematics, and metallicities of the damped Lya systems. This data has also been used to argue for variations in the fine-structure constant. We present new chemical abundance measurements for 10 damped Lya systems and a summary table of high-resolution metallicity measurements (including values from the literature) for 153 damped Lya systems at z>1.6. We caution, however, that this metallicity sample (and all previous ones) is biased to higher N(HI) values than a random sample.

We present the results of a numerical study of a galactic wind model and its implications on the properties of damped Lyman-alpha absorbers (DLAs) using cosmological hydrodynamic simulations. We vary both the wind strength and the internal parameters of the the wind model in a series of cosmological SPH simulations that include radiative cooling and heating by a UV background, star formation, and feedback from supernovae and galactic winds. To test our simulations, we examine the DLA `rate-of-incidence' as a function of halo mass, galaxy apparent magnitude, and impact parameter. We find that the statistical distribution of DLAs does not depend on the exact values of internal numerical parameters that control the decoupling of hydrodynamic forces when the gas is ejected from starforming regions. The DLA rate-of-incidence in our simulations at z=3 is dominated by the faint galaxies with apparent magnitude R_AB < 25.5. However, interestingly in a `strong wind' run, the differential distribution of DLA sight-lines is peaked at Mhalo = 10^{12} Msun/h (R_AB~27), and the mean DLA halo mass is Mmean=10^{12.4} Msun/h (R_AB ~ 26). These mass-scales are much larger than those if we ignore winds, because galactic wind feedback suppresses the DLA cross section in low-mass halos and increases the relative contribution to the DLA incidence from more massive halos. The DLAs in our simulations are more compact than the present-day disk galaxies, and the impact parameter distribution is very narrow unless we limit the search for the host galaxy to only bright LBGs. The comoving number density of DLAs is higher than that of LBGs down to R_AB=30 mag if the physical radius of each DLA is smaller than 5 kpc/h_70. We discuss conflicts between current simulations and observations, and potential problems with simulations based on the CDM model.

(Abridged) The Kennicutt--Schmidt law relates the face-on star formation rate (SFR) per unit area with the face-on gaseous column density in nearby galaxies. Applying this relation to damped Lya absorption systems (DLAs) of neutral hydrogen column density N > 1.6x10^21 cm^-2 leads to an estimate that three percent of the sky should be covered with extended sources brighter than mu_V~28.4 mag arcsec^-2, if DLAs at redshift z=[2.5,3.5] undergo in situ star formation. We test this hypothesis by searching the Hubble Ultra Deep Field (UDF) F606W image for low surface-brightness features of angular sizes, ranging between theta_dla=0.25" and 4.0". Our search yields upper limits on the comoving SFR densities that are between factors of 30 and 100 lower than predictions, suggesting a reduction by more than a factor of 10 in star formation efficiency at z~3. We consider several mechanisms that could reduce star formation efficiency at high redshift. We find that the cosmological increase with redshift of the critical surface density for the Toomre instability may be sufficient to suppress star formation to the levels implied by the UDF observations. However, the uncertainties are such that Toomre instabilities may still exist. In that case star formation at column densities less than 10^22 cm^-2 may be suppressed by the low molecular content of the DLA gas. The upper limits on in situ star formation reduce the predicted metallicities at z~3 to be significantly lower than observed, and reduce the heat input in the gas to be substantially lower than the inferred cooling rates. In contrast, the radiative output from compact Lyman Break Galaxies (LBGs) with R<27 is sufficient to balance the comoving cooling rate.

We have identified a metal-strong (logN(Zn+) > 13.15 or logN(Si+) > 15.95) DLA (MSDLA) population from an automated quasar (QSO) absorber search in the Sloan Digital Sky Survey Data Release 3 (SDSS-DR3) quasar sample, and find that MSDLAs comprise ~5% of the entire DLA population with z_abs > 2.2 found in QSO sightlines with r < 19.5. We have also acquired 27 Keck ESI follow-up spectra of metal-strong candidates to evaluate our automated technique and examine the MSDLA candidates at higher resolution. We demonstrate that the rest equivalent widths of strong ZnII 2026 and SiII 1808 lines in low-resolution SDSS spectra are accurate metal-strong indicators for higher-resolution spectra, and predict the observed equivalent widths and signal-to-noise ratios needed to detect certain extremely weak lines with high-resolution instruments. We investigate how the MSDLAs may affect previous studies concerning a dust-obscuration bias and the N(HI)-weighted cosmic mean metallicity <Z(z)>. Finally, we include a brief discussion of abundance ratios in our ESI sample and find that underlying mostly Type II supernovae enrichment are differential depletion effects due to dust (and in a few cases quite strong); we present here a handful of new Ti and Mn measurements, both of which are useful probes of depletion in DLAs. Future papers will present detailed examinations of particularly metal-strong DLAs from high-resolution KeckI/HIRES and VLT/UVES spectra.

We use 211 galaxy spectra from our survey for Lyman break galaxies (LBGs) associated with 11 damped Lyman alpha systems (DLAs) to measure the three-dimensional LBG auto-correlation and DLA-LBG cross-correlation functions with the primary goal of inferring the mass of DLAs at z~3. Conventional binning of the data while varying both r_0 and gamma parameters of the fiducial model of the correlation function xi(r)=(r/r_0)^{-gamma} resulted in the best fit values and 1 sigma uncertainties of r_0=2.65+/-0.48, gamma=1.55+/-0.40 for the LBG auto-correlation and r_0=3.32+/-1.25, gamma=1.74+/-0.36 for DLA-LBG cross-correlation function. To circumvent shortcomings found in binning small datasets, we perform a maximum likelihood analysis based on Poisson statistics. The best fit values and 1 sigma confidence levels from this analysis were found to be r_0=2.91(+1.0,-1.0), gamma=1.21(+0.6,-0.3) for the LBG auto-correlation and r_0=2.81(+1.4,-2.0), gamma=2.11(+1.3,-1.4) for the DLA-LBG cross-correlation function. We report a redshift spike of five LBGs with Delta z = 0.015 of the z=2.936 DLA in the PSS0808+5215 field and find that the DLA-LBG clustering signal survives when omitting this field from the analysis. Using the correlation functions measurements and uncertainties, we compute the z~3 LBG galaxy bias b_LBG to be 1.5<b_LBG<3 corresponding to an average halo mass of 10^(9.7)<M_LBG<10^(11.6) M_odot and the z~3 DLA galaxy bias b_DLA to be 1.3<b_DLA<4 corresponding to an average halo mass of 10^(9)<M_LBG<10^(12) M_odot. Lastly, two of the six QSOs discovered in this survey were found to lie within Delta z = 0.0125 of two of the survey DLAs. We estimate the probability of this occurring by chance is 1 in 940 and may indicate a possible relationship between the distribution of QSOs and DLAs at z~3.

We discuss the detectability of high-redshift galaxies via [CII] 158 micron line emission by coupling an analytic model with cosmological Smoothed Particle Hydrodynamics (SPH) simulations that are based on the concordance Lambda cold dark matter (CDM) model. Our analytic model describes a multiphase interstellar medium irradiated by the far ultra-violet radiation from local star-forming regions, and it calculates thermal and ionization equilibrium between cooling and heating. The model allows us to predict the mass fraction of a cold neutral medium (CNM) embedded in a warm neutral medium (WNM). Our cosmological SPH simulations include a treatment of radiative cooling/heating, star formation, and feedback effects from supernovae and galactic winds. Using our method, we make predictions for the [CII] luminosity from high-redshift galaxies which can be directly compared with upcoming observations by the Atacama Large Millimeter Array (ALMA) and the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We find that the number density of high-redshift galaxies detectable by ALMA and SPICA via [CII] emission depends significantly on the amount of neutral gas which is highly uncertain. Our calculations suggest that, in a CDM universe, most [CII] sources at z=3 are faint objects with \Snu < 0.01 mJy. Lyman-break galaxies (LBGs) brighter than R_AB=23.5 mag are expected to have flux densities \Snu = 1-3 mJy depending on the strength of galactic wind feedback. The recommended observing strategy for ALMA and SPICA is to aim at very bright LBGs or star-forming DRG/BzK galaxies.

As large optical quasar surveys for damped Lya become a reality and the study of star forming gas in the early Universe achieves statistical robustness, it is now vital to identify and quantify the sources of systematic error. Because the nature of optically-selected quasar surveys makes them vulnerable to dust obscuration, we have undertaken a radio-selected quasar survey for damped Lya systems to address this bias. We present the definition and results of this survey. We then combine our sample with the CORALS dataset to investigate the HI column density distribution function f(N) of damped Lya systems toward radio-selected quasars. We find that f(N) is well fit by a power-law f(N) = k_1 N^alpha_1, with log k_1 = 22.90 and alpha_1 = -2.18. This power-law is in excellent agreement with that of optically-selected samples at low N(HI), an important yet expected result given that obscuration should have negligible effect at these gas columns. However, because of the relatively small size of the radio-selected sample, 26 damped Lya systems in 119 quasars, f(N) is not well constrained at large N(HI) and the first moment of the HI distribution function, Omega_g, is, strictly speaking, a lower limit. The power-law is steep enough, however, that extrapolating it to higher column densities implies only a modest, logarithmic increase in Omega_g. The radio-selected value of Omega_g = 1.15 x 10^-3, agrees well with the results of optically-selected surveys. While our results indicate that dust obscuration is likely not a major issue for surveys of damped Lya systems, we estimate that a radio-selected sample of approximately 100 damped Lya systems will be required to obtain the precision necessary to absolutely confirm an absence of dust bias.

We use high resolution Eulerian hydrodynamics simulations to study kinematic properties of the low ionization species in damped Ly-alpha systems at redshift z=3. Our adaptive mesh refinement simulations include most key ingredients relevant for modeling neutral gas in high-column density absorbers: hydrodynamics, gravitational collapse, continuum radiative transfer and gas chemistry, but no star formation. We model high-resolution Keck spectra with unsaturated low ion transitions in two Si II lines (1526 and 1808 A), and compare simulated line profiles to the data from the SDSS DLA survey. We find that with increasing grid resolution the models show a trend in convergence towards the observed distribution of HI column densities. While in our highest resolution model we recover the cumulative number of DLAs per unit absorption distance, none of our models predicts DLA velocity widths as high as indicated by the data, suggesting that feedback from star formation might be important. At z=3 a non-negligible fraction of DLAs with column densities below 10^21 cm^-2 is caused by tidal tails due to galaxy-galaxy interactions in more massive halo environments. Lower column density absorbers with N_HI < 10^21.4 cm^-2 are sensitive to changes in the UV background resulting in a 10% reduction of the cumulative number of DLAs for twice the quasar background relative to the fiducial value. We find that the mass cut-off below which a large fraction of dwarf galaxies cannot retain gas after reionization is 7*10^7 msun, lower than the previous estimates. Finally, we show that models with self-shielding commonly used in the literature produce significantly lower DLA velocity widths than the full radiative transfer runs.

We present the first spectroscopic measurement of the spatial cross-correlation function between damped Lyman alpha systems (DLAs) and Lyman break galaxies (LBGs). We obtained deep u'BVRI images of nine QSO fields with 11 known z ~ 3 DLAs and spectroscopically confirmed 211 R < 25.5 photometrically selected z > 2 LBGs. We find strong evidence for an overdensity of LBGs near DLAs versus random, the results of which are similar to that of LBGs near other LBGs. A maximum likelihood cross-correlation analysis found the best fit correlation length value of r_0 = 2.9^(+1.4)_(-1.5) h^(-1)Mpc using a fixed value of gamma = 1.6. The implications of the DLA-LBG clustering amplitude on the average dark matter halo mass of DLAs are discussed.

Observations of damped Lyman alpha systems offer a unique window on the neutral-gas reservoirs that gave rise to galaxies at high redshifts. This review focuses on critical properties such as the H I and metal content of the gas and on independent evidence for star formation. Together, these provide an emerging picture of gravitationally bound objects in which accretion of gas from the IGM replenishes gas consumed by star formation. Other properties such as dust content, molecular content, ionized-gas content, gas kinematics, and galaxy identifications are also reviewed. These properties point to a multiphase ISM in which radiative and hydrodynamic feedback processes are present. Numerical simulations and other types of models used to describe damped Lyman alpha systems within the context of galaxy formation are also discussed.

We present the results from a damped Lya survey of the Sloan Digital Sky Survey, Data Release 3 based on over 500 new damped Lya systems at z>2.2. We measure the HI column density distribution f(N) and its zeroth and first moments (the incidence l(X) and gas mass-density O_dla of damped Lya systems, respectively) as a function of redshift. The key results include: (1) the f(N) distribution is well fit by a Gamma-function with `break' column density log N_g=10^21.5 and `faint-end' slope alpha=-1.8; (2) the shape of the f(N) distributions do not show evolution with redshift; (3) l(X) and O_dla decrease by 35% and 50% during ~1Gyr between redshift z=[3.,3.5] to z=[2.2,2.5]; and (4) l(X) and O_dla in the lowest SDSS redshift bin (z=2.2) are consistent with the current values. We investigate systematic errors in damped Lya analysis and identify only one important effect: we measure 40 +/- 20% higher O_dla values toward a subset of brighter quasars than toward a faint subset. This effect runs contrary to the bias associated with dust obscuration and suggests that gravitational lensing may be important. Comparing the results against models of galaxy formation, we find all of the models significantly underpredict l(X) at z=3 and only SPH models with significant feedback may reproduce O_dla at high redshift. We argue that the Lyman limit systems contribute ~1/3 of the universe's HI atoms at all redshifts z=2 to 5 and that the f(N) distribution for N(HI)<10^20 has an inflection with slope >-1. We advocate a new mass density definition -- the mass density of predominantly neutral gas O_neut -- to be contrasted with the mass density of gas associated with HI atoms. We contend the damped Lya systems contribute >80% of O_neut at all redshifts and therefore are the main reservoirs for star formation. [abridged]

We discuss interstellar temperature determinations using the excitation equilibrium of the ^2P levels of Si II and C II. We show how observations of the ^2P_3/2 fine structure levels of Si II and C II (which have significantly different excitation energies, corresponding to ~413 and 92 K, respectively) can be used to limit gas kinetic temperatures. We apply this method to the z=4.224 damped Lyman-alpha system toward the quasar PSS1443+27. The lack of significant absorption out of the SiII ^2P_3/2 level and the presence of very strong C II ^2P_3/2 provides an upper limit to the temperature of the C II*-bearing gas in this system. Assuming a solar Si/C ratio, the observations imply a 2-sigma limit T<954 K for this absorber; a super-solar Si/C ratio gives stricter limits, T<524 K. The observations suggest the presence of a cold neutral medium; such cold gas may serve as the fuel for star formation in this young galaxy.

We present a survey for z~3 Lyman break galaxies (LBGs) associated with damped Lyman alpha systems (DLAs) with the primary purpose of determining the DLA-LBG cross-correlation. This paper describes the acquisition and analysis of imaging and spectroscopic data of 9 quasar fields having 11 known z~3 DLAs covering an area of 465 arcmin^2. Using deep u'BVRI images, 796 LBG candidates to an apparent R_AB magnitude of 25.5 were photometrically selected from 17,343 sources detected in the field. Spectroscopic observations of 529 LBG candidates using Keck LRIS yielded 339 redshifts. We have conservatively identified 211 z>2 objects with <z>=3.02+/-0.32. We discuss our method of z~3 LBG identification and present a model of the u'BVRI photometric selection function. We use the 339 spectra to evaluate our u'BVRI z~3 Lyman break photometric selection technique.

We present Keck/HIRES abundance measurements and metal-line kinematic profiles of the damped Lyman alpha systems (DLAs) towards the quasars Q0933+733 (z_abs=1.479) and Q0948+433 (z_abs=1.233). These two DLAs have among the five highest HI column densities at any redshift: N(HI)=4.2E21 cm^{-2}. The metal-line data, presented here for the first time, reveal that these DLAs are noteworthy for several other reasons as well. 1) The Q0933+733 DLA exhibits simple kinematic structure with unusually narrow velocity widths as measured from its unsaturated metal lines (delta v=16 km/s). At 2.6% solar, it has the second lowest metallicity at z<2. 2) The Q0948+433 DLA has among the strongest metal-line transitions of any known DLA. The saturated SiII1808 line implies a high metallicity ([Si/H]> -1) and a significant alpha-enhancement. The strong metal lines of this DLA have made possible the detection of TiII1910, CoII2012, and MgI2026. 3) We find that the relative gas-phase abundances of both DLAs follow the general trend seen at high redshift, e.g., enhanced Zn/Fe and Si/Fe, and sub-solar Mn/Fe, indicating that there is little evolution in the nucleosynthetic patterns of DLAs down to this epoch. 4) Their high HI column densities imply that these DLAs dominate the column density-weighted cosmic mean metallicity, <Z>, of the universe at z<1.5. Using the 15 DLAs with measured metallicities in the redshift interval 0.4<z<1.5, we find <Z> = -0.89^{+0.40}_{-0.33}, where the uncertainties are 95% confidence limits.

We investigate the heat source of the neutral gas comprising DLAs. Unlike the Lya forest, where the extragalactic background radiation field ionizes and heats the gas, we find that grain photoelectric heating by the FUV background is not sufficient to balance the C II 158um cooling rate inferred from DLAs. In these systems, a local energy source is required. We show that in the case of the z=1.919 DLA toward Q2206-19, the local source is FUV emission from the associated galaxy found by Moller et al (2002): the mean intensity inferred from photometry is in good agreement with the intensity Jnu required to explain the cooling rate. The FUV mean intensity predicted for a cold neutral medium (CNM) model, Jnu=(1.7+2.7-1.0)x10^(-18) cgs (95% c.l.), is the largest expected from our CII* study of 45 DLAs. This may explain why this is the only confirmed DLA yet detected in emission at z>1.9. We argue that in most DLAs with detected CII* absorption, Jnu is between 10^{-19} and 10^{-18} and heats the gas which is a CNM. By contrast, in most DLAs with upper limits on CII* absorption the gas is a warm neutral medium (WNM). Surprisingly, the upper limits are compatible with the same range of Jnu values suggesting the majority of DLAs are heated by radiation fields generated by a limited range of star formation rates per unit H I area, between 10^{-3} and 10^{-2} Msol/kpc^2. We also show that CII* absorption is unlikely to arise in gas that is ionized.

We present evidence for variations in the fine-structure constant from Keck/HIRES spectra of 143 quasar absorption systems over the redshift range 0.2 < z_abs < 4.2. This includes 15 new systems, mostly at high-z (z_abs > 1.8). Our most robust estimate is a weighted mean da/a=(-0.57 +/- 0.11)x10^{-5}. We respond to recent criticisms of the many-multiplet method used to extract these constraints. The most important potential systematic error at low-z is the possibility of very different Mg heavy isotope abundances in the absorption clouds and laboratory: {\it higher} abundances of {25,26}Mg in the absorbers may explain the low-z results. Approximately equal mixes of {24}Mg and {25,26}Mg are required. Observations of Galactic stars generally show {\it lower} {25,26}Mg isotope fractions at the low metallicities typifying the absorbers. Higher values can be achieved with an enhanced population of intermediate mass stars at high redshift, a possibility at odds with observed absorption system element abundances. At present, all observational evidence is consistent with the varying-alpha results. Another promising method to search for variation of fundamental constants involves comparing different atomic clocks. Here we calculate the dependence of nuclear magnetic moments on quark masses and obtain limits on the variation of alpha and m_q/Lambda_QCD from recent atomic clock experiments with hyperfine transitions in H, Rb, Cs, Hg+ and an optical transition in Hg+.

Utilizing a new, high mass resolution hydrodynamic simulation we compute the metallicity evolution of damped Lyman alpha systems (DLAs) and find a reasonable agreement with observations. In particular, the observed slow evolution of the DLA metallicity occurs naturally in the simulation due to the combined effects of physical and observational selection. The slow metallicity evolution is caused by the steady transformation, with increasing time, of the highest metallicity systems to "galaxies", thus depleting this category, while all the lower metallicity systems show, individually, an increase in metallicity. Although the trend of DLA metallicity with redshift is in good agreement with observations, it appears that the average metallicity of simulated DLAs is higher than observed by 0.3-0.5 dex in the probed redshift range (z=0-5). Our study indicates that this difference may be attributed to observational selection effects due to dust obscuration. If we allow for a dust obscuration effect, our model reproduces the observed metallicity evolution in both amplitude and slope. We find that DLAs are not a simple population but probe a range of different systems and the mix changes with redshift. About 50% of all metals in the gaseous phase is in DLAs at all times from z=5 to z=1, making a rapid downturn at z\le 1 to ~20% by z=0, as metals are swept into the hotter components of the IGM as well as locked up in stars. While not the primary focus of this study, we find that the model provides good matches to observations with respect to column density distribution and evolution of neutral gas content, if the same dust obscuration is taken into account. We find Omega_{DLA,comp}=(1-3)E-3, depending on the effect of dust obscuration.

We present accurate metallicity measurements for 121 damped Lya systems at 0.5<z<5 including ~50 new measurements from our recently published Echellette Spectrograph and Imager surveys. This dataset is analysed to determine the age-metallicity relation of neutral gas in the universe. Contrary to previous datasets this sample shows statistically significant evolution in the mean metallicity. The best linear fit rate to metallicity vs. redshift is -0.26 +/- 0.07 dex corresponding to approximately a factor of 2 every Gyr at z=3. The DLA continue to maintain a floor in metallicity of ~1/700 solar independent of observational effects. This metallicity threshold limits the prevalence of primordial gas in high redshift galaxies and stresses the correspondence between damped systems and star formation (i.e. galaxy formation). This floor is significantly offset from the metallicity of the Lya forest and therefore we consider it to be more related to active star formation within these galaxies than scenarios of enrichment in the very early universe. Finally, we comment on an apparent 'missing metals problem': the mean metallicity of the damped systems is ~10x lower than the value expected from their observed star formation history. This problem is evident in current theoretical treatments of chemical evolution and galaxy formation; it may indicate a serious flaw in our understanding of the interplay between star formation and metal production.

This paper presents chemical abundance measurements for 37 damped Lya systems at z>2.5 observed with the Echellette Spectrograph and Imager on the KeckII telescope. We measure the HI column densities of these damped systems through Voigt profile fits to their Lya profiles and we implement the apparent optical depth method to determine ionic column densities. Figures and tables of all relevant data are presented. A full analysis of the chemical enrichment history described by these observations will be presented in a future paper. This dataset is also valuable for efficiently planning future echelle observations and for rough abundance pattern analyses. We aim to make this entire data set public within three years of this publication.

The discovery of metal-poor stars (where metal is any element more massive than helium) has enabled astronomers to probe the chemical enrichment history of the Milky Way. More recently, element abundances in gas inside high-redshift galaxies has been probed through the absorption lines imprinted on the spectra of background quasars, but these have typically yielded measurements of only a few elements. Furthermore, interpretation of these abundances is complicated by the fact that differential incorporation of metals into dust can produce an abundance pattern similar to that expected from nucleosynthesis by massive stars. Here we report the observation of over 25 elements in a galaxy at z=2.626. With these data, we can examine nucleosynthetic processes independent of the uncertainty arising from depletion. We find that the galaxy was enriched mainly by massive stars (M > 15 solar masses) and propose that it is the progenitor of a massive, elliptical galaxy. The detailed abundance patterns suggest that boron is produced through processes that act independently of metallicity, and may require alternative mechanisms for the nucleosynthesis of germanium.

Starting from the SFR per unit physical area, determined for DLAs using the C II* method, we obtain the SFR per unit comoving volume at $z$ $\approx$ 3, and find that it agrees with that for the Lyman Break Galaxies (LBGs). Though the mass of produced stars indicated by the SFRs is consistent with the current densities of known stellar populations, the mass of metals produced by $z$=2.5 is 30 times larger than detected in absorption in DLAs. The most likely solutions to this ``missing metals'' problem is that star formation occurs in compact bulge regions. We search for evidence of feedback and find no correlations between the SFR per unit area and N(H I), but possible correlations between SFR per unit area and low-ion velocity width and SFR per unit area and metal abundance. We show that (a) the correlation between cooling rate and dust-to-gas ratio is positive evidence for grain photoelectric heating, (b) the CMB does not significantly populate the C II excited fine-structure states, and (c) the ratio of CII* to resonance-line optical depths is a sensitive probe of the multi-phase structure of the DLA gas. We address recent arguments that DLAs are comprised only of WNM gas, and show them to be inconclusive. Despite the rough agreement between SFR per unit comoving volume for DLAs and LBGs, current evidence indicates these are distinct populations.

We describe a technique that for the first time measures star formation rates (SFRs) in damped Lyman alpha systems(DLAs) directly. We assume that massive stars form in DLAs, and that the FUV radiation they emit heats the gas by the grain photoelectric mechanism. We infer the heating rate from the cooling rate measured by the strength of CII* 1335.7 absorption. Since the heating rate is proportional to the dust-to-gas ratio and the SFR per unit area, we deduce the SFR per unit area for DLAs in which both quantities have been measured. We consider models in which the the dust comprises carbonaceous or silicate grains. We present two-phase models where the cold neutral medium (CNM) and warm neutral medium (WNM) are in pressuer equilibrium. In the CNM model the the sightline goes throught the CNM and WNM, while in the WNM model it goes only through the WNM. Since the grain photoelectric heating efficiency is at least 10 times higher in the CNM than in the WNM, CII* absorption mainly arises in the CNM in the CNM model. But in the WNM model all of the CII* absorption arises in the WNM. We use CII* absorption lines to derive the SFR per unit area for a sample of ~ 30 DLAs in which the dust-to-gas ratio has been inferred from element depletion patterns. We show that the resulting SFR per unit area corresponds to an average over the star forming volume of galaxy hosting the DLA rather than to local star formation along the line of sight. We find the average SFR per unit area and equals10$^{-2.2}$ M$_{\odot}$yr$^{-1}$kpc$^{-2}$ for the CNM model and 10$^{-1.3}$ M$_{\odot}$yr$^{-1}$kpc$^{-2}$ for the WNM model. The SFR per unit area in the CNM solution is similar to that measured in the Milky Way ISM.

Detailed Monte Carlo Inversion analysis of the spectral lines from three Lyman limit systems (LLS) [N(HI) > 10^{17} cm^-2] and nine lower N(HI) systems [2 10^{14} < N(HI) < 2 10^{16} cm^-2] observed in the VLT/UVES spectra of Q0347-3819 (in the range 2.21 <= z <= 3.14) and of APM BR J0307-4945 (at z = 4.21 and 4.81) is presented. Combined with the results from a previous work, the analyzed LLSs show that they are a heterogeneous population originating in different environments. A functional dependence of the line-of-sight velocity dispersion sigma_v on the absorber size L is confirmed: the majority of the analyzed systems follow the scaling relation sigma_v ~ (N_H L)^{0.3} (with N_H being the total gas column density). This means that most absorbers may be related to virialized systems like galaxies or their halos. Previously noted enhancement of the metal content in small size systems is also confirmed: metallicities of Z ~ (1/3-1/2)Z_solar are found in systems with L < 0.4 kpc, whereas we observe much lower metal abundances in systems with larger linear sizes. For the first time in LLSs, a pronounced [alpha-element/iron-peak] enrichment is revealed: the absorber at z = 2.21 shows [O/Fe] = 0.65+/-0.11, [Si/Fe] = 0.51+/-0.11, and [Mg/Fe] = 0.38+/-0.11. Several absorption systems exhibit characteristics which are very similar to that observed in high-velocity clouds in the Milky Way and may be considered as high-redshift counterparts of Galactic HVCs.

We present 14 N^0 measurements from our HIRES/Keck database of damped Lya abundances. These data are combined with measurements from the recent and past literature to build an homogeneous, uniform set of observations. We examine photoionization diagnostics like Fe^++ and Ar^0 in the majority of the complete sample and assess the impact of ionization corrections on N/alpha and alpha/H values derived from observed ionic column densities of N^0, Si^+, H^0, and S^+. Our final sample of 19 N/alpha, alpha/H pairs appears bimodal; the majority of systems show N/alpha values consistent with metal-poor emission regions in the local universe but a small sub-sample exhibit significantly lower N/alpha ratios. Contrary to previous studies of N/alpha in the damped systems, our sample shows little scatter within each sub-sample. We consider various scenarios to explain the presence of the low N/alpha sightlines and account for the apparent bimodality. We favor a model where at least some galaxies undergo an initial burst of star formation with suppressed formation of intermediate-mass stars. We found a power-law IMF with slope 0.10 or a mass cut of ~5-8 Msolar would successfully reproduce the observed LN-DLA values. If the bimodal distribution is confirmed by a larger sample of measurements, this may present the first observational evidence for a top heavy initial mass function in some early stellar populations.

We investigate the ionization state of the damped Lya system at z=2.62 toward GB1759+7539 through an analysis of ionic ratios sensitive to photoionization: ArI/SII, FeIII/FeII, NII/NI, AlIII/AlII. Approximately half of the metals arise in a mostly neutral velocity component with HI/H > 0.9, based on FeIII/FeII < 0.013. In contrast, the remaining half exhibits FeIII/FeII~0.3 indicative of a partially ionized medium with HI/H~0.5. These conclusions are supported by the observed NII/NI, AlIII/AlII, and ArI/SII ratios. We assess ionization corrections for the observed column densities through photoionization models derived from the CLOUDY software package. In the neutral gas, the ionization corrections are negligible except for ArI. However for the partially ionized gas, element abundance ratios differ from the ionic ratios by 0.1-0.3 dex for (SiII, SII, NiII, AlII)/FeII ratios and more for (NI, ArI)/FeII. Independent of the shape of the photoionizing spectrum and assumptions on the number of ionization phases, these ionization corrections have minimal impact (<0.1dex) on the total metallicity inferred for this damped Lya system. Measurements on the relative elemental abundances of the partially ionized gas, however, have a greater than ~0.15 dex uncertainty which hides the effects of nucleosynthesis and dust depletion. We caution the reader that this damped system is unusual for a number of reasons (e.g. a very low ArI/SII ratio) and we believe its ionization properties are special but not unique. Nevertheless, it clearly shows the value of examining photoionization diagnostics like FeIII/FeII in a larger sample of systems.

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 new chemical abundance measurements of 16 damped Lya systems at z>1.5 and update our previous abundance analyses. The entire database presented here was derived from HIRES observations on the Keck I telescope, reduced with the same software package, and analysed with identical techniques. Altogether, we present a large, homogeneous database of chemical abundance measurements for protogalaxies in the early universe, ideal for studying a number of important aspects of galaxy formation. In addition, we have established an online directory for this database and will continuously update the results.

We present a comprehensive analysis of the damped Lya abundance database presented in the first paper of this series. This database provides a homogeneous set of abundance measurements for many elements including Si, Cr, Ni, Zn, Fe, Al, S, Co, O, and Ar from 38 damped Lya systems with z > 1.5. With little exception, these damped Llya systems exhibit very similar relative abundances. There is no significant correlation in X/Fe with [Fe/H] metallicity and the dispersion in X/Fe is small at all metallicity. We search the database for trends indicative of dust depletion and in a few cases find strong evidence. Specifically, we identify a correlation between [Si/Ti] and [Zn/Fe] which is unambiguous evidence for depletion. We present a discussion on the nucleosynthetic history of the damped Lya systems by focusing on abundance patterns which are minimally affected by dust depletion. We find [Si/Fe] -> +0.25 dex as [Zn/Fe] -> 0 and that the [Si/Fe] values exhibit a plateau of ~+0.3 dex at [Si/H] < -1.5 dex. Together these trends indicate significant alpha-enrichment in the damped Lya systems at low metallicity, an interpretation further supported by the observed O/Fe, S/Fe and Ar/Fe ratios. We also discuss Fe-peak nucleosynthesis and the odd-even effect. To assess the impact of dust obscuration, we present estimates of the dust-to-gas ratios for the damped Lya sightlines and crudely calculate dust extinction corrections. The distribution of extinction corrections suggests the effects of dust obscuration are minimal and that the population of 'missing' damped systems has physical characteristics similar to the observed sample. We update our investigation on the chemical evolution of the early universe in neutral gas. [significantly abridged]

We describe the results of a search for time variability of the fine structure constant, alpha, using absorption systems in the spectra of distant quasars. Three large optical datasets and two 21cm/mm absorption systems provide four independent samples, spanning 23% to 87% of the age of the universe. Each sample yields a smaller alpha in the past and the optical sample shows a 4-sigma deviation: da/a = -0.72 +/- 0.18 x 10^{-5} over the redshift range 0.5 < z < 3.5. We find no systematic effects which can explain our results. The only potentially significant systematic effects push da/a towards positive values, i.e. our results would become more significant were we to correct for them.

An experimental search for variation in the fundamental coupling constants is strongly motivated by modern high-energy physics theories. Comparison of quasar absorption line spectra with laboratory spectra provides a sensitive probe for variability of the fine structure constant, alpha, over cosmological time-scales. We have previously developed and applied a new method providing an order of magnitude gain in precision over previous optical astrophysical constraints. Here we extend that work by including new quasar spectra of damped Lyman-alpha absorption systems. We also re-analyse our previous lower redshift data and confirm our initial results. The constraints on alpha come from simultaneous fitting of absorption lines of subsets of the following species: Mg I, Mg II, Al II, Al III, Si II, Cr II, Fe II, Ni II and Zn II. We present a detailed description of our methods and results based on an analysis of 49 quasar absorption systems (towards 28 QSOs) covering the redshift range 0.5 < z < 3.5. There is statistical evidence for a smaller alpha at earlier epochs: da/a = (-0.72 +/- 0.18) * 10^{-5}. The new and original samples are independent but separately yield consistent and significant non-zero values of da/a. We summarise the results of a thorough investigation of systematic effects published in a companion paper. The value we quote above is the raw value, not corrected for any of these systematic effects. The only significant systematic effects so far identified, if removed from our data, would lead to a more significant deviation of da/a from zero.

We compare new observations on the kinematic characteristics of the damped Lya systems against results from numerical SPH simulations to test the predictions of hierarchical galaxy formation. This exercise is particularly motivated by recent numerical results on the cross-section of damped Lya systems. Our analysis focuses on the velocity widths of ~50 low-ion absorption profiles from our sample of z>1.5 damped Lya systems. The results indicate that current numerical simulations fail to match the damped Lya observations at high confidence levels (>99.9%). Although we do not believe that our results present an insurmountable challenge to the paradigm of hierarchical cosmology, the damped Lya observations suggest that current numerical SPH simulations overlook an integral aspect of galaxy formation.

Comparison of quasar absorption line spectra with laboratory spectra provides a precise probe for variability of the fine structure constant, alpha, over cosmological time-scales. We constrain variation in alpha in 21 Keck/HIRES Si IV absorption systems using the alkali doublet (AD) method in which changes in alpha are related to changes in the doublet spacing. The precision obtained with the AD method has been increased by a factor of 3: da/a = (-0.5 +/- 1.3) * 10^{-5}. We also analyse potential systematic errors in this result. Finally, we compare the AD method with the many-multiplet method which has achieved an order of magnitude greater precision and we discuss the future of the AD method.

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 report the measurement of the primordial D/H abundance ratio towards QSO \object. The column density of the hydrogen in the $z \simeq 2.536$ Lyman limit system is high, \lnhi $= 19.422 \pm 0.009$ \cmm, allowing for the deuterium to be seen in 5 Lyman series transitions. The measured value of the D/H ratio towards QSO \object is found to be D/H$ = 2.54 \pm 0.23 \times 10^{-5}$. The metallicity of the system showing D/H is found to be $\simeq 0.01$ solar, indicating that the measured D/H is the primordial D/H within the measurement errors. The gas which shows D/H is neutral, unlike previous D/H systems which were more highly ionized. Thus, the determination of the D/H ratio becomes more secure since we are measuring it in different astrophysical environments, but the error is larger because we now see more dispersion between measurements. Combined with prior measurements of D/H, the best D/H ratio is now D/H$ = 3.0 \pm 0.4 \times 10^{-5}$, which is 10% lower than the previous value. The new values for the baryon to photon ratio, and baryonic matter density derived from D/H are $\eta = 5.6 \pm 0.5 \times 10^{-10} $ and \ob $=0.0205 \pm 0.0018$ respectively.

We measure feedback (heating rates) in damped Lyman alpha systems from the cooling rate of the neutral gas. Since cooling occurs through [C II] 158 micron emission, we infer cooling from C II^{*} 1335.7 absorption lines detected with HIRES on the Keck I telescope. The inferred heating rates are about 30 times lower than for the Galaxy ISM. At z = 2.8, the implied star formation rate per unit area is 10^{-2.4+-0.3} solar masses per kpc^{2} per year, and the the star formation rate per unit comoving volume is 10^{-0.8+-0.2} solar masses per Mpc^{3} per year. This is the first measurement of star formation rates in objects likely to be the progenitors of current galaxies.

We present the first results from an ongoing survey to discover and measure the metallicity of z>3 damped Lya systems with the Echellette Spectrograph and Imager (ESI) on the Keck II telescope. Our motivation arises from a recent study on the damped Lya systems suggesting only mild evolution in the cosmic metallicity from z~2 to 4. The Echellette Spectrograph and Imager, which provides two complementary spectroscopic modes, is the ideal instrument for a z>3 damped Lya survey. We describe our observing strategy and report on the discovery and analysis of 5 new z>3 damped Lya systems acquired in a single night of observing. These observations further support the principal conclusions of the previous study: (1) the cosmic metallicity in neutral gas inferred from the damped Lya systems does not evolve significantly from z~2 to 4; (2) the unweighted metallicity exhibits a statistically significant decrease with increasing redshift; and (3) not a single damped Lya system has a metallicity below [Fe/H]=-3. We discuss the implications of these results and comment on recent theoretical studies which attempt to explain the observations.

We investigate the kinematics of ionized and neutral gas in a sample of 35 damped Lyman alpha protogalaxies (DLAs) with HIRES on the Keck I 10 m telescope. Velocity profiles with resolution of ~ 8 km/s are obtained for high ions such as C IV and Si IV, and for intermediate ions such as Al III. Combining these profiles with accurate low-ion (e.g., Fe II) profiles, we investigate the kinematic state of DLAs in the redshift range 1.8 < z < 4.4 by comparisons between data for various ion pairs. We find the DLAs comprise distinct kinematic subsystems: a low ion subsystem in which the low ions are physically associated with intermediate ions, and a high subsystem containing neither low nor intermediate ions. The evidence for two subsystems stems from (a) differences between the widths of the velocity profiles, (b) misalignment in velocity space between the narrow components comprising the profiles in each subsystem, and (c) significant dissimilarities between the mean velocities of the high ion and low ion velocity profiles. In every case we find that test statistics such as velocity width and various asymmetry parameters distribute differently for low ions than for high ions. Despite misalignment in velocity space, the low and high ion kinematic subsystems are interrelated. This is indicated by detection of a statistically significant C IV versus low-ion cross correlation function, and by a systematic effect where the C IV velocity widths are greater than or equal to the low ion velocity widths in 29 out of 32 systems. These phenomena are consistent with the location of the two subsystems in the same potential well.

We test semi-analytic models for galaxy formation with accurate kinematic data of damped Lyman alpha protogalaxies (DLAs) presented in the companion paper I. The models envisage centrifugally supported exponential disks at the centers of dark matter halos which are filled with ionized gas undergoing radial infall to the disks. The halo masses are drawn from cross-section weighted mass distributions predicted by CDM cosmogonies, or by the null hypothesis (TF model) that the dark matter mass distribution has not evolved since z ~ 3. In our models, C IV absorption lines detected in DLAs arise in infalling ionized clouds while the low-ion absorption lines arise from neutral gas in the disks. Using Monte Carlo methods we find: (a) The CDM models are incompatible with the low-ion statistics at more than 99% confidence whereas some TF models cannot be excluded at more than 88% confidence. (b) Both CDM and TF models agree with the observed distribution of C IV velocity widths. (c) The CDM models generate differences between the mean velocities of C IV and low ion profiles in agreement with the data, while the TF model produces differences in the means that are too large. (d) Both CDM and TF models produce ratios of C IV to low-ion velocity widths that are too large. (e) Both CDM and TF models generate C IV versus low-ion cross-correlation functions incompatible with the data. While it is possible to select model parameters resulting in consistency with the data, the disk-halo configuration assumed in both cosmogonies still does not produce significant overlap in velocity space between C IV low-ion velocity profiles. We conjecture that including angular momentum of the infalling clouds will increase the overlap between C IV and low-ion profiles.

Evidence is presented that the damped Lyman alpha absorption systems are the high-redshift (z > 3) progenitors of galaxy disks. I discuss kinematic evidence that the damped Lyman Alpha systems are rotating disks. I also discuss implications of the lack of metal-poor damped Lyman alpha systems with line width Delta v > 100 {\kms}. I then present new evidence stemming from correlations between element-abundance ratios and [Fe/H], which connects damped systems to the thick stellar disk of the Galaxy. I discuss the connections between damped Lyman alpha systems and Lyman break galaxies, and how [CII] 158 micron emission from damped Lyman alpha systems discriminates among competing theories of galaxy formation. ~

We present first results from a program to measure the chemical abundances of a large (N>30) sample of thick disk stars with the principal goal of investigating the formation history of the Galactic thick disk. Our analysis confirms previous studies of O and Mg in the thick disk stars which reported enhancements in excess of the thin disk population. Furthermore, the observations of Si, Ca, Ti, Mn, Co, V, Zn, Al, and Eu all argue that the thick disk population has a distinct chemical history from the thin disk. With the exception of V and Co, the thick disk abundance patterns match or tend towards the values observed for halo stars with [Fe/H]~-1. This suggests that the thick disk stars had a chemical enrichment history similar to the metal-rich halo stars. With the possible exception of Si, the thick disk abundance patterns are in excellent agreement with the chemical abundances observed in the metal-poor bulge stars suggesting the two populations formed from the same gas reservoir at a common epoch. We discuss the implications for the formation of the thick disk, new insights presented for theories of nucleosynthesis, and the impact on interpretations of the abundance patterns of the damped Lya systems at high redshift. [significantly abridged]

We have completed spectroscopic observations using LRIS on the Keck 1 telescope of 30 very high redshift quasars, 11 selected for the presence of damped Ly-alpha absorption systems and 19 with redshifts z > 3.5 not previously surveyed for absorption systems. We have surveyed an additional 10 QSOs with the Lick 120'' and the Anglo-Australian Telescope. We have combined these with previous data resulting in a statistical sample of 646 QSOs and 85 damped Ly-alpha absorbers with column densities N(HI) >= 2 x 10^20 atoms/cm^2 covering the redshift range 0.008 <= z <= 4.694. To make the data in our statistical sample more readily available for comparison with scenarios from various cosmological models, we provide tables that includes all 646 QSOs from our new survey and previously published surveys. They list the minimum and maximum redshift defining the redshift path along each line of sight, the QSO emission redshift, and when an absorber is detected, the absorption redshift and measured HI column density. [see the paper for the complete abstract]

Observations of the damped Lya systems provide direct measurements on the chemical enrichment history of neutral gas in the early universe. In this Letter, we present new measurements for four damped Lya systems at high redshift. Combining these data with [Fe/H] values culled from the literature, we investigate the metallicity evolution of the universe from z~1.5-4.5. Contrary to our expectations and the predictions of essentially every chemical evolution model, the N(HI)-weighted mean [Fe/H] metallicity exhibits minimal evolution over this epoch. For the individual systems, we report tentative evidence for an evolution in the unweighted [Fe/H] mean and the scatter in [Fe/H] with the higher redshift systems showing lower scatter and lower typical [Fe/H] values. We also note that no damped Lya system has [Fe/H] < -2.7 dex. Finally, we discuss the potential impact of small number statistics and dust on our conclusions and consider the implications of these results on chemical evolution in the early universe.

We present chemical abundance measurements for 19 damped lya systems observed with HIRES on the 10m W.M. Keck Telescope. Our principal goal is to investigate the abundance patterns of the damped systems and thereby determine the underlying physical processes which dominate their chemical evolution. We place particular emphasis on gauging the relative importance of two complementary effects often invoked to explain the damped lya abundances: (1) nucleosynthetic enrichment from Type II supernovae and (2) an ISM-like dust depletion pattern. Similar to the principal results of Lu et al. (1996), our observations lend support both for dust depletion and Type II SN enrichment. Specifically, the observed overabundance of Zn/Fe and underabundance of Ni/Fe relative to solar abundances suggest significant dust depletion within the damped lya systems. Meanwhile, the relative abundances of Al, Si, and Cr vs. Fe are consistent with both dust depletion and Type II supernova enrichment. Our measurements of Ti/Fe and the Mn/Fe measurements from Lu et al. (1996), however, cannot be explained by dust depletion and indicate an underlying Type II SN pattern. Finally, the observed values of [S/Fe] are inconsistent with the combined effects of dust depletion and the nucleosynthetic yields expected for Type II supernovae. This last result emphasizes the need for another physical process to explain the damped lya abundance patterns. We also examine the metallicity of the damped lya systems both with respect to Zn/H and Fe/H. Our results confirm previous surveys by Pettini and collaborators, i.e., [<Zn/H>] = -1.15 +/- 0.15 dex. [abridged]

We present constraints on the deuterium to hydrogen ratio (D/H) in the metal-poor gas cloud at redshift $z=0.701$ towards QSO 1718+4807. We use new Keck spectra in addition to Hubble Space Telescope (HST) and International Ultraviolet Explorer (IUE) spectra. We use an improved redshift and a lower \HI column density to model the absorption. The HST spectrum shows an asymmetric Lyman-$\alpha$ (\lya) feature which is produced by either \HI at a second velocity, or a high abundance of D. Three models with a single simple H+D component give $8 \times 10^{-5} < D/H < 57 \times 10^{-5}$ (95%), a much larger range than reported by Webb et al (1997a,b). A more sophisticated velocity distribution, or a second component is necessary for lower D/H. With two components, which could be a part of one absorbing structure, or separate clouds in a galaxy halo, we find $D/H < 50 \times 10^{-5}$. We do not know if this second component is present, but it is reasonable because 40 -- 100% of absorption systems with similar redshifts and \HI column densities have more than one component.

We present a catalogue of absorption lines obtained from the analysis of the ultra-violet spectra of 66 quasars. The data were acquired with the Faint Object Spectrograph of the HST as part of the Quasar Absorption Line Survey, a Key Project for the first four cycles of HST observations. This is the third of a series of catalogues of absorption lines produced from the survey and increases the number of quasars whose higher resolution (R=1300) spectra we have published from 17 to 83. The general properties and execution of the survey are reviewed, including descriptions of the final sample of observed objects and the algorithmic processes used to construct the catalogue. The detection of a single damped Ly-a system in a path length of Delta_z=49 yields an observed number of damped systems per unit redshift of (dN/dz)_{damp}(z=0.58)=0.020 with 95% confidence boundaries of 0.001 to 0.096 systems per unit redshift. We include notes on our analysis of each of the observed quasars and the absorption systems detected in each spectrum. Some especially interesting systems include low redshift Ly-a absorbers suitable for extensive follow-up observations (in the spectra of TON28 and PG1216+069), possibly physically associated pairs of extensive metal line absorption systems (e.g., in the spectrum of PG0117+213), and systems known to be associated with galaxies (e.g., in the spectrum of 3C232).

We present an analysis of the evolution of the Ly alpha absorption lines in the redshift interval 0.0 to \sim 1.5 based upon a sample of 987 lines identified in the spectra of 63 QSOs obtained with the Faint Object Spectrograph on the Hubble Space Telescope. Fits to the number of absorbers per unit redshift dN/dz = A \times (1+z)^\gamma yield values of \gamma in the range 0.1 to 0.3, decidedly flatter than results from groundbased data pertaining to the redshift range z > 1.7. The combination of the space and groundbased data suggest a marked transition in the rate of evolution of the Ly alpha lines at a redshift of about 1.7. The distribution in redshift is well represented by a power law and the distribution of equivalent widths above a rest equivalent width of 0.1 \AA is fit well by an exponential. The sample of Ly alpha systems which has ions from heavy elements has a much larger value of \gamma than a sample of lines of comparable strength without heavy elements. We argue that this is not likely due to either line misidentification or incomplete spectral coverage and suggest that it is due to rapid evolution of the very highest column density systems, rather than differences in metallicity. We also find evidence that the rate of evolution increases with increasing equivalent width. We compare our results for the variation of line density with redshift to recent numerical simulations of Ly alpha absorbers, and find fairly good agreement. Our results for the stronger Ly alpha lines in our sample are compatible with the recent analysis of the absorber--galaxy correlation by Chen et al. 1998 but for the weaker lines our results suggest that the association between absorbers and galaxies is different from that for the stronger lines.

We present new observational results on the kinematics of the damped lya systems. Our full sample is now comprised of 31 low-ion profiles and exhibits similar characteristics to the sample from Paper I. The primary exception is that the new distribution of velocity widths includes values out to a maximum of nearly 300 km/s, approx 100 km/s greater than the previous maximum. These high velocity width systems will significantly leverage models introduced to explain the damped lya systems. Comparing the characteristics from low-redshift and high-redshift sub-samples, we find no evidence for significant evolution in the kinematic properties of protogalaxies from z = 2.0 - 3.3. The new observations give greater statistical significance to the main conclusions of our first paper. In particular, those models inconsistent with the damped lya observations in Paper I are ruled out at even higher levels of confidence. At the same time, the observations are consistent with a population of rapidly rotating, thick disks (the TRD model) at high redshift. Buoyed by the success of the TRD model, we investigate it more closely by considering more realistic disk properties. Our goal is to demonstrate the statistical power of the damped lya observations by investigating the robustness of the TRD model. In particular, we study the effects of warping, realistic rotation curves, and photoionization on the kinematics of disks in the TRD model. The principal results are: (1) disk warping has only minimal effect on the kinematic results, primarily influencing the effective disk thickness, (2) the TRD model is robust to more realistic rotation curves; (3) the effects of photoionization require thicker disks to give consistent velocity width distributions. [abridged]

The mean metallicity of the thick disk of the Galaxy is 0.5 dex higher than that of the damped Lyman alpha systems. This has been interpreted to argue that stars in the former do not arise out of gas in the latter. Using new metallicity and H I column-density data we show the metal-rich damped systems do contain sufficient baryons at the thick-disk metallicity to account for the stellar masses of thick disks. Comparing our kinematic data with the metallicities we show that damped Lyman alpha systems exhibiting the largest profile velocity widths span a narrow range of high metallicities, while systems with small velocity widths span a wider range of metallicities. This is naturally explained by passage of the damped Lyman alpha sightlines through rapidly rotating disks with negative radial gradients in metallicity. The systematically lower N(H I) of systems with high velocity widths indicates (a) the gaseous disks have centrally located holes, and (b) an apparent inconsistency with the protogalactic clump model for damped Lyman alpha systems. The higher metallicity of systems with low N(H I) further implies that stars rather than gas dominate the baryonic content of the most metal-rich damped systems.

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.

We present the first results of an ongoing program to investigate the kinematic characteristics of high redshift damped lya systems. Because damped lya systems are widely believed to be the progenitors of current massive galaxies, an analysis of their kinematic history allows a direct test of galaxy formation scenarios. We have collected a kinematically unbiased sample of 17 high S/N ratio, high resolution damped lya spectra taken with HIRES on the 10m W.M. Keck Telescope. Our study focuses on the unsaturated, low-ion transitions of these systems which reveal their kinematic traits. The profiles exhibit a nearly uniform distribution of velocity widths ranging from 20 - 200 km/s and a relatively high degree of asymmetry. In an attempt to explain these characteristics, we introduce several physical models, which have previously been attributed to damped lya systems, including rapidly rotating cold disks, slowly rotating hot disks, massive isothermal halos, and a hydrodynamic spherical accretion model. Using standard Monte Carlo techniques, we run sightlines through these model systems to derive simulated low-ion profiles. Comparing statistical measures of the simulated profiles with the observed profiles, we determine that the rapidly rotating cold disk model is the only tested model consistent with the data at high confidence levels. A Relative Likelihood Test of the rapidly rotating cold disk model indicates the disks must have large rotation speeds; v > 180 km/s at the 99% c.l. In turn, we demonstrate that the Cold Dark Matter Model, as developed by Kauffmann (1996), is inconsistent with the damped lya data at very high c.l. This is because the CDM Model does not predict a large enough fraction of rapidly rotating disks at z approx 2.5.

Radio spectroscopy at 323 MHz using the Arecibo Telescope and the VLA has produced a tentative detection of the 21cm line of neutral hydrogen at Z_abs=3.38716+-0.00007 in absorption against the radio continuum of the QSO MG0201+113 (Z_em = 3.61). This redshift roughly agrees with one determined by the Westerbork Synthesis Radio Telescope and reported earlier; however the observations at different telescopes produce different results for width and optical depth of the line. If the detection holds, it provides evidence for a high column density of neutral gas that is confined to a dynamically cold layer with velocity dispersion ~10 km/s. Although the interpretation is uncertain due to a lack of detailed knowledge of the extended radio structure of the background quasar and the relative quantities of neutral gas in the cold and turbulent components, the observations specify high spin temperatures, T_s > 1000 K for both this 21cm line absorbing cloud and the turbulent component, which together have N_HI ~ 10^21.4 cm^-2. Results of optical spectroscopy require the additional presence of metal-enriched clouds of still broader velocity dispersion than the 21cm line, creating a picture which is consistent with this system being a young disk galaxy that is sufficiently evolved by z=3.4 to have collapsed to a flattened system and produced a population of stars that have aged to pollute a turbulent halo. The observations constrain the neutral gas mass of a possible 21cm line emitter associated with the intervening absorber to be less than ~(Delta V/200 km/s) 10^13 Msolar for velocity widths Delta V (FWHM) in the range 200 to 1200 km/s (H_o=100 km/s/Mpc, Omega_o=1).

We report the discovery of a high-ionization broad absorption line system at a redshift of z_abs=0.695 in the spectrum of the z_em=1.052 radio-quiet quasar PG 2302+039. Broad absorption with FWHM from 3,000 to 5,000 km/s is detected from C IV, N V, and O VI in HST Faint Object Spectrograph spectra of the quasar. A narrow line system (FWHM approximately 250 km/s) at z_abs=0.7016 is resolved from the broad blend and includes absorption by Ly-alpha and the C IV, N V, and O VI doublets. No absorption by low-ionization metal species (e.g. Si II and Mg II) is detected in the HST or ground-based spectra for either the broad or the narrow system. The centroids of the broad system lines are displaced by approximately 56,000 km/s to the blue of the quasar's broad emission lines. The reddest extent of the broad line absorption is more than 50,000 km/s from the quasar. The properties of this system are unprecedented whether it is an intervening or an ejected system.

We present high resolution, high SNR spectra of the QSO Q2206-199 obtained with HIRES on the 10m W.M. Keck Telescope. Our analysis focuses on the two previously identified damped \lya systems found at $z=1.920$ and $z=2.076$. For each system, we measure accurate abundances. The $z=1.920$ system exhibits the highest metallicity we have measured for a damped \lya system. We report the first confident ($>5 \sigma$) detection of Ti in a QSO absorption line system. By contrast the $z=2.076$ system is the most metal poor we have analyzed, showing absorption features for only the strongest transitions. We find no positive evidence for the presence of dust in either system. The two damped systems exhibit significantly different kinematic characteristics, yet we contend the two systems are consistent with one physical description: that of a thick, rotating disk. We investigate a very strong Mg II system at $z=0.752$ which is very likely yet a third damped \lya system. The very weak Mn II and Ti II transitions have been positively measured and imply $\log \N{HI} > 19.0$. We analyze the abundance ratios [Mn/Fe] and [Ti/Fe] and their values are inconsistent with dust depletion, yet consistent with the abundance pattern detected for halo stars in the Galaxy (see Lu et al. 1996a). Finally, we identify a C IV system at $z=2.014$ that shows a very narrow feature in Si IV and C IV absorption. The corresponding $b$ values (5.5 \kms and 8.9 \kms for Si IV and C IV) for this component suggest a temperature of $4.7 \sci{4} \rm K$. Because collisional ionization can explain the observed abundances only for $T > 8 \sci{4} \rm K$, we contend these ions must have formed through a different physical process (e.g. photoionization).

We present high resolution ($\approx 8$ \kms) spectra of the QSO Q0201+365 obtained with HIRES, the echelle spectrograph on the 10m W.M. Keck Telescope. Although we identify over $80\%$ of the absorption features and analyze several of the more complex metal-line systems, we focus our analysis on the damped \Lya system at $z=2.462$. Ionization simulations suggest the hydrogen in this system is significantly neutral and all of the observed metals are predominantly singly ionized. We measure accurate abundances for Fe, Cr, Si, Ni and place a lower limit on the abundance of Zn: [Fe/H] = $-0.830 \pm 0.051$, [Cr/H] = $-0.902 \pm 0.064$, [Si/H] = $-0.376 \pm 0.052$, [Ni/H] = $-1.002 \pm 0.054$ and [Zn/H] $> -0.562 \pm 0.064$. We give evidence suggesting the actual Zn abundance is [Zn/H] $\approx -0.262$, implying the highest metallicity observed at a redshift $z \geq 2$. The relative abundances of these elements remains constant over essentially the entire system ($\approx 150$ \kms in velocity space), suggesting it is well mixed. Furthermore, we use the lack of abundance variations to infer properties of the dust responsible for element depletion. Finally, we discuss the kinematic characteristics of this damped \Lya system, comparing and contrasting it with other systems. The low-ion line profiles span $\approx 200$ \kms in velocity space and have an asymmetric shape with the strongest feature on the red edge. These kinematic characteristics are consistent with a rotating disk model.

We present evidence that clumps of \lya lines are physically associated with about half of the extensive metal-line systems (absorption systems with four or more observed metal-line species) found in this paper, demonstrate that all four \lylimit systems discussed here correspond to extensive metal-line absorption systems, and present an extraordinary pair of extensive metal-line absorption systems within 2000~km/s of each other at $z~=~0.95$ that are probably an early manifestation of large scale structure. These results are obtained using ultraviolet spectra, taken with the higher-resolution gratings of the Faint Object Spectrograph of the Hubble Space Telescope, for four quasars with emission-line redshifts between 1.0 and 1.3. We also determine the evolution of \lya absorption lines at redshifts less than 1.3 by combining the results for 13 smaller redshift quasars discussed in Paper I of this series with the 4 moderate redshift quasars analyzed in the present paper. For $z_{\rm abs}~\leq~1.3$, the density of \lya lines with equivalent widths greater than 0.24~\AA\ is adequately fit by $\left(dN/dz\right) =\left(dN/dz\right)_0 \cdot(1 + z)^\gamma$ with $(dN/dz)_0 = 24.3 \pm 6.6$ \lya lines per unit redshift, and $\gamma = 0.58 \pm 0.50$ (1-$\sigma$ uncertainties). This rate of evolution at low redshifts is less than the evolutionary rate inferred from several different ground-based data samples that pertain to high redshifts. The observed gaseous structures at redshifts of $0.5$ to $1.0$ with velocity dispersions of $6 \times 10^2$ \kms to $1.4 \times 10^3$ \kms (or velocity spans of $1.2 \times 10^3$ \kms to $3 \times 10^3$ \kms) constitute a constraint on cosmological models of structure formation.