Source author record

J. Xavier Prochaska

J. Xavier Prochaska appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

170works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

170 published item(s)

preprint2026arXiv

Featureless stars: Flux Calibration for Extremely Large Telescopes

The spectrophotometric flux calibration of recent spectroscopic surveys has reached a limiting systematic precision of approximately 1-3 percent, and is often biased near the wavelengths associated with H I Balmer absorption. As we prepare for the next generation of imaging and spectroscopic surveys, and high-precision cosmology experiments, we must find a way to address this systematic. Towards this goal, we have identified a global network of 29 bright (G < 17.5) featureless white dwarf stars that have a spectral energy distribution consistent with an almost pure blackbody form over the entire optical and near-infrared wavelength range. Based on this sample, we have computed the systematic uncertainty and AB magnitude offsets associated with Gaia, SDSS, SMSS, PanSTARRS, DES, and 2MASS, and we have also checked the consistency of our objects with both GALEX and WISE. The magnitude range of the featureless stars reported here are ideally suited to observations taken with the forthcoming generation of extremely large telescopes, as well as calibrating the survey data acquired by the Rubin, Euclid and Roman observatories. Finally, all of the high-precision spectrophotometric standard stars reported here have been included in the latest release of the PypeIt data reduction pipeline.

preprint2026arXiv

The Second CHIME/FRB Catalog of Fast Radio Bursts

We present a catalog of 4539 fast radio bursts (FRBs) observed with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope between 25 July 2018 and 15 September 2023. These bursts originate from 3641 unique sources, including 981 bursts from 83 known repeating sources. For each FRB, the catalog provides a $O(10')$ estimate of sky location along with corresponding measurements of cumulative exposure time and survey sensitivity over the observing period. It includes a total-intensity dynamic spectrum between 400 and 800 MHz at 0.983 ms resolution. From this spectrum, we constrain a model of the burst morphology and measure key parameters such as arrival time, intrinsic temporal width, dispersion measure, scattering time, and flux density. This second catalog includes all FRBs from the first catalog, with every event reprocessed using a uniform and improved analysis framework. We show that previously published inferences remain valid under the updated measurements. We assess consistency of the detection rate across observational parameters, present initial distributions of burst properties, and outline ongoing and future studies that will use this catalog to investigate the nature of FRBs and their utility as astrophysical and cosmological probes.

preprint2025arXiv

Localisation and host galaxy identification of new Fast Radio Bursts with MeerKAT

Accurately localising fast radio bursts (FRBs) is essential for understanding their birth environments and for their use as cosmological probes. Recent advances in radio interferometry, particularly with MeerKAT, have enabled the localisation of individual bursts with arcsecond precision. In this work, we present the localisation of 15 apparently non-repeating FRBs detected with MeerKAT. Two of the FRBs, discovered in 2022, were localised in 8 second images from the projects which MeerTRAP was commensal to, while eight were localised using the transient buffer (TB) pipeline, and another one through SeeKAT, all with arcsecond precision. Four additional FRBs lacked TB triggers and sufficient signal, limiting their localisation only to arcminute precision. For eight of the FRBs in our sample, we identify host galaxies with greater than 90% confidence, and one with 80% confidence, while two FRBs have ambiguous associations. We measured spectroscopic redshifts for six host galaxies, ranging from 0.33 to 0.85, demonstrating MeerKAT's sensitivity to high redshift FRBs. We modelled the spectral energy distributions of host galaxies with sufficient photometric coverage to derive their stellar population and star formation properties. This work represents one of the largest uniform samples of well-localised distant FRBs to date, laying the groundwork for using MeerKAT FRBs as cosmological probes and understand how FRB hosts evolve at high redshift.

preprint2023arXiv

SDSS DR17: The Cosmic Slime Value Added Catalog

The "cosmic web", the filamentary large-scale structure in a cold dark matter Universe, is readily apparent via galaxy tracers in spectroscopic surveys. However, the underlying dark matter structure is as of yet unobservable and mapping the diffuse gas permeating it lies beyond practical observational capabilities. A recently developed technique, inspired by the growth and movement of Physarum polycephalum "slime mold", has been used to map the cosmic web of a low redshift sub-sample of the SDSS spectroscopic galaxy catalog. This model, the Monte Carlo Physarum Machine (MCPM) was shown to promisingly reconstruct the cosmic web. Here, we improve the formalism used in calibrating the MCPM to better recreate the Bolshoi-Planck cosmological simulation's density distributions and apply them to a significantly larger cosmological volume than previous works using the Sloan Digital Sky Survey (SDSS, $z < 0.1$) and the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy (LRG, $z \lesssim 0.5$) spectroscopic catalogs. We present the "Cosmic Slime Value Added Catalog" which provides estimates for the cosmic overdensity for the sample of galaxies probed spectroscopically by the above SDSS surveys. In addition, we provide the fully reconstructed 3D density cubes of these volumes. These data products were released as part of Sloan Digital Sky Survey Data Release 17 and are publicly available. We present the input catalogs and the methodology for constructing these data products. We also highlight exciting potential applications to galaxy evolution, cosmology, the intergalactic and circumgalactic medium, and transient phenomenon localization.

preprint2023arXiv

Second Data Release of the COSMOS Lyman-alpha Mapping and Tomographic Observation: The First 3D Maps of the Detailed Cosmic Web at 2.05<z<2.55

We present the second data release of the COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) Survey conducted with the LRIS spectrograph on the Keck-I telescope. This project used Lyman-alpha forest absorption in the spectra of faint star forming galaxies and quasars at z ~ 2-3 to trace neutral hydrogen in the intergalactic medium. In particular, we use 320 objects over a footprint of ~0.2 deg^2 to reconstruct the absorption field at 2.05 < z < 2.55 at ~2 h^{-1}Mpc resolution. We apply a Wiener filtering technique to the observed data to reconstruct three dimensional maps of the field over a volume of 4.1 x 10^5 comoving cubic Mpc. In addition to the filtered flux maps, for the first time we infer the underlying dark matter field through a forward modeling framework from a joint likelihood of galaxy and Lyman-alpha forest data, finding clear examples of the detailed cosmic web consisting of cosmic voids, sheets, filaments, and nodes. In addition to traditional figures, we present a number of interactive three dimensional models to allow exploration of the data and qualitative comparisons to known galaxy surveys. We find that our inferred over-densities are consistent with those found from galaxy fields. Our reduced spectra, extracted Lyman-alpha forest pixel data, and reconstructed tomographic maps are available publicly at https://doi.org/10.5281/zenodo.7524313

preprint2022arXiv

A fast radio burst progenitor born in a galaxy merger

We report a Giant Metrewave Radio Telescope HI 21cm mapping study of the neutral atomic hydrogen (HI) in the host galaxy of the fast radio burst (FRB) FRB20180916B at $z \approx 0.03399$. We find that the FRB host has an HI mass of $\rm M_{HI} = (2.74 \pm 0.33) \times 10^9 \ M_\odot$ and a high HI-to-stellar mass ratio, $\approx 1.3$. The FRB host is thus a gas-rich but near-quiescent galaxy, that is likely to have acquired a significant mass of HI in the recent past. The HI distribution is disturbed, with extended HI 21cm emission detected in a north-eastern tail, a counter-tail towards the south, an HI hole between the galaxy centre and the FRB location, and a high HI column density measured close to the FRB position. The FRB host is part of a group with four companions detected in their HI 21cm emission, the nearest of which is only 22~kpc from the FRB location. The gas-richness and disturbed HI distribution indicate that the FRB host has recently undergone a minor merger, which increased its HI mass, disturbed the HI in the galaxy disk, and compressed the HI near the FRB location to increase its surface density. We propose that this merger caused the burst of star-formation in the outskirts of the galaxy that gave rise to the FRB progenitor. The evidence for a minor merger is consistent with scenarios in which the FRB progenitor is a massive star, formed due to the merger event.

preprint2022arXiv

A Massive, Dusty HI-Absorption-Selected Galaxy at $z \approx 2.46$ Identified in a CO Emission Survey

We report a NOrthern Extended Millimeter Array (NOEMA) and Atacama Large Millimeter/submillimeter Array (ALMA) search for redshifted CO emission from the galaxies associated with seven high-metallicity ([M/H] $\geq -1.03$) damped Lyman-$α$ absorbers (DLAs) at $z\approx1.64-2.51$. Our observations yielded one new detection of CO(3-2) emission from a galaxy at $z=2.4604$ using NOEMA, associated with the $z=2.4628$ DLA towards QSO B0201+365. Including previous searches, our search results in detection rates of CO emission of $\approx56^{+38}_{-24}$ % and $\approx11^{+26}_{-9}$ %, respectively, in the fields of DLAs with ${\rm [M/H]}>-0.3$ and ${\rm [M/H]}<-0.3$. Further, the HI-selected galaxies associated with five DLAs with [M/H] $>-0.3$ all have high molecular gas masses, $\gtrsim5\times10^{10}\ {\rm M}_\odot$. This indicates that the highest-metallicity DLAs at $z\approx2$ are associated with the most massive galaxies. The newly-identified $z\approx2.4604$ HI-selected galaxy, DLA0201+365g, has an impact parameter of $\approx7$ kpc to the QSO sightline, and an implied molecular gas mass of $(5.04\pm0.78) \times10^{10}\times(α_{\rm CO}/4.36)\times(r_{31}/0.55)\ {\rm M}_\odot$. Archival Hubble Space Telescope Wide Field and Planetary Camera 2 imaging covering the rest-frame near-ultraviolet (NUV) and far-ultraviolet (FUV) emission from this galaxy yield non-detections of rest-frame NUV and FUV emission, and a $5σ$ upper limit of 2.3 M$_\odot$ yr$^{-1}$ on the unobscured star formation rate (SFR). The low NUV-based SFR estimate, despite the very high molecular gas mass, indicates that DLA0201+365g either is a very dusty galaxy, or has a molecular gas depletion time that is around two orders of magnitude larger than that of star-forming galaxies at similar redshifts.

preprint2022arXiv

Constraining the Cosmic Baryon Distribution with Fast Radio Burst Foreground Mapping

The dispersion measures (DM) of fast radio bursts (FRBs) encode the integrated electron density along the line-of-sight, which is typically dominated by the intergalactic medium (IGM) contribution in the case of extragalactic FRBs. In this paper, we show that incorporating wide-field spectroscopic galaxy survey data in the foreground of localized FRBs can significantly improve constraints on the partition of diffuse cosmic baryons. Using mock DMs and realistic lightcone galaxy catalogs derived from the Millennium simulation, we define spectroscopic surveys that can be carried out with 4m and 8m-class wide field spectroscopic facilities. On these simulated surveys, we carry out Bayesian density reconstructions in order to estimate the foreground matter density field. In comparison with the `true' matter density field, we show that these can help reduce the uncertainties in the foreground structures by $\sim 2-3\times$ compared to cosmic variance. We calculate the Fisher matrix to forecast that $N=30\: (96)$ localized FRBs should be able to constrain the diffuse cosmic baryon fraction to $\sim 10\%\: (\sim 5\%) $, and parameters governing the size and baryon fraction of galaxy circumgalactic halos to within $\sim 20-25\%\: (\sim 8-12\%)$. From the Fisher analysis, we show that the foreground data increases the sensitivity of localized FRBs toward our parameters of interest by $\sim 25\times$. We briefly introduce FLIMFLAM, an ongoing galaxy redshift survey that aims to obtain foreground data on $\sim 30$ localized FRB fields.

preprint2022arXiv

Deep Learning of DESI Mock Spectra to Find Damped Lyα Systems

We have updated and applied a convolutional neural network (CNN) machine learning model to discover and characterize damped Ly$α$ systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99$\%$ for spectra which have signal-to-noise (S/N) above 5 per pixel. Classification accuracy is the rate of correct classifications. This accuracy remains above 97$\%$ for lower signal-to-noise (S/N) $\approx1$ spectra. This CNN model provides estimations for redshift and HI column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 per pixel. Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of Baryon Acoustic Oscillation (BAO) is investigated. The cosmological fitting parameter result for BAO has less than $0.61\%$ difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above $1.7\%$. We also compared the performance of CNN and Gaussian Process (GP) model. Our improved CNN model has moderately 14$\%$ higher purity and 7$\%$ higher completeness than an older version of GP code, for S/N $>$ 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by $24\%$ less standard deviation. A credible DLA catalog for DESI main survey can be provided by combining these two algorithms.

preprint2022arXiv

First discoveries and localisations of Fast Radio Bursts with MeerTRAP: a real-time, commensal MeerKAT survey

We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1'$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy's measured properties are consistent with other FRB hosts with secure associations.

preprint2022arXiv

HeII Ly$α$ Transmission Spikes and Absorption Troughs in Eight High-resolution Spectra Probing the End of HeII Reionization

We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=λ/Δλ\simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the HeII reionization epoch. The measured lengths of the troughs decrease dramatically from L>100cMpc at z>3 to L~30cMpc at z~2.7, signaling a significant progression of HeII reionization at these redshifts. Furthermore, unexpectedly long L~65cMpc troughs detected at z~2.9 suggest that the UV background fluctuates at larger scales than predicted by current models. By comparing the measured incidence of transmission spikes to predictions from forward-modeled mock spectra created from the outputs of a (146cMpc)^3 optically thin Nyx hydrodynamical simulation employing different UV background models, we infer the redshift evolution of the HeII photoionization rate $Γ_\mathrm{He\,II}(z)$. The photoionization rate decreases with increasing redshift from $\simeq 4.6\times 10^{-15}\mathrm{\,s^{-1}}$ at z~2.6 to $\simeq 1.2 \times 10^{-15}\mathrm{\,s^{-1}}$ at z~3.2, in agreement with previous inferences from the HeII effective optical depth, and following expected trends of current models of a fluctuating HeII-ionizing background.

preprint2022arXiv

Jansky Very Large Array detections of CO(1-0) emission in HI-absorption-selected galaxies at $z \gtrsim 2$

We report a Karl G. Jansky Very Large Array search for redshifted CO(1-0) emission from three HI-absorption-selected galaxies at $z \approx 2$, identified earlier in their CO(3-2) or CO(4-3) emission. We detect CO(1-0) emission from DLA B1228-113 at $z\approx2.1933$ and DLA J0918+1636 at $z\approx2.5848$; these are the first detections of CO(1-0) emission in high-$z$ HI-selected galaxies. We obtain high molecular gas masses, $\rm M_{mol}\approx10^{11}\times(α_{\rm CO}/4.36)\ M_\odot$, for the two objects with CO(1-0) detections, which are a factor of $\approx1.5-2$ lower than earlier estimates. We determine the excitation of the mid$-J$ CO rotational levels relative to the $J=1$ level, r$_{ J1}$, in HI-selected galaxies for the first time, obtaining r$_{\rm 31}=1.00\pm0.20$ and r$_{\rm 41}=1.03\pm0.23$ for DLA J0918+1636, and r$_{\rm 31}=0.86\pm0.21$ for DLA B1228-113. These values are consistent with thermal excitation of the $J=3,4$ levels. The excitation of the $J=3$ level in the HI-selected galaxies is similar to that seen in massive main-sequence and sub-mm galaxies at $z\gtrsim2$, but higher than that in main-sequence galaxies at $z\approx1.5$; the higher excitation of the galaxies at $z\gtrsim2$ is likely to be due to their higher star-formation rate (SFR) surface density. We use Hubble Space Telescope Wide Field Camera 3 imaging to detect the rest-frame near-ultraviolet emission of DLA B1228-113, obtaining an NUV SFR of $4.44\pm0.47$ M$_{\odot}$ yr$^{-1}$, significantly lower than that obtained from the total infrared luminosity, indicating significant dust extinction in the $z\approx2.1933$ galaxy.

preprint2022arXiv

Monte Carlo Physarum Machine: Characteristics of Pattern Formation in Continuous Stochastic Transport Networks

We present Monte Carlo Physarum Machine: a computational model suitable for reconstructing continuous transport networks from sparse 2D and 3D data. MCPM is a probabilistic generalization of Jones's 2010 agent-based model for simulating the growth of Physarum polycephalum slime mold. We compare MCPM to Jones's work on theoretical grounds, and describe a task-specific variant designed for reconstructing the large-scale distribution of gas and dark matter in the Universe known as the Cosmic web. To analyze the new model, we first explore MCPM's self-patterning behavior, showing a wide range of continuous network-like morphologies -- called "polyphorms" -- that the model produces from geometrically intuitive parameters. Applying MCPM to both simulated and observational cosmological datasets, we then evaluate its ability to produce consistent 3D density maps of the Cosmic web. Finally, we examine other possible tasks where MCPM could be useful, along with several examples of fitting to domain-specific data as proofs of concept.

preprint2021arXiv

A Long Stream of Metal-Poor Cool Gas around a Massive Starburst Galaxy at z = 2.67

We present the first detailed dissection of the circumgalactic medium (CGM) of massive starburst galaxies at z > 2. Our target is a submillimeter galaxy (SMG) at z = 2.674 that has a star formation rate of 1200 $M_\odot$/yr and a molecular gas reservoir of $1.3\times10^{11} M_\odot$. We characterize its CGM with two background QSOs at impact parameters of 93 kpc and 176 kpc. We detect strong HI and metal-line absorption near the redshift of the SMG towards both QSOs, each consisting of three main subsystems spanning over 1500 km/s. The absorbers show remarkable kinematic and metallicity coherence across a separation of 86 kpc. In particular, the cool gas in the CGM of the SMG exhibits high HI column densities ($\log N_{\rm HI}/{\rm cm}^{-2} = 20.2, 18.6$), low metallicities ([M/H] $\approx$ -2.0), and similar radial velocities ($\approx$ -300 km/s). While the HI column densities match previous results on the CGM around QSOs at z > 2, the metallicities are lower by more than an order of magnitude, making it an outlier in the line width$-$metallicity relation of damped Ly$α$ absorbers. The large physical extent, the velocity coherence, the high surface density, and the low metallicity are all consistent with the cool, inflowing, and near-pristine gas streams predicted to penetrate hot massive halos at z > 1.5. We estimate a total gas accretion rate of ~100 $M_\odot$/yr from three such streams, which falls short of the star formation rate but is consistent with simulations. At this rate, it takes about a gigayear to acquire the molecular gas reservoir of the central starburst.

preprint2021arXiv

A Multiwavelength Study of ELAN Environments (AMUSE$^2$): Detection of a dusty star-forming galaxy within the enormous Lyman $α$ nebula at $z=2.3$ sheds light on its origin

We present ALMA observations on and around the radio-quiet quasar UM287 at $z=2.28$. Together with a companion quasar, UM287 is believed to play a major role in powering the surrounding enormous Ly$α$ nebula (ELAN), dubbed the Slug ELAN, that has an end-to-end size of 450 physical kpc. In addition to the quasars, we detect a new dusty star-forming galaxy (DSFG), dubbed the Slug-DSFG, in 2 mm continuum with a single emission line consistent with CO(4-3). The Slug-DSFG sits at a projected distance of 100 kpc south-east from UM287, with a systemic velocity difference of $-360\pm30$ km s$^{-1}$ with respect to UM287, suggesting it being a possible contributor to the powering of the Slug ELAN. With careful modeling of SED and dynamical analyses it is found that the Slug-DSFG and UM287 appear low in both gas fraction and gas-to-dust ratio, suggesting environmental effects due to the host massive halo. In addition, our Keck long-slit spectra reveal significant Ly$α$ emissions from the Slug-DSFG, as well as a Ly$α$ tail that starts at the location and velocity of the Slug-DSFG and extends towards the south, with a projected length of about 100 kpc. Supported by various analytical estimates we propose that the Ly$α$ tail is a result of the Slug-DSFG experiencing ram pressure stripping. The gas mass stripped is estimated to be about 10$^9$ M$_\odot$, contributing to the dense warm/cool gas reservoir that is believed to help power the exceptional Ly$α$ luminosity.

preprint2021arXiv

A Multiwavelength Study of ELAN Environments (AMUSE$^2$): Ubiquitous dusty star-forming galaxies associated with enormous Ly$α$ nebulae on megaparsec scales

We have been undertaking a systematic survey at 850 $μ$m based on a sample of four prototypical $z\sim2-3$ enormous Ly$α$ nebulae (ELANe) as well as their megaparsec-scale (Mpc-scale) environments to study the physical connections between ELANe and their coeval dusty submillimeter galaxies (SMGs). By analysing the SCUBA-2 data with self-consistent Monte Carlo simulations to construct the number counts, here, we report on the overabundance of 850 $μ$m-selected submillimeter sources around all the four ELANe, by a factor of 3.6$\pm$0.6 (weighted average) compared to the blank fields. This suggests that the excessive number of submillimeter sources are likely to be part of the Mpc-scale environment around the ELANe, corroborating the co-evolution scenario for SMGs and quasars; this is a process which may be more commonly observed in the ELAN fields. If the current form of the underlying count models continues toward the fainter end, our results would suggest an excess of the 850 $μ$m extragalactic background light by a factor of between 2-10, an indication of significant background light fluctuations on the survey scales. Finally, by assuming that all the excessive submillimeter sources are associated with their corresponding ELAN environments, we estimate the SFR densities of each ELAN field, as well as a weighted average of $Σ$SFR=1200$\pm$300 $M_\odot$ yr$^{-1}$ Mpc$^{-3}$, consistent with that found in the vicinity of other quasar systems or proto-clusters at similar redshifts; in addition, it is a factor of about 300 greater than the cosmic mean.

preprint2021arXiv

A Multiwavelength Study of ELAN Environments (AMUSE$^2$). Mass budget, satellites spin alignment and gas infall in a massive $z\sim3$ quasar host halo

The systematic targeting of extended Ly$α$ emission around high-redshift quasars resulted in the discovery of rare and bright Enormous Ly$α$ Nebulae (ELANe) associated with multiple active galactic nuclei (AGN). We here initiate "a multiwavelength study of ELAN environments" (AMUSE$^2$) focusing on the ELAN around the $z\sim3$ quasar SDSS J1040+1020, a.k.a. the Fabulous ELAN. We report on VLT/HAWK-I, APEX/LABOCA, JCMT/SCUBA-2, SMA/850$μ$m, ALMA/CO(5-4) and 2mm observations and compare them to previously published VLT/MUSE data. The continuum and line detections enable a first estimate of the star-formation rates, dust, stellar and molecular gas masses in four objects associated with the ELAN (three AGNs and one Ly$α$ emitter), confirming that the quasar host is the most star-forming (${\rm SFR}\sim500$ M$_\odot$ yr$^{-1}$) and massive galaxy ($M_{\rm star}\sim10^{11}$ M$_{\odot}$) in the system, and thus can be assumed as central. All four embedded objects have similar molecular gas reservoirs ($M_{\rm H_2}\sim10^{10}$ M$_{\odot}$), resulting in short depletion time scales. This fact together with the estimated total dark-matter halo mass, $M_{\rm DM}=(0.8-2)\times10^{13}$ M$_{\odot}$, implies that this ELAN will evolve into a giant elliptical galaxy. Consistently, the constraint on the baryonic mass budget for the whole system indicates that the majority of baryons should reside in a massive warm-hot reservoir (up to $10^{12}$ M$_{\odot}$), needed to complete the baryons count. Additionally, we discuss signatures of gas infall on the compact objects as traced by Ly$α$ radiative transfer effects and the evidence for the alignment between the satellites' spins and their directions to the central.

preprint2021arXiv

Characterizing the FRB host galaxy population and its connection to transients in the local and extragalactic Universe

We present the localization and host galaxies of one repeating and two apparently non-repeating Fast Radio Bursts. FRB20180301A was detected and localized with the Karl G. Jansky Very Large Array to a star-forming galaxy at $z=0.3304$. FRB20191228A, and FRB20200906A were detected and localized by the Australian Square Kilometre Array Pathfinder to host galaxies at $z=0.2430$ and $z=0.3688$, respectively. We combine these with 13 other well-localized FRBs in the literature, and analyze the host galaxy properties. We find no significant differences in the host properties of repeating and apparently non-repeating FRBs. FRB hosts are moderately star-forming, with masses slightly offset from the star-forming main-sequence. Star formation and low-ionization nuclear emission-line region (LINER) emission are major sources of ionization in FRB host galaxies, with the former dominant in repeating FRB hosts. FRB hosts do not track stellar mass and star formation as seen in field galaxies (more than 95% confidence). FRBs are rare in massive red galaxies, suggesting that progenitor formation channels are not solely dominated by delayed channels which lag star formation by Gigayears. The global properties of FRB hosts are indistinguishable from core-collapse supernovae (CCSNe) and short gamma-ray bursts (SGRBs) hosts, and the spatial offset (from galaxy centers) of FRBs is mostly inconsistent with that of the Galactic neutron star population (95% confidence). The spatial offsets of FRBs (normalized to the galaxy effective radius) also differ from those of globular clusters (GCs) in late- and early-type galaxies with 95% confidence.

preprint2021arXiv

On the Environments of Giant Radio Galaxies

We test the hypothesis that environments play a key role in enabling the growth of enormous radio structures spanning more than 700 kpc, an extreme population of radio galaxies called giant radio galaxies (GRGs). To achieve this, we explore (1) the relationships between the occurrence of GRGs and the surface number density of surrounding galaxies, including satellite galaxies and galaxies from neighboring halos, as well as (2) the GRG locations towards large-scale structures. The analysis is done by making use of a homogeneous sample of 110 GRGs detected from the LOFAR Two-metre Sky Survey in combination with photometric galaxies from the DESI Legacy Imaging Surveys and a large-scale filament catalog from the Sloan Digital Sky Survey. Our results show that the properties of galaxies around GRGs are similar with that around the two control samples, consisting of galaxies with optical colors and luminosity matched to the properties of the GRG host galaxies. Additionally, the properties of surrounding galaxies depend on neither their relative positions to the radio jet/lobe structures nor the sizes of GRGs. We also find that the locations of GRGs and the control samples with respect to the nearby large-scale structures are consistent with each other. These results demonstrate that there is no correlation between the GRG properties and their environments traced by stars, indicating that external galaxy environments are not the primary cause of the large sizes of the radio structures. Finally, regarding radio feedback, we show that the fraction of blue satellites does not correlate with the GRG properties, suggesting that the current epoch of radio jets have minimal influence on the nature of their surrounding galaxies.

preprint2020arXiv

3D Distribution Map of HI Gas and Galaxies Around an Enormous Ly$α$ Nebula and Three QSOs at $z=2.3$ Revealed by the HI Tomographic Mapping Technique

We present an IGM HI tomography map in a survey volume of $16 \times 19 \times 131 \ h^{-3} {\rm comoving \ Mpc}^{3}$ (cMpc$^3$) centered at MAMMOTH-1 nebula and three neighbouring quasars at $z=2.3$. MAMMOTH-1 nebula is an enormous Ly$α$ nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over $1\ h^{-1}$ cMpc with not fully clear physical origin. Here we investigate the HI-gas distribution around MAMMOTH1-QSO with the ELAN and three neighbouring type-I quasars, making the IGM HI tomography map with a spatial resolution of $2.6\ h^{-1}$ cMpc. Our HI tomography map is reconstructed with HI Ly$α$ forest absorption of bright background objects at $z=2.4-2.9$: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of HI flux overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with significantly weak HI absorption. This suggests that MAMMOTH1-QSO has a proximity zone where quasar illuminates and photo-ionizes the surrounding HI gas and suppresses HI absorption, and that the ELAN is probably a photo-ionized cloud embedded in the cosmic web. The HI radial profile of MAMMOTH1-QSO is very similar to those of three neighbouring type-I quasars at $z=2.3$, which is compatible with the AGN unification model. We compare the distributions of the HI absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and HI gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

preprint2020arXiv

A Cold, Massive, Rotating Disk Galaxy 1.5 Billion Years after the Big Bang

Massive disk galaxies like the Milky Way are expected to form at late times in traditional models of galaxy formation, but recent numerical simulations suggest that such galaxies could form as early as a billion years after the Big Bang through the accretion of cold material and mergers. Observationally, it has been difficult to identify disk galaxies in emission at high redshift, in order to discern between competing models of galaxy formation. Here we report imaging, with a resolution of about 1.3 kiloparsecs, of the 158-micrometre emission line from singly ionized carbon, the far-infrared dust continuum and the near-ultraviolet continuum emission from a galaxy at a redshift of 4.2603, identified by detecting its absorption of quasar light. These observations show that the emission arises from gas inside a cold, dusty, rotating disk with a rotational velocity of 272 kilometres per second. The detection of emission from carbon monoxide in the galaxy yields a molecular mass that is consistent with the estimate from the ionized carbon emission of about 72 billion solar masses. The existence of such a massive, rotationally supported, cold disk galaxy when the Universe was only 1.5 billion years old favours formation through either cold-mode accretion or mergers, although its large rotational velocity and large content of cold gas remain challenging to reproduce with most numerical simulations.

preprint2020arXiv

A Data-Driven Technique Using Millisecond Transients to Measure the Milky Way Halo

We introduce a new technique to constrain the line-of-sight integrated electron density of our Galactic halo $\text{DM}_\text{MW,halo}$ through analysis of the observed dispersion measure distributions of pulsars $\text{DM}_\text{pulsar}$ and fast radio bursts $\text{DM}_\text{FRB}$. We model these distributions, correcting for the Galactic interstellar medium, with kernel density estimation---well-suited to the small data regime---to find lower/upper bounds to the corrected $\text{DM}_\text{pulsar}$/$\text{DM}_\text{FRB}$ distributions: $\max[\text{DM}_\text{pulsar}] \approx 7\pm2 \text{ (stat)} \pm 9 \text{ (sys) pc cm}^{-3}$ and $\min[\text{DM}_\text{FRB}] \approx 63^{+27}_{-21} \text{ (stat)} \pm 9 \text{ (sys) pc cm}^{-3}$. Using bootstrap resampling to estimate uncertainties, we set conservative limits on the Galactic halo dispersion measure $-2 < \text{DM}_\text{MW,halo} < 123 \text{pc cm}^{-3}$ (95\% c.l.). The upper limit is especially conservative because it may include a non-negligible contribution from the FRB host galaxies and a non-zero contribution from the cosmic web. It strongly disfavors models where the Galaxy has retained the majority of its baryons with a density profile tracking the presumed dark matter density profile. Last, we perform Monte Carlo simulations of larger FRB samples to validate our technique and assess the sensitivity of ongoing and future surveys. We recover bounds of several tens $\text{pc cm}^{-3}$ which may be sufficient to test whether the Galaxy has retained a majority of its baryonic mass. We estimate that a sample of several thousand FRBs will significantly tighten constraints on $\text{DM}_\text{MW,halo}$ and offer a valuable complement to other analyses.

preprint2020arXiv

A metal-poor damped Ly-alpha system at redshift 6.4

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

preprint2020arXiv

A search for supernova-like optical counterparts to ASKAP-localised Fast Radio Bursts

Fast radio bursts (FRBs) are millisecond-scale radio pulses, which originate in distant galaxies and are produced by unknown sources. The mystery remains partially because of the typical difficulty in localising FRBs to host galaxies. Accurate localisations delivered by the Commensal Real-time ASKAP Fast Transients (CRAFT) survey now provide an opportunity to study the host galaxies and potential transient counterparts of FRBs at a large range of wavelengths. In this work, we investigate whether the first three FRBs accurately localised by CRAFT have supernova-like transient counterparts. We obtained two sets of imaging epochs with the Very Large Telescope for three host galaxies, one soon after the burst detection and one several months later. After subtracting these images no optical counterparts were identified in the associated FRB host galaxies, so we instead place limits on the brightness of any potential optical transients. A Monte Carlo approach, in which supernova light curves were modelled and their base properties randomised, was used to estimate the probability of a supernova associated with each FRB going undetected. We conclude that Type Ia and IIn supernovae are unlikely to accompany every apparently non-repeating FRB.

preprint2020arXiv

ALMA [CII] 158um imaging of an HI-Selected Major Merger at z~4

We present high spatial-resolution (~2kpc) Atacama Large Millimeter/submillimeter Array (ALMA) observations of [CII] 158um and dust-continuum emission from a galaxy at z=3.7978 selected by its strong HI absorption (a damped Lya absorber, DLA) against a background QSO. Our ALMA images reveal a pair of star-forming galaxies separated by ~6kpc (projected) undergoing a major merger. Between these galaxies is a third emission component with highly elevated (2x) [CII] 158um emission relative to the dust continuum, which is likely to arise from stripped gas associated with the merger. This merger of two otherwise-normal galaxies is not accompanied by enhanced star-formation, contrary to mergers detected in most luminosity-selected samples. The DLA associated with the merger exhibits extreme kinematics, with a velocity width for the low-ionization metal lines of Dv90~470km/s that spans the velocity spread revealed in the [CII] 158um emission. We propose that DLAs with high Dv90 values are a signpost of major mergers in normal galaxies at high redshifts, and use the distribution of the velocity widths of metal lines in high-z DLAs to provide a rough estimate the fraction of z>3 galaxies that are undergoing a major merger.

preprint2020arXiv

Confronting the Magnetar Interpretation of Fast Radio Bursts Through Their Host Galaxy Demographics

We explore the prompt magnetar progenitor scenario in the context of fast radio burst (FRB) host galaxies demographics and offset distributions. Magnetars are neutron stars with strong magnetic fields on the order of $10^{15}$ G with a short decay lifetime of less than $10^4$ years. Due to their extremely short lifetimes, magnetars should follow the demographics of galaxies according to their current star-formation rate (SFR). Moreover, we hypothesize that magnetars should follow the SFR profile within galaxies, which we assume to follow an exponential profile. We construct a simple model for the host galaxies of magnetars assuming these events track SFR in all galaxies and compare it to observed properties from a sample of \nsecure\ secure FRB hosts. We find the distribution of observed SFRs is inconsistent with the model at $>95\%$ c.l. The offset distribution is consistent with this scenario; however, this could be due to the limited sample size and the seeing limited estimates for the effective radii of the FRB host galaxies. Despite the recent association of an FRB with a magnetar in the Milky Way, magnetars may not be the only source of FRBs in the universe, yet any other successful model must account for the demographics of the FRB host in SFR and their observed galactocentric offsets.

preprint2020arXiv

Constraining Magnetic Fields in the Circumgalactic Medium

We study the properties of magnetic fields in the circumgalactic medium (CGM) of $z<1$ galaxies by correlating Faraday rotation measures of $\sim 1,000$ high-redshift radio sources with the foreground galaxy number density estimated from the DESI Legacy Imaging Surveys. This method enables us to extract signals of rotation measures contributed by intervening gas around multiple galaxies. Our results show that there is no detectable correlation between the distribution of rotation measures and the number of foreground galaxies, contrary to several previous results. Utilizing the non-detection signals, we estimate $3σ$ upper limits to the rotation measures from the CGM of $\sim 20 \rm \ rad/m^{2}$ within~50 kpc and $\sim 10 \rm \ rad / m^{2}$ at separations of $100$ kpc. By adopting a column density distribution of ionized gas obtained from absorption line measurements, we further estimate the strengths of coherent magnetic fields parallel to the line of sight of $<\rm 2 \ μG$ in the CGM. We show that the estimated upper limits of rotation measures and magnetic field strengths are sufficient to constrain outputs of recent galaxy magneto-hydrodynamic simulations. Finally, we discuss possible causes for the inconsistency between our results and previous works.

preprint2020arXiv

Deep Hubble Space Telescope Imaging on the Extended Ly$α$ Emission of a QSO at $z = 2.19$ with Damped Lyman Alpha System as a Natural Coronagraph

Recent observations suggest that proximate damped Ly$α$ (PDLA) systems can be used to study the host galaxies of Quasi-stellar objects (QSOs), because the PDLAs can block the bright point-spread-function (PSF) from central QSOs. Using six-orbits of narrowband imaging with $\it{HST}$/WFC3, we present the first high resolution narrowband image of the Ly$α$ emission in the PDLA trough of the QSO SDSSJ115432.67-021537.9 . We detect one major component and one minor component in the narrowband imaging. Combining the $\it{HST}$/WFC3 imaging with deep Magellan/MagE spectra, we measure that the Ly$α$ flux F$_{\rm{Lyα}}$ = 1.56$\pm0.10 \times10^{-16}$ $\rm{erg\ s^{-1}\ cm^{-2}}$, which is among the luminous ($\approx$ 2.7 L$^{\star}_{\rm{Lyα}}$) Ly$α$ emitters at $\it{z} =$ 2.19. The $\rm{Lyα}$-based star formation rate (SFR) is $\sim$ 7 $\textrm{M_{\sun} \ yr^{-1}}$. These observational results favor that the star formation from the host galaxy could be the main mechanism to power the Ly$α$ emission. This new method sheds new light on the study of the kinematic structure and the spatial distribution of the extended Ly$α$ emitting regions around the QSO host.

preprint2020arXiv

Effective Opacity of the Intergalactic Medium from Galaxy Spectra Analysis

We measure the effective opacity ($τ_{eff}$) of the Intergalactic Medium (IGM) from the composite spectra of 281 Lyman-Break Galaxies (LBGs) in the redshift range $2 \lesssim z \lesssim 3$. Our spectra are taken from the COSMOS Lyman-Alpha Mapping And Tomographic Observations (CLAMATO) survey derived from the Low Resolution Imaging Spectrometer (LRIS) on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Ly$α$ forest, we measure the effective Ly$α$ opacity ($τ_{eff}$) in the $2.02 \leq z \leq 2.44$ range. At $z = 2.22$, we estimate $τ_{eff} = 0.159 \pm 0.001$ from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at $z<2.5$ indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian Processes model of our results and previous $τ_{eff}$ measurements that describes the steep redshift evolution in $τ_{eff}$ from $z = 1.5 - 4$.

preprint2020arXiv

First Constraints on Compact Dark Matter from Fast Radio Burst Microstructure

Despite existing constraints, it remains possible that up to $35\%$ of all dark matter is comprised of compact objects, such as the black holes in the 10-100\,M$_\odot$ range whose existence has been confirmed by LIGO. The strong gravitational lensing of transients such as FRBs and GRBs has been suggested as a more sensitive probe for compact dark matter than intensity fluctuations observed in microlensing experiments. Recently ASKAP has reported burst substructure down to $15μ$s timescales in FRBs in the redshift range $0.3-0.5$. We investigate here the implications of this for the detectability of compact dark matter by FRBs. We find that a sample size of $\sim130$ FRBs would be required to constrain compact dark matter to less than the existing 35$\%$ limit with 95$\%$ confidence, if it were distributed along $\gtrsim 1\,$Gpc-long FRB sightlines through the cosmic web. Conversely, existing constraints on the fraction of compact dark matter permit as many as 1 in $\approx 40$ of all $z \lesssim 0.4$ FRBs to exhibit micro-lensed burst structure. Approximately $170$ FRBs intercepting halos within $\sim 50\,$kpc would be required to exclude the fraction of compact dark matter in each intercepted halo to a similar level. Furthermore, we consider the cumulative effects of lensing of the FRB signal by a macroscopic dark matter distribution. We conclude that lensing from a uniform distribution of compact objects is likely not observable, but suggest that FRBs may set meaningful limits on power-law distributions of dark matter.

preprint2020arXiv

Fully probabilistic quasar continua predictions near Lyman-α with conditional neural spline flows

Measurement of the red damping wing of neutral hydrogen in quasar spectra provides a probe of the epoch of reionization in the early Universe. Such quantification requires precise and unbiased estimates of the intrinsic continua near Lyman-$α$ (Ly$α$), a challenging task given the highly variable Ly$α$ emission profiles of quasars. Here, we introduce a fully probabilistic approach to intrinsic continua prediction. We frame the problem as a conditional density estimation task and explicitly model the distribution over plausible blue-side continua ($1190\ \unicode{xC5} \leq λ_{\text{rest}} < 1290\ \unicode{xC5}$) conditional on the red-side spectrum ($1290\ \unicode{xC5} \leq λ_{\text{rest}} < 2900\ \unicode{xC5}$) using normalizing flows. Our approach achieves state-of-the-art precision and accuracy, allows for sampling one thousand plausible continua in less than a tenth of a second, and can natively provide confidence intervals on the blue-side continua via Monte Carlo sampling. We measure the damping wing effect in two $z>7$ quasars and estimate the volume-averaged neutral fraction of hydrogen from each, finding $\bar{x}_\text{HI}=0.304 \pm 0.042$ for ULAS J1120+0641 ($z=7.09$) and $\bar{x}_\text{HI}=0.384 \pm 0.133$ for ULAS J1342+0928 ($z=7.54$).

preprint2020arXiv

High time resolution and polarisation properties of ASKAP-localised fast radio bursts

Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescope has provided the first ensemble of bursts with this information. In this paper, we present the high time and spectral resolution, full polarisation observations of five localised FRBs to complement the results published for the previously studied ASKAP FRB~181112. We find that every FRB is highly polarised, with polarisation fractions ranging from 80 -- 100\%, and that they are generally dominated by linear polarisation. While some FRBs in our sample exhibit properties associated with an emerging archetype (i.e., repeating or apparently non-repeating), others exhibit characteristic features of both, implying the existence of a continuum of FRB properties. When examined at high time resolution, we find that all FRBs in our sample have evidence for multiple sub-components and for scattering at a level greater than expected from the Milky Way. We find no correlation between the diverse range of FRB properties (e.g., scattering time, intrinsic width, and rotation measure) and any global property of their host galaxy. The most heavily scattered bursts reside in the outskirts of their host galaxies, suggesting that the source-local environment rather than the host interstellar medium is likely the dominant origin of the scattering in our sample.

preprint2020arXiv

Modelling the prompt optical emission of GRB 180325A: the evolution of a spike from the optical to gamma-rays

The transition from prompt to the afterglow emission is one of the most exciting and least understood phases in gamma-ray bursts (GRBs). Correlations among optical, X-ray and gamma-ray emission in GRBs have been explored, to attempt to answer whether the earliest optical emission comes from internal and/or external shocks. We present optical photometric observations of GRB 180325A collected with the TAROT and RATIR ground-based telescopes. These observations show two strong optical flashes with separate peaks at $\sim50\;$s and $\sim120\;$s, followed by a temporally extended optical emission. We also present X-rays and gamma-ray observations of GRB 180325A, detected by the Burst Alert Telescope (BAT) and X-ray Telescope (XRT), on the Neil Gehrels Swift observatory, which both observed a narrow flash at $\sim80\;$s. We show that the prompt gamma- and X-ray early emission shares similar temporal and spectral features consistent with internal dissipation within the relativistic outflow (e.g. by internal shocks or magnetic reconnection), while the early optical flashes are likely generated by the reverse shock that decelerates the ejecta as it sweeps up the external medium.

preprint2020arXiv

Multi-phase outflows in post starburst E+A galaxies -- II. A direct connection between the neutral and ionized outflow phases

Post starburst E+A galaxies are systems that hosted a powerful starburst that was quenched abruptly. Simulations suggest that these systems provide the missing link between major merger ULIRGs and red and dead ellipticals, where AGN feedback is responsible for the expulsion or destruction of the molecular gas. However, many details remain unresolved and little is known about AGN-driven winds in this short-lived phase. We present spatially-resolved IFU spectroscopy with MUSE/VLT of SDSS J124754.95-033738.6, a post starburst E+A galaxy with a recent starburst that started 70 Myrs ago and ended 30 Myrs ago, with a peak SFR of $\sim 150\,\mathrm{M_{\odot}/yr}$. We detect disturbed gas throughout the entire field of view, suggesting triggering by a minor merger. We detect fast-moving multi-phased gas clouds, embedded in a double-cone face-on outflow, which are traced by ionized emission lines and neutral NaID \emph{emission} and absorption lines. We find remarkable similarities between the kinematics, spatial extents, and line luminosities of the ionized and neutral gas phases, and propose a model in which they are part of the same outflowing clouds, which are exposed to both stellar and AGN radiation. Our photoionization model provides consistent ionized line ratios, NaID absorption optical depths and EWs, and dust reddening. Using the model, we estimate, for the first time, the neutral-to-ionized gas mass ratio (about 20), the sodium neutral fraction, and the size of the outflowing clouds. This is one of the best ever observed direct connections between the neutral and ionized outflow phases in AGN.

preprint2020arXiv

Polyphorm: Structural Analysis of Cosmological Datasets via Interactive Physarum Polycephalum Visualization

This paper introduces Polyphorm, an interactive visualization and model fitting tool that provides a novel approach for investigating cosmological datasets. Through a fast computational simulation method inspired by the behavior of Physarum polycephalum, an unicellular slime mold organism that efficiently forages for nutrients, astrophysicists are able to extrapolate from sparse datasets, such as galaxy maps archived in the Sloan Digital Sky Survey, and then use these extrapolations to inform analyses of a wide range of other data, such as spectroscopic observations captured by the Hubble Space Telescope. Researchers can interactively update the simulation by adjusting model parameters, and then investigate the resulting visual output to form hypotheses about the data. We describe details of Polyphorm's simulation model and its interaction and visualization modalities, and we evaluate Polyphorm through three scientific use cases that demonstrate the effectiveness of our approach.

preprint2020arXiv

Project AMIGA: The Circumgalactic Medium of Andromeda

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.

preprint2020arXiv

PypeIt: The Python Spectroscopic Data Reduction Pipeline

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

preprint2020arXiv

Revealing the Dark Threads of the Cosmic Web

Modern cosmology predicts that matter in our Universe has assembled today into a vast network of filamentary structures colloquially termed the Cosmic Web. Because this matter is either electromagnetically invisible (i.e., dark) or too diffuse to image in emission, tests of this cosmic web paradigm are limited. Wide-field surveys do reveal web-like structures in the galaxy distribution, but these luminous galaxies represent less than 10% of baryonic matter. Statistics of absorption by the intergalactic medium (IGM) via spectroscopy of distant quasars support the model yet have not conclusively tied the diffuse IGM to the web. Here, we report on a new method inspired by the Physarum polycephalum slime mold that is able to infer the density field of the Cosmic Web from galaxy surveys. Applying our technique to galaxy and absorption-line surveys of the local Universe, we demonstrate that the bulk of the IGM indeed resides in the Cosmic Web. From the outskirts of Cosmic Web filaments, at approximately the cosmic mean matter density (rho_m) and approx. 5 virial radii from nearby galaxies, we detect an increasing H I absorption signature towards higher densities and the circumgalactic medium, to approx. 200 rho_m. However, the absorption is suppressed within the densest environments, suggesting shock-heating and ionization deep within filaments and/or feedback processes within galaxies.

preprint2020arXiv

Spectropolarimetric analysis of FRB 181112 at microsecond resolution: Implications for Fast Radio Burst emission mechanism

We have developed a new coherent dedispersion mode to study the emission of Fast Radio Bursts that trigger the voltage capture capability of the Australian SKA Pathfinder (ASKAP) interferometer. In principle the mode can probe emission timescales down to 3 ns with full polarimetric information preserved. Enabled by the new capability, here we present a spectropolarimetric analysis of FRB 181112 detected by ASKAP, localized to a galaxy at redshift 0.47. At microsecond time resolution the burst is resolved into four narrow pulses with a rise time of just $15 μ$s for the brightest. The pulses have a diversity of morphology, but do not show evidence for temporal broadening by turbulent plasma along the line of sight, nor is there any evidence for periodicity in their arrival times. The pulses are highly polarized (up to 95%), with the polarization position angle varying both between and within pulses. The pulses have apparent rotation measures that vary by $15\pm 2\, {\rm rad \,m^{-2}}$ and apparent dispersion measures that vary by $0.041\pm 0.004\,{\rm pc\,cm^{-3}}$. Conversion between linear and circular polarization is observed across the brightest pulse. We conclude that the FRB 181112 pulses are most consistent with being a direct manifestation of the emission process or the result of propagation through a relativistic plasma close to the source. This demonstrates that our method, which facilitates high-time-resolution polarimetric observations of FRBs, can be used to study not only burst emission processes, but also a diversity of propagation effects present on the gigaparsec paths they traverse.

preprint2020arXiv

Statistical correlation between the distribution of Ly$α$ emitters and IGM HI at $z\sim2.2$ mapped by Subaru/Hyper Suprime-Cam

The correlation between neutral Hydrogen (HI) in the intergalactic medium (IGM) and galaxies now attracts great interests. We select four fields which include several coherently strong Ly$α$ absorption systems at $z\sim2.2$ detected by using background quasars from the whole SDSS/(e)BOSS database. Deep narrow-band and $g$-band imaging are performed using the Hyper Suprime-Cam on the Subaru Telescope. We select out 2,642 Ly$α$ emitter (LAE) candidates at $z=2.177\pm0.023$ down to the Ly$α$ luminosity of $L_{\text{Ly}α}\approx 2 \times 10^{42} {\rm erg~s}^{-1}$ to construct the galaxy overdensity maps, covering an effective area of 5.39 deg$^2$. Combining the sample with the Ly$α$ absorption estimated from 64 (e)BOSS quasar spectra, we find a moderate to strong correlation between the LAE overdensity $δ_{\rm LAE}$ and the effective optical depth $τ_{\rm LoS}$ in line-of-sights, with $P$-value$=0.09\%$ ($<0.01\%$) when the field that contains a significant quasar overdensity is in(ex)cluded. The cross-correlation analysis also clearly suggests that up to $4\pm1$ pMpc, LAEs tend to cluster in the regions rich in HI gas, indicated by the high $τ_{\rm LoS}$, and avoid the low $τ_{\rm LoS}$ region where the HI gas is deficient. By averaging the $τ_{\rm LoS}$ as a function of the projected distance ($d$) to LAEs, we find a $30\%$ excess signal at $2σ$ level at $d<200$ pkpc, indicating the dense HI in circumgalactic medium, and a tentative excess at $400<d<600$ pkpc in IGM regime, corroborating the cross-correlation signal detected at about $0.5$ pMpc. These statistical analyses indicate that galaxy$-$IGM HI correlations exist on scales ranging from several hundred pkpc to several pMpc at $z\sim2.2$.

preprint2020arXiv

Testing Galaxy Formation Simulations with Damped Lyman-$α$ Abundance and Metallicity Evolution

We examine the properties of damped Lyman-$α$ absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: {\sc Simba} and {\sc Technicolor Dawn}. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multi-frequency radiative transfer to model an inhomogeneous ultraviolet background (UVB) self-consistently and is calibrated to match the Thomson scattering optical depth, UVB amplitude, and Ly-$α$ forest mean transmission at $z>5$. Both simulations are in reasonable agreement with the measured stellar mass and star formation rate functions at $z\geq 3$, and both reproduce the observed neutral hydrogen cosmological mass density, $Ω_{\rm HI}(z)$. However, the DLA abundance and metallicity distribution are sensitive to the galactic outflows' feedback and the UVB amplitude. Adopting a strong UVB and/or slow outflows under-produces the observed DLA abundance, but yields broad agreement with the observed DLA metallicity distribution. By contrast, faster outflows eject metals to larger distances, yielding more metal-rich DLAs whose observational selection may be more sensitive to dust bias. The DLA metallicity distribution in models adopting an ${\rm H}_2$-regulated star formation recipe includes a tail extending to $[M/H] \ll -3$, lower than any DLA observed to date, owing to curtailed star formation in low-metallicity galaxies. Our results show that DLA observations play an imporant role in constraining key physical ingredients in galaxy formation models, complementing traditional ensemble statistics such as the stellar mass and star formation rate functions.

preprint2020arXiv

The host galaxies and progenitors of Fast Radio Bursts localized with the Australian Square Kilometre Array Pathfinder

The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11<z<0.48$. All four are massive galaxies (log( $M_{*}/ M_{\odot}$) $\sim 9.4 -10.4$) with modest star-formation rates of up to $2M_{\odot}$yr$^{-1}$ -- very different to the host galaxy of the first repeating FRB 121102, which is a dwarf galaxy with a high specific star-formation rate. The FRBs localized by ASKAP typically lie in the outskirts of their host galaxies, which appears to rule out FRB progenitor models that invoke active galactic nuclei (AGN) or free-floating cosmic strings. The stellar population seen in these host galaxies also disfavors models in which all FRBs arise from young magnetars produced by superluminous supernovae (SLSNe), as proposed for the progenitor of FRB 121102. A range of other progenitor models (including compact-object mergers and magnetars arising from normal core-collapse supernovae) remain plausible.

preprint2020arXiv

The PHLEK Survey: A New Determination of the Primordial Helium Abundance

We present Keck NIRSPEC and Keck NIRES spectroscopy of sixteen metal-poor galaxies that have pre-existing optical observations. The near-infrared (NIR) spectroscopy specifically targets the HeI 10830A emission line, due to its sensitivity to the physical conditions of the gas in HII regions. We use these NIR observations, combined with optical spectroscopy, to determine the helium abundance of sixteen galaxies across a metallicity range 12 + log (O/H) = 7.13 - 8.00. This data set is combined with two other samples where metallicity and helium abundance measurements can be secured: star-forming galaxies selected from the Sloan Digital Sky Survey spectroscopic database and existing low-metallicity systems in the literature. We calculate a linear fit to these measurements, accounting for intrinsic scatter, and report a new determination of the primordial helium number abundance, y_P = 0.0805+/-0.0017, which corresponds to a primordial helium mass fraction Y_P = 0.2436+/-0.0040. Using our determination of the primordial helium abundance in combination with the latest primordial deuterium measurement, (D/H)_P x 10^5 =2.527+/-0.030, we place a bound on the baryon density Omega_bh^2 = 0.0215+/-0.0005 and the effective number of neutrino species Neff = 2.85+/-0.28. These values are in 1.3-sigma agreement with those deduced from the Planck satellite observations of the temperature fluctuations imprinted on the Cosmic Microwave Background.

preprint2020arXiv

The Third Data Release of the KODIAQ Survey

We present and make publicly available the third data release (DR3) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR3 consists of a fully-reduced sample of 727 quasars at 0.1 < z < 6.4 observed with ESI at moderate resolution (4000 < R < 10000). DR3 contains 872 spectra available in flux calibrated form, representing a sum total exposure time of approximately 2.8 megaseconds. These co-added spectra arise from a total of 2753 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a data reduction package made available through the XIDL distribution. DR3 is publicly available to the community, housed as a higher level science product at the KOA and in the igmspec database.

preprint2019arXiv

Discovery of a Rare Late-Type, Low-Mass Wolf-Rayet Star in the LMC

We report the serendipitous discovery of an object, UVQS J060819.93-715737.4, with a spectrum dominated by extremely intense, narrow C II emission lines. The spectrum is similar to those of the very rare, late-type [WC11] low-mass Wolf-Rayet stars. Despite the recognition of these stars as a distinct class decades ago, there remains barely a handful of Galactic members, all of which are also planetary-nebula central stars. Although no obvious surrounding nebulosity is present in J0608, [O II], [N II], and [S II] emission suggest the presence of an inconspicuous, low-excitation nebula. There is low-amplitude incoherent photometric variability on timescales of days to years, as well as numerous prominent P Cygni profiles, implying mass loss. There are indications of a binary companion. The star is located on the outskirts of the LMC, and the observed radial velocity (~250 km/s) and proper motion strongly suggest membership. If indeed an LMC member, this is the first extragalactic late [WC] star, and the first with an accurately determined luminosity, as the Galactic examples are too distant for precise parallax determinations. A high-quality, broad-coverage spectrum of the prototype of the late [WC] class, CPD -56 8032, is also presented. We discuss different excitation mechanisms capable of producing the great strength of the C II emission. Numerous autoionizing levels of C II are definitely populated by processes other than dielectronic recombination. Despite the spectacular emission spectra, observational selection makes objects such as these difficult to discover. Members of the [WC11] class may in fact be considerably more common than the handful of previously known late [WC] stars.

preprint2019arXiv

Evolution of the Cool Gas in the Circumgalactic Medium (CGM) of Massive Halos -- A Keck Cosmic Web Imager (KCWI) Survey of Ly$α$ Emission around QSOs at $z\approx2$

Motivated by the recent discovery of the near-ubiquity of Ly$α$ emission around $z \gtrsim 3$ QSOs, we performed a systematic study of QSO circumgalactic Ly$α$ emission at $z\approx2$, utilizing the unique capability of the Keck Cosmic Web Imager (KCWI) -- a new wide-field, blue sensitive integral-field spectrograph (IFU). In this paper, we present KCWI observations on a sample of 16 ultraluminous Type-I QSOs at $z=2.1-2.3$ with ionizing luminosities of $L_{\rm{ν_{\rm{LL}}}}=10^{31.1-32.3}$ erg s$^{-1}$ Hz$^{-1}$. We found that 14 out of 16 QSOs are associated with Ly$α$ nebulae with projected linear-sizes larger than 50 physical kpc (pkpc). Among them, four nebulae have enormous Ly$α$ emission with the Ly$α$ surface brightness $SB_{\rm{Lyα}} >10^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ on the $>100$ kpc scale, extending beyond the field of view of KCWI. Our KCWI observations reveal that most $z\approx2$ QSO nebulae have a more irregular morphology compared to those at $z\gtrsim3$. In turn, we measure that the circularly-averaged surface brightness (SB) at $z\approx2$ is 0.4 dex fainter than the redshift-corrected, median SB at $z\gtrsim 3$. The Ly$α$ SB profile (SB$_{\rm{Lyα}}$) of QSOs at $z\approx 2$ can be described by a power law of SB$_{\rm{Lyα},z\approx2.3}= 3.7\times10^{-17}\times(r/40)^{-1.8}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$, with the slope similar to that at $z\gtrsim3$. The observed lower redshift-corrected, circularly-averaged SB may be mainly due to the lower covering factor of cool gas clouds in massive halos at $z\approx2$.

preprint2019arXiv

GRB 180620A: Evidence for late-time energy injection

The early optical emission of gamma-ray bursts gives an opportunity to understand the central engine and first stages of these events. About 30\% of GRBs present flares whose origin is still a subject of discussion. We present optical photometry of GRB 180620A with the COATLI telescope and RATIR instrument. COATLI started to observe from the end of prompt emission at $T+39.3$~s and RATIR from $T+121.4$~s. We supplement the optical data with the X-ray light curve from \emph{Swift}/XRT. %The optical and X-ray light curves show very unusual behavior with features clearly beyond the standard fireball model. We observe an optical flare from $T+110$ to $T+550$~s, with a temporal index decay $α_\mathrm{O,decay}=1.32\pm 0.01$, and a $Δt/t=1.63$, which we interpret as the signature of a reverse shock component. After the initial normal decay the light curves show a long plateau from $T+500$ to $T+7800$~s both in X-rays and the optical before decaying again after an achromatic jet break at $T+7800$~s. Fluctuations are seen during the plateau phase in the optical. Adding to the complexity of GRB afterglows, the plateau phase (typically associated with the coasting phase of the jet) is seen in this object after the ``normal'' decay phase (emitted during the deceleration phase of the jet) and the jet break phase occurs directly after the plateau. We suggest that this sequence of events can be explained by a rapid deceleration of the jet with $t_d\lesssim 40$ s due to the high density of the environment ($\approx 100$ cm$^{-3}$) followed by reactivation of the central engine which causes the flare and powers the plateau phase.

preprint2019arXiv

MUSE Analysis of Gas around Galaxies (MAGG) -- I: Survey design and the environment of a near pristine gas cloud at z~3.5

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.

preprint2019arXiv

The low density and magnetization of a massive galaxy halo exposed by a fast radio burst

Present-day galaxies are surrounded by cool and enriched halo gas extending to hundreds of kiloparsecs. This halo gas is thought to be the dominant reservoir of material available to fuel future star formation, but direct constraints on its mass and physical properties have been difficult to obtain. We report the detection of a fast radio burst (FRB 181112) with arcsecond precision, which passes through the halo of a foreground galaxy. Analysis of the burst shows the halo gas has low net magnetization and turbulence. Our results imply predominantly diffuse gas in massive galactic halos, even those hosting active supermassive black holes, contrary to some previous results.

preprint2016arXiv

A Deep Search For Faint Galaxies Associated With Very Low-redshift C IV Absorbers: III. The Mass- and Environment-dependent Circumgalactic Medium

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.

preprint2016arXiv

A Universal Density Structure for Circum-Galactic Gas

We develop a new method to constrain the physical conditions in the cool (~10^4 K) circumgalactic medium (CGM) from measurements of ionic column densities, by assuming that the cool CGM spans a large range of gas densities and that small high-density clouds are hierarchically embedded in large low-density clouds. The new method combines the information available from different sightlines during the photoionization modeling, thus yielding tighter constraints on CGM properties compared to traditional methods which model each sightline individually. Applying this new technique to the COS-Halos survey of low-redshift ~L* galaxies, we find that we can reproduce all observed ion columns in all 44 galaxies in the sample, from the low-ions to OVI, with a single universal density structure for the cool CGM. The gas densities span the range 50 < ρ/ρ_mean < 5x10^5 (ρ_mean is the cosmic mean), while the physical size of individual clouds scales as ~ρ^-1, from ~35 kpc of the low density OVI clouds to ~6 pc of the highest density low-ion clouds. The deduced cloud sizes are too small for this density structure to be driven by self-gravity, thus its physical origin is unclear. The implied cool CGM mass within the virial radius is 1.3x10^10 M_sun (~1% of the halo mass), distributed rather uniformly over the four decades in density. The mean cool gas density profile scales as R^-1.0, where R is the distance from the galaxy center. We construct a 3D model of the cool CGM based on our results, which we argue provides a benchmark for the CGM structure in hydrodynamic simulations. Our results can be tested by measuring the coherence scales of different ions.

preprint2016arXiv

An Ultraviolet Spectrum of the Tidal Disruption Flare ASASSN-14li

We present a Hubble Space Telescope STIS spectrum of ASASSN-14li, the first rest-frame UV spectrum of a tidal disruption flare (TDF). The underlying continuum is well fit by a blackbody with $T_{\mathrm{UV}} = 3.5 \times 10^{4}$ K, an order of magnitude smaller than the temperature inferred from X-ray spectra (and significantly more precise than previous efforts based on optical and near-UV photometry). Super-imposed on this blue continuum, we detect three classes of features: narrow absorption from the Milky Way (probably a high-velocity cloud), and narrow absorption and broad (FWHM $\approx 2000$-8000 km s$^{-1}$) emission lines at/near the systemic host velocity. The absorption lines are blueshifted with respect to the emission lines by $Δv = -(250$-400) km s$^{-1}$. Due both to this velocity offset and the lack of common low-ionization features (Mg II, Fe II), we argue these arise from the same absorbing material responsible for the low-velocity outflow discovered at X-ray wavelengths. The broad nuclear emission lines display a remarkable abundance pattern: N III], N IV], He II are quite prominent, while the common quasar emission lines of C III] and Mg II are weak or entirely absent. Detailed modeling of this spectrum will help elucidate fundamental questions regarding the nature of the emission process(es) at work in TDFs, while future UV spectroscopy of ASASSN-14li would help to confirm (or refute) the previously proposed connection between TDFs and "N-rich" quasars.

preprint2016arXiv

Constraining the lifetime and opening angle of quasars using fluorescent Lyman $α$ emission: the case of Q0420-388

A toy model is developed to understand how the spatial distribution of fluorescent emitters in the vicinity of bright quasars could be affected by the geometry of the quasar bi-conical radiation field and by its lifetime. The model is then applied to the distribution of high equivalent width Lyman $α$ emitters (with rest-frame equivalent widths above 100 A, threshold used in e.g. Trainor & Steidel, 2013) identified in a deep narrow-band 36x36 arcmin$^2$ image centered on the luminous quasar Q0420-388. These emitters are found to the edge of the field and show some evidence of an azimuthal asymmetry on the sky of the type expected if the quasar is radiating in a bipolar cone. If these sources are being fluorescently illuminated by the quasar, the two most distant objects require a lifetime of at least 15 Myr for an opening angle of 60 degrees or more, increasing to more than 40 Myr if the opening angle is reduced to a minimum 30 degrees. However, few of the other expected signatures of boosted fluorescence are seen at the current survey limits, e.g. a fall off in Lyman $α$ brightness, or equivalent width, with distance. Furthermore, to have most of the Lyman $α$ emission of the two distant sources to be fluorescently boosted would require the quasar to have been significantly brighter in the past. This suggests that these particular sources may not be fluorescent, invalidating the above lifetime constraints. This would cast doubt on the use of this relatively low equivalent width threshold and thus on the lifetime analysis in Trainor & Steidel (2013).

preprint2016arXiv

Low-metallicity Absorbers Account for Half of the Dense Circumgalactic Gas at z < 1

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.

preprint2016arXiv

MUSE searches for galaxies near very metal-poor gas clouds at z~3: new constraints for cold accretion models

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α 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α 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α 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α nebula at the quasar redshift and three Lyα emitters reminiscent of a "dark galaxy" population.

preprint2016arXiv

Nature and statistical properties of quasar associated absorption systems in the XQ-100 Legacy Survey

We statistically study the physical properties of a sample of narrow absorption line (NAL) systems looking for empirical evidences to distinguish between intrinsic and intervening NALs without taking into account any a priori definition or velocity cut-off. We analyze the spectra of 100 quasars with 3.5 < z$\rm_{em}$ < 4.5, observed with X-shooter/VLT in the context of the XQ-100 Legacy Survey. We detect a $\sim$ 8 $σ$ excess in the number density of absorbers within 10,000 km/s of the quasar emission redshift with respect to the random occurrence of NALs. This excess does not show a dependence on the quasar bolometric luminosity and it is not due to the redshift evolution of NALs. It extends far beyond the standard 5000 km/s cut-off traditionally defined for associated absorption lines. We propose to modify this definition, extending the threshold to 10,000 km/s when also weak absorbers (equivalent width < 0.2 Å) are considered. We infer NV is the ion that better traces the effects of the quasar ionization field, offering the best statistical tool to identify intrinsic systems. Following this criterion we estimate that the fraction of quasars in our sample hosting an intrinsic NAL system is 33 percent. Lastly, we compare the properties of the material along the quasar line of sight, derived from our sample, with results based on close quasar pairs investigating the transverse direction. We find a deficiency of cool gas (traced by CII) along the line of sight associated with the quasar host galaxy, in contrast with what is observed in the transverse direction.

preprint2016arXiv

Precise limits on cosmological variability of the fine-structure constant with zinc and chromium quasar absorption lines

The strongest transitions of Zn and CrII are the most sensitive to relative variations in the fine-structure constant ($Δα/α$) among the transitions commonly observed in quasar absorption spectra. They also lie within just 40Å of each other (rest frame), so they are resistant to the main systematic error affecting most previous measurements of $Δα/α$: long-range distortions of the wavelength calibration. While Zn and CrII absorption is normally very weak in quasar spectra, we obtained high signal-to-noise, high-resolution echelle spectra from the Keck and Very Large Telescopes of 9 rare systems where it is strong enough to constrain $Δα/α$ from these species alone. These provide 12 independent measurements (3 quasars were observed with both telescopes) at redshifts 1.0--2.4, 11 of which pass stringent reliability criteria. These 11 are all consistent with $Δα/α=0$ within their individual uncertainties of 3.5--13 parts per million (ppm), with a weighted mean $Δα/α= 1.2\pm1.7_{\rm stat}\pm0.9_{\rm sys}$ ppm (1$σ$ statistical and systematic uncertainties), indicating no significant cosmological variations in $α$. This is the first statistical sample of absorbers that is resistant to long-range calibration distortions (at the $<$1 ppm level), with a precision comparable to previous large samples of $\sim$150 (distortion-affected) absorbers. Our systematic error budget is instead dominated by much shorter-range distortions repeated across echelle orders of individual spectra.

preprint2016arXiv

Quasars Probing Quasars VIII. The Physical Properties of the Cool Circumgalactic Medium Surrounding z ~ 2-3 Massive Galaxies Hosting Quasars

We characterize the physical properties of the cool T ~10^4 K circumgalactic medium surrounding z ~2-3 quasar host galaxies, which are predicted to evolve into present day massive ellipticals. Using a statistical sample of 14 quasar pairs with projected separation < 300 kpc and high dispersion, high S/N spectra, we find extreme kinematics with low metal ion lines typically spanning ~ 500 km/s, exceeding any previously studied galactic population. The CGM is significantly enriched, even beyond the virial radius, with a median metallicity [M/H] ~ -0.6. The alpha/Fe abundance ratio is enhanced, suggesting that halo gas is primarily enriched by core-collapse supernovae. The projected cool gas mass within the virial radius is estimated to be 1.9*10^11 M_sun (R_\perp/160 kpc)^2, accounting for ~ 1/3 of the galaxy halo baryonic budget. The ionization state of CGM gas increases with projected distance from the foreground quasars, contrary to expectation if the quasar dominates the ionizing radiation flux. However, we also found peculiarities not exhibited in the CGM of other galaxy populations. In one absorption system, we may be detecting unresolved fluorescent Ly-alpha emission, and another system shows strong NV lines. Taken together these anomalies suggest that transverse sightlines are at least in some cases possibly illuminated. We also discovered a peculiar case where detection of the CII* fine structure line implies an electron density > 100 cm^-3 and subparsec scale gas clumps.

preprint2016arXiv

The COS-Halos Survey: Origins of the Highly Ionized Circumgalactic Medium of Star-Forming Galaxies

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by OVI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of NV/OVI for 39 fitted Voigt profile components of OVI. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to NV/OVI and measurements of N$_{\rm OVI}$ imply unphysically large path lengths > 100 kpc. Furthermore, very broad OVI absorption (b > 40 km s$^{-1}$) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed OVI: (1) highly structured gas clouds photoionized primarily by local high energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic OVI does not align with any low-ionization state gas within $\pm$50 km s$^{-1}$ and is found only in halos with M$_{\rm halo}$ < 10$^{12}$ M$_{\odot}$. We suggest that this type of unmatched OVI absorption traces the hot corona itself at a characteristic temperature of 10$^{5.5}$ K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos.

preprint2016arXiv

The Cosmic Evolution of the Metallicity Distribution of Ionized Gas Traced by Lyman Limit Systems

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.

preprint2016arXiv

The HI Content of the Universe over the Past 10 Gyrs

We use the Hubble Space Telescope (HST) archive of ultraviolet (UV) quasar spectroscopy to conduct the first blind survey for damped Ly-alpha absorbers (DLAs) at low redshift (z < 1.6). Our statistical sample includes 463 quasars with spectral coverage spanning a total redshift path, dz = 123.3 or an absorption path, dX = 229.7. Within this survey path, we identify 4 DLAs, defined as absorbers with HI column density N(HI) >= 10^20.3cm-2, which implies an incidence per absorption length, l(X)= 0.017(+0.014-0.008) at a median survey path redshift of z=0.623. While our estimate of l(X) is lower than earlier estimates at z ~ 0 from HI 21cm emission studies, the results are consistent within the measurement uncertainties. Our dataset is too small to properly sample the N(HI) frequency distribution function f(N(HI),X), but the observed distribution agrees with previous estimates at z > 2. Adopting the z > 2 shape of f(N(HI),X), we infer an HI mass density at z ~ 0.6 of rho_HI = 0.25(+0.20-0.12) x 10^8 Msol Mpc-3. This is significantly lower than previous estimates from targeted DLA surveys with the HST, but consistent with results from low-z HI 21cm observations, and suggests that the neutral gas density of the universe has been decreasing over the past 10 Gyrs.

preprint2016arXiv

The Optical Variability of SDSS Quasars from Multi-epoch Spectroscopy. III. A Sudden UV Cutoff in Quasar SDSS J2317+0005

We have collected near-infrared to X-ray data of 20 multi-epoch heavily reddened SDSS quasars to investigate the physical mechanism of reddening. Of these, J2317+0005 is found to be a UV cutoff quasar. Its continuum, which usually appears normal, decreases by a factor 3.5 at 3000Å, compared to its more typical bright state during an interval of 23 days. During this sudden continuum cut-off, the broad emission line fluxes do not change, perhaps due to the large size of the Broad Line Region (BLR), r > 23 / (1+z) days. The UV continuum may have suffered a dramatic drop out. However, there are some difficulties with this explanation. Another possibility is that the intrinsic continuum did not change, but was temporarily blocked out, at least towards our line of sight. As indicated by X-ray observations, the continuum rapidly recovers after 42 days. A comparison of the bright state and dim states would imply an eclipse by a dusty cloud with a reddening curve having a remarkably sharp rise shortward of 3500Å. Under the assumption of being eclipsed by a Keplerian dusty cloud, we characterized the cloud size with our observations, however, which is a little smaller than the 3000Å continuum-emitting size inferred from accretion disk models. Therefore, we speculate this is due to a rapid outflow or inflow with a dusty cloud passing through our line-of-sight to the center.

preprint2016arXiv

The Stacked Lyman-Alpha Emission Profile from the Circum-Galactic Medium of z~2 Quasars

In the context of the FLASHLIGHT survey, we obtained deep narrow band images of 15 $z\sim2$ quasars with GMOS on Gemini-South in an effort to measure Ly$α$ emission from circum- and inter-galactic gas on scales of hundreds of kpc from the central quasar. We do not detect bright giant Ly$α$ nebulae (SB~10$^{-17}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at distances >50 kpc) around any of our sources, although we routinely ($\simeq47$%) detect smaller scale <50 kpc Ly$α$ emission at this SB level emerging from either the extended narrow emission line regions powered by the quasars or by star-formation in their host galaxies. We stack our 15 deep images to study the average extended Ly$α$ surface brightness profile around $z\sim2$ quasars, carefully PSF-subtracting the unresolved emission component and paying close attention to sources of systematic error. Our analysis, which achieves an unprecedented depth, reveals a surface brightness of SB$_{\rm Lyα}\sim10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at $\sim200$ kpc, with a $2.3σ$ detection of Ly$α$ emission at SB$_{\rm Lyα}=(5.5\pm3.1)\times10^{-20}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ within an annulus spanning 50 kpc <R< 500 kpc from the quasars. Assuming this Ly$α$ emission is powered by fluorescence from highly ionized gas illuminated by the bright central quasar, we deduce an average volume density of $n_{\rm H}=0.6\times10^{-2}$ cm$^{-3}$ on these large scales. Our results are in broad agreement with the densities suggested by cosmological hydrodynamical simulations of massive ($M\simeq10^{12.5}M_\odot$) quasar hosts, however they indicate that the typical quasars at these redshifts are surrounded by gas that is a factor of ~100 times less dense than the (~1 cm$^{-3}$) gas responsible for the giant bright Ly$α$ nebulae around quasars recently discovered by our group.

preprint2016arXiv

The UV-bright Quasar Survey (UVQS): DR1

We present the first data release (DR1) from our UV-bright Quasar Survey (UVQS) for new $z \sim 1$ active galactic nuclei (AGN) across the sky. Using simple GALEX UV and WISE near-IR color selection criteria, we generated a list of 1450 primary candidates with $FUV < 18.5$ mag. We obtained discovery spectra, primarily on 3m-class telescopes, for 1040 of these candidates and confirmed 86% as AGN with redshifts generally at $z>0.5$. Including a small set of observed secondary candidates, we report the discovery of 217 AGN with $FUV < 18$ mag that had no previously reported spectroscopic redshift. These are excellent potential targets for UV spectroscopy before the end of the {\it Hubble Space Telescope} mission. The main data products are publicly released through the Mikulski Archive for Space Telescopes.

preprint2015arXiv

A search for H$α$ emission in high-metallicity damped Lyman-$α$ systems at $z \sim 2.4$

We report on a sensitive search for redshifted H$α$ line-emission from three high-metallicity damped Ly$α$ absorbers (DLAs) at $z \approx 2.4$ with the Near-infrared Integral Field Spectrometer (NIFS) on the Gemini-North telescope, assisted by the ALTtitude conjugate Adaptive optics for the InfraRed (ALTAIR) system with a laser guide star. Within the NIFS field-of-view, $\approx 3.22" \times 2.92"$ corresponding to $\approx 25$ kpc $ \times 23$ kpc at $z=2.4$, we detect no statistically significant line-emission at the expected redshifted H$α$ wavelengths. The measured root-mean-square noise fluctuations in $0.4"$ apertures are $1-3\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$. Our analysis of simulated, compact, line-emitting sources yields stringent limits on the star-formation rates (SFRs) of the three DLAs, $< 2.2$~M$_{\odot}$ yr$^{-1}$ ($3σ$) for two absorbers, and $< 11$~M$_{\odot}$ yr$^{-1}$ ($3σ$) for the third, at all impact parameters within $\approx 12.5$~kpc to the quasar sightline at the DLA redshift. For the third absorber, the SFR limit is $< 4.4$~M$_\odot$ yr$^{-1}$ for locations away from the quasar sightline. These results demonstrate the potential of adaptive optics-assisted, integral field unit searches for galaxies associated with high-$z$ DLAs.

preprint2015arXiv

Cold gas and a Milky Way-type 2175 Å bump in a metal-rich and highly depleted absorption system

We report the detection of a strong Milky Way-type 2175 Å$ $ extinction bump at $z$ = 2.1166 in the quasar spectrum towards SDSS J121143.42+083349.7 from the Sloan Digital Sky Survey (SDSS) Data Release 10. We conduct follow up observations with the Echelle Spectrograph and Imager (ESI) onboard the Keck-II telescope and the Ultraviolet and Visual Echelle Spectrograph (UVES) on the VLT. This 2175 Å$ $ absorber is remarkable in that we simultaneously detect neutral carbon (C I), neutral chlorine (Cl I), and carbon monoxide (CO). It also qualifies as a damped Lyman alpha system. The J1211+0833 absorber is found to be metal-rich and has a dust depletion pattern resembling that of the Milky Way disk clouds. We use the column densities of the C I fine structure states and the C II/C I ratio (under the assumption of ionization equilibrium) to derive the temperature and volume density in the absorbing gas. A Cloudy photoionization model is constructed, which utilizes additional atoms/ions to constrain the physical conditions. The inferred physical conditions are consistent with a canonical cold (T $\sim$ 100 K) neutral medium with a high density ($n$(H I) $\sim$ 100 cm$^{-3}$) and a slightly higher pressure than the local interstellar medium. Given the simultaneous presence of C I, CO, and the 2175 Å$ $ bump, combined with the high metallicity, high dust depletion level and overall low ionization state of the gas, the absorber towards J1211+0833 supports the scenario that the presence of the bump requires an evolved stellar population.

preprint2015arXiv

Data Reduction with the MIKE Spectrometer

This manuscript describes the design, usage, and data-reduction pipeline developed for the Magellan Inamori Kyocera Echelle (MIKE) spectrometer used with the Magellan telescope at the Las Campanas Observatory. We summarize the basic characteristics of the instrument and discuss observational procedures recommended for calibrating the standard data products. We detail the design and implementation of an IDL based data-reduction pipeline for MIKE data (since generalized to other echelle spectrometers, e.g. Keck/HIRES, VLT/UVES). This includes novel techniques for flat-fielding, wavelength calibration, and the extraction of echelle spectroscopy. Sufficient detail is provided in this manuscript to enable inexperienced observers to understand the strengths and weaknesses of the instrument and software package and an assessment of the related systematics.

preprint2015arXiv

Deep HeII and CIV Spectroscopy of a Giant Lyman alpha Nebula: Dense Compact Gas Clumps in the Circumgalactic Medium of a z~2 Quasar

The recent discovery by Cantalupo et al. (2014) of the largest (~500 kpc) and luminous Ly-alpha nebula associated with the quasar UM287 (z=2.279) poses a great challenge to our current understanding of the astrophysics of the halos hosting massive z~2 galaxies. Either an enormous reservoir of cool gas is required $M\simeq10^{12}$ $M_{\odot}$, exceeding the expected baryonic mass available, or one must invoke extreme gas clumping factors not present in high-resolution cosmological simulations. However, observations of Ly-alpha emission alone cannot distinguish between these two scenarios. We have obtained the deepest ever spectroscopic integrations in the HeII and CIV lines with the goal of detecting extended line emission, but detect neither line to a 3$σ$ limiting SB $\simeq10^{-18}$ erg/s/cm$^2$/arcsec$^2$. We construct models of the expected emission spectrum in the highly probable scenario that the nebula is powered by photoionization from the central hyper-luminous quasar. The non-detection of HeII implies that the nebular emission arises from a mass $M_{\rm c}\lesssim6.4\times10^{10}$ $M_{\odot}$ of cool gas on ~200 kpc scales, distributed in a population of remarkably dense ($n_{\rm H}\gtrsim3$ cm$^{-3}$) and compact ($R\lesssim20$ pc) clouds, which would clearly be unresolved by current cosmological simulations. Given the large gas motions suggested by the Ly-alpha line ($v\simeq$ 500 km/s), it is unclear how these clouds survive without being disrupted by hydrodynamic instabilities. Our study serves as a benchmark for future deep integrations with current and planned wide-field IFU such as MUSE, KCWI, and KMOS. Our work suggest that a $\simeq$ 10 hr exposure would likely detect ~10 rest-frame UV/optical emission lines, opening up the possibility of conducting detailed photoionization modeling to infer the physical state of gas in the CGM.

preprint2015arXiv

iPTF14yb: The First Discovery of a GRB Afterglow Independent of a High-Energy Trigger

We report here the discovery by the Intermediate Palomar Transient Factory (iPTF) of iPTF14yb, a luminous ($M_{r}\approx-27.8$ mag), cosmological (redshift 1.9733), rapidly fading optical transient. We demonstrate, based on probabilistic arguments and a comparison with the broader population, that iPTF14yb is the optical afterglow of the long-duration gamma-ray burst GRB 140226A. This marks the first unambiguous discovery of a GRB afterglow prior to (and thus entirely independent of) an associated high-energy trigger. We estimate the rate of iPTF14yb-like sources (i.e., cosmologically distant relativistic explosions) based on iPTF observations, inferring an all-sky value of $\Re_{\mathrm{rel}}=610$ yr$^{-1}$ (68% confidence interval of 110-2000 yr$^{-1}$). Our derived rate is consistent (within the large uncertainty) with the all-sky rate of on-axis GRBs derived by the Swift satellite. Finally, we briefly discuss the implications of the nondetection to date of bona fide "orphan" afterglows (i.e., those lacking detectable high-energy emission) on GRB beaming and the degree of baryon loading in these relativistic jets.

preprint2015arXiv

MApping the Most Massive Overdensities Through Hydrogen (MAMMOTH) I: Methodology

Modern cosmology predicts that a galaxy overdensity is associated to a large reservoir of the intergalactic gas, which can be traced by the Ly$α$ forest absorption. We have undertaken a systematic study of the relation between Coherently Strong intergalactic Ly$α$ Absorption systems (CoSLAs), which have highest optical depth ($τ$) in $τ$ distribution, and mass overdensities on the scales of $\sim$ 10 - 20 $h^{-1}$ comoving Mpc. On such large scales, our cosmological simulations show a strong correlation between the effective optical depth ($τ_{\rm{eff}}$) of the CoSLAs and the 3-D mass overdensities. In moderate signal-to-noise spectra, however, the profiles of CoSLAs can be confused with high column density absorbers. For $z>2.6$, where the corresponding Ly$β$ is redshifted to the optical, we have developed the technique to differentiate between these two alternatives. We have applied this technique to SDSS-III quasar survey at $z = 2.6$ - 3.3, and we present a sample of five CoSLA candidates with $τ_{\rm{eff}}$ on 15 $h^{-1}$ Mpc greater than $4.5\times$ the mean optical depth. At lower redshifts of $z < 2.6$, where the background quasar density is higher, the overdensity can be traced by intergalactic absorption groups using multiple sight lines. Our overdensity searches fully utilize the current and next generation of Ly$α$ forest surveys which cover a survey volume of $> (1\ h^{-1}$ Gpc)$^3$. In addition, systems traced by CoSLAs will build a uniform sample of the most massive overdensities at $z > 2$ to constrain the models of structure formation, and offer a unique laboratory to study the interactions between galaxy overdensities and the intergalactic medium.

preprint2015arXiv

Metal-enriched, sub-kiloparsec gas clumps in the circumgalactic medium of a faint z = 2.5 galaxy

We report the serendipitous detection of a 0.2 L$^*$, Lyman-$α$ emitting galaxy at redshift 2.5 at an impact parameter of 50 kpc from a bright background QSO sightline. A high-resolution spectrum of the QSO reveals a partial Lyman-limit absorption system ($N_\mathrm{HI}=10^{16.94\pm0.10}$ cm$^{-2}$) with many associated metal absorption lines at the same redshift as the foreground galaxy. Using photoionization models that carefully treat measurement errors and marginalise over uncertainties in the shape and normalisation of the ionizing radiation spectrum, we derive the total hydrogen column density $N_\mathrm{H}=10^{19.4\pm0.3}$ cm$^{-2}$, and show that all the absorbing clouds are metal enriched, with $Z=0.1$-$0.6 Z_\odot$. These metallicities and the system's large velocity width ($436$ km$\,$s$^{-1}$) suggest the gas is produced by an outflowing wind. Using an expanding shell model we estimate a mass outflow rate of $\sim5 M_\odot\,$yr$^{-1}$. Our photoionization model yields extremely small sizes ($<$100-500 pc) for the absorbing clouds, which we argue are typical of high column density absorbers in the circumgalactic medium (CGM). Given these small sizes and extreme kinematics, it is unclear how the clumps survive in the CGM without being destroyed by hydrodynamic instabilities. The small cloud sizes imply that even state-of-the-art cosmological simulations require more than a $1000$-fold improvement in mass resolution to resolve the hydrodynamics relevant for cool gas in the CGM.

preprint2015arXiv

Quasar Quartet Embedded in Giant Nebula Reveals Rare Massive Structure in Distant Universe

All galaxies once passed through a hyperluminous quasar phase powered by accretion onto a supermassive black hole. But because these episodes are brief, quasars are rare objects typically separated by cosmological distances. In a survey for Lyman-alpha emission at redshift z ~ 2, we discovered a physical association of four quasars embedded in a giant nebula. Located within a substantial overdensity of galaxies, this system is probably the progenitor of a massive galaxy cluster. The chance probability of finding a quadruple quasar is estimated to be ~10^-7, implying a physical connection between Lyman-alpha nebulae and the locations of rare protoclusters. Our findings imply that the most massive structures in the distant universe have a tremendous supply (~ 10^11 solar masses) of cool dense (volume density ~1 cm^-3) gas, which is in conflict with current cosmological simulations.

preprint2015arXiv

Shadow of a Colossus: A z=2.45 Galaxy Protocluster Detected in 3D Ly-a Forest Tomographic Mapping of the COSMOS Field

Using moderate-resolution optical spectra from 58 background Lyman-break galaxies and quasars at $z\sim 2.3-3$ within a $11.5'\times13.5'$ area of the COSMOS field ($\sim 1200\,\mathrm{deg}^2$ projected area density or $\sim 2.4\,h^{-1}\,\mathrm{Mpc}$ mean transverse separation), we reconstruct a 3D tomographic map of the foreground Ly$α$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $σ_{3d}\approx3.5\,h^{-1}\,\mathrm{Mpc}$ comoving. Comparing with 61 coeval galaxies with spectroscopic redshifts in the same volume, we find that the galaxy positions are clearly biased towards regions with enhanced IGM absorption in the tomographic map. We find an extended IGM overdensity with deep absorption troughs at $z=2.45$ associated with a recently-discovered galaxy protocluster at the same redshift. Based on simulations matched to our data, we estimate the enclosed dark matter mass within this IGM overdensity to be $M_{\rm dm} (z=2.45) = (9\pm4)\times 10^{13}\,h^{-1}\,\mathrm{M_\odot}$, and argue based on this mass and absorption strength that it will form at least one $z\sim0$ galaxy cluster with $M(z=0) = (3\pm 2) \times 10^{14}\,h^{-1}\mathrm{M_\odot}$, although its elongated nature suggests that it will likely collapse into two separate clusters. We also point out a compact overdensity of six MOSDEF galaxies at $z=2.30$ within a $r\sim 1\,h^{-1}\,\mathrm{Mpc}$ radius and $Δz\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$α$ forest tomography on larger volumes to study galaxy properties as a function of environment, as well as revealing the large-scale IGM overdensities associated with protoclusters and other features of large-scale structure.

preprint2015arXiv

The Bright Symbiotic Mira EF Aquilae

An incidental spectrum of the poorly studied long period variable EF Aquilae shows [O III] emission indicative of a symbiotic star. Strong GALEX detections in the UV reinforce this classification, providing overt evidence for the presence of the hot subluminous companion. Recent compilations of the photometric behavior strongly suggest that the cool component is a Mira variable. Thus EF Aql appears to be a member of the rare symbiotic Mira subgroup.

preprint2015arXiv

The chemistry of the most metal-rich damped Lyman $α$ systems at z$\sim2$ II. Context with the Local Group

Using our sample of the most metal-rich damped Lyman $α$ systems (DLAs) at z$\sim2$, and two literature compilations of chemical abundances in 341 DLAs and 2818 stars, we present an analysis of the chemical composition of DLAs in the context of the Local Group. The metal-rich sample of DLAs at z$\sim2$ probes metallicities as high as the Galactic disc and the most metal-rich dwarf spheroidals (dSphs), permitting an analysis of many elements typically observed in DLAs (Fe, Zn, Cr, Mn, Si, and S) in comparison to stellar abundances observed in the Galaxy and its satellites (in particular dSphs). Our main conclusions are: (1) non-solar [Zn/Fe] abundances in metal-poor Galactic stars and in dSphs over the full metallicity range probed by DLAs, suggest that Zn is not a simple proxy for Fe in DLAs and therefore not a suitable indicator of dust depletion. After correcting for dust depletion, the majority of DLAs have subsolar [Zn/Fe] similar to dSphs; (2) at [Fe/H]$\sim-0.5$, a constant [Mn/Fe]$\sim-0.5$ and near-solar [$α$/Fe] (requiring an assumption about dust depletion) are in better agreement with dwarf galaxies than Galactic disc stars; (3) [$α$/Zn] is usually solar or subsolar in DLAs. However, although low ratios of [$α$/Fe] are usually considered more `dwarf-like' than `Milky Way-like', subsolar [Zn/Fe] in Local Group dwarfs leads to supersolar [$α$/Zn] in the dSphs, in contrast with the DLAs. Therefore, whilst DLAs exhibit some similarities with the Local Group dwarf population, there are also notable differences.

preprint2015arXiv

The Keck+Magellan Survey for Lyman Limit Absorption III: Sample Definition and Column Density Measurements

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.

preprint2015arXiv

The Neutral Hydrogen Cosmological Mass Density at z=5

We present the largest homogeneous survey of $z>4.4$ damped Lyman-$α$ 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, $Ω_{\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 $Ω_{\rm HI}$. We find $Ω_{\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 $Ω_{\rm HI}$ can be described by the functional form $Ω_{\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.

preprint2015arXiv

The physical properties of z>2 Lyman limit systems: new constraints for feedback and accretion models

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.

preprint2015arXiv

Towards the statistical detection of the warm-hot intergalactic medium in inter-cluster filaments of the cosmic web

Modern analyses of structure formation predict a universe tangled in a 'cosmic web' of dark matter and diffuse baryons. These theories further predict that at low-z, a significant fraction of the baryons will be shock-heated to $T \sim 10^{5}-10^{7}$K yielding a warm-hot intergalactic medium (WHIM), but whose actual existence has eluded a firm observational confirmation. We present a novel experiment to detect the WHIM, by targeting the putative filaments connecting galaxy clusters. We use HST/COS to observe a remarkable QSO sightline that passes within $Δd = 3$ Mpc from the 7 inter-cluster axes connecting 7 independent cluster-pairs at redshifts $0.1 \le z \le 0.5$. We find tentative excesses of total HI, narrow HI (NLA; Doppler parameters $b<50$ km/s), broad HI (BLA; $b \ge 50$ km/s) and OVI absorption lines within rest-frame velocities of $Δv \lesssim 1000$ km/s from the cluster-pairs redshifts, corresponding to $\sim 2$, $\sim 1.7$, $\sim 6$ and $\sim 4$ times their field expectations, respectively. Although the excess of OVI likely comes from gas close to individual galaxies, we conclude that most of the excesses of NLAs and BLAs are truly intergalactic. We find that the covering fractions, $f_c$, of BLAs close to cluster-pairs are $\sim 4-7$ times higher than the random expectation (at the $\sim 2 σ$ c.l.), whereas the $f_c$ of NLAs and OVI are not significantly enhanced. We argue that a larger relative excess of BLAs compared to those of NLAs close to cluster-pairs may be a signature of the WHIM in inter-cluster filaments. By extending the present analysis to tens of sightlines our experiment offers a promising route to detect the WHIM.

preprint2015arXiv

Using machine learning to classify the diffuse interstellar bands

Using over a million and a half extragalactic spectra we study the correlations of the Diffuse Interstellar Bands (DIBs) in the Milky Way. We measure the correlation between DIB strength and dust extinction for 142 DIBs using 24 stacked spectra in the reddening range E(B-V) < 0.2, many more lines than ever studied before. Most of the DIBs do not correlate with dust extinction. However, we find 10 weak and barely studied DIBs with correlations that are higher than 0.7 with dust extinction and confirm the high correlation of additional 5 strong DIBs. Furthermore, we find a pair of DIBs, 5925.9A and 5927.5A which exhibits significant negative correlation with dust extinction, indicating that their carrier may be depleted on dust. We use Machine Learning algorithms to divide the DIBs to spectroscopic families based on 250 stacked spectra. By removing the dust dependency we study how DIBs follow their local environment. We thus obtain 6 groups of weak DIBs, 4 of which are tightly associated with C2 or CN absorption lines.

preprint2014arXiv

A Budget and Accounting of Metals at z~0: Results from the COS-Halos Survey

We present a budget and accounting of metals in and around star-forming galaxies at $z\sim 0$. We combine empirically derived star formation histories with updated supernova and AGB yields and rates to estimate the total mass of metals produced by galaxies with present-day stellar mass of $10^{9.3}$--$10^{11.6} M_{\odot}$. On the accounting side of the ledger, we show that a surprisingly constant 20--25% mass fraction of produced metals remain in galaxies' stars, interstellar gas and interstellar dust, with little dependence of this fraction on the galaxy stellar mass (omitting those metals immediately locked up in remnants). Thus, the bulk of metals are outside of galaxies, produced in the progenitors of today's $L^*$ galaxies. The COS-Halos survey is uniquely able to measure the mass of metals in the circumgalactic medium (to impact parameters of $< 150$ kpc) of low-redshift $\sim L^*$ galaxies. Using these data, we map the distribution of CGM metals as traced by both the highly ionized OVI ion and a suite of low-ionization species; combined with constraints on circumgalactic dust and hotter X-ray emitting gas out to similar impact parameters, we show that $\sim 40$% of metals produced by $M_{\star}\sim 10^{10}M_{\odot}$ galaxies can be easily accounted for out to 150 kpc. With the current data, we cannot rule out a constant mass of metals within this fixed physical radius. This census provides a crucial boundary condition for the eventual fate of metals in galaxy evolution models.

preprint2014arXiv

A cosmic web filament revealed in Lyman-alpha emission around a luminous high-redshift quasar

Simulations of structure formation in the Universe predict that galaxies are embedded in a "cosmic web", where the majority of baryons reside as rarefied and highly ionized gas. This material has been studied for decades in absorption against background sources, but the sparseness of these inherently one-dimensional probes preclude direct constraints on the three-dimensional morphology of the underlying web. Here we report observations of a cosmic web filament in Lyman-alpha emission, discovered during a survey for cosmic gas fluorescently "illuminated" by bright quasars at z=2.3. With a projected size of approximately 460 physical kpc, the Lyman-alpha emission surrounding the radio-quiet quasar UM287 extends well beyond the virial radius of any plausible associated dark matter halo. The estimated cold gas mass of the nebula from the observed emission is at least ten times larger than what is typically found by cosmological simulations, suggesting that a population of intergalactic gas clumps with sub-kpc sizes may be missing within current numerical models.

preprint2014arXiv

A Deep Narrowband Imaging Search for CIV and He II Emission from Ly$α$ Blobs

We conduct a deep narrow-band imaging survey of 13 Ly$α$ blobs (LABs) located in the SSA22 proto-cluster at z~3.1 in the CIV and HeII emission lines in an effort to constrain the physical process powering the Ly$α$ emission in LABs. Our observations probe down to unprecedented surface brightness limits of 2.1 $-$ 3.4 $\times$ 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ per 1 arcsec$^2$ aperture (5$σ$) for the HeII$λ$1640 and CIV$λ$1549 lines, respectively. We do not detect extended HeII and CIV emission in any of the LABs, placing strong upper limits on the HeII/Ly$α$ and CIV/Ly$α$ line ratios, of 0.11 and 0.16, for the brightest two LABs in the field. We conduct detailed photoionization modeling of the expected line ratios and find that, although our data constitute the deepest ever observations of these lines, they are still not deep enough to rule out a scenario where the Ly$α$ emission is powered by the ionizing luminosity of an obscured AGN. Our models can accommodate HeII/Ly$α$ and CIV/Ly$α$ ratios as low as $\simeq$0.05 and $\simeq$0.07 respectively, implying that one needs to reach surface brightness as low as 1 $-$ 1.5 $\times$ 10$^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ (at 5$σ$) in order to rule out a photoionization scenario. These depths will be achievable with the new generation of image-slicing integral field units such as VLT/MUSE or Keck/KCWI. We also model the expected HeII/Ly$α$ and CIV/Ly$α$ in a different scenario, where Ly$α$ emission is powered by shocks generated in a large-scale superwind, but find that our observational constraints can only be met for shock velocities $v_{\rm s} \gtrsim$ 250 km s$^{-1}$, which appear to be in conflict with recent observations of quiescent kinematics in LABs.

preprint2014arXiv

A Trio of GRB-SNe: GRB 120729A, GRB 130215A / SN 2013ez and GRB 130831A / SN 2013fu

We present optical and near-infrared (NIR) photometry for three gamma-ray burst supernovae (GRB-SNe): GRB 120729A, GRB 130215A / SN 2013ez and GRB 130831A / SN 2013fu. In the case of GRB 130215A / SN 2013ez, we also present optical spectroscopy at t-t0=16.1 d, which covers rest-frame 3000-6250 Angstroms. Based on Fe II (5169) and Si (II) (6355), our spectrum indicates an unusually low expansion velocity of 4000-6350 km/s, the lowest ever measured for a GRB-SN. Additionally, we determined the brightness and shape of each accompanying SN relative to a template supernova (SN 1998bw), which were used to estimate the amount of nickel produced via nucleosynthesis during each explosion. We find that our derived nickel masses are typical of other GRB-SNe, and greater than those of SNe Ibc that are not associated with GRBs. For GRB 130831A / SN 2013fu, we use our well-sampled R-band light curve (LC) to estimate the amount of ejecta mass and the kinetic energy of the SN, finding that these too are similar to other GRB-SNe. For GRB 130215A, we take advantage of contemporaneous optical/NIR observations to construct an optical/NIR bolometric LC of the afterglow. We fit the bolometric LC with the millisecond magnetar model of Zhang & Meszaros (2001), which considers dipole radiation as a source of energy injection to the forward shock powering the optical/NIR afterglow. Using this model we derive an initial spin period of P=12 ms and a magnetic field of B=1.1 x 10^15 G, which are commensurate with those found for proposed magnetar central engines of other long-duration GRBs.

preprint2014arXiv

Constraints on the gas masses of low-z damped Lyman-$α$ systems

We report a deep search for redshifted HI 21 cm emission from three damped and sub-damped Lyman-$α$ absorbers (DLAs) at $z \approx 0.1$ with the Green Bank Telescope (GBT). No evidence for a redshifted HI 21 cm emission signal was obtained in the GBT spectra of two absorbers, with the data on the third rendered unusable by terrestrial interference. The non-detections of HI 21 cm emission yield strong constraints on the HI masses of the associated galaxies, M$_{\rm HI} < 2.3 \times 10^9 \times (ΔV/100)^{1/2}$ M$_\odot$ for the sub-DLA at $z = 0.0830$ towards J1553+3548, and M$_{\rm HI} < 2.7 \times 10^9 \times (ΔV/100)^{1/2}$ M$_\odot$ for the DLA at $z = 0.0963$ towards J1619+3342, where $ΔV$ is the HI 21 cm line width, in km s$^{-1}$. This continues the trend of low HI masses found in all low-$z$ DLAs and sub-DLAs that have been searched for redshifted HI 21 cm emission. Low-redshift absorbers with relatively low HI column densities, $\lesssim few \times 10^{20}$ cm$^{-2}$, thus do not typically arise in massive gas-rich galaxies.

preprint2014arXiv

Directly imaging damped Ly-alpha galaxies at z>2. II: Imaging and spectroscopic observations of 32 quasar fields

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.

preprint2014arXiv

Directly imaging damped Ly-alpha galaxies at z>2. III: The star formation rates of neutral gas reservoirs at z~2.7

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.

preprint2014arXiv

Dissecting the Gaseous Halos of z~2 Damped Ly$α$ Systems with Close Quasar Pairs

We use spectroscopy of close pairs of quasars to study diffuse gas in the circumgalactic medium (CGM) surrounding a sample of 40 Damped Lya systems (DLAs). The primary sightline in each quasar pair probes an intervening DLA in the redshift range 1.6 < z_DLA < 3.6, such that the second quasar sightline then probes Lya, CII, SiII, and CIV absorption in the CGM transverse to the DLA to projected distances $R_{\perp} < 300$ kpc. Analysis of the Lya profiles in these CGM sightlines constrains the covering fraction (f_C) of optically thick HI (having column density N_HI > 10^17.2 cm^-2) to be greater than ~30% within $R_{\perp} < 200$ kpc of DLAs. Strong SiII 1526 absorption with equivalent width W_1526 > 0.2 Ang occurs with an incidence f_C (W_1526 > 0.2 Ang) = 20(+12/-8)% within $R_{\perp}<100$ kpc, indicating that low-ionization metal absorption associated with DLAs probes material at a physical distance R_3D < 30 kpc. However, we find that strong CIV 1548 absorption is ubiquitous in these environments (f_C (W_1548 > 0.2 Ang) = 57(+12/-13)% within $R_{\perp} < 100$ kpc), and in addition exhibits a high degree of kinematic coherence on scales up to ~175 kpc. We infer that this high-ionization material arises predominantly in large, quiescent structures extending beyond the scale of the DLA host dark matter halos rather than in ongoing galactic winds. The Lya equivalent width in the DLA-CGM is anticorrelated with $R_{\perp}$ at >98% confidence, suggesting that DLAs arise close to the centers of their host halos rather than on their outskirts. Finally, the average Lya, CII and CIV equivalent widths are consistent with those measured around z~2 Lyman Break Galaxies. Assuming that DLAs trace a galaxy population with lower masses and luminosities, this finding implies that the absorption strength of cool circumgalactic material has a weak dependence on dark matter halo mass for M_h < 10^12 M_sun.

preprint2014arXiv

Dusting off the diffuse interstellar bands: DIBs and dust in extragalactic SDSS spectra

Using over a million and a half extragalactic spectra we study the properties of the mysterious Diffuse Interstellar Bands (DIBs) in the Milky Way. These data provide us with an unprecedented sampling of the skies at high Galactic-latitude and low dust-column-density. We present our method, study the correlation of the equivalent width of 8 DIBs with dust extinction and with a few atomic species, and the distribution of four DIBs - 5780.6A, 5797.1A, 6204.3A, and 6613.6A - over nearly 15000 squared degrees. As previously found, DIBs strengths correlate with extinction and therefore inevitably with each other. However, we show that DIBs can exist even in dust free areas. Furthermore, we find that the DIBs correlation with dust varies significantly over the sky. DIB under- or over-densities, relative to the expectation from dust, are often spread over hundreds of square degrees. These patches are different for the four DIBs, showing that they are unlikely to originate from the same carrier, as previously suggested.

preprint2014arXiv

Galactic and Circumgalactic OVI and its Impact on the Cosmological Metal and Baryon Budgets at 2<z<3.5

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 (τ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 α, β forests. We found 20 absorbers that satisfy these rules: 7 Lyman limit systems (LLSs), 8 super-LLSs (SLLSs) and 5 damped Lyα (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 τ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 τLL > 1 absorbers could contain as much as 3-14% of the cosmic baryon budget at z~2-3, only second to the Lyα 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.

preprint2014arXiv

Lyman-alpha Forest Tomography from Background Galaxies: The First Megaparsec-Resolution Large-Scale Structure Map at z>2

We present the first observations of foreground Lyman-$α$ forest absorption from high-redshift galaxies, targeting 24 star-forming galaxies (SFGs) with $z\sim 2.3-2.8$ within a $5' \times 15'$ region of the COSMOS field. The transverse sightline separation is $\sim 2\,h^{-1}\mathrm{Mpc}$ comoving, allowing us to create a tomographic reconstruction of the 3D Ly$α$ forest absorption field over the redshift range $2.20\leq z\leq 2.45$. The resulting map covers $6\,h^{-1}\mathrm{Mpc} \times 14\,h^{-1}\mathrm{Mpc}$ in the transverse plane and $230\,h^{-1}\mathrm{Mpc}$ along the line-of-sight with a spatial resolution of $\approx 3.5\,h^{-1}\mathrm{Mpc}$, and is the first high-fidelity map of large-scale structure on $\sim\mathrm{Mpc}$ scales at $z>2$. Our map reveals significant structures with $\gtrsim 10\,h^{-1}\mathrm{Mpc}$ extent, including several spanning the entire transverse breadth, providing qualitative evidence for the filamentary structures predicted to exist in the high-redshift cosmic web. Simulated reconstructions with the same sightline sampling, spectral resolution, and signal-to-noise ratio recover the salient structures present in the underlying 3D absorption fields. Using data from other surveys, we identified 18 galaxies with known redshifts coeval with our map volume enabling a direct comparison to our tomographic map. This shows that galaxies preferentially occupy high-density regions, in qualitative agreement with the same comparison applied to simulations. Our results establishes the feasibility of the CLAMATO survey, which aims to obtain Ly$α$ forest spectra for $\sim 1000$ SFGs over $\sim 1 \,\mathrm{deg}^2$ of the COSMOS field, in order to map out IGM large-scale structure at $\langle z \rangle \sim 2.3$ over a large volume $(100\,h^{-1}\mathrm{Mpc})^3$.

preprint2014arXiv

Precious Metals in SDSS Quasar Spectra II: Tracking the Evolution of Strong, 0.4 < z < 2.3 MgII Absorbers with Thousands of Systems

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

preprint2014arXiv

Probing the Physical Conditions of Atomic Gas at High Redshift

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

preprint2014arXiv

Quasars Probing Quasars VII. The Pinnacle of the Cool Circumgalactic Medium Surrounds Massive z~2 Galaxies

We survey the incidence and absorption strength of the metal-line transitions CII 1334 and CIV from the circumgalactic medium (CGM) surrounding z~2 quasars, which act as signposts for massive dark matter halos M_halo~10^12.5 Msun. On scales of the virial radius (Mvir~160kpc), we measure a high covering fraction fC=0.73+/-0.10 to strong CII absorption (rest equivalent width W1334>0.2A), implying a massive reservoir of cool (T~10^4K) metal enriched gas. We conservatively estimate a metal mass exceeding 10^8 Msun. We propose these metals trace enrichment of the incipient intragroup/intracluster medium that these halos eventually inhabit. This cool CGM around quasars is the pinnacle amongst galaxies observed at all epochs, as regards covering fraction and average equivalent width of HI Lya and low-ion metal absorption. We argue that the properties of this cool CGM primarily reflect the halo mass, and that other factors such as feedback, star-formation rate, and accretion from the intergalactic medium are secondary. We further estimate, that the CGM of massive, z~2 galaxies accounts for the majority of strong MgII absorption along random quasar sightlines. Lastly, we detect an excess of strong CIV absorption (W1548>0.3A) over random incidence to 1Mpc physical impact parameter and measure the quasar-CIV cross-correlation function: xi(r)=(r/r0)^-g with r0 = 7.5Mpc and g=1.7. Consistent with previous work on larger scales, we infer that this highly ionized CIV gas traces massive (10^12 Msun) halos.

preprint2014arXiv

The COS-Dwarfs Survey: The Carbon Reservoir Around sub-L* Galaxies

We report new observations of circumgalactic gas from the COS-Dwarfs survey, a systematic investigation of the gaseous halos around 43 low-mass z $\leq$ 0.1 galaxies using background QSOs observed with the Cosmic Origins Spectrograph. From the projected 1D and 2D distribution of C IV absorption, we find that C IV absorption is detected out to ~ 0.5 R$_{vir}$ of the host galaxies. The C IV absorption strength falls off radially as a power law and beyond 0.5 R$_{vir}$, no C IV absorption is detected above our sensitivity limit of ~ 50-100 m$Å$. We find a tentative correlation between detected C IV absorption strength and star formation, paralleling the strong correlation seen in highly ionized oxygen for L~L* galaxies by the COS-Halos survey. The data imply a large carbon reservoir in the CGM of these galaxies, corresponding to a minimum carbon mass of $\gtrsim$ 1.2$\times 10^6$ $M_\odot$ out to ~ 110 kpc. This mass is comparable to the carbon mass in the ISM and more than the carbon mass currently in stars of these galaxies. The C IV absorption seen around these sub-L* galaxies can account for almost two-thirds of all $W_r$> 100 m$Å$ C IV absorption detected at low z. Comparing the C IV covering fraction with hydrodynamical simulations, we find that an energy-driven wind model is consistent with the observations whereas a wind model of constant velocity fails to reproduce the CGM or the galaxy properties.

preprint2014arXiv

The COS-Halos Survey: Physical Conditions and Baryonic Mass in the Low-Redshift Circumgalactic Medium

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.

preprint2014arXiv

The Giant Gemini GMOS survey of z>4.4 quasars - I. Measuring the mean free path across cosmic time

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 ($λ_\mathrm{r}<912$Å) to assess the mean free path $λ_\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 $λ_\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 $λ_\mathrm{mfp}^{912}$ measurements from lower redshifts, the data are well modelled by a simple power-law $λ_\mathrm{mfp}^{912}=A[(1+z)/5]^η$ with $A=(37\pm 2)h_{70}^{-1}$ Mpc and $η= -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 $λ_\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.

preprint2014arXiv

The most metal-rich damped Lyman alpha systems at z>1.5 I: The Data

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

preprint2014arXiv

The Rapid Decline in Metallicity of Damped Ly-alpha Systems at z~5

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

preprint2014arXiv

The Role of Stellar Feedback in the Dynamics of HII Regions

Stellar feedback is often cited as the biggest uncertainty in galaxy formation models today. This uncertainty stems from a dearth of observational constraints as well as the great dynamic range between the small scales (<1 pc) where the feedback occurs and the large scales of galaxies (>1 kpc) that are shaped by this feedback. To bridge this divide, in this paper we aim to assess observationally the role of stellar feedback at the intermediate scales of HII regions. In particular, we employ multiwavelength data to examine several stellar feedback mechanisms in a sample of 32 HII regions in the Large and Small Magellanic Clouds (LMC and SMC, respectively). Using optical, infrared, radio, and X-ray images, we measure the pressures exerted on the shells from the direct stellar radiation, the dust-processed radiation, the warm ionized gas, and the hot X-ray emitting gas. We find that the warm ionized gas dominates over the other terms in all of the sources, although two have comparable dust-processed radiation pressures to their warm gas pressures. The hot gas pressures are comparatively weak, while the direct radiation pressures are 1-2 orders of magnitude below the other terms. We discuss the implications of these results, particularly highlighting evidence for hot gas leakage from the HII shells and regarding the momentum deposition from the dust-processed radiation to the warm gas. Furthermore, we emphasize that similar observational work should be done on very young HII regions to test whether direct radiation pressure and hot gas can drive the dynamics at early times.

preprint2014arXiv

Towards a Unified Description of the Intergalactic Medium at Redshift z~2.5

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.

preprint2013arXiv

A Deep Search for Faint Galaxies Associated with Very Low-Redshift C IV Absorbers: A Case with Cold-Accretion Characteristics

Studies of QSO absorber-galaxy connections are often hindered by inadequate information on whether faint/dwarf galaxies are located near the QSO sight lines. To investigate the contribution of faint galaxies to QSO absorber populations, we are conducting a deep galaxy redshift survey near low-z C IV absorbers. Here we report a blindly-detected C IV absorption system (z(abs) = 0.00348) in the spectrum of PG1148+549 that appears to be associated either with an edge-on dwarf galaxy with an obvious disk (UGC 6894, z(gal) = 0.00283) at an impact parameter of rho = 190 kpc or with a very faint dwarf irregular galaxy at rho = 23 kpc, which is closer to the sightline but has a larger redshift difference (z(gal) = 0.00107, i.e., dv = 724 km/s). We consider various gas/galaxy associations, including infall and outflows. Based on current theoretical models, we conclude that the absorber is most likely tracing (1) the remnants of an outflow from a previous epoch, a so-called 'ancient outflow', or (2) intergalactic gas accreting onto UGC 6894, 'cold mode' accretion. The latter scenario is supported by H I synthesis imaging data that shows the rotation curve of the disk being codirectional with the velocity offset between UGC 6894 and the absorber, which is located almost directly along the major axis of the edge-on disk.

preprint2013arXiv

A high molecular fraction in a sub-damped absorber at z=0.56

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.

preprint2013arXiv

A Refined Measurement of the Mean Transmitted Flux in the Ly-alpha Forest over 2 < z < 5 Using Composite Quasar Spectra

We present new measurements of the mean transmitted flux in the Ly-alpha forest over 2 < z < 5 made using 6065 quasar spectra from the Sloan Digital Sky Survey DR7. We exploit the general lack of evolution in the mean quasar continuum to avoid the bias introduced by continuum fitting over the Ly-alpha forest at high redshifts, which has been the primary systematic uncertainty in previous measurements of the mean Ly-alpha transmission. The individual spectra are first combined into twenty-six composites with mean redshifts spanning 2.25 < z_comp < 5.08. The flux ratios of separate composites at the same rest wavelength are then used, without continuum fitting, to infer the mean transmitted flux, F(z), as a fraction of its value at z~2. Absolute values for F(z) are found by scaling our relative values to measurements made from high-resolution data by Faucher-Giguere et al. (2008) at z < 2.5, where continuum uncertainties are minimal. We find that F(z) evolves smoothly with redshift, with no evidence of a previously reported feature at z~3.2. This trend is consistent with a gradual evolution of the ionization and thermal state of the intergalactic medium over 2 < z < 5. Our results generally agree with the most careful measurements to date made from high-resolution data, but offer much greater precision and extend to higher redshifts. This work also improves upon previous efforts using SDSS spectra by significantly reducing the level of systematic error.

preprint2013arXiv

A Search for Boron in Damped Lyα Systems

We present the first systematic study of boron beyond the Local Group. This analysis is performed on a sample of 30 damped Lyα systems (DLAs) with strong metal-lines, which are expected to trace the interstellar medium of high z galaxies. We report on two boron detections at > 3σ significance; one new detection and one confirmation. The ratios of B/O and, for the first time, B/S are compared with previous stellar and interstellar measurements in the Milky Way and Small Magellanic Cloud. The novel comparison with sulphur, which tracks oxygen's abundance, alleviates the uncertainty associated with stellar oxygen measurements. For both detections, the inferred B/S ratio is in excess of the prediction of primary boron production from spallation processes. Possible sources of contamination are discussed, as well as physical effects that could impact the observed ratios. However taken at face value, the implication of these measurements suggest potentially higher cosmic ray fluxes in DLAs. The prospects for future boron detections in other high redshift DLAs to confirm our results is also discussed.

preprint2013arXiv

An Explanation for the Different X-ray to Optical Column Densities in the Environments of Gamma Ray Bursts: A Progenitor Embedded in a Dense Medium

We study the > 10 ratios in the X-ray to optical column densities inferred from afterglow spectra of Gamma Ray Bursts due to gas surrounding their progenitors. We present time-evolving photoionization calculations for these afterglows and explore different conditions for their environment. We find that homogenous models of the environment (constant density) predict X-ray columns similar to those found in the optical spectra, with the bulk of the opacity being produced by neutral material at large distances from the burst. This result is independent of gas density or metallicity. Only models assuming a progenitor immersed in a dense (10^(2-4) cm-3) cloud of gas (with radius ~10 pc), with a strong, declining gradient of density for the surrounding interstellar medium are able to account for the large X-ray to optical column density ratios. However, to avoid an unphysical correlation between the size of this cloud, and the size of the ionization front produced by the GRB, the models also require that the circumburst medium is already ionized prior to the burst. The inferred cloud masses are <10^6 M_solar, even if low metallicities in the medium are assumed (Z~0.1 Z_solar). These cloud properties are consistent with those found in giant molecular clouds and our results support a scenario in which the progenitors reside within intense star formation regions of galaxies. Finally, we show that modeling over large samples of GRB afterglows may offer strong constraints on the range of properties in these clouds, and the host galaxy ISM.

preprint2013arXiv

Confronting Simulations of Optically Thick Gas in Massive Halos with Observations at z=2-3

Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high resolution "zoom-in" simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M_vir~ 2x10^11 - 4x10^12 M_sun at z~2. We find that both the mass fraction of cool (T <= 3x10^4 K) gas and the covering fraction of optically-thick Lyman limit systems (LLSs) depend only weakly on halo mass, even around the critical value for the formation of stable virial shocks. The covering fraction of LLSs interior to the virial radius varies between f_c~0.05 - 0.2, with significant scatter among halos. Our simulations of massive halos (M_vir >= 10^12 M_sun) underpredict the covering fraction of optically-thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M_vir >= 5x10^11 M_sun) hosting star-forming galaxies, the predicted covering factor agrees with observations, however current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale auto-correlation function of optically-thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z~2-3, as well as the characteristic size and covering factor of the CGM.

preprint2013arXiv

Discovery of a transparent sightline at ρ < 20 kpc from an interacting pair of galaxies

We report the discovery of a transparent sightline at projected distances of ρ < 20 kpc to an interacting pair of mature galaxies at z = 0.12. The sightline of the UV-bright quasar PG1522+101 at z = 1.328 passes at ρ = 11.5 kpc from the higher-mass galaxy (M_* = 10^10.6 M_Sun) and ρ = 20.4 kpc from the lower-mass one (M_* = 10^10.0 M_Sun). The two galaxies are separated by 9 kpc in projected distance and 30 km/s in line-of-sight velocity. Deep optical images reveal tidal features indicative of close interactions. Despite the small projected distances, the quasar sightline shows little absorption associated with the galaxy pair with a total HI column density it no greater than log N(HI) = 13.65. This limiting HI column density is already two orders-of-magnitude less than what is expected from previous halo gas studies. In addition, we detect no heavy-element absorption features associated with the galaxy pair with 3-σ limits of log N(MgII) < 12.2 and log N(OVI) < 13.7. The probability of seeing such little absorption in a sightline passing at a small projected distance from two non-interacting galaxies is 0.2%. The absence of strong absorbers near the close galaxy pair suggests that the cool gas reservoirs of the galaxies have been significantly depleted by the galaxy interaction. These observations therefore underscore the potential impact of galaxy interactions on the gaseous halos around galaxies.

preprint2013arXiv

Dissecting the properties of optically-thick hydrogen at the peak of cosmic star formation history

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.

preprint2013arXiv

Evidence for Ubiquitous Collimated Galactic-Scale Outflows along the Star-Forming Sequence at z~0.5

We present an analysis of the MgII 2796, 2803 and FeII 2586, 2600 absorption line profiles in individual spectra of 105 galaxies at 0.3<z<1.4. The galaxies, drawn from redshift surveys of the GOODS fields and the Extended Groth Strip, fully sample the range in star formation rates (SFRs) occupied by the star-forming sequence with stellar masses log M_*/M_sun > 9.5 at 0.3<z<0.7. Using the Doppler shifts of the MgII and FeII absorption lines as tracers of cool gas kinematics, we detect large-scale winds in 66+/-5% of the galaxies. HST/ACS imaging and our spectral analysis indicate that the outflow detection rate depends primarily on galaxy orientation: winds are detected in ~89% of galaxies having inclinations (i) <30 degrees (face-on), while the wind detection rate is only ~45% in objects having i>50 degrees (edge-on). Combined with the comparatively weak dependence of the wind detection rate on intrinsic galaxy properties, this suggests that biconical outflows are ubiquitous in normal, star-forming galaxies at z~0.5. We find that the wind velocity is correlated with host galaxy M_* at 3.4-sigma significance, while the equivalent width of the flow is correlated with host galaxy SFR at 3.5-sigma significance, suggesting that hosts with higher SFR may launch more material into outflows and/or generate a larger velocity spread for the absorbing clouds. Assuming that the gas is launched into dark matter halos with simple, isothermal density profiles, the wind velocities measured for the bulk of the cool material (~200-400 km/s) are sufficient to enable escape from the halo potentials only for the lowest-M_* systems in the sample. However, the outflows typically carry sufficient energy to reach distances of >50 kpc, and may therefore be a viable source of cool material for the massive circumgalactic medium observed around bright galaxies at z~0. [abridged]

preprint2013arXiv

Metal-Poor, Cool Gas in the Circumgalactic Medium of a z = 2.4 Star-Forming Galaxy: Direct Evidence for Cold Accretion?

In our current galaxy formation paradigm, high-redshift galaxies are predominantly fuelled by accretion of cool, metal-poor gas from the intergalactic medium. Hydrodynamical simulations predict that this material should be observable in absorption against background sightlines within a galaxy's virial radius, as optically thick Lyman-limit systems (LLSs) with low metallicities. Here we report the discovery of exactly such a strong metal-poor absorber at an impact parameter R_perp = 58 kpc from a star-forming galaxy at z = 2.44. Besides strong neutral hydrogen [N(HI) = 10^(19.50 +/- 0.16) cm^-2] we detect neutral deuterium and oxygen, allowing a precise measurement of the metallicity: log10(Z / Zsolar) = -2.0 +/- 0.17, or (7-15) x 10^-3 solar. Furthermore, the narrow deuterium linewidth requires a cool temperature < 20,000 K. Given the striking similarities between this system and the predictions of simulations, we argue that it represents the direct detection of a high redshift cold-accretion stream. The low-metallicity gas cloud is a single component of an absorption system exhibiting a complex velocity, ionization, and enrichment structure. Two other components have metallicities > 0.1 solar, ten times larger than the metal-poor component. We conclude that the photoionized circumgalactic medium (CGM) of this galaxy is highly inhomogeneous: the majority of the gas is in a cool, metal-poor and predominantly neutral phase, but the majority of the metals are in a highly-ionized phase exhibiting weak neutral hydrogen absorption but strong metal absorption. If such inhomogeneity is common, then high-resolution spectra and detailed ionization modelling are critical to accurately appraise the distribution of metals in the high-redshift CGM.

preprint2013arXiv

Precious Metals in SDSS Quasar Spectra I: Tracking the Evolution of Strong, 1.5 < z < 4.5 CIV Absorbers with Thousands of Systems

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 Å. 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 Å. This suggests a monotonic and significant increase in the enrichment of gas outside galaxies over the 12 Gyr lifetime of the universe.

preprint2013arXiv

QSO Absorption Systems Detected in Ne VIII: High-Metallicity Clouds with a Large Effective Cross Section

Using high resolution, high signal-to-noise ultraviolet spectra of the z = 0.9754 quasar PG1148+549 obtained with the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope, we study the physical conditions and abundances of NeVIII+OVI absorption line systems at z(abs) =0.68381, 0.70152, 0.72478. In addition to NeVIII and OVI, absorption lines from multiple ionization stages of oxygen (OII, OIII, OIV) are detected and are well-aligned with the more highly ionized species. We show that these absorbers are multiphase systems including hot gas (T ~ 10^{5.7} K) that produces NeVIII and OVI, and the gas metallicity of the cool phase ranges from Z = 0.3 Z_{solar} to supersolar. The cool (~10^{4} K) phases have densities n_{H} ~ 10^{-4} cm^{-3} and small sizes (< 4kpc); these cool clouds are likely to expand and dissipate, and the NeVIII may be within a transition layer between the cool gas and a surrounding, much hotter medium. The NeVIII redshift density, dN/dz = 7^{+7}_{-3}, requires a large number of these clouds for every L > 0.1L* galaxy and a large effective absorption cross section (>~ 100 kpc), and indeed, we find a star forming ~L* galaxy at the redshift of the z(abs)=0.72478 system, at an impact parameter of 217 kpc. Multiphase absorbers like these NeVIII systems are likely to be an important reservoir of baryons and metals in the circumgalactic media of galaxies.

preprint2013arXiv

Quasars Probing Quasars IV: Joint Constraints on the Circumgalactic Medium from Absorption and Emission

We have constructed a sample of 29 close projected quasar pairs where the background quasar spectrum reveals absorption from optically thick HI gas associated with the foreground quasar. These unique sightlines allow us to study the quasar circumgalactic medium (CGM) in absorption and emission simultaneously, because the background quasar pinpoints large concentrations of gas where Ly-a emission, resulting from quasar-powered fluorescence, resonant Ly-a scattering, and/or cooling radiation, is expected. A sensitive slit-spectroscopic search (1-sigma limits of SB_Lya ~= 3e-18 erg/s/cm^2/arcsec^2) for diffuse Ly-a emission in the environments of the foreground quasars is conducted. We fail to detect large-scale ~ 100 kpc Ly-a emission, either at the location of the optically thick absorbers or in the foreground quasar halos, in all cases except a single system. We interpret these non-detections as evidence that the gas detected in absorption is shadowed from the quasar UV radiation due to obscuration effects, which are frequently invoked in unified models of AGN. Small-scale R_perp <~ 50 kpc extended Ly-a nebulosities are detected in 34% of our sample, which are likely the high-redshift analogs of the extended emission-line regions commonly observed around low-redshift (z < 0.5) quasars. We also detect a compact high rest-frame equivalent width (W_Lya > 50 A) Ly-alpha-emitter with luminosity L_Lya =2.1+-0.32e41 erg/s at small impact parameter R_perp=134 kpc from one foreground quasar, and argue that it is more likely to result from quasar-powered fluorescence, than simply be a star-forming galaxy clustered around the quasar. Our observations imply that much deeper integrations with upcoming integral-field spectrometers such as MUSE and KCWI will be able to routinely detect a diffuse Ly-a glow around bright quasars on scales R ~ 100 kpc and thus directly image the CGM. [abridged]

preprint2013arXiv

Quasars Probing Quasars VI. Excess HI Absorption Within One Proper Mpc of z~2 Quasars

With close pairs of quasars at different redshifts, a background quasar sightline can be used to study a foreground quasar's environment in absorption. We use a sample of 650 projected quasar pairs to study the HI Lya absorption transverse to luminous, z~2 quasars at proper separations of 30kpc < R < 1Mpc. In contrast to measurements along the line-of-sight, regions transverse to quasars exhibit enhanced HI Lya absorption and a larger variance than the ambient intergalactic medium, with increasing absorption and variance toward smaller scales. Analysis of composite spectra reveals excess absorption characterized by a Lya equivalent width profile W = 2.3A (R/100kpc)^-0.46. We also observe a high (~60%) covering factor of strong, optically thick HI absorbers (HI column log NHI > 17.3) at separations R<200kpc, which decreases to ~20% at R~1Mpc, but still represents a significant excess over the cosmic average. This excess of optically thick absorption can be described by a quasar-absorber cross-correlation function xi_QA(r) = (r/r_0)^gamma with a large correlation length r_0 = 12.5+2.7-1.4 Mpc/h (comoving) and gamma = 1.68+0.14-0.30. The HI absorption measured around quasars exceeds that of any previously studied population, consistent with quasars being hosted by massive dark matter halos Mhalo~10^12.5 Msun at z~2.5. The environments of these massive halos are highly biased towards producing optically thick gas, and may even dominate the cosmic abundance of Lyman limit systems and hence the intergalactic opacity to ionizing photons at z~2.5. The anisotropic absorption around quasars implies the transverse direction is much less likely to be illuminated by ionizing radiation than the line-of-sight, which we interpret in terms of the same obscuration effects frequently invoked in unified models of active galactic nuclei.

preprint2013arXiv

The circumgalactic medium of massive galaxies at z~3: a test for stellar feedback, galactic outflows, and cold streams

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift 3, using the "Eris" suite of cosmological "zoom-in" simulations. The reference run adopts a blastwave scheme for supernova feedback that produces galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. Synthetic spectra through the multiphase CGM produce interstellar absorption line strengths of Lya, CII, CIV, SiII, and SiIV as a function of galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high-redshift by Steidel et al. (2010). Only about one third of all the gas within R_vir is outflowing. The fraction of sightlines within one virial radius that intercept optically thick material is 27%, in agreement with recent observations by Rudie et al. (2012). Such optically thick absorption is shown to trace inflowing "cold" streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar, give origin to strong (log N > 13) CII absorption with a covering factor of 22% within R_vir and 10% within 2 R_vir. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large OVI halo, with a typical column density log N>14 and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. (2011). Our zoom-in simulations of this single system appear to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings. (Abridged)

preprint2013arXiv

The COS-Halos Survey: Rationale, Design, and A Census of Circumgalactic Neutral Hydrogen

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.

preprint2013arXiv

The Fundamental Plane of Damped Lyman Alpha Systems

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

preprint2013arXiv

The High-Ion Content and Kinematics of Low-Redshift Lyman Limit Systems

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.

preprint2013arXiv

The HST/ACS+WFC3 Survey for Lyman Limit Systems II: Science

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

preprint2012arXiv

A search of CO emission lines in blazars: the low molecular gas content of BL Lac objects compared to quasars

BL Lacertae (Lac) objects that are detected at very-high energies (VHE) are of fundamental importance to study multiple astrophysical processes, including the physics of jets, the properties of the extragalactic background light and the strength of the intergalactic magnetic field. Unfortunately, since most blazars have featureless optical spectra that preclude a redshift determination, a substantial fraction of these VHE extragalactic sources cannot be used for cosmological studies. To assess whether molecular lines are a viable way to establish distances, we have undertaken a pilot program at the IRAM 30m telescope to search for CO lines in three BL Lac objects with known redshifts. We report a positive detection of M_H2 ~ 3x10^8 Msun toward 1ES 1959+650, but due to the poor quality of the baseline, this value is affected by a large systematic uncertainty. For the remaining two sources, W Comae and RGB J0710+591, we derive 3sigma upper limits at, respectively, M_H2 < 8.0x10^8 Msun and M_H2 < 1.6x10^9 Msun, assuming a line width of 150 km/s and a standard conversion factor alpha=4 M_sun/(K km/s pc^2). If these low molecular gas masses are typical for blazars, blind redshift searches in molecular lines are currently unfeasible. However, deep observations are still a promising way to obtain precise redshifts for sources whose approximate distances are known via indirect methods. Our observations further reveal a deficiency of molecular gas in BL Lac objects compared to quasars, suggesting that the host galaxies of these two types of active galactic nuclei (AGN) are not drawn from the same parent population. Future observations are needed to assess whether this discrepancy is statistically significant, but our pilot program shows how studies of the interstellar medium in AGN can provide key information to explore the connection between the active nuclei and the host galaxies.

preprint2012arXiv

A Substantial Mass of Cool, Metal-Enriched Gas Surrounding the Progenitors of Modern-Day Ellipticals

The hosts of luminous z~2 quasars evolve into today's massive elliptical galaxies. Current theories predict that the circum-galactic medium (CGM) of these massive, dark-matter halos (M~10^12.5 Msun) should be dominated by a T~10^7 K virialized plasma. We test this hypothesis with observations of 74 close-projected quasar pairs, using spectra of the background QSO to characterize the CGM of the foreground one. Surprisingly, our measurements reveal a cool (T~10^4 K), massive (M_CGM > 10^10 Msun), and metal-enriched (Z > ~0.1 Zsun) medium extending to at least the expected virial radius (r_vir = 160 kpc). The average equivalent widths of HI Lya (<W_lya> = 2.1 pm 0.15Ang for impact parameters R<200 kpc) and CII 1334 (<W_1334> = 0.7 pm 0.1Ang) exceed the corresponding CGM measurements of these transitions from all galaxy populations studied previously. Furthermore, we conservatively estimate that the quasar CGM has a 64% covering fraction of optically thick gas (N_HI>10^17.2) within r_vir; this covering factor is twice that of the contemporaneous Lyman Break Galaxy population. This unexpected reservoir of cool gas is rarely detected "down-the-barrel" to quasars, and hence it is likely that our background sightlines intercept gas which is shadowed from the quasar ionizing radiation by the same obscuring medium often invoked in models of AGN unification. Because the high-z halos inhabited by quasars predate modern groups and clusters, these observations are also relevant to the formation and enrichment history of the intragroup/intracluster medium.

preprint2012arXiv

An Empirical Relation between Sodium Absorption and Dust Extinction

Dust extinction and reddening are ubiquitous in astronomical observations and are often a major source of systematic uncertainty. We present here a study of the correlation between extinction in the Milky Way and the equivalent width of the NaI D absorption doublet. Our sample includes more than 100 high resolution spectra from the KECK telescopes and nearly a million low resolution spectra from the Sloan Digital Sky Survey (SDSS). We measure the correlation to unprecedented precision, constrain its shape, and derive an empirical relation between these quantities with a dispersion of order 0.15 magnitude in E(B-V). From the shape of the curve of growth we further show that a typical sight line through the Galaxy, as seen within the SDSS footprint, crosses about three dust clouds. We provide a brief guide on how to best estimate extinction to extragalactic sources such as supernovae, using the NaI D absorption feature, under a variety of circumstances.

preprint2012arXiv

HI content, metallicities and spin temperatures of damped and sub-damped Lyman alpha systems in the redshift desert (0.6 < z < 1.7)

The HI 21cm absorption optical depth and the N(HI) derived from Lya absorption can be combined to yield the spin temperature (Ts) of DLAs. Although Ts measurements exist for samples of DLAs with z <0.6 and z >1.7, the intermediate redshift regime currently contains only 2 HI 21cm detections, leading to a `redshift desert' that spans 4 Gyrs of cosmic time. To connect the low and high z regimes, we present observations of the Lya line of six 0.6<z<1.7 HI 21cm absorbers. The dataset is complemented by both VLBA observations (to derive the absorber covering factor, f), and optical echelle spectra from which metal abundances are determined. Our dataset therefore not only offers the largest statistical study of HI 21cm absorbers to date, and bridges the redshift desert, but is also the first to use a fully f-corrected dataset to look for metallicity-based trends. In agreement with trends found in Galactic sightlines, we find that the lowest N(HI) absorbers tend to be dominated by warm gas. In the DLA regime, spin temperatures show a wider range of values than Galactic data, as may be expected in a heterogenous galactic population. However, we find that low metallicity DLAs are dominated by small cold gas fractions and only absorbers with relatively high metallicities exhibit significant fractions of cold gas. Using a compilation of HI 21cm absorbers which are selected to have f-corrected spin temperatures, we confirm an anti-correlation between metallicity and Ts at 3.4 sigma significance. Finally, one of the DLAs in our sample is a newly-discovered HI 21cm absorber (at z=0.602 towards J1431+3952), which we find to have the lowest f-corrected spin temperature yet reported in the literature: Ts=90+-23 K. The observed distribution of Ts and metallicities in DLAs and the implications for understanding the characteristics of the interstellar medium in high redshift galaxies are discussed.

preprint2012arXiv

MAGICC haloes: confronting simulations with observations of the circumgalactic medium at z=0

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Lyαabsorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

preprint2012arXiv

Metallicity Evolution of Damped Lyman-alpha Systems out to z~5

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

preprint2012arXiv

Not Dead Yet: Cool Circumgalactic Gas in the Halos of Early Type 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.

preprint2012arXiv

On the redshift of the blazar PKS0447-439

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.

preprint2011arXiv

Absorption line systems in simulated galaxies fed by cold streams

Hydro cosmological simulations reveal that massive galaxies at high redshift are fed by long narrow streams of merging galaxies and a smoother component of cold gas. We post-process seven high-resolution simulated galaxies with radiative transfer to study the absorption characteristics of the gas in galaxies and streams, in comparison with the statistics of observed absorption-line systems. We find that much of the stream gas is ionized by UV radiation from background and local stellar sources, but still optically thick (N_HI > 10^17 cm^-2) so that the streams appear as Lyman-limit systems (LLSs). At z>3, the fraction of neutral gas in streams becomes non-negligible, giving rise to damped Lyman-alpha (DLA) absorbers as well. The gas in the central and incoming galaxies remains mostly neutral, responsible for DLAs. Within one (two) virial radii, the covering factor of optically thick gas is <25% (10%) for LLSs and <5% (1%) for DLAs, slowly declining with time following the universal expansion. Nevertheless, galaxies and their cold streams in the studied mass range, M_vir = 10^10 - 10^12 Msun, account for >30% of the observed absorbers in the foreground of quasars, the rest possibly arising from smaller galaxies or the intergalactic medium. The mean metallicity in the streams is ~1% solar, much lower than in the galaxies. The simulated galaxies reproduce the Lyalpha-absorption equivalent widths observed around Lyman-break galaxies, but they severely underpredict the equivalent widths in metal lines, suggesting that the latter may arise from outflows. We conclude that the observed metal-poor LLSs are likely detections of the predicted cold streams. Revised analysis of the observed LLSs kinematics and simulations with more massive outflows in conjunction with the inflows may enable a clearer distinction between the signatures of the various gas modes.

preprint2011arXiv

Detection of Pristine Gas Two Billion Years after the Big Bang

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.

preprint2011arXiv

Evidence for Cold Accretion: Primitive Gas Flowing onto a Galaxy at z~0.274

We present UV and optical observations from the Cosmic Origins Spectrograph on the Hubble Space Telescope and Keck of a z= 0.27395 Lyman limit system (LLS) seen in absorption against the QSO PG1630+377. We detect H I absorption with log N(HI)=17.06\pm0.05 as well as Mg II, C III, Si III, and O VI in this system. The column densities are readily explained if this is a multi-phase system, with the intermediate and low ions arising in a very low metallicity ([Mg/ H] =-1.71 \pm 0.06) photoionized gas. We identify via Keck spectroscopy and Large Binocular Telescope imaging a 0.3 L_* star-forming galaxy projected 37 kpc from the QSO at nearly identical redshift (z=0.27406, Δv = -26 \kms) with near solar metallicity ([O/ H]=-0.20 \pm 0.15). The presence of very low metallicity gas in the proximity of a near-solar metallicity, sub-L_* galaxy strongly suggests that the LLS probes gas infalling onto the galaxy. A search of the literature reveals that such low metallicity LLSs are not uncommon. We found that 50% (4/8) of the well-studied z < 1 LLSs have metallicities similar to the present system and show sub-L_* galaxies with rho < 100 kpc in those fields where redshifts have been surveyed. We argue that the properties of these primitive LLSs and their host galaxies are consistent with those of cold mode accretion streams seen in galaxy simulations.

preprint2011arXiv

Probing the IGM/Galaxy Connection IV: The LCO/WFCCD Galaxy Survey of 20 Fields Surrounding UV Bright Quasars

We publish the survey for galaxies in 20 fields containing ultraviolet bright quasars (with z_em 0.1 to 0.5) that can be used to study the association between galaxies and absorption systems from the low-z intergalactic medium (IGM). The survey is magnitude limited (R~19.5 mag) and highly complete out to 10' from the quasar in each field. It was designed to detect dwarf galaxies (L ~ 0.1 L*) at an impact parameter rho 1Mpc (z=0.1) from a quasar. The complete sample (all 20 fields) includes R-band photometry for 84718 sources and confirmed redshifts for 2800 sources. This includes 1198 galaxies with 0.005 < z < (z_em - 0.01) at a median redshift of 0.18, which may associated with IGM absorption lines. All of the imaging was acquired with cameras on the Swope 40" telescope and the spectra were obtained via slitmask observations using the WFCCD spectrograph on the Dupont 100" telescope at Las Campanas Observatory (LCO). This paper describes the data reduction, imaging analysis, photometry, and spectral analysis of the survey. We tabulate the principal measurements for all sources in each field and provide the spectroscopic dataset online.

preprint2011arXiv

Probing the IGM/Galaxy Connection V: On the Origin of Lya and OVI Absorption at z<0.2

We analyze the association of galaxies to Lya and OVI absorption, the most commonly detected transitions in the low-z intergalactic medium (IGM), in the fields of 14 quasars with z_em = 0.06-0.57. Confirming previous studies, we observe a high covering fraction for Lya absorption to impact parameter rho = 300kpc: 33/37 of our L>0.01L* galaxies show Lya equivalent width W_Lya>50mA. Galaxies of all luminosity L>0.01L* and spectral type are surrounded by a diffuse and ionized circumgalactic medium (CGM), whose baryonic mass is estimated at ~10^(10.5 +/- 0.3) Msun for a constant N_H. The virialized halos and extended CGM of present-day galaxies are responsible for most strong Lya absorbers (W_Lya > 300mA) but cannot reproduce the majority of observed lines in the Lya forest. We conclude that the majority of Lya absorption with W_Lya=30-300mA occurs in the cosmic web predicted by cosmological simulations and estimate a characteristic width for these filaments of ~400kpc. Regarding OVI, we observe a near unity covering fraction to rho=200kpc for L>0.1L* galaxies and to rho = 300kpc for sub-L* (0.1 L*<L<L*) galaxies. Similar to our Lya results, stronger OVI systems (W_OVI > 70mA) arise in the virialized halos of L>0.1L* galaxies. Unlike Lya, the weaker OVI systems (W_OVI~30mA) arise in the extended CGM of sub-L* galaxies. The majority of OVI gas observed in the low-z IGM is associated with a diffuse medium surrounding individual galaxies with L~0.3L*, and rarely originates in the so-called warm-hot IGM (WHIM) predicted by cosmological simulations.

preprint2011arXiv

Shining Light on Merging Galaxies I: The Ongoing Merger of a Quasar with a `Green Valley' Galaxy

Serendipitous observations of a pair z = 0.37 interacting galaxies (one hosting a quasar) show a massive gaseous bridge of material connecting the two objects. This bridge is photoionized by the quasar (QSO) revealing gas along the entire projected 38 kpc sightline connecting the two galaxies. The emission lines that result give an unprecedented opportunity to study the merger process at this redshift. We determine the kinematics, ionization parameter (log U ~ -2.5 +- 0.03), column density (N_H ~ 10^{21} cm^{-2}), metallicity ([M/H] ~ -0.20 +- 0.15), and mass (~ 10^8 Msun) of the gaseous bridge. We simultaneously constrain properties of the QSO-host (M_DM>8.8x 10^{11} Msun) and its companion galaxy (M_DM>2.1 x 10^{11} Msun; M_star ~ 2 x 10^{10} Msun; stellar burst age=300-800 Myr; SFR~6 Msun/yr; and metallicity 12+log (O/H)= 8.64 +- 0.2). The general properties of this system match the standard paradigm of a galaxy-galaxy merger caught between first and second passage while one of the galaxies hosts an active quasar. The companion galaxy lies in the so-called `green valley', with a stellar population consistent with a recent starburst triggered during the first passage of the merger and has no detectable AGN activity. In addition to providing case-studies of quasars associated with galaxy mergers, quasar/galaxy pairs with QSO-photoionized tidal bridges such as this one offer unique insights into the galaxy properties while also distinguishing an important and inadequately understood phase of galaxy evolution.

preprint2011arXiv

Simple Models of Metal-Line Absorption and Emission from Cool Gas Outflows

We analyze the absorption and emission-line profiles produced by a set of simple, cool gas wind models motivated by galactic-scale outflow observations. We implement monte carlo radiative transfer techniques that track the propagation of scattered and fluorescent photons to generate 1D spectra and 2D spectral images. We focus on the MgII 2796,28303 doublet and FeII UV1 multiplet at ~2600A, but the results are applicable to other transitions that trace outflows (e.g. NaI, Lya, SiII). By design, the resonance transitions show blue-shifted absorption but one also predicts strong resonance and fine-structure line-emission at roughly the systemic velocity. This line-emission `fills-in' the absorption reducing the equivalent width by up to 50%, shift the absorption-lin centroid by tens of km/s, and reduce the effective opacity near systemic. Analysis of cool gas outflows that ignores this line-emission may incorrectly infer that the gas is partially covered, measure asignificantly lower peak optical depth, and/or conclude that gas at systemic velocity is absent. Because the FeII lines are connected by optically-thin transitions to fine-structure levels, their profiles more closely reproduce the intrinsic opacity of the wind. Together these results naturally explain the absorption and emission-line characteristics observed for star-forming galaxies at z<1. We also study a scenario promoted to describe the outflows of z~3 Lyman break galaxies and find prfiles inconsistent with the observations due to scattered photon emission. Although line-emission complicates the analysis of absorption-line profiles, the surface brightness profiles offer a unique means of assessing the morphology and size of galactic-scale winds. Furthermore, the kinematics and line-ratios offer powerful diagnostics of outflows, motivating deep, spatially-extended spectroscopic observations.

preprint2011arXiv

Spectral Polarization of the Redshifted 21 cm Absorption Line Toward 3C 286

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

preprint2011arXiv

The COS-Halos Survey: Keck LRIS and Magellan MagE Optical Spectroscopy

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.

preprint2011arXiv

The Direct Detection of Cool, Metal-Enriched Gas Accretion onto Galaxies at z ~ 0.5

We report on the discovery of cool gas inflow toward six star-forming galaxies with redshifts z ~ 0.35 - 1. Analysis of MgII and FeII resonance-line absorption in Keck/LRIS spectroscopy of this sample reveals velocity shifts of 80 - 200 km/s and equivalent widths for inflowing gas of >~ 0.6 Ang in five of the six objects. The host galaxies exhibit a wide range in star formation rates (SFR ~ 1 - 40 M_sun/yr) and have stellar masses similar to that of the Milky Way (log M_*/M_sun ~ 9.6 - 10.5). Imaging from the Hubble Space Telescope Advanced Camera for Surveys indicates that five of the six galaxies have highly inclined (i > 55 deg), disk-like morphologies. These data represent the first unambiguous detection of inflow into isolated, star-forming galaxies in the distant universe. We suggest that the inflow is due to the infall of enriched material from dwarf satellites and/or a galactic fountain within the galaxies. Assuming that the material has been enriched to 0.1Z_sun and has a physical extent approximately equal to that of the galaxies, we infer mass inflow rates of dM_in/dt >~ 0.2 - 3 M_sun/yr for four of these systems. Finally, from comparison of these absorption lines to the profiles of MgII and FeII absorption in a larger spectroscopic sample of ~100 objects, we measure a covering fraction of cool inflow of at least 6%, but cannot rule out the presence of enriched infall onto as many as ~40 of these galaxies.

preprint2011arXiv

The End of Helium Reionization at z~2.7 Inferred from Cosmic Variance in HST/COS HeII Lyman Alpha Absorption Spectra

We report on the detection of strongly varying intergalactic HeII absorption in HST/COS spectra of two z~3 quasars. From our homogeneous analysis of the HeII absorption in these and three archival sightlines, we find a marked increase in the mean HeII effective optical depth from tau~1 at z~2.3 to tau>5 at z~3.2, but with a large scatter of 2< tau <5 at 2.7< z <3 on scales of ~10 proper Mpc. This scatter is primarily due to fluctuations in the HeII fraction and the HeII-ionizing background, rather than density variations that are probed by the co-eval HI forest. Semianalytic models of HeII absorption require a strong decrease in the HeII-ionizing background to explain the strong increase of the absorption at z>2.7, probably indicating HeII reionization was incomplete at z>2.7. Likewise, recent three-dimensional numerical simulations of HeII reionization qualitatively agree with the observed trend only if HeII reionization completes at z=2.7 or even below, as suggested by a large tau>3 in two of our five sightlines at z<2.8. By doubling the sample size at 2.7< z <3, our newly discovered HeII sightlines for the first time probe the diversity of the second epoch of reionization when helium became fully ionized.

preprint2011arXiv

The Hidden Mass and Large Spatial Extent of a Poststarburst Galaxy Outflow

Outflowing winds of multiphase plasma have been proposed to regulate the buildup of galaxies, but key aspects of these outflows have not been probed with observations. Using ultraviolet absorption spectroscopy, we show that "warm-hot" plasma at 10^{5.5} K contains 10-150 times more mass than the cold gas in a poststarburst galaxy wind. This wind extends to distances >68 kiloparsecs, and at least some portion of it will escape. Moreover, the kinematical correlation of the cold and warm-hot phases indicates that the warm-hot plasma is related to the interaction of the cold matter with a hotter (unseen) phase at >>10^{6} K. Such multiphase winds can remove substantial masses and alter the evolution of poststarburst galaxies.

preprint2011arXiv

The Large, Oxygen-Rich Halos of Star-Forming Galaxies Are A Major Reservoir of Galactic Metals

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.

preprint2011arXiv

The Last Eight-Billion Years of Intergalactic SiIV Evolution

We identified 24 SiIV absorption systems with z <~ 1 from a blind survey of 49 low-redshift quasars with archival Hubble Space Telescope ultraviolet spectra. We relied solely on the characteristic wavelength separation of the doublet to automatically detect candidates. After visual inspection, we defined a sample of 20 definite (group G = 1) and 4 "highly-likely" (G = 2) doublets with rest equivalent widths W_r for both lines detected at > 3 sigma. The absorber line density of the G = 1 doublets was dN_SiIV/dX = 1.4+0.4/-0.3 for log N(Si+3) > 12.9. The best-fit power law to the G = 1 frequency distribution of column densities f(N(Si+3)) had normalization k = (1.2+0.5/-0.4) x 10^-14 cm2 and slope alpha = -1.6+0.3/-0.3. Using the power-law model of f(N(Si+3)), we measured the Si+3 mass density relative to the critical density: Omega(Si+3) = (3.7+2.8/-1.7) x 10^-8 for 13 < log N(Si+3) < 15. From Monte Carlo sampling of the distributions, we estimated our value to be a factor of 4.8+3.0/-1.9 higher than the 2 < z < 4.5 <Omega(Si+3)>. From a simple linear fit to Omega(Si+3) over the age of the Universe, we estimated a slow and steady increase from z = 5.5 --> 0 with dOmega/dt_age = (0.61+/-0.23) x 10^-8 Gyr^-1. We compared our ionic ratios N(Si+3)/N(C+3) to a 2 < z < 4.5 sample and concluded, from survival analysis, that the two populations are similar, with median <N(Si+3)/N(C+3)> = 0.16.

preprint2011arXiv

What Drives the Expansion of Giant HII Regions?: A Study of Stellar Feedback in 30 Doradus

Observations show that star formation is an inefficient and slow process. This result can be attributed to the injection of energy and momentum by stars that prevents free-fall collapse of molecular clouds. The mechanism of this stellar feedback is debated theoretically: possible sources of pressure include the classical warm HII gas, the hot gas generated by shock-heating from stellar winds and supernovae, direct radiation of stars, and the dust-processed radiation field trapped inside the HII shell. In this paper, we measure observationally the pressures associated with each component listed above across the giant HII region 30 Doradus in the Large Magellanic Cloud. We exploit high-resolution, multi-wavelengh images (radio, infrared, optical, and X-ray) to map these pressures as a function of position. We find that radiation pressure dominates within 75 pc of the central star cluster, R136, while the HII gas pressure dominates at larger radii. By contrast, the dust-processed radiation pressure and hot gas pressure are generally weak and not dynamically important, although the hot gas pressure may have played a more significant role at early times. Based on the low X-ray gas pressures, we demonstrate that the hot gas is only partially confined and must be leaking out the HII shell. Additionally, we consider the implications of a dominant radiation pressure on the early dynamics of 30 Doradus.

preprint2010arXiv

Directly imaging damped Lyman-alpha galaxies at z>2. I: Methodology and First Results

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.

preprint2010arXiv

GALEX far-UV color selection of UV-bright high-redshift quasars

We study the small population of z>2.7 quasars detected by GALEX, whose far-UV emission is not extinguished by intervening HI Lyman limit systems. These quasars are of particular importance to detect intergalactic HeII absorption along their sightlines. We correlate verified z>2.7 quasars to the GALEX GR4 source catalog, yielding 304 S/N>3 sources. However, ~50% of these are only detected in the GALEX NUV band, signaling the truncation of the FUV flux by low-redshift Lyman limit systems. We exploit the GALEX UV color to cull the most promising targets for follow-up studies, with blue (red) GALEX colors indicating transparent (opaque) sightlines. Monte Carlo simulations indicate a HeII detection rate of ~60% for quasars with FUV-NUV<1 at z<3.5, a ~50% increase over GALEX searches that do not include color information. We regard 52 quasars detected at S/N>3 to be most promising for HST follow-up, with an additional 114 quasars if we consider S/N>2 detections in the FUV. SDSS provides just half of the NUV-bright quasars that should have been detected by SDSS & GALEX. We revise the SDSS quasar selection function, finding that SDSS systematically misses quasars with blue u-g<2 colors at 3<z<3.5 due to overlap with the stellar locus in color space. Our color-dependent SDSS selection function naturally explains the inhomogeneous u-g color distribution of SDSS quasars with redshift and the color difference between color-selected and radio-selected SDSS quasars. Moreover, it yields excellent agreement between the observed and the predicted number of UV-bright SDSS quasars. We confirm our previous claims that SDSS preferentially selects 3<z<3.5 quasars with intervening HI Lyman limit systems. Our results imply that broadband optical color surveys for 3<z<3.5 quasars have likely underestimated their space density by selecting IGM sightlines with an excess of strong HI absorbers.

preprint2010arXiv

H I Column Densities, Metallicities, and Dust Extinction of Metal-Strong Damped Lya Systems

With the Blue Channel Spectrograph (BCS) on the MMT telescope, we have obtained spectra to the atmospheric cutoff of quasars previously known to show at least one absorption system at z>1.6 with very strong metal lines (candidate metal-strong damped Lya systems; cMSDLAs). The BCS/MMT spectra yield precise estimates of the HI column densities (NHI) of the systems through Voigt profile analysis of their Lya transitions. Nearly all of the cMSDLAs (41/43) satisfy the NHI criterion of DLAs, 10^20.3. As a population, these systems have systematically higher NHI values than DLAs chosen randomly from quasar sightlines. Combining our NHI measurements with previously measured metal column densities, we estimate metallicities for the MSDLAs. These systems have significantly higher values than randomly selected DLAs; at z~2, the MSDLAs show a median metallicity [M/H] ~ -0.67 that is 0.6dex higher than a corresponding control sample. This establishes MSDLAs as having amongst the most metal-rich gas in the high z universe. Our measurements extend the observed correlation between SiII 1526 equivalent width and the gas metallicity to higher values. If interpreted as a mass-metallicity relation, this implies the MSDLAs are the high mass subset of the DLA population. We demonstrate that dust in the MSDLAs reddens their background quasars, with a median shift in the spectral slope of Da = 0.29. Assuming an SMC extinction law, this implies a median reddening E(B-V)=0.025mag and visual extinction A_V=0.076mag. Future studies of MSDLAs offer the opportunity to study the extinction, nucleosynthesis, and kinematics of the most chemically evolved, gas-rich galaxies at high z. [abridged]

preprint2010arXiv

High Dust Depletion in Two Intervening Quasar Absorption Line Systems with the 2175-Å Extinction Bump at $z\sim1.4$

We present the column densities of heavy-elements and dust depletion studies in twostrong Mg~II absorption systems at $z\sim1.4$ displaying the 2175-Å dust extinction feature. Column densities are measured from low-ionization absorption lines using Apparent Optical Depth Method on the Keck/ESI spectra. We find the dust depletion patterns resemble to that of cold diffuse clouds in the Milky Way (MW). The values, [Fe/Zn]$\approx -1.5$ and [Si/Zn]$<-0.67$, are among the highest dust depletion measured for quasar absorption line systems. In another 2175-Å absorber at $z$=1.64 toward the quasar SDSS J160457.50+220300.5, Noterdaeme et al. (2009) reported a similar dust depletion measurement ([Fe/Zn]=$-1.47$ and [Si/Zn]=$-1.07$) and detected C~I and CO absorption lineson its VLT/UVES spectrum. We conclude that heavy dust depletion (i.e. a characteristic of cold dense clouds in MW) is required to produce a pronounced 2175-Å extinction bump.

preprint2010arXiv

Ionization corrections in a multi-phase interstellar medium: Lessons from a z~2 sub-DLA

We present a high resolution (FWHM=2.7 km/s), high S/N echelle spectrum for the z = 2.26 QSO J2123-0050 and determine elemental abundances for the z = 2.06 sub-DLA in its line of sight. This high redshift sub-DLA has a complex kinematic structure and harbours detections of neutral (SI, CI), singly (e.g. CII, SII) and multiply ionized (e.g. CIV, SiIV) species as well as molecular H and HD. The plethora of detected transitions in various ionization stages is indicative of a complex multi-phase structure present in this high redshift galaxy. We demonstrate that the ionization corrections in this sub-DLA are significant (up to ~0.7 dex). For example, if no ionization correction is applied, a super-solar metallicity is derived ([S/H] = +0.36), whereas a single phase ionization correction reduces this to [S/H] = -0.19. The theoretical impact of a multi-phase medium is investigated through Cloudy modelling and it is found that the abundances of Si, S and Fe are always over-estimated (by up to 0.15 dex in our experiments) if a single-phase is assumed. Therefore, although Cloudy models improve estimates of metal column densities, the simplification of a single phase medium leaves a systematic error in the result, so that even ionization-corrected abundances may still be too high. Without ionization corrections the properties of this sub-DLA appear to require extreme scenarios of nucleosynthetic origins. After ionization corrections are applied the ISM of this galaxy appears to be similar to some of the sightlines through the Milky Way.

preprint2010arXiv

Keck ESI Observations of Metal-Poor Damped Lyman-α Systems

We present the first results from a survey of SDSS quasars selected for strong H I damped Lyman-α (DLA) absorption with corresponding low equivalent width absorption from strong low-ion transitions (e.g. C II λ1334 and Si II λ1260). These metal-poor DLA candidates were selected from the SDSS DR5 quasar spectroscopic database, and comprise a large new sample for probing low metallicity galaxies. Medium-resolution echellette spectra from the Keck ESI spectrograph for an initial sample of 35 systems were obtained to explore the metal-poor tail of the DLA distribution and to investigate the nucleosynthetic patterns at these metallicities. We have estimated saturation corrections for the moderately under-resolved spectra, and systems with very narrow Doppler parameter (b \le 5 km s-1) will likely have underestimated abundances. For those systems with Doppler parameters b > 5 km s-1, we have measured low metallicity DLA gas with [X/H] < -2.4 for at least one of C, O, Si, or Fe. Assuming non-saturated components, we estimate that several DLA systems have [X/H] < -2.8, including five DLA systems with both low equivalent widths and low metallicity in transitions of both C II and O I. All of the measured DLA metallicities, however, exceed or are consistent with a metallicity of at least 1/1000 of solar, regardless of the effects of saturation in our spectra. Our results indicate that the metal-poor tail of galaxies at z \sim 3 drops exponentially at [X/H]< -3. The observed ratio of [C/O] for values of [O/H] < -2.5 exceeds values seen in moderate metallicity DLA systems, and also exceeds theoretical nucleosynthesis predictions for higher mass Population III stars.

preprint2010arXiv

Keck telescope constraint on cosmological variation of the proton-to-electron mass ratio

Molecular transitions recently discovered at redshift z_abs=2.059 toward the bright background quasar J2123-0050 are analysed to limit cosmological variation in the proton-to-electron mass ratio, mu=m_p/m_e. Observed with the Keck telescope, the optical echelle spectrum has the highest resolving power and largest number (86) of H_2 transitions in such analyses so far. Also, (seven) HD transitions are used for the first time to constrain mu-variation. These factors, and an analysis employing the fewest possible free parameters, strongly constrain mu's relative deviation from the current laboratory value: dmu/mu =(+5.6+/-5.5_stat+/-2.9_sys)x10^{-6}, indicating an insignificantly larger mu in the absorber. This is the first Keck result to complement recent null constraints from three systems at z_abs>2.5 observed with the Very Large Telescope. The main possible systematic errors stem from wavelength calibration uncertainties. In particular, distortions in the wavelength solution on echelle order scales are estimated to contribute approximately half the total systematic error component, but our estimate is model dependent and may therefore under or overestimate the real effect, if present. To assist future mu-variation analyses of this kind, and other astrophysical studies of H_2 in general, we provide a compilation of the most precise laboratory wavelengths and calculated parameters important for absorption-line work with H_2 transitions redwards of the hydrogen Lyman limit.

preprint2010arXiv

Low-ionization Line Emission from Starburst Galaxies: A New Probe of Galactic-Scale Outflows

We study the kinematically narrow, low-ionization line emission from a bright, starburst galaxy at z = 0.69 using slit spectroscopy obtained with Keck/LRIS. The spectrum reveals strong absorption in MgII and FeII resonance transitions with Doppler shifts of -200 to -300 km/s, indicating a cool gas outflow. Emission in MgII near and redward of systemic velocity, in concert with the observed absorption, yields a P Cygni-like line profile similar to those observed in the Ly alpha transition in Lyman Break Galaxies. Further, the MgII emission is spatially resolved, and extends significantly beyond the emission from stars and HII regions within the galaxy. Assuming the emission has a simple, symmetric surface brightness profile, we find that the gas extends to distances > ~7 kpc. We also detect several narrow FeII* fine-structure lines in emission near the systemic velocity, arising from energy levels which are radiatively excited directly from the ground state. We suggest that the MgII and FeII* emission is generated by photon scattering in the observed outflow, and emphasize that this emission is a generic prediction of outflows. These observations provide the first direct constraints on the minimum spatial extent and morphology of the wind from a distant galaxy. Estimates of these parameters are crucial for understanding the impact of outflows in driving galaxy evolution.

preprint2010arXiv

Metallicities and dust content of proximate damped Lyman alpha systems in the Sloan Digital Sky Survey

Composite spectra of 85 proximate absorbers (log N(HI)>20 and velocity difference between the absorption and emission redshift, dv<10,000 km/s) in the SDSS are used to investigate the trends of metal line strengths with velocity separation from the QSO. We construct composites in 3 velocity bins: dv<3000 km/s, 3000<dv<6000 km/s and dv>6000 km/s, with further sub-samples to investigate the metal line dependence on N(HI) and QSO luminosity. Low (e.g. SiII and FeII) and high ionization (e.g. SiIV and CIV) species alike have equivalent widths (EWs) that are larger by factors of 1.5 -- 3 in the dv<3000 km/s composite, compared to the dv>6000 km/s spectrum. The EWs show an even stronger dependence on dv if only the highest neutral hydrogen column density (log N(HI)>20.7) absorbers are considered. We conclude that PDLAs generally have higher metallicities than intervening absorbers, with the enhancement being a function of both dv and N(HI). It is also found that absorbers near QSOs with lower rest-frame UV luminosities have significantly stronger metal lines. We speculate that absorbers near to high luminosity QSOs may have had their star formation prematurely quenched. Finally, we search for the signature of dust reddening by the PDLAs, based on an analysis of the QSO continuum slopes relative to a control sample and determine a limit of E(B-V)<0.014 for an SMC extinction curve. This work provides an empirical motivation for distinguishing between proximate and intervening DLAs, and establishes a connection between the QSO environment and galaxy properties at high redshifts.

preprint2010arXiv

New limit on a varying proton-to-electron mass ratio from high-resolution optical quasar spectra

Molecular transitions recently discovered at redshift z_abs=2.059 toward the bright background quasar J2123-0050 are analysed to limit cosmological variation in the proton-to-electron mass ratio, mu=m_p/m_e. Observed with the Keck telescope, the optical spectrum has the highest resolving power and largest number (86) of H_2 transitions in such analyses so far. Also, (7) HD transitions are used for the first time to constrain mu-variation. These factors, and an analysis employing the fewest possible free parameters, strongly constrain mu's relative deviation from the current laboratory value: dmu/mu =(+5.6+/-5.5_stat+/-2.7_sys)x10^{-6}. This is the first Keck result to complement recent constraints from three systems at z_abs>2.5 observed with the Very Large Telescope.

preprint2010arXiv

On the Analysis of DLA Kinematics

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

preprint2010arXiv

The nature of proximate damped Lyman alpha systems

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

preprint2010arXiv

Understanding Physical Conditions in High Redshift Galaxies through C I Fine Structure Lines: Data and Methodology

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

preprint2009arXiv

A Definitive Survey for Lyman Limit Systems at z~3.5 with the Sloan Digital Sky Survey

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.

preprint2009arXiv

A Direct Measurement of the IGM Opacity to HI Ionizing Photons

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.

preprint2009arXiv

A Survey of UV Bright Sources Behind the Halo of M31

We have performed a wide-area ultraviolet (UV) imaging survey using the GALaxy Evolution eXplorer (GALEX) to search for bright, point-like UV sources behind M31's extended halo. Our survey consisted of 46 pointings covering an effective area of ~50 deg^2, in both the far-UV and near-UV channels. We combined these data with optical R-band observations acquired with the WIYN Mosaic-1 imager on the Kitt Peak National Observatory 0.9m WIYN telescope. An analysis of the brightness and colors of sources matched between our photometric catalogs yielded ~100 UV-bright quasar candidates. We have obtained discovery spectra for 76 of these targets with the Kast spectrometer on the Lick 3m telescope and confirm 30 active galactic nuclei and quasars, 29 galaxies at z > 0.02 including several early-type systems, 16 Galactic stars (hot main-sequence stars), and one featureless source previously identified as a BL Lac object. Future UV spectroscopy of the brightest targets with the Cosmic Origins Spectrograph on the Hubble Space Telescope will enable a systematic search for diffuse gas in the extended halo of M31.

preprint2009arXiv

An Observational Determination of the Proton to Electron Mass Ratio in the Early Universe

In an effort to resolve the discrepancy between two measurements of the fundamental constant mu, the proton to electron mass ratio, at early times in the universe we reanalyze the same data used in the earlier studies. Our analysis of the molecular hydrogen absorption lines in archival VLT/UVES spectra of the damped Lyman alpha systems in the QSOs Q0347-383 and Q0405-443 yields a combined measurement of a (Delta mu)/mu value of (-7 +/- 8) x 10^{-6}, consistent with no change in the value of mu over a time span of 11.5 gigayears. Here we define (Delta mu) as (mu_z - mu_0) where mu_z is the value of mu at a redshift of z and mu_0 is the present day value. Our null result is consistent with the recent measurements of King et al. 2009, (Delta mu)/u = (2.6 +/- 3.0) x 10^{-6}, and inconsistent with the positive detection of a change in mu by Reinhold et al. 2006. Both of the previous studies and this study are based on the same data but with differing analysis methods. Improvements in the wavelength calibration over the UVES pipeline calibration is a key element in both of the null results. This leads to the conclusion that the fundamental constant mu is unchanged to an accuracy of 10^{-5} over the last 80% of the age of the universe, well into the matter dominated epoch. This limit provides constraints on models of dark energy that invoke rolling scalar fields and also limits the parameter space of Super Symmetric or string theory models of physics. New instruments, both planned and under construction, will provide opportunities to greatly improve the accuracy of these measurements.

preprint2009arXiv

Direct Evidence of Cold Gas in DLA 0812+32B

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

preprint2009arXiv

Galaxies probing galaxies: cool halo gas from a z = 0.47 post-starburst galaxy

We study the cool gas around a galaxy at z = 0.4729 using Keck/LRIS spectroscopy of a bright (B = 21.7) background galaxy at z = 0.6942 at a transverse distance of 16.5/h_70 kpc. The background galaxy spectrum reveals strong FeII, MgII, MgI, and CaII absorption at the redshift of the foreground galaxy, with a MgII 2796 rest equivalent width of 3.93 +/- 0.08 Angstroms, indicative of a velocity width exceeding 400 km/s. Because the background galaxy is large (> 4/h_70 kpc), the high covering fraction of the absorbing gas suggests that it arises in a spatially extended complex of cool clouds with large velocity dispersion. Spectroscopy of the massive (log M_*/M_sun = 11.15 +/- 0.08) host galaxy reveals that it experienced a burst of star formation about 1 Gyr ago and that it harbors a weak AGN. We discuss the possible origins of the cool gas in its halo, including multiphase cooling of hot halo gas, cold inflow, tidal interactions, and galactic winds. We conclude the absorbing gas was most likely ejected or tidally stripped from the interstellar medium of the host galaxy or its progenitors during the past starburst event. Adopting the latter interpretation, these results place one of only a few constraints on the radial extent of cool gas driven or stripped from a galaxy in the distant Universe. Future studies with integral field unit spectroscopy of spatially extended background galaxies will provide multiple sightlines through foreground absorbers and permit analysis of the morphology and kinematics of the gas surrounding galaxies with a diverse set of properties and environments.

preprint2009arXiv

GRB 090426: The Environment of a Rest-Frame 0.35-second Gamma-Ray Burst at Redshift z=2.609

We present the discovery of an absorption-line redshift of z = 2.609 for GRB 090426, establishing the first firm lower limit to a redshift for a gamma-ray burst with an observed duration of <2 s. With a rest-frame burst duration of T_90z = 0.35 s and a detailed examination of the peak energy of the event, we suggest that this is likely (at >90% confidence) a member of the short/hard phenomenological class of GRBs. From analysis of the optical-afterglow spectrum we find that the burst originated along a very low HI column density sightline, with N_HI < 3.2 x 10^19 cm^-2. Our GRB 090426 afterglow spectrum also appears to have weaker low-ionisation absorption (Si II, C II) than ~95% of previous afterglow spectra. Finally, we also report the discovery of a blue, very luminous, star-forming putative host galaxy (~2 L*) at a small angular offset from the location of the optical afterglow. We consider the implications of this unique GRB in the context of burst duration classification and our understanding of GRB progenitor scenarios.

preprint2009arXiv

MASE: A New Data--Reduction Pipeline for the Magellan Echellette Spectrograph

We introduce a data reduction package written in Interactive Data Language (IDL) for the Magellan Echellete Spectrograph (MAGE). MAGE is a medium-resolution (R ~4100), cross-dispersed, optical spectrograph, with coverage from ~3000-10000 Angstroms. The MAGE Spectral Extractor (MASE) incorporates the entire image reduction and calibration process, including bias subtraction, flat fielding, wavelength calibration, sky subtraction, object extraction and flux calibration of point sources. We include examples of the user interface and reduced spectra. We show that the wavelength calibration is sufficient to achieve ~5 km/s RMS accuracy and relative flux calibrations better than 10%. A light-weight version of the full reduction pipeline has been included for real-time source extraction and signal-to-noise estimation at the telescope.

preprint2009arXiv

Metal-Enriched Plasma in Protogalactic Halos: A Survey of N V Absorption in High-z Damped & Sub-Damped Lyman-alpha Systems

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

preprint2009arXiv

Molecular hydrogen deficiency in HI-poor galaxies and its implications for star formation

We use a sample of 47 homogeneous and high sensitivity CO images taken from the Nobeyama and BIMA surveys to demonstrate that, contrary to common belief, a significant number (~40%) of HI-deficient nearby spiral galaxies are also depleted in molecular hydrogen. While HI-deficiency by itself is not a sufficient condition for molecular gas depletion, we find that H2 reduction is associated with the removal of HI inside the galaxy optical disk. Those HI-deficient galaxies with normal H2 content have lost HI mainly from outside their optical disks, where the H2 content is low in all galaxies. This finding is consistent with theoretical models in which the molecular fraction in a galaxy is determined primarily by its gas column density. Our result is supported by indirect evidence that molecular deficient galaxies form stars at a lower rate or have dimmer far infrared fluxes than gas rich galaxies, as expected if the star formation rate is determined by the molecular hydrogen content. Our result is consistent with a scenario in which, when the atomic gas column density is lowered inside the optical disk below the critical value required to form molecular hydrogen and stars, spirals become quiescent and passive evolving systems. We speculate that this process would act on the time-scale set by the gas depletion rate and might be a first step for the transition between the blue and red sequence observed in the color-magnitude diagram.

preprint2009arXiv

The Keck+Magellan Survey for Lyman Limit Absorption II: A Case Study on Metallicity Variations

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.

preprint2009arXiv

The Last Eight-Billion Years of Intergalactic CIV Evolution

We surveyed the HST UV spectra of 49 low-redshift quasars for z < 1 CIV candidates, relying solely on the characteristic wavelength separation of the doublet. After consideration of the defining traits of CIV doublets (e.g., consistent line profiles, other associated transitions, etc.), we defined a sample of 38 definite (group G = 1) and five likely (G = 2) doublets with rest equivalent widths W_r for both lines detected at >= 3 sigma. We conducted Monte-Carlo completeness tests to measure the unblocked redshift (dz) and co-moving pathlength (dX) over which we were sensitive to CIV doublets of a range of equivalent widths and column densities. The absorber line density of (G = 1+2) doublets is dN/dX = 4.1+0.7/-0.6 for log N(C^+3) >= 13.2, and dN/dX has not evolved significantly since z = 5. The best-fit power-law to the G = 1 frequency distribution of column densities f(N(C^+3)) = k(N(C^+3)/N_0)^alpha_N has coefficient k = (0.67+0.18/-0.16) x 10^-14 cm^2 and exponent alpha_N = -1.50+0.17/-0.19, where N_0 = 10^14 cm^-2. Using the power-law model of f(N(C^+3)), we measured the C^+3 mass density relative to the critical density: Omega(C^+3) = (6.20+1.82/-1.52) x 10^-8 for 13 <= log N(C^+3) <= 15. This value is a 2.8+/-0.7 increase in Omega(C^+3) compared to the error-weighted mean from several 1 < z < 5 surveys for CIV absorbers. A simple linear regression to Omega(C^+3) over the age of the Universe indicates that Omega(C^+3) has slowly but steadily increased from z = 5 --> 0, with dOmega(C^+3)/dt_age = (0.42+/-0.2) x 10^-8 Gyr^-1.

preprint2009arXiv

The Persistence of Cool Galactic Winds in High Stellar Mass Galaxies Between z~1.4 and ~1

We present an analysis of the MgII 2796, 2803 and FeII 2586, 2600 absorption line profiles in coadded spectra of 468 galaxies at 0.7 < z < 1.5. The galaxy sample, drawn from the Team Keck Treasury Redshift Survey of the GOODS-N field, has a range in stellar mass (M_*) comparable to that of the sample at z~1.4 analyzed in a similar manner by Weiner et al. (2009; W09), but extends to lower redshifts and has specific star formation rates which are lower by ~0.6 dex. We identify outflows of cool gas from the Doppler shift of the MgII absorption lines and find that the equivalent width (EW) of absorption due to outflowing gas increases on average with M_* and star formation rate (SFR). We attribute the large EWs measured in spectra of the more massive, higher-SFR galaxies to optically thick absorbing clouds having large velocity widths. The outflows have hydrogen column densities N(H) > 10^19.3 cm^-2, and extend to velocities of ~500 km/s. While galaxies with SFR > 10 Msun/yr host strong outflows in both this and the W09 sample, we do not detect outflows in lower-SFR (i.e., log M_*/Msun < 10.5) galaxies at lower redshifts. Using a simple galaxy evolution model which assumes exponentially declining SFRs, we infer that strong outflows persist in galaxies with log M_*/Msun > 10.5 as they age between z=1.4 and z~1, presumably because of their high absolute SFRs. Finally, using high resolution HST/ACS imaging in tandem with our spectral analysis, we find evidence for a weak trend (at 1 sigma significance) of increasing outflow absorption strength with increasing galaxy SFR surface density.

preprint2008arXiv

Quasars Probing Quasars III: New Clues to Feedback, Quenching, and the Physics of Massive Galaxy Formation

Galaxies hosting z~2 quasars are the high-$z$ progenitors of today's massive `red-and-dead' galaxies. With close pairs of quasars at different redshifts, a background quasar can be used to study a foreground quasar's halo gas in absorption, providing a wealth of information about feedback, quenching, and the physics of massive galaxy formation. We present a Keck/HIRES spectrum of the bright background quasar in a projected pair with angular separation 13.3" corresponding to 108kpc at the redshift of the foreground quasar z_fg=2.4360 +/- 0.0005, precisely determined from Gemini/GNIRS near-IR spectroscopy. Our echelle spectrum reveals optically thick gas (NHI~10^19.7), coincident with the foreground quasar redshift. The ionic transitions of associated metal-lines reveal the following properties of the foreground quasar's halo: (1) the kinematics are extreme with absorption extending to +780km/s relative to z_fg; (2) the metallicity is nearly solar; (3) the temperature of the predominantly ionized gas is T<~20,000K; (4) the electron density is n_e~1 cm^-3 indicating a characteristic size ~10 - 100pc for the absorbing `clouds'; (7) there is a negligible amount of warm gas 10^5K < T < 10^6K; (8) the gas is unlikely illuminated by the foreground quasar, implying anisotropic or intermittent emission. The mass of cold T~10^4K gas implied by our observations is significant, amounting to a few percent of the total expected baryonic mass density of the foreground quasar's dark halo at r~100kpc. The origin of this material is still unclear, and we discuss several possibilities in the context of current models of feedback and massive galaxy formation.