Project AMIGA (Absorption Maps In the Gas of Andromeda) is a large ultraviolet Hubble Space Telescope program, which has assembled a sample of 43 QSOs that pierce the circumgalactic medium (CGM) of Andromeda (M31) from R=25 to 569 kpc (25 of them probing gas from 25 kpc to about the virial radius-Rvir = 300 kpc-of M31). Our large sample provides an unparalleled look at the physical conditions and distribution of metals in the CGM of a single galaxy using ions that probe a wide range of gas phases (Si II, Si III, Si IV, C II, C IV, and O VI, the latter being from the Far Ultraviolet Spectroscopic Explorer). We find that Si III and O VI have near unity covering factor maintained all the way out to 1.2Rvir and 1.9Rvir, respectively. We show that Si III is the dominant ion over Si II and Si IV at any R. While we do not find that the properties of the CGM of M31 depend strongly on the azimuth, we show that they change remarkably around 0.3-0.5Rvir, conveying that the inner regions of the CGM of M31 are more dynamic and have more complicated multi-phase gas-structures than at R>0.5Rvir. We estimate the metal mass of the CGM within Rvir as probed by Si II, Si III, and Si IV is 2x10^7 Msun and by O VI is >8x10^7 Msun, while the baryon mass of the 10^4-10^5.5 K gas is ~4x10^10 (Z/0.3 Zsun)^(-1) Msun within Rvir. We show that different zoom-in cosmological simulations of L* galaxies better reproduce the column density profile of O VI with R than Si III or the other studied ions. We find that observations of the M31 CGM and zoom-in simulations of L* galaxies have both lower ions showing higher column density dispersion and dependence on R than higher ions, indicating that the higher ionization structures are larger and/or more broadly distributed.

We address the spatial scale, ionization structure, mass and metal content of gas at the Milky Way disk-halo interface detected as absorption in the foreground of seven closely-spaced, high-latitude halo blue horizontal branch stars (BHBs) with heights z = 3 - 14 kpc. We detect transitions that trace multiple ionization states (e.g. CaII, FeII, SiIV, CIV) with column densities that remain constant with height from the disk, indicating that the gas most likely lies within z < 3.4 kpc. The intermediate ionization state gas traced by CIV and SiIV is strongly correlated over the full range of transverse separations probed by our sightlines, indicating large, coherent structures greater than 1 kpc in size. The low ionization state material traced by CaII and FeII does not exhibit a correlation with either N$_{\rm HI}$ or transverse separation, implying cloudlets or clumpiness on scales less than 10 pc. We find that the observed ratio log(N_SiIV/ N_CIV), with a median value of -0.69+/-0.04, is sensitive to the total carbon content of the ionized gas under the assumption of either photoionization or collisional ionization. The only self-consistent solution for photoionized gas requires that Si be depleted onto dust by 0.35 dex relative to the solar Si/C ratio, similar to the level of Si depletion in DLAs and in the Milky Way ISM. The allowed range of values for the areal mass infall rate of warm, ionized gas at the disk-halo interface is 0.0003 < dM_gas / dtdA [M_sun kpc^-2 yr^-] < 0.006. Our data support a physical scenario in which the Milky Way is fed by complex, multiphase processes at its disk-halo interface that involve kpc-scale ionized envelopes or streams containing pc-scale, cool clumps.

We present the design, methods, and first results of the MUSE Analysis of Gas around Galaxies (MAGG) survey, a large programme on the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope (VLT) which targets 28 z > 3.2 quasars to investigate the connection between optically-thick gas and galaxies at z~3-4. MAGG maps the environment of 52 strong absorption line systems at z > 3, providing the first statistical sample of galaxies associated with gas-rich structures in the early Universe. In this paper, we study the galaxy population around a very metal poor gas cloud at z~3.5 towards the quasar J124957.23-015928.8. We detect three Lyman alpha emitters within <200km/s of the cloud redshift, at projected separations <185 kpc (physical). The presence of star-forming galaxies near a very metal-poor cloud indicates that metal enrichment is still spatially inhomogeneous at this redshift. Based on its very low metallicity and the presence of nearby galaxies, we propose that the most likely scenario for this LLS is that it lies within a filament which may be accreting onto a nearby galaxy. Taken together with the small number of other LLSs studied with MUSE, the observations to date show a range of different environments near strong absorption systems. The full MAGG survey will significantly expand this sample and enable a statistical analysis of the link between gas and galaxies to pin down the origin of these diverse environments at z~3-4.

We characterize the metallicities and physical properties of cool, photoionized gas in a sample of 152 z<1 strong Lya forest systems (SLFSs, absorbers with 15<log N(HI)<16.2). The sample is drawn from our COS circumgalactic medium (CGM) compendium (CCC), an ultraviolet survey of HI-selected circumgalactic gas around z<1 galaxies that targets 262 absorbers with 15<log N(HI)<19. We show that the metallicity probability distribution function of the SLFSs at z<1 is unimodal, skewed to low metallicities with a mean and median of [X/H]=-1.47$ and -1.18 dex. Very metal-poor gas with [X/H]<-1.4 represents about half of the population of absorbers with 15<log N(HI)<18. Thus, there are important reservoirs of primitive (though not pristine) gas around z<1 galaxies. The photoionized gas around z<1 galaxies is highly inhomogeneous based on the wide range of metallicities observed (-3<[X/H]<+0.4) and that there are large metallicity variations (factors of 2 to 25) for most of the closely-spaced absorbers (Dv<300 km/s) along the same sightlines. These absorbers show a complex evolution with redshift and HI column density, and we identify subtle cosmic evolution effects that affect the interpretation of metallicity distributions and comparison with other of absorbers samples. We discuss the physical conditions and cosmic baryon and metal budgets of the CCC absorbers. Finally, we compare the CCC results to recent cosmological zoom simulations and explore the origins of the 15<log N(HI)<19 absorbers within the EAGLE high-resolution simulations.

We present a search for HI in the circumgalactic medium (CGM) of 21 massive ($\langle \log M_\star \rangle \sim 11.4$), luminous red galaxies (LRGs) at $z\sim0.5$. Using UV spectroscopy of QSO sightlines projected within 500 kpc ($\sim R_{vir}$) of these galaxies, we detect HI absorption in 11/21 sightlines, including two partial Lyman limit systems and two Lyman limit systems. The covering factor of $\log N(HI) \ge 16.0$ gas within the virial radius of these LRGs is $f_c(\rho \le R_{vir}) = 0.27^{+0.11}_{-0.10}$, while for optically-thick gas ($\log N(HI) \ge 17.2$) it is $f_c(\rho \le R_{vir}) = 0.15^{+0.10}_{-0.07}$. Combining this sample of massive galaxies with previous galaxy-selected CGM studies, we find no strong dependence of the HI covering factor on galaxy mass, although star-forming galaxies show marginally higher covering factors. There is no evidence for a critical mass above which dense, cold ($T \sim 10^4$ K) gas is suppressed in the CGM of galaxies (spanning stellar masses $9.5 \lesssim \log M_\star \lesssim 11.8$). The metallicity distribution in LRGs is indistinguishable from those found about lower-mass star-forming galaxies, and we find low-metallicity gas with $[{\rm X/H}] \approx -1.8$ (1.5% solar) and below about massive galaxies. About half the cases show super-solar [FeII/MgII] abundances as seen previously in cool gas near massive galaxies. While the high-metallicity cold gas seen in LRGs could plausibly result from condensation from a corona, the low-metallicity gas is inconsistent with this interpretation.

We present a novel absorption line survey using 54 blue horizontal branch stars (BHBs) in the Milky Way halo as background sources for detecting gas flows at the disk-halo interface. Distance measurements to high-latitude ($b$ > 60{\deg}) background stars at 3.1-13.4 kpc, combined with unprecedented spatial sampling and spectral resolution, allow us to examine the 3-dimensional spatial distribution and kinematics of gas flows near the disk. We detect absorption signatures of extraplanar CaII and NaI in Keck HIRES spectra and find that their column densities exhibit no trend with distance to the background sources, indicating that these clouds lie within 3.1 kpc of the disk. We calculate covering fractions of $f_{CaII}$ = 63%, $f_{NaI}$ = 26%, and $f_{HI}$ = 52%, consistent with a picture of the CGM that includes multi-phase clouds containing small clumps of cool gas within hotter, more diffuse gas. Our measurements constrain the scale of any substructure within these cool clouds to <0.5 kpc. CaII and NaI absorption features exhibit an intermediate-velocity (IV) component inflowing at velocities of -75 km/s < v < -25 km/s relative to the local standard of rest, consistent with previously-studied HI structures in this region. We report the new detection of an inflow velocity gradient ${\Delta}v_z$ ~ 6-9 km/s/kpc across the Galactic plane. These findings place constraints on the physical and kinematic properties of CGM gas flows through the disk-halo interface, and support a galactic fountain model in which cold gas rains back onto the disk.

High-multiplex and deep spectroscopic follow-up of upcoming panoramic deep-imaging surveys like LSST, Euclid, and WFIRST is a widely recognized and increasingly urgent necessity. No current or planned facility at a U.S. observatory meets the sensitivity, multiplex, and rapid-response time needed to exploit these future datasets. FOBOS, the Fiber-Optic Broadband Optical Spectrograph, is a near-term fiber-based facility that addresses these spectroscopic needs by optimizing depth over area and exploiting the aperture advantage of the existing 10m Keck II Telescope. The result is an instrument with a uniquely blue-sensitive wavelength range (0.31-1.0 um) at R~3500, high-multiplex (1800 fibers), and a factor 1.7 greater survey speed and order-of-magnitude greater sampling density than Subaru's Prime Focus Spectrograph (PFS). In the era of panoramic deep imaging, FOBOS will excel at building the deep, spectroscopic reference data sets needed to interpret vast imaging data. At the same time, its flexible focal plane, including a mode with 25 deployable integral-field units (IFUs) across a 20 arcmin diameter field, enables an expansive range of scientific investigations. Its key programmatic areas include (1) nested stellar-parameter training sets that enable studies of the Milky Way and M31 halo sub-structure, as well as local group dwarf galaxies, (2) a comprehensive picture of galaxy formation thanks to detailed mapping of the baryonic environment at z~2 and statistical linking of evolving populations to the present day, and (3) dramatic enhancements in cosmological constraints via precise photometric redshifts and determined redshift distributions. In combination with Keck I instrumentation, FOBOS also provides instant access to medium-resolution spectroscopy for transient sources with full coverage from the UV to the K-band.

Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.

We survey the highly ionized circumgalactic media (CGM) of 29 blindly selected galaxies at 0.49 < z_(gal) < 1.44 based on high-S/N ultraviolet spectra of z > 1 QSOs and the galaxy database from the COS Absorption Survey of Baryon Harbors (CASBaH). We detect the Ne VIII doublet in nine of the galaxies, and for gas with N(Ne VIII) > 10^13.3 cm^-2 (> 10^13.5 cm^-2), we derive a Ne VIII covering fraction f_c = 75 +15/-25% (44 +22/-20%) within impact parameter (rho) < 200 kpc of M_* = 10^(9.5-11.5) Msol galaxies and f_c = 70 +16/-22% (f_c = 42 +20/-17%) within rho < 1.5 virial radii. We estimate the mass in Ne VIII-traced gas to be M_gas(Ne VIII) > 10^9.5 Msol (Z/Zsol)^-1, or 6-20% of the expected baryonic mass if the Ne VIII absorbers have solar metallicity. Ionizing Ne VII to Ne VIII requires 207 eV, and photons with this energy are scarce in the CGM. However, for the median halo mass and redshift of our sample, the virial temperature is close to the peak temperature for the Ne VIII ion, and the Ne VIII-bearing gas is plausibly collisionally ionized near this temperature. Moreover, we find that photoionized Ne VIII requires cool and low-density clouds that would be highly underpressured (by approximately two orders of magnitude) relative to the putative, ambient virialized medium, complicating scenarios where such clouds could survive. Thus, more complex (e.g., non-equilibrium) models may be required; this first statistical sample of Ne VIII absorber/galaxy systems will provide stringent constraints for future CGM studies.

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

The ionizing photons produced by massive stars are key actors in galaxy evolution. Ionizing photon production and escape is poorly understood. Improved space-based, spatially-resolved, multiplexed spectroscopic capabilities covering observed wavelengths of 1000 to 3000 Angstrom, in concert with spectroscopy from the ELTs and JWST, would lead to definitive answers as to how ionizing photons are produced and leaked, what populations of galaxies are responsible for ionizing photon leakage, what determines whether escape is possible, and how ionizing galaxy populations evolve over cosmic time.

We present results from a survey of galaxies in the fields of six z>3 Damped Lyman alpha systems (DLAs) using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We report a high detection rate of up to ~80% of galaxies within 1000 km/s from DLAs and with impact parameters between 25 and 280 kpc. In particular, we discovered 5 high-confidence Lyman alpha emitters associated with three DLAs, plus up to 9 additional detections across five of the six fields. The majority of the detections are at relatively large impact parameters (>50 kpc) with two detections being plausible host galaxies. Among our detections, we report four galaxies associated with the most metal-poor DLA in our sample (Z/Z_sun = -2.33), which trace an overdense structure resembling a filament. By comparing our detections with predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological simulations and a semi-analytic model designed to reproduce the observed bias of DLAs at z>2, we conclude that our observations are consistent with a scenario in which a significant fraction of DLAs trace the neutral regions within halos with a characteristic mass of 10^11-10^12 M_sun, in agreement with the inference made from the large-scale clustering of DLAs. We finally show how larger surveys targeting ~25 absorbers have the potential of constraining the characteristic masses of halos hosting high-redshift DLAs with sufficient accuracy to discriminate between different models.

Galaxies are surrounded by extended atmospheres, which are often called the circumgalactic medium (CGM) and are the least understood part of galactic ecosystems. The CGM serves as a reservoir of both diffuse, metal-poor gas accreted from the intergalactic medium, and metal-rich gas that is either ejected from galaxies by energetic feedback or stripped from infalling satellites. As such, the CGM is empirically multi-phased and complex in dynamics. Significant progress has been made in the past decade or so in observing the cosmic-ray/B-field, as well as various phases of the CGM. But basic questions remain to be answered. First, what are the energy, mass, and metal contents of the CGM? More specifically, how are they spatially distributed and partitioned in the different components? Moreover, how are they linked to properties of host galaxies and their global clustering and intergalactic medium environments? Lastly, what are the origin, state, and life-cycle of the CGM? This question explores the dynamics of the CGM. Here we illustrate how these questions may be addressed with multi-wavelength observations of the CGM.

Two groundbreaking new facilities will commence operations early in the 2020s and thereafter define much of the broad landscape of US optical-infrared astronomy in the remaining decade. The Large Synoptic Survey Telescope (LSST), perched atop Cerro Pachon in the Chilean Andes, will revolutionize the young field of Time Domain Astronomy through its wide-field, multi-band optical imaging survey. At the same time, the James Webb Space Telescope (JWST), orbiting at the Sun-Earth L2 Lagrange point, will provide stunningly high-resolution views of selected targets from the red end of the optical spectrum to the mid-infrared. However, the spatial resolution of the LSST observations will be limited by atmospheric seeing, while JWST will be limited in its time-domain capabilities. This paper highlights the scientific opportunities lying between these two landmark missions, i.e., science enabled by systems capable of astronomical observations with both high cadence in the time domain and high resolution in the spatial domain. The opportunities range from constraining the late phases of stellar evolution in nearby resolved populations to constraining dark matter distributions and cosmology using lensed transient sources. We describe a system that can deliver the required capabilities.

The circumgalactic medium (CGM) of galaxies serves as a record of the influences of outflows and accretion that drive the evolution of galaxies. Feedback from star formation drives outflows that carry mass and metals away from galaxies to the CGM, while infall from the intergalactic medium (IGM) is thought to bring in fresh gas to fuel star formation. Such exchanges of matter between IGM-CGM-galaxies have proven critical to producing galaxy scaling relations in cosmological simulations that match observations. However, the nature of these processes, of the physics that drives outflows and accretion, and their evolution with cosmic time are not fully characterized. One approach to constraining these processes is to characterize the metal enrichment of gas around and beyond galaxies. Measurements of the metallicity distribution functions of CGM/IGM gas over cosmic time provide independent tests of cosmological simulations. We have made great progress over the last decade as direct result of a very sensitive, high-resolution space-based UV spectrograph and the rise of ground-based spectroscopic archives. We argue the next transformative leap to track CGM/IGM metals during the epoch of galaxy formation and transformation into quiescent galaxies will require 1) a larger space telescope with an even more sensitive high-resolution spectrograph covering both the far- and near-UV (1,000-3,000 \AA); and 2) ground-based archives housing science-ready data.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

We present simulations from the new "Figuring Out Gas & Galaxies in Enzo" (FOGGIE) project. In contrast to most extant simulations of galaxy formation, which concentrate computational resources on galactic disks and spheroids with fluid and particle elements of fixed mass, the FOGGIE simulations focus on extreme spatial and mass resolution in the circumgalactic medium (CGM) surrounding galaxies. Using the Enzo code and a new refinement scheme, FOGGIE reaches spatial resolutions of 381 comoving $h^{-1}$ pc and resolves extremely low masses ($\lesssim 1$--$100$ Msun out to 100 comoving $h^{-1}$ kpc from the central halo. At these resolutions, cloud and filament-like structures giving rise to simulated absorption are smaller, and better resolved, than the same structures simulated with standard density-dependent refinement. Most of the simulated absorption arises in identifiable and well-resolved structures with masses $\lesssim 10^4$ Msun, well below the mass resolution of typical zoom simulations. However, integrated quantities such as mass surface density and ionic covering fractions change at only the $\lesssim 30$% level as resolution is varied. This relatively small changes in projected quantities---even when the sizes and distribution of absorbing clouds change dramatically---indicate that commonly used observables provide only weak constraints on the physical structure of the underlying gas. Comparing the simulated absorption features to the KODIAQ (Keck Observatory Database of Ionized Absorption toward Quasars) survey of $z \sim2$--$3.5$ Lyman limit systems, we show that high-resolution FOGGIE runs better resolve the internal kinematic structure of detected absorption, and better match the observed distribution of absorber properties. These results indicate that CGM resolution is key in properly testing simulations of galaxy evolution with circumgalactic observations.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

We report the discovery and analysis of only the third Lyman-limit system in which a high-quality resolution, echelle spectrum reveals no metal absorption lines, implying a metallicity $\lesssim$1/10000 solar. Our HIRES spectrum of the background quasar, PSS1723$+$2243, provides a neutral hydrogen column density range for LLS1723 of $N_\textrm{HI}=10^{\text{17.9--18.3}}$ cm$^{-2}$ at redshift $z_\textrm{abs}\approx4.391$. The lower bound on this range, and the lack of detectable absorption from the strongest low-ionisation metal lines, are combined in photoionisation models to infer a robust, conservative upper limit on the metallicity: $\log(Z/Z_\odot)<-4.14$ at 95% confidence. Such a low metallicity raises the question of LLS1723's origin and enrichment history. Previous simulations of the circumgalactic medium imply that LLS1723 is a natural candidate for a cold gas stream accreting towards a galaxy. Alternatively, LLS1723 may represent a high-density portion of the intergalactic medium containing either pristine gas -- unpolluted by stellar debris for 1.4 Gyr after the Big Bang -- or the remnants of low-energy supernovae from (likely low-mass) Population III stars. Evidence for the circumgalactic scenario could be obtained by mapping the environment around LLS1723 with optical integral-field spectroscopy. The intergalactic possibilities highlight the need for -- and opportunity to test -- simulations of the frequency with which such high-density, very low-metallicity systems arise in the intergalactic medium.

We present the results from our COS circumgalactic medium (CGM) compendium (CCC), a survey of the CGM at z<1 using HI-selected absorbers with 15<log N(HI) <19. We focus here on 82 partial Lyman limit systems (pLLSs, 16.2<log N(HI) <17.2) and 29 LLSs (17.2<log N(HI) <19). Using Bayesian techniques and Markov-chain Monte Carlo sampling of a grid of photoionization models, we derive the posterior probability distribution functions (PDFs) for the metallicity of each absorber in CCC. We show that the combined pLLS metallicity PDF at z<1 has two main peaks at [X/H]=-1.7 and -0.4, with a strong dip at [X/H]=-1. The metallicity PDF of the LLSs might be more complicated than an unimodal or bimodal distribution. The pLLSs and LLSs probe a similar range of metallicities -3<[X/H]<+0.4, but the fraction of very metal-poor absorbers with [X/H]<-1.4 is much larger for the pLLSs than the LLSs. In contrast, absorbers with log N(HI)>19 have mostly -1<[X/H]<0 at z<1. The metal-enriched gas probed by pLLSs and LLSs confirms that galaxies that have been enriching their CGM over billions of years. Surprisingly, despite this enrichment, there is also abundant metal-poor CGM gas (41-59% of the pLLSs have [X/H]<-1.4), representing a reservoir of near-pristine gas around z<1 galaxies. We compare our empirical results to recent cosmological zoom simulations, finding some discrepancies, including an overabundance of metal-enriched CGM gas in simulations.

Observing the circumgalactic medium (CGM) in emission provides 3D maps of the spatial and kinematic extent of the gas that fuels galaxies and receives their feedback. We present mock emission-line maps of highly resolved CGM gas from the FOGGIE project (Figuring Out Gas & Galaxies in Enzo) and link these maps back to physical and spatial properties of the gas. By increasing the spatial resolution alone, the total luminosity of the line emission increases by an order of magnitude. This increase arises in the abundance of dense small-scale structure resolved when the CGM gas is simulated to < 100 pc scales. Current integral field unit instruments like KCWI and MUSE should be able to detect the brightest knots and filaments of such emission, and from this to infer the bulk kinematics of the CGM gas with respect to the galaxy. We conclude that accounting for small-scale structure well below the level of instrument spatial resolution is necessary to properly interpret such observations in terms of the underlying gas structure driving observable emission.

Lyman Limit systems (LLSs) trace the low-density circumgalactic medium and the most dense regions of the intergalactic medium, so their number density and evolution at high redshift, just after reionisation, are important to constrain. We present a survey for LLSs at high redshifts, $z_{\rm LLS} =3.5$--5.4, in the homogeneous dataset of 153 optical quasar spectra at $z \sim 5$ from the Giant Gemini GMOS survey. Our analysis includes detailed investigation of survey biases using mock spectra which provide important corrections to the raw measurements. We estimate the incidence of LLSs per unit redshift at $z \approx 4.4$ to be $\ell(z) = 2.6 \pm 0.4$. Combining our results with previous surveys at $z_{\rm LLS} <4$, the best-fit power-law evolution is $\ell(z) = \ell_* [(1+z)/4]^\alpha$ with $\ell_* = 1.46 \pm 0.11$ and $\alpha = 1.70 \pm 0.22$ (68\% confidence intervals). Despite hints in previous $z_{\rm LLS} <4$ results, there is no indication for a deviation from this single power-law soon after reionization. Finally, we integrate our new results with previous surveys of the intergalactic and circumgalactic media to constrain the hydrogen column density distribution function, $f(N_{\rm HI},X)$, over 10 orders of magnitude. The data at $z \sim 5$ are not well described by the $f(N_{\rm HI},X)$ model previously reported for $z \sim 2$--3 (after re-scaling) and a 7-pivot model fitting the full $z \sim 2$--5 dataset is statistically unacceptable. We conclude that there is significant evolution in the shape of $f(N_{\rm HI},X)$ over this $\sim$2 billion year period.

The Giant Southern Stream (GSS) of M31, a keystone signature of a major accretion event, yields crucial constraints on M31 formation and evolution models. Currently, our understanding of the GSS, in terms of both its geometry and its chemistry, results from either wide-field imaging probing only a few magnitudes below the red giant branch tip, or deep imaging or spectroscopy of isolated regions. Here, we take an alternative approach, using Hubble Space Telescope (HST) imaging to characterize the horizontal branch red clump (RC) using unbinned maximum likelihood fits to luminosity functions (LFs) from observed color-magnitude diagrams (CMDs). Comparing the RC mean magnitude across three fields at projected distances of 21, 52 and 80 kpc from M31, we find a line of sight distance gradient identical to recent literature measurements in fields along the core. We also find tentative evidence that the line of sight distance dispersion increases with projected distance from M31. Meanwhile, the metallicity in the 52 kpc field westward of the GSS core is at least as high as that in the 21 kpc GSS core field, and the peak colors of the RC in these two fields imply identical metallicities to within 0.2 dex. We discuss implications for distance and metallicity gradients both along and perpendicular to the GSS in the context of recent ground-based photometric and spectroscopic results, including evidence for a dropoff in metallicity moving westward from the GSS, as well as prospects for further constraining stellar populations in the vicinity of the GSS.

We present a study of cold gas absorption from a damped Lyman-$\alpha$ absorber (DLA) at redshift $z_{\rm abs}=1.946$ towards two lensed images of the quasar J144254.78+405535.5 at redshift $z_{\rm QSO} = 2.590$. The physical separation of the two lines of sight at the absorber redshift is $d_{\rm abs}=0.7$~kpc based on our lens model. We observe absorption lines from neutral carbon and H$_2$ along both lines of sight indicating that cold gas is present on scales larger than $d_{\rm abs}$. We measure column densities of HI to be $\log N(\rm H\,\i) = 20.27\pm0.02$ and $20.34\pm0.05$ and of H$_2$ to be $\log N(\rm H_2) = 19.7\pm0.1$ and $19.9\pm0.2$. The metallicity inferred from sulphur is consistent with Solar metallicity for both sightlines: $[{\rm S/H}]_A = 0.0\pm0.1$ and $[{\rm S/H}]_B = -0.1\pm0.1$. Based on the excitation of low rotational levels of H$_2$, we constrain the temperature of the cold gas phase to be $T=109\pm20$ and $T=89\pm25$ K for the two lines of sight. From the relative excitation of fine-structure levels of CI, we constrain the hydrogen volumetric densities in the range of $40-110$ cm$^{-3}$. Based on the ratio of observed column density and volumetric density, we infer the average individual `cloud' size along the line of sight to be $l\approx0.1$ pc. Using the transverse line-of-sight separation of 0.7 kpc together with the individual cloud size, we are able to put an upper limit to the volume filling factor of cold gas of $f_{\rm vol} < 0.2$ %. Nonetheless, the projected covering fraction of cold gas must be large (close to unity) over scales of a few kpc in order to explain the presence of cold gas in both lines of sight. Compared to the typical extent of DLAs (~10-30 kpc), this is consistent with the relative incidence rate of CI absorbers and DLAs.

We decompose the Lyman-{\alpha} (Ly{\alpha}) forest of an extensive sample of 74 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in the peculiar velocity field have the smallest Doppler parameters, resulting in a low-$b$ cutoff in the $b$-$N_{\text{HI}}$ set by the thermal state of intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range $2.0\leq z \leq 3.4$, which allows us to measure the evolution of the IGM temperature-density ($T= T_0 (\rho/ \rho_0)^{\gamma-1}$) relation parameters $T_0$ and $\gamma$. We calibrate our measurements against Ly$\alpha$ forest simulations, using 21 different thermal models of the IGM at each redshift, also allowing for different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of $T_0$ shows a suggestive peak at $z=2.8$, while our evolution of $\gamma$ is consistent with $\gamma\simeq 1.4$ and disfavors inverted temperature-density relations. Our measured evolution of $T_0$ and $\gamma$ are generally in good agreement with previous determinations in the literature. Both the peak in the evolution of $T_0$ at $z = 2.8$, as well as the high temperatures $T_0\simeq 15000-20000\,$K that we observe at $2.4 < z < 3.4$, strongly suggest that a significant episode of heating occurred after the end of HI reionization, which was most likely the cosmic reionization of HeII.

We present a neutral hydrogen-selected absorption-line survey of gas with HI column densities 15<log N(HI)<19 at z<1 using the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our main aim is to determine the metallicity distribution of these absorbers. Our sample consists of 224 absorbers selected on the basis of their HI absorption strength. Here we discuss the properties of our survey and the immediate empirical results. We find singly and doubly ionized metal species and HI typically have similar velocity profiles, implying they probe gas in the same or similar environments. The column density ionic ratios (e.g., CII/CIII, OI/CII) indicate the gas in these absorbers is largely ionized, and the ionization conditions are quite comparable across the sampled N(HI) range. The Doppler parameters of the HI imply T<50,000 K on average, consistent with the gas being photoionized. The MgII column densities span >2 orders of magnitude at any given N(HI), indicating a wide range of metallicities (from solar to <1/100 solar). In the range 16.2<log N(HI)<17, there is a gap in the N(MgII) distribution corresponding to gas with ~10% solar metallicity, consistent with the gap seen in the previously identified bimodal metallicity distribution in this column density regime. Less than 3% of the absorbers in our sample show no detectable metal absorption, implying truly-pristine gas at z<1 is uncommon. We find [FeII/MgII] = -0.4+/-0.3, and since alpha-enhancement can affect this ratio, dust depletion is extremely mild.

We present the discovery of a gravitationally lensed dust-reddened QSO at $z=2.517$ discovered in a survey for red QSOs by infrared selection. $Hubble~Space~Telescope$ imaging in the WFC3/IR F160W and F125W filters reveals a quadruply lensed system in a cusp configuration. We find that compared to the central image of the cusp, the nearby, brightest image is anomalous by a factor of $\sim7-11$. Although the source is extremely bright in the mid-infrared, a magnification by a factor of $\sim50-120$ places an upper limit of 1.35 mJy on the intrinsic mid-infrared brightness, well below the $WISE~W4$ detection limit of 6 mJy. We find that this QSO is moderately reddened, with $E(B-V)=0.7$ and that $\sim1\%$ of the intrinsic spectrum is leaked back into the line of sight resulting in an upturn in its UV spectrum. We conclude that the QSO's reddening is intrinsic and not due to the lens. Consistent with previous red quasar samples, this source exhibits outflows in its spectrum as well as morphological properties suggestive of it being in a merger-driven transitional phase. Depending on how $L_{\rm bol}$ is computed, the quasar's accretion rate may be as high as $0.26~L_{\rm Edd}$. We detect two Lyman limit systems, at $z=2.102$ and $z=2.431$, with absorption by metal lines likely at small impact parameter to the QSO, and a putative lens redshift of $z=0.599$. Given the rarity of quad lenses, the discovery of this source allows detailed study of a less luminous, more typical infrared-selected quasar at high redshift.

We present the first paper of the series Origin of Metals around Galaxies (OMG) aimed to study the origin of the metals observed in the circumgalactic and intergalactic media. In this work we extract and build the catalogs of metal absorbers that will be used in future analyses, and make our results publicly available to the community. We design a fully automatic algorithm to search for absorption metal-line doublets of the species CIV, NV, SiIV and MgII in high-resolution ($R\gtrsim30\,000$) quasar spectra without human intervention, and apply it to the high-resolution and signal-to-noise ratio spectra of 690 quasars, observed with the UVES and HIRES instruments. We obtain $5\,656$ CIV doublets, $7\,919$ doublets of MgII, $2\,258$ of SiIV, and 239 of NV, constituting the largest high-resolution metal-doublet samples to date, and estimate the dependence of their completeness and purity on various doublet parameters such as equivalent width and redshift, using real and artificial quasar spectra. The catalogs include doublets with rest-frame line equivalent widths down to a few ${\rm m\AA}$, all detected at a significance above 3$\sigma$, and covering the redshifts between $1<z \lesssim 5$, properties that make them useful for a wide range of chemical evolution studies.

We present the first results from our survey of intervening and proximate Lyman limit systems (LLSs) at $z$$\sim$2.0-2.5 using the Wide Field Camera 3 on-board the Hubble Space Telescope. The quasars in our sample are projected pairs with proper transverse separations $R_\perp$$\leq$150 kpc and line of sight velocity separations $\lesssim$11,000 km/s. We construct a stacked ultraviolet (rest-frame wavelengths 700-2000\AA) spectrum of pairs corrected for the intervening Lyman forest and Lyman continuum absorption. The observed spectral composite presents a moderate flux excess for the most prominent broad emission lines, a $\sim$30% decrease in flux at $\lambda$=800-900\AA\ compared to a stack of brighter quasars not in pairs at similar redshifts, and lower values of the mean free path of the HI ionizing radiation for pairs ($\lambda_{\rm mfp}^{912}=140.7\pm20.2~h_{70}^{-1}$Mpc) compared to single quasars ($\lambda_{\rm mfp}^{912}=213.8\pm28~h_{70}^{-1}$Mpc) at the average redshift $z\simeq2.44$. From the modelling of LLS absorption in these pairs, we find a higher ($\sim$20%) incidence of proximate LLSs with $\log N_{\rm HI}\geq17.2$ at $\delta v$$<$5,000 km/s compared to single quasars ($\sim$6%). These two rates are different at the 5$\sigma$ level. Moreover, we find that optically-thick absorbers are equally shared between foreground and background quasars. Based on these pieces of evidence, we conclude that there is a moderate excess of gas absorbing Lyman continuum photons in our closely-projected quasar pairs compared to single quasars. We argue that this gas arises mostly within large-scale structures or partially neutral regions inside the dark matter haloes where these close pairs reside.

We investigate the limitations of statistical absorption measurements with the SDSS optical spectroscopic surveys. We show that changes in the data reduction strategy throughout different data releases have led to a better accuracy at long wavelengths, in particular for sky line subtraction, but a degradation at short wavelengths with the emergence of systematic spectral features with an amplitude of about one percent. We show that these features originate from inaccuracy in the fitting of modeled F-star spectra used for flux calibration. The best-fit models for those stars are found to systematically over-estimate the strength of metal lines and under-estimate that of Lithium. We also identify the existence of artifacts due to masking and interpolation procedures at the wavelengths of the hydrogen Balmer series leading to the existence of artificial Balmer $\alpha$ absorption in all SDSS optical spectra. All these effects occur in the rest-frame of the standard stars and therefore present Galactic longitude variations due to the rotation of the Galaxy. We demonstrate that the detection of certain weak absorption lines reported in the literature are solely due to calibration effects. Finally, we discuss new strategies to mitigate these issues.

We present a photometric detection of the first brightness dips of the unique variable star KIC 8462852 since the end of the Kepler space mission in 2013 May. Our regular photometric surveillance started in October 2015, and a sequence of dipping began in 2017 May continuing on through the end of 2017, when the star was no longer visible from Earth. We distinguish four main 1-2.5% dips, named "Elsie," "Celeste," "Skara Brae," and "Angkor", which persist on timescales from several days to weeks. Our main results so far are: (i) there are no apparent changes of the stellar spectrum or polarization during the dips; (ii) the multiband photometry of the dips shows differential reddening favoring non-grey extinction. Therefore, our data are inconsistent with dip models that invoke optically thick material, but rather they are in-line with predictions for an occulter consisting primarily of ordinary dust, where much of the material must be optically thin with a size scale <<1um, and may also be consistent with models invoking variations intrinsic to the stellar photosphere. Notably, our data do not place constraints on the color of the longer-term "secular" dimming, which may be caused by independent processes, or probe different regimes of a single process.

We present a new measurement of the Ly${\alpha}$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise, and corrects for the bias induced by masking missing data, DLAs, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes $k > 0.02\,\mathrm{s\,km^{-1}}$, unaccessible to previous SDSS/BOSS analyses. It is well known that these high-$k$ modes are highly sensitive to the thermal state of the intergalactic medium, however contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power, and find a modest effect on modes with $k < 0.1\,\mathrm{s\,km^{-1}} $. As a result, by masking metals and restricting to $k < 0.1\,\mathrm{s\,km^{-1}}$ their impact is completely mitigated. We present an end-to-end Bayesian forward modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Ly${\alpha}$ forest simulations. These mocks are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform MCMC fits. Our results agree well with BOSS on scales $k < 0.02\,\mathrm{s\,km^{-1}}$ where the measurements overlap. The combination of BOSS' percent level low-$k$ precision with our $5-15\%$ high-$k$ measurements, results in a powerful new dataset for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables.

The age and evolutionary status of MWC349A, the unique emission line star with maser and laser radiation in hydrogen recombination lines, remain unknown because the spectrum of the star is veiled by bright emission from the ionized disk and wind. The major argument for this massive ($>$10 M_sun) star being evolved is its association with a close-by (2.4 arcsec) companion, MWC349B, whose B0III spectrum implies an age of a few Myrs. However, newly obtained high-resolution spectra of MWC349B reveal a difference $\approx$ 35 km/s in the radial velocities of the two stars, which makes their being gravitationally bound highly improbable. An estimate of the relative proper motion of the two stars seems to confirm this conclusion. This reopens the previously suggested possibility that MWC349A is a young massive star in a region of active star formation close to the Cyg OB2 association. MWC349B, which moves with a speed $\geq$ 35 km/s relative to Cyg OB2, may be a runaway star from this association.

The Large Ultraviolet / Optical / Infrared Surveyor (LUVOIR) is one of four large mission concepts currently undergoing community study for consideration by the 2020 Astronomy and Astrophysics Decadal Survey. The LUVOIR Ultraviolet Multi-Object Spectrograph, LUMOS, is being designed to support all of the UV science requirements of LUVOIR, from exoplanet host star characterization to tomography of circumgalactic halos to water plumes on outer solar system satellites. LUMOS offers point source and multi-object spectroscopy across the UV bandpass, with multiple resolution modes to support different science goals. The instrument will provide low (R = 8,000-18,000) and medium (R = 30,000-65,000) resolution modes across the far-ultraviolet (FUV: 100-200 nm) and near-ultraviolet (NUV: 200-400 nm) windows, and a very low resolution mode (R = 500) for spectroscopic investigations of extremely faint objects in the FUV. Imaging spectroscopy will be accomplished over a 3 x 1.6 arcminute field-of-view by employing holographically-ruled diffraction gratings to control optical aberrations, microshutter arrays (MSA), advanced optical coatings for high-throughput in the FUV, and next generation large-format photon-counting detectors. The spectroscopic capabilities of LUMOS are augmented by an FUV imaging channel (100-200nm, 13 milliarcsecond angular resolution, 2 x 2 arcminute field-of-view) that will employ a complement of narrow and medium-band filters. We present an overview of LUMOS' observing modes and estimated performance curves for effective area, spectral resolution, and imaging performance. Example "LUMOS 100-hour Highlights" observing programs are presented to demonstrate the potential power of LUVOIR's ultraviolet spectroscopic capabilities.

Flux estimates for ionizing radiation escaping from star-forming galaxies with characteristic UV luminosities ($L^{*}_{1500(1+z)}$) derived from GALEX and the VIMOS-VLT Deep Survey, are presented as a function of redshift and assumed escape fraction. These estimates offer guidance to the design of instrumentation and observing strategies, be they spectroscopic or photometric, attempting to detect LyC escaping star-forming galaxies for redshifts $ z < 3$. Examples are given that relate the absolution escape fraction ($f^{e}_{LyC}$) of ionizing photons, integrated over the entire extreme UV (EUV) bandpass, to the relative escape fraction ($f^{e}_{900}$) observed just shortward of the ionization edge at 911.8 \AA\ as a function of HI, HeI, and HeII column densities. We find that for $\log{N_{HI}(cm^{-2})} \gtrsim 17.0$ $f^{e}_{LyC}$ is significantly greater than $f^{e}_{900}$. Detection of LyC "drop-ins" in the rest-frame EUV will provide enhanced fidelity to determinations of the integrated fraction of ionizing photons $f^{e}_{LyC}$ that escape star-forming galaxies and contribute to the meta-galactic ionizing background.

We present a deep search for HI 21-cm emission from the gaseous halo of Messier 31 as part of Project AMIGA, a large program Hubble Space Telescope program to study the circumgalactic medium of the Andromeda galaxy. Our observations with the Robert C. Byrd Green Bank Telesope target sight lines to 48 background AGNs, more than half of which have been observed in the ultraviolet with the Cosmic Origins Spectrograph, with impact parameters $25 \lesssim \rho \lesssim 330$ kpc ($0.1 \lesssim \rho / R_{\rm vir} \lesssim 1.1$). We do not detect any 21-cm emission toward these AGNs to limits of $N({\rm HI}) \approx 4 \times10^{17}$ cm$^{-2}$ ($5\sigma$, per 2 kpc diameter beam). This column density corresponds to an optical depth of $\sim2.5$ at the Lyman limit, thus our observations overlap with absorption line studies of Lyman limit systems at higher redshift. Our non-detections place a limit on the covering factor of such optically-thick gas around M31 to $f_c < 0.051$ (at 90\% confidence) for $\rho \leq R_{\rm vir}$. While individual clouds have previously been found in the region between M31 and M33, the covering factor of strongly optically-thick gas is quite small. Our upper limits on the covering factor are consistent with expectations from recent cosmological "zoom" simulations. Recent COS-Halos ultraviolet measurements of \HI\ absorption about an ensemble of galaxies at $z \approx 0.2$ show significantly higher covering factors within $\rho \lesssim 0.5 R_{\rm vir}$ at the same $N({\rm H I})$, although the metal ion-to-H I ratios appear to be consistent with those seen in M31.

We present Keck Cosmic Web Imager spectroscopy of the four putative images of the lensed quasar candidate J014709+463037 recently discovered by Berghea et al. (2017). The data verify the source as a quadruply lensed, broad absorption-line quasar having z_S = 2.377 +/- 0.007. We detect intervening absorption in the FeII 2586, 2600, MgII 2796, 2803, and/or CIV 1548, 1550 transitions in eight foreground systems, three of which have redshifts consistent with the photometric-redshift estimate reported for the lensing galaxy (z_L ~ 0.57). By virtue of their positions on the sky, the source images probe these absorbers over transverse physical scales of ~0.3-21 kpc, permitting assessment of the variation in metal-line equivalent width W_r as a function of sight-line separation. We measure differences in W_r,2796 of <40% across all sight-line pairs subtending 7-21 kpc, suggestive of a high degree of spatial coherence for MgII-absorbing material. W_r,2600 is observed to vary by >50% over the same scales across the majority of sight-line pairs, while CIV absorption exhibits a wide range in W_r,1548 differences of ~5-80% within transverse distances less than ~3 kpc. J014709+463037 is one of only a handful of z > 2 quadruply lensed systems for which all four source images are very bright (r = 15.4-17.7 mag) and are easily separated in ground-based seeing conditions. As such, it is an ideal candidate for higher-resolution spectroscopy probing the spatial variation in the kinematic structure and physical state of intervening absorbers.

We present and make publicly available the second data release (DR2) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR2 consists of a fully-reduced sample of 300 quasars at 0.07 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000). DR2 contains 831 spectra available in continuum normalized form, representing a sum total exposure time of ~4.9 megaseconds on source. These co-added spectra arise from a total of 1577 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed. DR2 extends DR1 by adding 130 new quasars to the sample, including additional observations of QSOs in DR1. All new data in DR2 were obtained with the single-chip Tektronix TK2048 CCD configuration of HIRES in operation between 1995 and 2004. DR2 is publicly available to the community, housed as a higher level science product at the KOA and in the igmspec database (v03).

We present new observations acquired with the Multi Unit Spectroscopic Explorer instrument on the Very Large Telescope in a quasar field that hosts a high column-density damped Ly{\alpha} absorber (DLA) at z~3.25. We detect Ly{\alpha} emission from a nebula at the redshift of the DLA with line luminosity (27+/-1)x1e41 erg/s, which extends over 37+/-1 kpc above a surface brightness limit of 6x1e-19 erg/s/cm2/arcsec2 at a projected distance of 30.5+/-0.5 kpc from the quasar sightline. Two clumps lie inside this nebula, both with Ly{\alpha} rest-frame equivalent width > 50 A and with relative line-of-sight velocities aligned with two main absorption components seen in the DLA spectrum. In addition, we identify a compact galaxy at a projected distance of 19.1+/-0.5 kpc from the quasar sightline. The galaxy spectrum is noisy but consistent with that of a star-forming galaxy at the DLA redshift. We argue that the Ly{\alpha} nebula is ionized by radiation from star formation inside the two clumps, or by radiation from the compact galaxy. In either case, these data imply the presence of a structure with size >>50 kpc inside which galaxies are assembling, a picture consistent with galaxy formation in groups and filaments as predicted by cosmological simulations such as the EAGLE simulations.

We measured metallicities for 33 z=3.4-4.2 absorption line systems drawn from a sample of H I-selected-Lyman limit systems (LLSs) identified in Sloan Digital Sky Survey (SDSS) quasar spectra and stratified based on metal line features. We obtained higher-resolution spectra with the Keck Echellette Spectrograph and Imager, selecting targets according to our stratification scheme in an effort to fully sample the LLS population metallicity distribution. We established a plausible range of H I column densities and measured column densities (or limits) for ions of carbon, silicon, and aluminum, finding ionization-corrected metallicities or upper limits. Interestingly, our ionization models were better constrained with enhanced $\alpha$-to-aluminum abundances, with a median abundance ratio of [$\alpha$/Al]=0.3. Measured metallicities were generally low, ranging from [M/H]=-3 to -1.68, with even lower metallicities likely for some systems with upper limits. Using survival statistics to incorporate limits, we constructed the cumulative distribution function (CDF) for LLS metallicities. Recent models of galaxy evolution propose that galaxies replenish their gas from the low-metallicity intergalactic medium (IGM) via high-density H I "flows" and eject enriched interstellar gas via outflows. Thus, there has been some expectation that LLSs at the peak of cosmic star formation ($z\approx3$) might have a bimodal metallicity distribution. We modeled our CDF as a mix of two Gaussian distributions, one reflecting the metallicity of the IGM and the other representative of the interstellar medium of star-forming galaxies. This bimodal distribution yielded a poor fit. A single Gaussian distribution better represented the sample with a low mean metallicity of [M/H] $\approx -2.5$.

We present the first results from our KODIAQ Z survey aimed to determine the metallicity distribution and physical properties of the z>2 partial and full Lyman limit systems (pLLSs and LLSs; 16.2<log N(HI)<19), which are probed of the interface regions between the intergalactic medium (IGM) and galaxies. We study 31 HI-selected pLLSs and LLSs at 2.3<z<3.3 observed with Keck/HIRES in absorption against background QSOs. We compare the column densities of metal-ions to HI and use photoionization models to assess the metallicity. The metallicity distribution of the pLLSs/LLSs at 2.3<z<3.3 is consistent with a unimodal distribution peaking at [X/H]=-2. The metallicity distribution of these absorbers therefore evolves markedly with z since at z<1 it is bimodal with peaks at [X/H]=-1.8 and -0.3. There is a substantial fraction (25-41%) of pLLSs/LLSs with metallicities well below those of damped Lya absorbers (DLAs) at any studied z from z<1 to z~2-4, implying reservoirs of metal-poor cool, dense gas in the IGM/galaxy interface at all z. However, the gas probed by pLLSs and LLSs is rarely pristine, with a fraction 3-18% for pLLSs/LLSs with [X/H]<-3. We find C/alpha enhancement in several pLLSs and LLSs in the metallicity range -2<[X/H]<-0.5, where C/alpha is 2-5 times larger than observed in Galactic metal-poor stars or high redshift DLAs at similar metallicities. This is likely caused by preferential ejection of carbon from metal-poor galaxies into their surroundings.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

Using Hubble Space Telescope Cosmic Origins Spectrograph observations of 89 QSO sightlines through the Sloan Digital Sky Survey footprint, we study the relationships between C IV absorption systems and the properties of nearby galaxies as well as large-scale environment. To maintain sensitivity to very faint galaxies, we restrict our sample to 0.0015 < z < 0.015, which defines a complete galaxy survey to L > 0.01 L* or stellar mass log M_* > 8 Msun. We report two principal findings. First, for galaxies with impact parameter rho < 1 rvir, C IV detection strongly depends on the luminosity/stellar mass of the nearby galaxy. C IV is preferentially associated with galaxies with log M_* > 9.5 Msun; lower mass galaxies rarely exhibit significant C IV absorption (covering fraction f = 9 +12-6% for 11 galaxies with log M_* < 9.5 Msun). Second, C IV detection within the log M_* > 9.5 Msun population depends on environment. Using a fixed-aperture environmental density metric for galaxies with rho < 160 kpc at z < 0.055, we find that 57+/-12% (8/14) of galaxies in low-density regions (regions with fewer than seven L > 0.15 L* galaxies within 1.5 Mpc) have affiliated C IV absorption; however, none (0/7) of the galaxies in denser regions show C IV. Similarly, the C IV detection rate is lower for galaxies residing in groups with dark-matter halo masses of log Mhalo > 12.5 Msun. In contrast to C IV, H I is pervasive in the CGM without regard to mass or environment. These results indicate that C IV absorbers with log N(C IV) > 13.5 cm^-2 trace the halos of log M_* > 9.5 Msun galaxies but also reflect larger scale environmental conditions.

We present an analysis of the metallicity distribution of the dense circumgalactic medium (CGM) of galaxies at 0.1 < z < 1.1 as probed by partial Lyman limit systems (pLLSs, 16.1 < log N(H I) < 17.2) and LLSs (17.2 < log N(H I) < 17.7 in our sample). The new H I-selected sample, drawn from our HST COS G140L snapshot survey of 61 QSOs, has 20 pLLSs and 10 LLSs. Combined with our previous survey, we have a total of 44 pLLSs and 11 LLSs. We find that the metallicity distribution of the pLLSs is bimodal at z < 1, with a minimum at [X/H] = -1. The low-metallicity peak comprises (57 +/- 8)% of the pLLSs and is centered at [X/H] ~ -1.87 (1.3% solar metallicity), while the high-metallicity peak is centered at [X/H] ~ -0.32 (48% solar metallicity). Although the sample of LLSs is still small, there is some evidence that the metallicity distributions of the LLSs and pLLSs are different, with a far lower fraction of very metal-poor ([X/H] < -1.4) LLSs than pLLSs. The fraction of LLSs with [X/H] < -1 is similar to that found in pLLSs (~56%). However, higher H I column density absorbers (log N(H I) > 19.0) show a much lower fraction of metal-poor gas; therefore, the metallicity distribution of gas in and around galaxies depends sensitively on N(H I) at z < 1. We interpret the high-metallicity ([X/H] > -1) pLLSs and LLSs as arising in outflows, recycling winds, and tidally-stripped gas around galaxies. The low-metallicity pLLSs and LLSs imply that the CGM of z < 1 galaxies is also host to a substantial mass of cool, dense, low-metallicity gas that may ultimately accrete onto the galaxies.

We report on the search for galaxies in the proximity of two very metal-poor gas clouds at z~3 towards the quasar Q0956+122. With a 5-hour MUSE integration in a ~500x500 kpc^2 region centred at the quasar position, we achieve a >80% complete spectroscopic survey of continuum-detected galaxies with m<25 mag and Ly{\alpha} emitters with luminosity L>3e41 erg/s. We do not identify galaxies at the redshift of a z~3.2 Lyman limit system (LLS) with log Z/Zsun = -3.35 +/- 0.05, placing this gas cloud in the intergalactic medium or circumgalactic medium of a galaxy below our sensitivity limits. Conversely, we detect five Ly{\alpha} emitters at the redshift of a pristine z~3.1 LLS with log Z/Zsun < -3.8, while ~0.4 sources were expected given the z~3 Ly{\alpha} luminosity function. Both this high detection rate and the fact that at least three emitters appear aligned in projection with the LLS suggest that this pristine cloud is tracing a gas filament that is feeding one or multiple galaxies. Our observations uncover two different environments for metal-poor LLSs, implying a complex link between these absorbers and galaxy halos, which ongoing MUSE surveys will soon explore in detail. Moreover, in agreement with recent MUSE observations, we detected a ~90 kpc Ly{\alpha} nebula at the quasar redshift and three Ly{\alpha} emitters reminiscent of a "dark galaxy" population.

We study the properties of HI gas in the outer regions (~2r_25) of a spiral galaxy, UGC 00439 (z = 0.01769), using HI 21-cm absorption towards different components of an extended background radio source, J0041$-$0043 (z = 1.679). The radio source exhibits a compact core coincident with the optical quasar and two lobes separated by ~7 kpc, all at an impact parameter ~25 kpc. The HI 21-cm absorption detected towards the southern lobe is found to extend over ~2 kpc^2. The absorbing gas shows sub-kpc-scale structures with the line-of-sight velocities dominated by turbulent motions. Much larger optical depth variations over 4-7 kpc-scale are revealed by the non-detection of HI 21-cm absorption towards the radio core and the northern lobe, and the detection of NaI and CaII absorption towards the quasar. This could reflect a patchy distribution of cold gas in the extended HI disc. We also detect HI 21-cm emission from UGC 00439 and two other galaxies within ~150 kpc to it, that probably form an interacting group. However, no HI 21-cm emission from the absorbing gas is detected. Assuming a linear extent of ~4 kpc, as required to cover both the core and the southern lobe, we constrain the spin temperature <~300 K for the absorbing gas. The kinematics of the gas and the lack of signatures of any ongoing in situ star formation are consistent with the absorbing gas being at the kinematical minor axis and corotating with the galaxy. Deeper HI 21-cm observations would help to map in greater detail both the large- and small-scale structures in the HI gas associated with UGC 00439.

The first stars, known as Population III (PopIII), produced the first heavy elements, thereby enriching their surrounding pristine gas. Previous detections of metals in intergalactic gas clouds, however, find a heavy element enrichment larger than $1/1000$ times that of the solar environment, higher than expected for PopIII remnants. In this letter we report the discovery of a Lyman limit system (LLS) at $z=3.53$ with the lowest metallicity seen in gas with discernable metals, $10^{-3.41\pm0.26}$ times the solar value, at a level expected for PopIII remnants. We make the first relative abundance measurement in such low metallicity gas: the carbon-to-silicon ratio is $10^{-0.26\pm0.17}$ times the solar value. This is consistent with models of gas enrichment by a PopIII star formation event early in the Universe, but also consistent with later, Population II enrichment. The metals in all three components comprising the LLS, which has a velocity width of 400 km s$^{-1}$, are offset in velocity by $\sim+6$ km s$^{-1}$ from the bulk of the hydrogen, suggesting the LLS was enriched by a single event. Relative abundance measurements in this near-pristine regime open a new avenue for testing models of early gas enrichment and metal mixing.

We study the physical properties of a homogeneous sample of 157 optically-thick absorption line systems at redshifts ~1.8-4.4, selected from a high-dispersion spectroscopic survey of Lyman limit systems (LLSs). By means of multiple ionisation models and Bayesian techniques, we derive the posterior probability distribution functions for the density, metallicity, temperature, and dust content of the absorbing gas. We find that z>2 LLSs are highly ionised with ionisation parameters between -3<log U<-2, depending on the HI column density. LLSs are characterised by low temperatures (T<5x10^4 K) and reside in dust-poor environments. Between z~2.5-3.5, ~80% of the LLSs have physical densities between n(H)~10^-3.5-10^-2 cm^-3 for the assumed UV background, but we caution that a degeneracy between the ionisation parameter and the intensity of the radiation field prevents robust inference on the density and sizes of LLSs. Conversely, metallicity estimates are less sensitive to the assumptions behind ionisation corrections. LLSs at z>2 are characterised by a broad unimodal distribution over >4 orders of magnitude, with a peak at log Z/Zsun~-2. LLSs are metal poor, significantly less enriched than DLAs, with ~70% of the metallicity PDF below log Z/Zsun<-1.5. The median metallicity of super LLSs with log N(HI)>19 rapidly evolves with redshift, with a ten-fold increase between z~2.1-3.6 (~1.5 Gyr). Based on this sample, we find that LLSs at z=2.5-3.5 account for ~15% of all the metals produced by UV-selected galaxies. The implications for theories of cold gas accretion and metal ejection from galaxies are also discussed.

The formation and evolution of galaxies require large reservoirs of cold, neutral gas. The damped Lya systems (DLAs), seen in absorption towards distant quasars and gamma-ray bursts, are predicted to be the dominant reservoirs for this gas. Detailed properties of DLAs have been studied extensively for decades with great success. However, their size, fundamental in understanding their nature, has remained elusive, as quasar and gamma-ray-burst sightlines only probe comparatively tiny areas of the foreground DLAs. Here, we introduce a new approach to measure the full extent of DLAs in the sightlines toward extended background sources. We present the discovery of a high-column-density (log N(HI) = 21.1 +/-0.4 cm^-2) DLA at z ~ 2.4 covering 90-100% of the luminous extent of a line-of-sight background galaxy. Estimates of the size of the background galaxy range from a minimum of a few kpc^2, to ~100 kpc^2, and demonstrate that high-column density neutral gas can span continuous areas 10^8 - 10^10 times larger than previously explored in quasar or gamma-ray burst sightlines. The DLA presented here is the first from a sample of DLAs in our pilot survey that searches Lyman break and Lyman continuum galaxies at high redshift. The low luminosities, large sizes, and mass contents (>~10^6 - 10^9 M_solar) implied by this DLA and the early data suggest that DLAs contain the necessary fuel for galaxies, with many systems consistent with relatively massive, low-luminosity primeval galaxies.

To investigate the evolution of metal-enriched gas over recent cosmic epochs as well as to characterize the diffuse, ionized, metal-enriched circumgalactic medium (CGM), we have conducted a blind survey for C IV absorption systems in 89 QSO sightlines observed with the Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS). We have identified 42 absorbers at z < 0.16, comprising the largest uniform blind sample size to date in this redshift range. Our measurements indicate an increasing C IV absorber number density per comoving path length (dN/dX = 7.5 +/- 1.1) and modestly increasing mass density relative to the critical density of the Universe (Omega(C IV) = 10.0 +/- 1.5 x 10^-8 ) from z ~ 1.5 to the present epoch, consistent with predictions from cosmological hydrodynamical simulations. Furthermore, the data support a functional form for the column density distribution function that deviates from a single power-law, also consistent with independent theoretical predictions. As the data also probe heavy element ions in addition to C IV at the same redshifts, we identify, measure, and search for correlations between column densities of these species where components appear aligned in velocity. Among these ion-ion correlations, we find evidence for tight correlations between C II and Si II, C II and Si III, and C IV and Si IV, suggesting that these pairs of species arise in similar ionization conditions. However, the evidence for correlations decreases as the difference in ionization potential increases. Finally, when controlling for observational bias, we find only marginal evidence for a correlation (86.8% likelihood) between the Doppler line width b(C IV) and column density N(C IV).

We present an absorption-line survey of optically thick gas clouds -- Lyman Limit Systems (LLSs) -- observed at high dispersion with spectrometers on the Keck and Magellan telescopes. We measure column densities of neutral hydrogen NHI and associated metal-line transitions for 157 LLSs at z=1.76-4.39 restricted to 10^17.3 < NHI < 10^20.3. An empirical analysis of ionic ratios indicates an increasing ionization state of the gas with decreasing NHI and that the majority of LLSs are highly ionized, confirming previous expectations. The Si^+/H^0 ratio spans nearly four orders-of-magnitude, implying a large dispersion in the gas metallicity. Fewer than 5% of these LLSs have no positive detection of a metal transition; by z~3, nearly all gas that is dense enough to exhibit a very high Lyman limit opacity has previously been polluted by heavy elements. We add new measurements to the small subset of LLS (~5-10) that may have super-solar abundances. High Si^+/Fe^+ ratios suggest an alpha-enhanced medium whereas the Si^+/C^+ ratios do not exhibit the super-solar enhancement inferred previously for the Lya forest.

We present the largest homogeneous survey of $z>4.4$ damped Lyman-$\alpha$ systems (DLAs) using the spectra of 163 QSOs that comprise the Giant Gemini GMOS (GGG) survey. With this survey we make the most precise high-redshift measurement of the cosmological mass density of neutral hydrogen, $\Omega_{\rm HI}$. At such high redshift important systematic uncertainties in the identification of DLAs are produced by strong intergalactic medium absorption and QSO continuum placement. These can cause spurious DLA detections, result in real DLAs being missed, or bias the inferred DLA column density distribution. We correct for these effects using a combination of mock and higher-resolution spectra, and show that for the GGG DLA sample the uncertainties introduced are smaller than the statistical errors on $\Omega_{\rm HI}$. We find $\Omega_{\rm HI}=0.98^{+0.20}_{-0.18}\times10^{-3}$ at $\langle z\rangle=4.9$, assuming a 20% contribution from lower column density systems below the DLA threshold. By comparing to literature measurements at lower redshifts, we show that $\Omega_{\rm HI}$ can be described by the functional form $\Omega_{\rm HI}(z)\propto(1+z)^{0.4}$. This gradual decrease from $z=5$ to $0$ is consistent with the bulk of HI gas being a transitory phase fuelling star formation, which is continually replenished by more highly-ionized gas from the intergalactic medium, and from recycled galactic winds.

Cold accretion is a primary growth mechanism of simulated galaxies, yet observational evidence of "cold flows" at redshifts where they should be most efficient ($z=2$-4) is scarce. In simulations, cold streams manifest as Lyman-limit absorption systems (LLSs) with low heavy-element abundances similar to those of the diffuse IGM. Here we report on an abundance survey of 17 H I-selected LLSs at $z=3.2$-4.4 which exhibit no metal absorption in SDSS spectra. Using medium-resolution spectra obtained at Magellan, we derive ionization-corrected metallicities (or limits) with a Markov-Chain Monte Carlo sampling that accounts for the large uncertainty in $N_{\rm HI}$ measurements typical of LLSs. The metal-poor LLS sample overlaps with the IGM in metallicity and is best described by a model where $71^{+13}_{-11}\%$ are drawn from the IGM chemical abundance distribution. These represent roughly half of all LLSs at these redshifts, suggesting that 28-40$\%$ of the general LLS population at $z\sim3.7$ could trace unprocessed gas. An ancillary sample of ten LLSs without any a priori metal-line selection is best fit with $48^{+14}_{-12}\%$ of metallicities drawn from the IGM. We compare these results with regions of a moving-mesh simulation; the simulation finds only half as many baryons in IGM-metallicity LLSs, and most of these lie beyond the virial radius of the nearest galaxy halo. A statistically significant fraction of all LLSs have low metallicity and therefore represent candidates for accreting gas; large-volume simulations can establish what fraction of these candidates actually lie near galaxies and the observational prospects for detecting the presumed hosts in emission.

We present and make publicly available the first data release (DR1) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. The KODIAQ survey is aimed at studying galactic and circumgalactic gas in absorption at high-redshift, with a focus on highly-ionized gas traced by OVI, using the HIRES spectrograph on the Keck-I telescope. KODIAQ DR1 consists of a fully-reduced sample of 170 quasars at 0.29 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000) between 2004 and 2012. DR1 contains 247 spectra available in continuum normalized form, representing a sum total exposure time of ~1.6 megaseconds. These co-added spectra arise from a total of 567 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a HIRES data reduction package made available through the XIDL distribution. DR1 is publicly available to the community, housed as a higher level science product at the KOA. We will provide future data releases that make further QSOs, including those with pre-2004 observations taken with the previous-generation HIRES detectors.

The ionising continuum from active galactic nuclei (AGN) is fundamental for interpreting their broad emission lines and understanding their impact on the surrounding gas. Furthermore, it provides hints on how matter accretes onto supermassive black holes. Using HST's Wide Field Camera 3 we have constructed the first stacked ultraviolet (rest-frame wavelengths 600-2500\AA) spectrum of 53 luminous quasars at z=2.4, with a state-of-the-art correction for the intervening Lyman forest and Lyman continuum absorption. The continuum slope ($f_\nu \propto \nu^{\alpha_\nu}$) of the full sample shows a break at ~912\AA\ with spectral index $\alpha_\nu=-0.61\pm0.01$ at $\lambda>912$\AA\ and a softening at shorter wavelengths ($\alpha_\nu=-1.70 \pm 0.61$ at $\lambda\leq 912$\AA). Our analysis proves that a proper intergalactic medium absorption correction is required to establish the intrinsic continuum emission of quasars. We interpret our average ultraviolet spectrum in the context of photoionisation, accretion disk models, and quasar contribution to the ultraviolet background. We find that observed broad line ratios are consistent with those predicted assuming an ionising slope of $\alpha_\mathrm{ion}=$-2.0, similar to the observed ionising spectrum in the same wavelength range. The continuum break and softening are consistent with accretion disk plus X-ray corona models when black hole spin is taken into account. Our spectral energy distribution yields a 30% increase to previous estimates of the specific quasar emissivity, such that quasars may contribute significantly to the total specific Lyman limit emissivity estimated from the Ly$\alpha$ forest at z<3.2.

We present results from a survey designed to probe the star formation properties of 32 damped Ly-alpha systems (DLAs) at z~2.7. By using the "double-DLA" technique that eliminates the glare of the bright background quasars, we directly measure the rest-frame FUV flux from DLAs and their neighbouring galaxies. At the position of the absorbing gas, we place stringent constraints on the unobscured star formation rates (SFRs) of DLAs to 2-sigma limits of <0.09-0.27 M/yr, corresponding to SFR surface densities <10^(-2.6)-10^(-1.5) M/yr/kpc^2. The implications of these limits for the star formation law, metal enrichment, and cooling rates of DLAs are examined. By studying the distribution of impact parameters as a function of SFRs for all the galaxies detected around these DLAs, we place new direct constraints on the bright end of the UV luminosity function of DLA hosts. We find that <13% of the hosts have SFR>2 M/yr at impact parameters b < (SFR/(M/yr))^(0.8)+6 kpc, differently from current samples of confirmed DLA galaxies. Our observations also disfavor a scenario in which the majority of DLAs arise from bright LBGs at distances 20 < b < 100 kpc. These new findings corroborate a picture in which DLAs do not originate from highly star forming systems that are coincident with the absorbers, and instead suggest that DLAs are associated with faint, possibly isolated, star-forming galaxies. Potential shortcomings of this scenario and future strategies for further investigation are discussed.

We analyze the physical conditions of the cool, photoionized (T $\sim 10^4$ K) circumgalactic medium (CGM) using the COS-Halos suite of gas column density measurements for 44 gaseous halos within 160 kpc of $L \sim L^*$ galaxies at $z \sim 0.2$. These data are well described by simple photoionization models, with the gas highly ionized (n$_{\rm HII}$/n$_{\rm H} \gtrsim 99\%$) by the extragalactic ultraviolet background (EUVB). Scaling by estimates for the virial radius, R$_{\rm vir}$, we show that the ionization state (tracked by the dimensionless ionization parameter, U) increases with distance from the host galaxy. The ionization parameters imply a decreasing volume density profile n$_{\rm H}$ = (10$^{-4.2 \pm 0.25}$)(R/R$_{\rm vir})^{-0.8\pm0.3}$. Our derived gas volume densities are several orders of magnitude lower than predictions from standard two-phase models with a cool medium in pressure equilibrium with a hot, coronal medium expected in virialized halos at this mass scale. Applying the ionization corrections to the HI column densities, we estimate a lower limit to the cool gas mass M$_{\rm CGM}^{\rm cool} > 6.5 \times 10^{10}$ M$_{\odot}$ for the volume within R $<$ R$_{\rm vir}$. Allowing for an additional warm-hot, OVI-traced phase, the CGM accounts for at least half of the baryons purported to be missing from dark matter halos at the 10$^{12}$ M$_{\odot}$ scale.

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 the first results from our NASA Keck Observatory Database of Ionized Absorbers toward Quasars (KODIAQ) survey which aims to characterize the properties of the highly ionized gas of high redshift galaxies and their circumgalactic medium (CGM) at 2<z<4. We select absorbers optically thick at the Lyman limit ({\tau}LL > 1, log N(HI) > 17.3) as probes of these galaxies and their CGM where both transitions of the O VI doublet have little contamination from the Ly {\alpha}, {\beta} forests. We found 20 absorbers that satisfy these rules: 7 Lyman limit systems (LLSs), 8 super-LLSs (SLLSs) and 5 damped Ly{\alpha} (DLAs). The O VI detection rate is 100% for the DLAs, 71% for the LLSs, and 63% for the SLLSs. When O VI is detected, log N(OVI)=14.9+/-0.3, an average O VI column density substantially larger and with a smaller dispersion than found in blind O VI surveys at similar redshifts. Strong O VI absorption is therefore nearly ubiquitous in the CGM of z~2-3 galaxies. The total velocity widths of the O VI profiles are also large (200<Dv(OVI)<400 km/s). These properties are quite similar to those seen for O VI in low z star-forming galaxies, and therefore we hypothesize that these strong CGM O VI absorbers (with {\tau}LL > 1) at 2<z<3.5 also probe outflows of star-forming galaxies. The LLSs and SLLSs with no O VI absorption have properties consistent with those seen in cosmological simulations tracing cold streams feeding galaxies. When the highly ionized (Si IV and O VI) gas is taken into account, we determine that the {\tau}LL > 1 absorbers could contain as much as 3-14% of the cosmic baryon budget at z~2-3, only second to the Ly{\alpha} forest. We conservatively show that 5-20% of the metals ever produced at z~2-3 are in form of highly ionized metals ejected in the CGM of galaxies.

We have obtained spectra of 163 quasars at $z_\mathrm{em}>4.4$ with the Gemini Multi Object Spectrometers on the Gemini North and South telescopes, the largest publicly available sample of high-quality, low-resolution spectra at these redshifts. From this homogeneous data set, we generated stacked quasar spectra in three redshift intervals at $z\sim 5$. We have modelled the flux below the rest-frame Lyman limit ($\lambda_\mathrm{r}<912$\AA) to assess the mean free path $\lambda_\mathrm{mfp}^{912}$ of the intergalactic medium to HI-ionizing radiation. At mean redshifts $z_\mathrm{q}=4.56$, 4.86 and 5.16, we measure $\lambda_\mathrm{mfp}^{912}=(22.2\pm 2.3, 15.1\pm 1.8, 10.3\pm 1.6)h_{70}^{-1}$ proper Mpc with uncertainties dominated by sample variance. Combining our results with $\lambda_\mathrm{mfp}^{912}$ measurements from lower redshifts, the data are well modelled by a simple power-law $\lambda_\mathrm{mfp}^{912}=A[(1+z)/5]^\eta$ with $A=(37\pm 2)h_{70}^{-1}$ Mpc and $\eta = -5.4\pm 0.4$ between $z=2.3$ and $z=5.5$. This rapid evolution requires a physical mechanism -- beyond cosmological expansion -- which reduces the cosmic effective Lyman limit opacity. We speculate that the majority of HI Lyman limit opacity manifests in gas outside galactic dark matter haloes, tracing large-scale structures (e.g. filaments) whose average density (and consequently neutral fraction) decreases with cosmic time. Our measurements of the strongly redshift-dependent mean free path shortly after the completion of HI reionization serve as a valuable boundary condition for numerical models thereof. Having measured $\lambda_\mathrm{mfp}^{912}\approx 10$ Mpc at $z=5.2$, we confirm that the intergalactic medium is highly ionized by that epoch and that the redshift evolution of the mean free path does not show a break that would indicate a recent end to HI reionization.

We examine recent measurements of the z~2.5 intergalactic medium (IGM) which constrain the HI frequency distribution f(N) and the mean free path MFP to ionizing radiation. We argue that line-blending and the clustering of strong absorption-line systems have led previous authors to systematically overestimate the effective Lyman limit opacity, yielding too small of a MFP for the IGM. We further show that recently published measurements of f(N) at log NHI ~ 16 lie in strong disagreement, implying underestimated uncertainty from sample variance and/or systematics like line-saturation. Allowing for a larger uncertainty in the f(N) measurements, we provide a new cubic Hermite spline model for f(N) which reasonably satisfies all of the observational constraints under the assumption of randomly distributed absorption systems. We caution, however, that this formalism is invalid in light of absorber clustering and use a toy model to estimate the effects. Future work must properly account for the non-Poissonian nature of the IGM.

We have performed an analysis of over 34,000 MgII doublets at 0.36 < z < 2.29 in Sloan Digital Sky Survey (SDSS) Data-Release 7 quasar spectra; the catalog, advanced data products, and tools for analysis are publicly available. The catalog was divided into 14 small redshift bins with roughly 2500 doublets in each, and from Monte-Carlo simulations, we estimate 50% completeness at rest equivalent width W_r ~ 0.8 Angstrom. The equivalent-width frequency distribution is described well by an exponential model at all redshifts, and the distribution becomes flatter with increasing redshift, i.e., there are more strong systems relative to weak ones. Direct comparison with previous SDSS MgII surveys reveal that we recover at least 70% of the doublets in these other catalogs, in addition to detecting thousands of new systems. We discuss how these surveys come by their different results, which qualitatively agree but, due to the very small uncertainties, differ by a statistically significant amount. The estimated physical cross-section of MgII-absorbing galaxy halos increased three-fold, approximately, from z = 0.4 --> 2.3, while the W_r >= 1 Angstrom absorber line density grew, dN_MgII/dX, by roughly 45%. Finally, we explore the different evolution of various absorber populations - damped Lyman-alpha absorbers, Lyman-limit systems, strong CIV absorbers, and strong and weaker MgII systems - across cosmic time (0 < z < 6).

We report on the detection of deuterated molecular hydrogen, HD, at $z = 0.18$. HD and H$_{\rm 2}$ are detected in HST/COS data of a low metallicity ($Z \sim 0.07Z_\odot$) damped Ly$\alpha$ system at $z = 0.18562$ toward QSO B0120$-$28, with log $N$(H I) = 20.50 $\pm$ 0.10. Four absorption components are clearly resolved in H$_{\rm 2}$ while two components are resolved in HD; the bulk of the molecular hydrogen is associated with the components traced by HD. We find total column densities log $N$(HD) = 14.82 $\pm$ 0.15 and log $N$(H$_{\rm 2}$) = 20.00 $\pm$ 0.10. This system has a high molecular fraction, $f$(H$_{\rm 2}$) = 0.39 $\pm$ 0.10 and a low HD to H$_{\rm 2}$ ratio, log (HD/2H$_{\rm 2}$) $= -5.5 \pm 0.2$ dex. The excitation temperature, $T_{01} = 65 \pm 2$ K, in the component containing the bulk of the molecular gas is lower than in other DLAs. These properties are unlike those in other higher redshift DLA systems known to contain HD, but are consistent with what is observed in dense clouds in the Milky Way.

We present the cross-correlation function of MgII absorbers with respect to a volume-limited sample of luminous red galaxies (LRGs) at z=0.45-0.60 using the largest MgII absorber sample and a new LRG sample from SDSS DR7. We present the clustering signal of absorbers on projected scales r_p = 0.3-35 Mpc/h in four Wr(2796) bins spanning Wr(2796)=0.4-5.6A. We found that on average MgII absorbers reside in halos < log M_h > \approx 12.1, similar to the halo mass of an L_* galaxy. We report that the weakest absorbers in our sample with W_r(2796)=0.4-1.1A reside in relatively massive halos with < log M_h > \approx 12.5^{+0.6}_{-1.3}, while stronger absorbers reside in halos of similar or lower masses < log M_h > \approx 11.6^{+0.9}. We compared our bias data points, b, and the frequency distribution function of absorbers, f_{W_r}, with a simple model incorporating an isothermal density profile to mimic the distribution of absorbing gas in halos. We also compared the bias data points with Tinker & Chen (2008) who developed halo occupation distribution models of MgII absorbers that are constrained by b and f_{W_r}. The simple isothermal model can be ruled at a \approx 2.8\sigma level mostly because of its inability to reproduce f_{W_r}. However, b values are consistent with both models, including TC08. In addition, we show that the mean b of absorbers does not decrease beyond W_r(2796) \approx 1.6A. The flat or potential upturn of b for Wr(2796) \gtrsim 1.6A absorbers suggests the presence of additional cool gas in massive halos.

We study the high-ionization phase and kinematics of the circumgalactic medium around low-redshift galaxies using a sample of 23 Lyman Limit Systems (LLSs) at 0.08<z<0.93 observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. In Lehner et al. (2013), we recently showed that low-z LLSs have a bimodal metallicity distribution. Here we extend that analysis to search for differences between the high-ion and kinematic properties of the metal-poor and metal-rich branches. We find that metal-rich LLSs tend to show higher O VI columns and broader O VI profiles than metal-poor LLSs. The total H I line width (dv90 statistic) in LLSs is not correlated with metallicity, indicating that the H I kinematics alone cannot be used to distinguish inflow from outflow and gas recycling. Among the 17 LLSs with O VI detections, all but two show evidence of kinematic sub-structure, in the form of O VI-H I centroid offsets, multiple components, or both. Using various scenarios for how the metallicity in the high-ion and low-ion phases of each LLS compare, we constrain the ionized hydrogen column in the O VI phase to lie in the range log N(H II)~17.6-20. The O VI phase of LLSs is a substantial baryon reservoir, with M(high-ion)~10^{8.5-10.9}(r/150 kpc)^2 solar masses, similar to the mass in the low-ion phase. Accounting for the O VI phase approximately doubles the contribution of low-z LLSs to the cosmic baryon budget.

We present the design and methods of the COS-Halos survey, a systematic investigation of the gaseous halos of 44 z = 0.15-0.35 galaxies using background QSOs observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. This survey has yielded 39 spectra of z_em ~ 0.5 QSOs with S/N ~ 10-15 per resolution element. The QSO sightlines pass within 150 physical kpc of the galaxies, which span early and late types over stellar mass log M* / Msun= 9.5 - 11.5. We find that the CGM exhibits strong HI, averaging 1 Ang in Lya equivalent width out to 150 kpc, with 100% covering fraction for star-forming galaxies and 75% covering for passive galaxies. We find good agreement in column densities between this survey and previous studies over similar range of impact parameter. There is weak evidence for a difference between early- and late-type galaxies in the strength and distribution of HI. Kinematics indicate that the detected material is bound to the host galaxy, such that >~90% of the detected column density is confined within +/-200 km s^-1 of the galaxies. This material generally exists well below the halo virial temperatures at T<~ 10^5 K. We evaluate a number of possible origin scenarios for the detected material, and in the end favor a simple model in which the bulk of the detected HI arises in a bound, cool, low-density photoionized diffuse medium that is generic to all L* galaxies and may harbor a total gaseous mass comparable to galactic stellar masses.

We present results of a blind survey of Lyman limit systems (LLSs) detected in absorption against 105 quasars at z~3 using the blue sensitive MagE spectrograph at the Magellan Clay telescope. By searching for Lyman limit absorption in the wavelength range ~3000-4000A, we measure the number of LLSs per unit redshift l(z)=1.21+/-0.28 at z~2.8. Using a stacking analysis, we further estimate the mean free path of ionizing photons in the z~3 Universe lambda^912 = 100+/-29 Mpc/h_70.4. Combined with our LLS survey, we conclude that systems with log N_HI >= 17.5 cm^-2 contribute only ~40% to the observed mean free path at these redshifts. Further, with the aid of photo-ionization modeling, we infer that a population of ionized and metal poor systems is likely required to reproduce the metal line strengths observed in a composite spectrum of 20 LLSs with N_HI = 17.5-19 cm^-2 at z~2.6-3.0. Finally, with a simple toy model, we deduce that gas in the halos of galaxies can alone account for the totality of LLSs at z<~3, but a progressively higher contribution from the intergalactic medium is required beyond z~3.5. We also show how the weakly evolving number of LLSs per unit redshift at z<~3 can be modeled either by requiring that the spatial extent of the circumgalactic medium is redshift invariant in the last ~10 Gyr of cosmic evolution or by postulating that LLSs arise in halos that are rare fluctuations in the density field at each redshift.

We assess the metal content of the cool (10^4 K) circumgalactic medium (CGM) about galaxies at z<1 using an H I-selected sample of 28 Lyman limit systems (LLS, defined here as absorbers with 16.2<log N(H I)<18.5) observed in absorption against background QSOs by the Cosmic Origins Spectrograph on-board the Hubble Space Telescope. The N(H I) selection avoids metallicity biases inherent in many previous studies of the low-redshift CGM. We compare the column densities of weakly ionized metal species (e.g., O II, Si II, Mg II) to N(H I) in the strongest H I component of each absorber. We find that the metallicity distribution of the LLS (and hence the cool CGM) is bimodal with metal-poor and metal-rich branches peaking at [X/H]=-1.6 and -0.3 (or about 2.5% and 50% solar metallicities). The cool CGM probed by these LLS is predominantly ionized. The metal-rich branch of the population likely traces winds, recycled outflows, and tidally stripped gas; the metal-poor branch has properties consistent with cold accretion streams thought to be a major source of fresh gas for star forming galaxies. Both branches have a nearly equal number of absorbers. Our results thus demonstrate there is a significant mass of previously-undiscovered cold metal-poor gas and confirm the presence of metal enriched gas in the CGM of z<1 galaxies.

Measuring rest-frame ultraviolet rotational transitions from the Lyman and Werner bands in absorption against a bright background continuum is one of the few ways to directly measure molecular hydrogen (H2). Here we report the detection of Lyman-Werner absorption from H2 at z=0.56 in a sub-damped Ly-alpha system with neutral hydrogen column density N(HI) = 10^(19.5 +/- 0.2) cm^-2. This is the first H2 system analysed at a redshift < 1.5 beyond the Milky Way halo. It has a surprisingly high molecular fraction: log f(H2) > -1.93 +/- 0.36 based on modelling the line profiles, with a robust model-independent lower limit of f(H2) > 10^-3. This is higher than f(H2) values seen along sightlines with similar N(HI) through the Milky Way disk and the Magellanic clouds. The metallicity of the absorber is 0.19 +0.21 -0.10 solar, with a dust-to-gas ratio < 0.36 times the value in the solar neighbourhood. Absorption from associated low-ionisation metal transitions such as OI and FeII is observed in addition to OVI. Using Cloudy models we show that there are three phases present; a ~100 K phase giving rise to H2, a ~10^4 K phase where most of the low-ionisation metal absorption is produced; and a hotter phase associated with OVI. Based on similarities to high velocity clouds in the Milky Way halo showing H2 and the presence of two nearby galaxy candidates with impact parameters of ~10 kpc, we suggest that the absorber may be produced by a tidally-stripped structure similar to the Magellanic Stream.

We present the first science results from our Hubble Space Telescope Survey for Lyman limit absorption systems (LLS) using the low dispersion spectroscopic modes of the Advanced Camera for Surveys and the Wide Field Camera 3. Through an analysis of 71 quasars, we determine the incidence frequency of LLS per unit redshift and per unit path length, l(z) and l(x) respectively, over the redshift range 1 < z< 2.6, and find a weighted mean of l(x)=0.29 +/-0.05 for 2.0 < z < 2.5 through a joint analysis of our sample and that of Ribaudo et al. (2011). Through stacked spectrum analysis, we determine a median (mean) value of the mean free path to ionizing radiation at z=2.4 of lambda_mfp = 243(252)h^(-1) Mpc, with an error on the mean value of +/- 43h^(-1) Mpc. We also re-evaluate the estimates of lambda_mfp from Prochaska et al. (2009) and place constraints on the evolution of lambda_mfp with redshift, including an estimate of the "breakthrough" redshift of z = 1.6. Consistent with results at higher z, we find that a significant fraction of the opacity for absorption of ionizing photons comes from systems with N_HI <= 10^{17.5} cm^(-2) with a value for the total Lyman opacity of tau_lyman = 0.40 +/- 0.15. Finally, we determine that at minimum, a 5-parameter (4 power-law) model is needed to describe the column density distribution function f(N_HI, X) at z \sim 2.4, find that f(N_HI,X) undergoes no significant change in shape between z \sim 2.4 and z \sim 3.7, and provide our best fit model for f(N_HI,X).

We have vastly increased the CIV statistics at intermediate redshift by surveying the thousands of quasars in the Sloan Digital Sky Survey Data-Release 7. We visually verified over 16,000 CIV systems with 1.46 < z < 4.55---a sample size that renders Poisson error negligible. Detailed Monte Carlo simulations show we are approximately 50% complete down to rest equivalent widths W_r ~ 0.6 \AA. We analyzed the sample as a whole and in ten small redshift bins with approximately 1500 doublets each. The equivalent width frequency distributions f(W_r) were well modeled by an exponential, with little evolution in shape. In contrast with previous studies that modeled the frequency distribution as a single power law, the fitted exponential gives a finite mass density for the CIV ions. The co-moving line density dN_CIV/dX evolved smoothly with redshift, increasing by a factor of 2.37+/-0.09 from z = 4.55 to 1.96, then plateauing at dN_CIV/dX ~ 0.34 for z = 1.96 to 1.46. Comparing our SDSS sample with z < 1 (ultraviolet) and z > 5 (infrared) surveys, we see an approximately 10-fold increase in dN_CIV/dX over z ~ 6 --> 0, for W_r >= 0.6 \AA. This suggests a monotonic and significant increase in the enrichment of gas outside galaxies over the 12 Gyr lifetime of the universe.

We present the equivalent width and column density measurements for low and intermediate ionization states of the circumgalactic medium (CGM) surrounding 44 low-z, L ~ L* galaxies drawn from the COS-Halos survey. These measurements are derived from far-UV transitions observed in HST/COS and Keck/HIRES spectra of background quasars within an impact parameter R < 160 kpc to the targeted galaxies. The data show significant metal-line absorption for 33 of the 44 galaxies, including quiescent systems, revealing the common occurance of a cool (T ~ 10^{4 - 5} K), metal-enriched CGM. The detection rates and column densities derived for these metal lines decrease with increasing impact parameter, a trend we interpret as a declining metal surface density profile for the CGM. A comparison of the relative column densities of adjacent ionization states indicates the gas is predominantly ionized. The large surface density in metals demands a large reservoir of metals and gas in the cool CGM (very conservatively, M_ CGMcool > 10^9 MSun), which likely traces a distinct density and/or temperature regime from the highly-ionized CGM traced by OVI absorption. The large dispersion in absorption strengths (including non-detections) suggests the cool CGM traces a wide range of densities or a mix of local ionizing conditions. Lastly, the kinematics inferred from the metal-line profiles are consistent with the cool CGM being bound to the dark matters halos hosting the galaxies; this gas may serve as fuel for future star-formation. Future work will leverage this dataset to provide estimates on the mass, metallicity, dynamics, and origin of the cool CGM in low-z, L* galaxies.

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We find that detections of HI surrounding early type galaxies are typically as common and strong as around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in the two populations. For early type galaxies, the covering fraction for HI absorption above 10^16 cm^2 is ~40-50% within ~150 kpc. Line widths and kinematics of the detected material show it to be cold (T ~< 10^5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, photoionized CGM baryons may be up to 10^9 --- 10^11 Msun. Contrary to some theoretical expectations, strong halo HI absorbers do not disappear as part of the quenching of star-formation. Even passive galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently observed inside ETGs.

PKS0447-439 is one of the BL Lacertae objects that have been detected at very high energy. There has been a recent report of a lower limit of z>1.246 for the redshift of this blazar, challenging the current paradigm in which very high-energy gamma-rays cannot freely propagate in the z>1 universe. In this research note, we present a new MagE/Magellan spectrum of PKS0447-439 with exquisite signal-to-noise (S/N>150 at 6500A). Our analysis confirms the presence of the previously-reported absorption line at 6280A, which we identify, however, with a known telluric absorption, invalidating the claim that this blazar lies at z>1. Since no other extragalactic spectral features are detected, we cannot establish a redshift based on our spectrum.

The circumgalactic medium (CGM) is fed by galaxy outflows and accretion of intergalactic gas, but its mass, heavy element enrichment, and relation to galaxy properties are poorly constrained by observations. In a survey of the outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos of ionized oxygen surrounding star-forming galaxies, but we find much less ionized oxygen around galaxies with little or no star formation. This ionized CGM contains a substantial mass of heavy elements and gas, perhaps far exceeding the reservoirs of gas in the galaxies themselves. It is a basic component of nearly all star-forming galaxies that is removed or transformed during the quenching of star formation and the transition to passive evolution.

In the current cosmological model, only the three lightest elements were created in the first few minutes after the Big Bang; all other elements were produced later in stars. To date, however, heavy elements have been observed in all astrophysical environments. We report the detection of two gas clouds with no discernible elements heavier than hydrogen. These systems exhibit the lowest heavy-element abundance in the early universe and thus are potential fuel for the most metal poor halo stars. The detection of deuterium in one system at the level predicted by primordial nucleosynthesis provides a direct confirmation of the standard cosmological model. The composition of these clouds further implies that the transport of heavy elements from galaxies to their surroundings is highly inhomogeneous.

We present high signal-to-noise optical spectra for 67 low-redshift (0.1 < z < 0.4) galaxies that lie within close projected distances (5 kpc < rho < 150 kpc) of 38 background UV-bright QSOs. The Keck LRIS and Magellan MagE data presented here are part of a survey that aims to construct a statistically sampled map of the physical state and metallicity of gaseous galaxy halos using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). We provide a detailed description of the optical data reduction and subsequent spectral analysis that allow us to derive the physical properties of this uniquely data-rich sample of galaxies. The galaxy sample is divided into 38 pre-selected L ~ L*, z ~ 0.2 "target" galaxies and 29 "bonus" galaxies that lie in close proximity to the QSO sightlines. We report galaxy spectroscopic redshifts accurate to +/- 30 km s-1, impact parameters, rest-frame colors, stellar masses, total star formation rates, and gas-phase interstellar medium oxygen abundances. When we compare the distribution of these galaxy characteristics to those of the general low-redshift population, we find good agreement. The L ~ L* galaxies in this sample span a diverse range of color (1.0 < u-r < 3.0), stellar mass (10^9.5 < M/M_sun < 10^11.5), and SFRs (0.01 - 20 M_sun yr-1). These optical data, along with the COS UV spectroscopy, comprise the backbone of our efforts to understand how halo gas properties may correlate with their host galaxy properties, and ultimately to uncover the processes that drive gas outflow and/or are influenced by gas inflow.

We present a new determination of the intergalactic CIV mass density at 4.3 < z < 6.3. Our constraints are derived from high signal-to-noise spectra of seven quasars at z > 5.8 obtained with the newly commissioned FIRE spectrograph on the Magellan Baade telescope, coupled with six observations of northern objects taken from the literature. We confirm the presence of a downturn in the CIV abundance at <z>=5.66 by a factor of 4.1 relative to its value at <z>=4.96, as measured in the same sightlines. In the FIRE sample, a strong system previously reported in the literature as CIV at z=5.82 is re-identified as MgII at z=2.78, leading to a substantial downward revision in $\Omega_{CIV}$ for these prior studies. Additionally we confirm the presence of at least two systems with low-ionization CII, SiII, and OI absorption but relatively weak signal from CIV. The latter systems systems may be of interest if the downward trend in $\Omega_{CIV}$ at high redshift is driven in part by ionization effects.

We have serendipitously detected a strong O VI-bearing Lyman limit system at z_abs = 0.3558 toward the QSO J1009+0713 (z_em = 0.456) in our survey of low-redshift galaxy halos with the Hubble Space Telescope's Cosmic Origins Spectrograph. Its rest-frame equivalent width of W_r = 835 +/- 49 mA is the highest for an intervening absorber yet detected in any low-redshift QSO sightline, with absorption spanning 400 km s^-1 in its rest frame. HST/WFC3 images of the galaxy field show that the absorber is associated with two galaxies lying at 14 and 46 kpc from the QSO line of sight. The bulk of the absorbing gas traced by H I resides in two strong, blended component groups that possess a total logN(HI) = 18 - 18.8. The ion ratios and column densities of C, N, O, Mg, Si, S, and Fe, except the O VI, can be accommodated into a simple photoionization model in which diffuse, low-metallicity halo gas is exposed to a photoionizing field from stars in the nearby galaxies that propagates into the halo at 10% efficiency. We constrain the metallicity firmly within the range 0.1 - 1 Zsun, and photoionization modeling indirectly indicates a subsolar metallicity of 0.05 - 0.5 Zsun. The appearance of strong O VI and nine Mg II components and our review of similar systems in the literature support the "interface" picture of high-velocity O VI: the total strength of the O VI shows a positive correlation with the number of detected components in the low-ionization gas, however the total O VI column densities still far exceed the values expected from interface models for the number of detected clouds.

We report on the first Cosmic Origins Spectrograph (COS) observations of damped and sub-damped Lyman-alpha (DLA) systems discovered in a new survey of the gaseous halos of low-redshift galaxies. From observations of 37 sightlines, we have discovered three DLAs and four sub-DLAs. We measure the neutral gas density Omega(HI), and redshift density dN/dz, of DLA and sub-DLA systems at z<0.35. We find dN/dz=0.25 and Omega(HI)=1.4x10^-3 for DLAs, and dN/dz=0.08 with Omega(HI)=4.2x10^-5 for sub-DLAs over a redshift path delta z=11.9. To demonstrate the scientific potential of such systems, we present a detailed analysis of the DLA at z=0.1140 in the spectrum of SDSS J1009+0713. Profile fits to the absorption lines determine log N(H I)=20.68pm0.10 with a metallicity determined from the undepleted element Sulfur of [S/H]=-0.62pm0.18. The abundance pattern of this DLA is similar to that of higher z DLAs, showing mild depletion of the refractory elements Fe and Ti with [S/Fe]=+0.24pm0.22 and [S/Ti]=+0.28pm0.15. Nitrogen is underabundant in this system with [N/H]=-1.40pm0.14, placing this DLA below the plateau of the [N/alpha] measurements in the local Universe at similar metallicities. This DLA has a simple kinematic structure with only two components required to fit the profiles and a kinematic width of 52 km/s. Imaging of the QSO field with WFC3 reveals a spiral galaxy at very small impact parameter to the QSO and several galaxies within 10". Followup spectra with LRIS reveal that none of the nearby galaxies are at the redshift of the DLA. The spiral galaxy is identified as the host galaxy of the QSO based on the near perfect alignment of the nucleus and disk of the galaxy as well as spectra of an H II region showing emission lines at the QSO redshift. A small feature appears 0.70" from the nucleus of the QSO after PSF subtraction, providing another candidate for the host galaxy of the DLA. (abb)

We present the methodology for, and the first results from, a new imaging program aimed at identifying and characterizing the host galaxies of damped Lyman-alpha absorbers (DLAs) at z>2. We target quasar sightlines with multiple optically-thick HI absorbers and use the higher-redshift system as a "blocking filter" (via its Lyman-limit absorption) to eliminate all far-ultraviolet (FUV) emission from the quasar. This allows us to directly image the rest-frame FUV continuum emission of the lower-redshift DLA, without any quasar contamination and with no bias towards large impact parameters. We introduce a formalism based on galaxy number counts and Bayesian statistics with which we quantify the probability that a candidate is the DLA host galaxy. This method will allow the identification of a bona fide sample of DLAs that are too faint to be spectroscopically confirmed. The same formalism can be adopted to the study of other quasar absorption line systems (e.g. MgII absorbers). We have applied this imaging technique to two QSO sightlines. For the z~2.69 DLA towards J073149+285449, a galaxy with impact parameter b=1.54"=11.89 kpc and implied star formation rate (SFR) of ~5 M/yr is identified as the most reliable candidate. In the case of the z~2.92 DLA towards J211444-005533, no likely host is found down to a 3-sigma SFR limit of 1.4 M/yr. Studying the HI column density as a function of the impact parameter, including 6 DLAs with known hosts from the literature, we find evidence that the observed HI distribution is more extended than what is generally predicted from numerical simulation.

We perform a semi-automated survey for tau>=2 Lyman Limit systems (LLSs) in quasar spectra from the Sloan Digital Sky Survey, Data Release 7. From a starting sample of 2473 quasars with zem=3.6-4.4, we analyze 469 spectra meeting strict seletion criteria for a total redshift path Dz=93.8 and identify 192 intervening systems at z>3.3. The incidence of tau>=2 LLSs per unit redshift, l(z), is well described by a single-power law at these redshifts: l(z) = C_LLS [(1+z)/(1+z_*)]^gamma, with z_*=3.7, C_LLS = 1.9+/-0.2, and gamma = 5.2+/-1.5 (68% c.l.). These values are systematically lower than previous estimates (especially at z<4) but are consistent with recent measurements of the mean free path to ionizing radiation. Extrapolations of this power-law to z=0 are inconsistent with previous estimations of l(z) at z<1 and suggest a break at z~2, similar to that observed for the Lya forest. Our results also indicate that the systems giving rise to LLS absorption decrease by ~50% in comoving number density and/or physical size from z=4 to 3.3, perhaps due to an enhanced extragalactic ultraviolet background. The observations place an integral constraint on the HI frequency distribution f(N_HI,X) and indicate that the power-law slope beta= dln[f(N,X)]/dln[N] is likely shallower than beta = -1 at N_HI=10^18 cm^-2. Including other constraints on f(N_HI,X) from the literature, we infer that beta is steeper than beta = -1.7 at N_HI~10^15 cm^-2, implying at least two inflections in f(N_HI,X). We also perform a survey for proximate LLSs (PLLSs) and find that l(z)_PLLS is systematically lower ~25% than intervening systems. Finally, we estimate that systematic effects impose an uncertainty of 10-20% in the l(z) measurements; these effects may limit the precision of all future surveys.

We present an absorption line analysis of the Lyman limit system (LLS) at z=3.55 in our Magellan/MIKE spectrum of PKS2000-330. Our analysis of the Lyman limit and full HI Lyman series constrains the total HI column density of the LLS (N_HI = 10^[18.0 +/- 0.25] cm^{-2} for b_HI >= 20 km/s) and also the N_HI values of the velocity subsystems comprising the absorber. We measure ionic column densities for metal-line transitions associated with the subsystems and use these values to constrain the ionization state (>90% ionized) and relative abundances of the gas. We find an order of magnitude dispersion in the metallicities of the subsystems, marking the first detailed analysis of metallicity variations in an optically thick absorber. The results indicate that metals are not well mixed within the gas surrounding high $z$ galaxies. Assuming a single-phase photoionization model, we also derive an N_H-weighted metallicity, <[Si/H]> = -1.66 +/- 0.25, which matches the mean metallicity in the neutral ISM in high z damped Lya systems (DLAs). Because the line density of LLSs is ~10 times higher than the DLAs, we propose that the former dominate the metal mass-density at z~3 and that these metals reside in the galaxy/IGM interface. Considerations of a multi-phase model do not qualitatively change these conclusions. Finally, we comment on an anomalously large O^0/Si^+ ratio in the LLS that suggests an ionizing radiation field dominated by soft UV sources (e.g. a starburst galaxy). Additional abundance analysis is performed on the super-LLS systems at z=3.19.

We present a new method to directly measure the opacity from HI Lyman limit (LL) absorption k_LL along quasar sightlines by the intergalactic medium (IGM). The approach analyzes the average (``stacked'') spectrum of an ensemble of quasars at a common redshift to infer the mean free path (MFP) to ionizing radiation. We apply this technique to 1800 quasars at z=3.50-4.34 drawn from the Sloan Digital Sky Survey (SDSS), giving the most precise measurements on k_LL at any redshift. From z=3.6 to 4.3, the opacity increases steadily as expected and is well parameterized by MFP = (48.4 +/- 2.1) - (38.0 +/- 5.3)*(z-3.6) h^-1 Mpc (proper distance). The relatively high MFP values indicate that the incidence of systems which dominate k_LL evolves less strongly at z>3 than that of the Lya forest. We infer a mean free path three times higher than some previous estimates, a result which has important implications for the photo-ionization rate derived from the emissivity of star forming galaxies and quasars. Finally, our analysis reveals a previously unreported, systematic bias in the SDSS quasar sample related to the survey's color targeting criteria. This bias potentially affects all z~3 IGM studies using the SDSS database.

We present optical echelle spectra of four gamma-ray burst (GRB) afterglows (GRB 050730, GRB 050820, GRB 051111, and GRB 060418) discovered during the first 1.5 years of operation of the Swift satellite and localized by either the Swift telescope or follow-up ground-based imaging. We analyze the spectra to derive accurate column density measurements for the transitions arising in the interstellar medium (ISM) of the GRB host galaxies. These measurements can be used to constrain the physical properties of the ISM including the metallicity, dust-to-gas ratio, ionization state, and chemical abundances of the gas. We also present measurements of the strong MgII systems in the GRB afterglow spectra. With the publication of this paper, we provide the first data release of echelle afterglow spectra by the GRAASP collaboration to the general community.

We present the results of a survey for super Lyman limit systems (SLLS; defined to be absorbers with 19.0 <= log(NHI) <= 20.3 cm^-2) from a large sample of high resolution spectra acquired using the Keck and Magellan telescopes. Specifically, we present 47 new SLLS from 113 QSO sightlines. We focus on the neutral hydrogen frequency distribution f(N,X) of the SLLS and its moments, and compare these results with the Lyman-alpha forest and the damped Lyman alpha systems (DLA; absorbers with log(NHI) >= 20.3 cm^-2). We find that that f(N,X) of the SLLS can be reasonably described with a power-law of index alpha = -1.43^{+0.15}_{-0.16} or alpha = -1.19^{+0.20}_{-0.21} depending on whether we set the lower N(HI) bound for the analysis at 10^{19.0} cm^-2 or 10^{19.3}$ cm^-2, respectively. The results indicate a flattening in the slope of f(N,X) between the SLLS and DLA. We find little evidence for redshift evolution in the shape of f(N,X) for the SLLS over the redshift range of the sample 1.68 < z < 4.47 and only tentative evidence for evolution in the zeroth moment of f(N,X), the line density l_lls(X). We introduce the observable distribution function O(N,X) and its moment, which elucidates comparisons of HI absorbers from the Lyman-alpha through to the DLA. We find that a simple three parameter function can fit O(N,X) over the range 17.0 <= log(NHI) <=22.0. We use these results to predict that f(N,X) must show two additional inflections below the SLLS regime to match the observed f(N,X) distribution of the Lyman-alpha forest. Finally, we demonstrate that SLLS contribute a minor fraction (~15%) of the universe's hydrogen atoms and, therefore, an even small fraction of the mass in predominantly neutral gas.

We present a measurement of the D/H abundance ratio in a metal-poor damped Lyman alpha (DLA) system along the sightline of QSO SDSS1558-0031. The DLA system is at redshift z = 2.70262, has a neutral column density of log(NHI)=20.67+/-0.05 cm^2, and a gas-phase metallicity [O/H]= -1.49 which indicates that deuterium astration is negligible. Deuterium absorption is observed in multiple Lyman series with a column density of log(NDI)=16.19+/-0.04 cm^2, best constrained by the deuterium Lyman-11 line. We measure log(D/H) = -4.48+/-0.06, which when combined with previous measurements along QSO sightlines gives a best estimate of log(D/H) = -4.55+/-0.04, where the 1-sigma error estimate comes from a jackknife analysis of the weighted means. Using the framework of standard big bang nucleosynthesis, this value of D/H translates into a baryon density of Omega_b h^2 = 0.0213 +/- 0.0013 +/- 0.0004 where the error terms represent the 1-sigma errors from D/H and the uncertainties in the nuclear reaction rates respectively. Combining our new measurement with previous measurements of D/H, we no longer find compelling evidence for a trend of D/H with NHI.

We present abundance measurements for two super Lyman Limit systems (SLLS; quasar absorption line systems with 10^19 cm^-2 < N_HI < 10^20.3 cm^-2) selected from a set of metal-strong absorbers in the Sloan Digital Sky Survey quasar database. After applying estimate corrections for photoionization effects, we derive gas-phase metallicities of [M/H]=+0.7 +/- 0.2 dex for the SLLS at z=1.7749 toward SDSS0927+5621 and [M/H]=+0.05 +/- 0.1 dex for the SLLS at z=1.7678 toward SDSS0953+5230. The former exhibits among the highest gas metallicity of any astrophysical environment and its total metal surface density exceeds that of nearly every known damped Lya system. The properties of these absorbers -- high metallicity and large velocity width (> 300 km/s) -- resemble those of gas observed in absorption in the spectra of bright, star-forming galaxies at high redshift. We discuss the metal mass density of the SLLS based on these observations and our ongoing SLLS survey and argue that a conservative estimate to the total metal budget at z=2 is greater than 15% of the total, suggesting that the metal-rich LLS may represent the dominant metal reservoir in the young universe.

We present initial results from a Magellan telescope program to image galaxies that give rise to damped Lyman absorbers (DLAs) at 1.63 < z_DLA < 2.37. Our program differs from previous efforts in that we target quasars with intervening Lyman limit systems (LLS) along the line of sight at redshift z_LLS > 3.5. The higher-redshift LLS is applied as a blocking filter to remove the glare of the background quasar at the rest-frame ultraviolet wavelengths of the foreground galaxy. The complete absence of quasar light offers an unimpeded view along the sightline to the redshift of the LLS, allowing an exhaustive search for the DLA galaxy to the sensitivity limit of the imaging data (at or better than 0.25L*). In both of our pilot fields (PKS2000-330, z_DLA=2.033 and SDSS0322-0558, z_DLA=1.69), we identify an L* galaxy within 5" of the sightline which has optical colors consistent with star-forming galaxies at z~2. We examine the correlation between absorption-line properties and galaxy luminosity and impact distance, and compare the high-redshift galaxy and absorber pairs with those known at z<1.

Using H I absorption alone, we attempt to separate H I absorption lines in quasar spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density Lines), and we discuss the difference in their physical properties. We deblend and fit all H I lines with Voigt profiles, and make an unbiased sample of H I lines covering a wide column density range (12 < log N_HI < 19 cm^-2). To reduce the influence of line blending, we simultaneously fit several Lyman series lines. As a result of a two-point correlation analysis, we found that higher column density H I lines are clustering at $\Delta$ v < 200 km/s, while lower ones at $\Delta$ v < 100 km/s. We define HDLs as H I lines with 15 < log N_HI < 19 cm^-2 and all H I lines within +/- 200 km/s of a line with log N_HI > 15 cm^-2, and LDLs as others with 12 < log N_HI < 15 cm^-2. We found that the HDLs have smaller minimum b-values for a given column density than the LDLs. This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler dense gas near to galaxies.

We present a moderate resolution (~1.15 Angstroms/pixel) survey of 79 quasars obtained using the Kast spectrograph on the Shane 3m telescope at Lick observatory. The spectra span the wavelength range of 3175-5880 Angstroms, and have typical signal to noise of 6-20 in the regions of the spectra showing Lyman alpha forest absorption. The quasars have a mean emission redshift of z=2.17, and nearly all cover the entire Lyman alpha forest between Lyman alpha and Lyman beta. Although the quasars were selected to avoid BAL, two quasars in the survey are BAL, one of which is a new discovery. We list the HI and metal ions observed in a total of 140 absorption systems. We also identify 526 emission lines, and list their observed wavelengths, along with new redshifts of the quasars. We determine the rest wavelengths of 3 emission lines or line blends in the forest to be 1070.95 +/- 1.00, 1123.13 +/- 0.51, and 1175.88 +/- 0.30 Angstroms.

We identify a concordant model for the intergalactic medium (IGM) at redshift z=1.9 that uses popular values for cosmological and astrophysical parameters and accounts for all baryons with an uncertainty of 6%. We have measured the amount of absorption, DA, in the Ly-alpha forest at redshift 1.9 in spectra of 77 QSO from the Kast spectrograph. We calibrated the continuum fits with realistic artificial spectra, and we found that averaged over all 77 QSOs the mean continuum level is within 1-2% of the correct value. Absorption from all lines in the Ly-alpha forest at z=1.9 removes DA=15.1 +/- 0.7% of the flux between 1070 and 1170 (rest) Angstroms. This is the first measurement using many QSOs at this z, and the first calibrated measurement at any redshift. Metal lines absorb 2.3 +/- 0.5%, and LLS absorb 1.0 +/- 0.4% leaving 11.8 +/- 1.0% from the lower density bulk of the IGM. Averaging over Delta z=0.1 or 154 Mpc, the dispersion is 6.1 +/- 0.3% including LLS and metal lines, or 3.9 (+0.5, -0.7)% for the lower density IGM alone, consistent with the usual description of large scale structure. LLS and metal lines are major contributors to the variation in the mean flux, and they make the flux field significantly non-Gaussian. We find that a hydrodynamic simulation on a 1024 cubed grid in a 75.7 Mpc box reproduces the observed DA from the low density IGM with parameters values H_o=71 km/s/Mpc, Omega_Lambda=0.73, Omega_m=0.27, Omega_b=0.044, sigma_8=0.9 and a UV background that has an ionization rate that is 1.08 +/- 0.08 times the prediction by Madau, Haardt & Rees (1999).

We present a method to make predictions with sets of correlated data values, in this case QSO flux spectra. We predict the continuum in the Lyman-Alpha forest of a QSO, from 1020 -- 1216 A, using the spectrum of that QSO from 1216 -- 1600 A . We find correlations between the unabsorbed flux in these two wavelengths regions in the HST spectra of 50 QSOs. We use principal component analysis (PCA) to summarize the variety of these spectra and we relate the weights of the principal components for 1020 -- 1600 A to the weights for 1216 -- 1600 A, and we apply this relation to make predictions. We test the method on the HST spectra, and we find an average absolute flux error of 9%, with a range 3 -- 30%, where individual predictions are systematically too low or too high. We mention several ways in which the predictions might be improved.

We describe a new method to calibrate the relative flux levels in spectra from the HIRES echelle spectrograph on the Keck-I telescope. Standard data reduction techniques that transfer the instrument response between HIRES integrations leave errors in the flux of 5 - 10%, because the effective response varies. The flux errors are most severe near the ends of each spectral order, where there can be discontinuous jumps. The source of these errors is uncertain, but may include changes in the vignetting connected to the optical alignment. Our new flux calibration method uses a calibrated reference spectrum of each target to calibrate individual HIRES integrations. We determine the instrument response independently for each integration, and hence we avoid the need to transfer the instrument response between HIRES integrations. The procedure can be applied to any HIRES spectrum, or any other spectrum. While the accuracy of the method depends upon many factors, we have been able to flux calibrate a HIRES spectrum to 1% over scales of 200 A that include order joins. We illustrate the method with spectra of Q1243+3047 towards which we have measured the deuterium to hydrogen abundance ratio.

We report the detection of Deuterium absorption at redshift 2.525659 towards Q1243+3047. We describe improved methods to estimate the Deuterium to Hydrogen abundance ratio (D/H) in absorption systems, including improved modeling of the continuum level, the Ly-alpha forest and the velocity structure of the absorption. Together with improved relative flux calibration, these methods give D/H = 2.42^+0.35_-0.25 x 10^-5 cm^-2 from our Keck-I HIRES spectra of Q1243+3047, where the error is from the uncertainty in the shape of the continuum level and the amount of D absorption in a minor second component. The measured D/H is likely the primordial value because the [O/H] = -2.79 +/- 0.05. This absorption system has a neutral Hydrogen column density of 19.73 +/- 0.04 cm^-2, it shows five D lines and is mostly ionized. The best estimate of the primordial D/H is 2.78^+0.44_-0.38 x 10^-5, from the log D/H values towards five QSOs. The dispersion in the five values is larger than we expect from their individual measurement errors and we suspect this is because some of these errors were underestimated. We observe a trend in D/H with neutral H column density that we also suspect is spurious. The D/H corresponds to a baryon-to-photon ratio ETA = 5.9 +/- 0.5 x 10^-10 and a cosmological baryon density Omega_b h^2 = 0.0214 +/- 0.0020 (9.3%) that agrees with values from the pre-MAP measurements of the anisotropy of the Cosmic Microwave Background.

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.

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.

The z ~ 0.701 absorption system towards the QSO PG1718+4807 is the only example of a QSO absorption system which might have a deuterium/hydrogen ratio approximately ten times the value found towards other QSOs. We have obtained new STIS spectra from the Hubble Space Telescope of the Lyman alpha and Lyman limit regions of the system. These spectra give the redshift and velocity dispersion of the neutral hydrogen which produces most of the observed absorption. The Lyman alpha line is too narrow to account for all of the observed absorption. It was previously known that extra absorption is needed on the blue side of the main H I near to the expected position of deuterium. The current data suggests with a 98% confidence level that the extra absorption is not deuterium. Some uncertainty persists because we have a low signal to noise ratio and the extra absorption - be it deuterium or hydrogen - is heavily blended with the Lyman alpha absorption from the main hydrogen component.

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.

Big Bang Nucleosynthesis (BBN) is the synthesis of the light nuclei, Deuterium, He3, He4 and Li7, during the first few minutes of the universe. This review concentrates on recent improvements in the measurement of the primordial (after BBN, and prior to modification) abundances of these nuclei. We mention improvement in the standard theory, and the non-standard extensions which are limited by the data. (abridged)

We have discovered a third quasar absorption system which is consistent with a low deuterium to hydrogen abundance ratio, D/H = 3.4 times 10^-5. The z ~ 2.8 partial Lyman limit system towards QSO 0130-4021 provides the strongest evidence to date against large D/H ratios because the H I absorption, which consists of a single high column density component with unsaturated high order Lyman series lines, is readily modeled -- a task which is more complex in other D/H systems. We have obtained twenty-two hours of spectra from the HIRES spectrograph on the W.M. Keck telescope, which allow a detailed description of the Hydrogen. We see excess absorption on the blue wing of the H I Lyman alpha line, near the expected position of Deuterium. However, we find that Deuterium cannot explain all of the excess absorption, and hence there must be contamination by additional absorption, probably H I. This extra H I can account for most or all of the absorption at the D position, and hence D/H = 0 is allowed. We find an upper limit of D/H < 6.7 times 10^-5 in this system, consistent with the value of D/H ~ 3.4 times 10^-5 deduced towards QSO 1009+2956 and QSO 1937-1009 by Burles and Tytler (1998a, 1998b). This absorption system shows only weak metal line absorption, and we estimate [Si/H] < -2.6 -- indicating that the D/H ratio of the system is likely primordial. All four of the known high redshift absorption line systems simple enough to provide useful limits on D are consistent with D/H = 3.4 +/- 0.25 times 10^-5. Conversely, this QSO provides the third case which is inconsistent with much larger values.