Source author record

Joseph F. Hennawi

Joseph F. Hennawi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.IM astro-ph.HE physics.data-an

Catalog footprint

What is connected

79works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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

(Nearly) Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at $z\sim 5-7$ Using Dark Pixel Fractions in Ly$α$ and Ly$β$ Forests

Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.73<z<7.09$. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of "dark" pixels identified in the Lyman $α$ (Ly$α$) and Lyman $β$ (Ly$β$) forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at $z\sim6.2-6.8$ using quasar absorption measurements. Our results give upper limits of $\bar{x}_{\rm HI}(z=6.3) < 0.79\pm0.04$ (1$σ$), $\bar{x}_{\rm HI} (z=6.5) < 0.87\pm0.03$ (1$σ$), and $\bar{x}_{\rm HI} (z=6.7) < 0.94^{+0.06}_{-0.09}$ (1$σ$). The dark pixel fractions at $z>6.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018.

preprint2022arXiv

Constraining IGM enrichment and metallicity with the C IV forest correlation function

The production and distribution of metals in the diffuse intergalactic medium (IGM) have implications for galaxy formation models and the baryon (re)cycling process. The relative abundance of metals in high versus low-ionization states has also been argued to be sensitive to the Universe's reionization history. However, measurements of the background metallicity of the IGM at z~4 are sparse and in poor agreement with one another, and reduced sensitivity in the near-IR renders detecting individual metal absorbers nearly impossible. We present a new clustering-based technique that enables the detection of these weak IGM absorbers by statistically averaging over all spectral pixels, here applied to the C IV forest. We simulate the z=4.5 IGM with different models of inhomogeneous metal distributions and investigate its two-point correlation function (2PCF) using mock skewers of the C IV forest. The 2PCF demonstrates a clear peak at the doublet separation of the C IV line. The peak amplitude scales quadratically with metallicity, while enrichment morphology affects both the shape and amplitude of the 2PCF. The effect of enrichment topology can also be framed in terms of the metal mass- and volume-filling factors, and we show their trends as a function of the enrichment topology. For models consistent with the distribution of metals at z~3, we find that we can constrain [C/H] to within 0.2 dex, log$\,M_{\rm{min}}$ to within 0.4 dex, and $R$ to within 15%. We show that strong absorbers arising from the circumgalactic medium of galaxies can be easily identified and masked, allowing one to recover the underlying IGM signal. The auto-correlation of the metal-line forest presents a new and compelling avenue to simultaneously constrain IGM metallicity and enrichment topology with high precision at z>4, thereby pushing such measurements into the Epoch of Reionization.

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

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$α$ systems. We present improved measurements of the mean Ly-$α$ transmission over $4.9<z<6.1$. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-$α$ transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to $z\leq5.2$ ($<1 σ$), and mild tension ($\sim 2.5 σ$) at $z=5.3$. Homogeneous UVB models are ruled out by excess Ly-$α$ transmission scatter at $z\geq5.4$ with high confidence ($>3.5 σ$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation.

preprint2022arXiv

Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference

We present a new approach to measure the power-law temperature density relationship $T=T_0 (ρ/ \barρ)^{γ-1}$ and the UV background photoionization rate $Γ_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$α$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape of the $b$-$N_{\rm HI}$ distribution is sensitive to the IGM thermal parameters $T_0$ and $γ$, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density d$N$/d$z$, and we demonstrate that precise constraints can also be obtained on $Γ_{\rm HI}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b$-$N_{\rm HI}$ distribution on IGM parameters trained on an ensemble of 624 Nyx hydrodynamical simulations at $z = 0.1$, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS datasets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly$α$ forest spectra at $z\simeq 0.1$, one can measure the thermal and ionization state of the IGM with very high precision ($σ_{\log T_0} \sim 0.08$ dex, $σ_γ\sim 0.06$, and $σ_{\log Γ_{\rm HI}} \sim 0.07$ dex).

preprint2021arXiv

A Luminous Quasar at Redshift 7.642

Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at $z\ge7.5$, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at $z=7.642$, J0313$-$1806, the most distant quasar yet known. This quasar has a bolometric luminosity of $3.6\times10^{13} L_\odot$. Deep spectroscopic observations reveal a SMBH with a mass of $(1.6\pm0.4) \times10^9M_\odot$ in this quasar. The existence of such a massive SMBH just $\sim$670 million years after the Big Bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in CIV and SiIV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted CIV emission line, indicate that there is a strong active galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum and [CII] emission from the quasar host galaxy, yielding an accurate redshift of $7.6423 \pm 0.0013$ and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of $\rm\sim 200 ~M_\odot ~yr^{-1}$ and a dust mass of $\sim7\times10^7~M_\odot$. Followup observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive galaxies.

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

ALMA reveals large overdensity and strong clustering of galaxies in quasar environments at $z\sim4$

We present an Atacama Large Millimeter/submillimeter Array (ALMA) survey of CO(4--3) line emitting galaxies in 17 quasar fields at $z\sim4$ aimed to perform the first systematic search of dusty galaxies in high$-z$ quasar environments. Our blind search of galaxies around the quasars results in 5 CO emitters with S/N$\geq5.6$ within a projected radius of $R\lesssim1.5\,h^{-1}\,$cMpc and a velocity range of $\rmΔ\,v=\pm1000\,$km\,s$^{-1}$ around the quasar. In blank fields, we expect to detect only 0.28 CO emitters within the same volume, implying a total overdensity of $17.6^{+11.9}_{-7.6}$ in our fields, and indicating that quasars trace massive structures in the early universe. We quantify this overdensity by measuring the small-scale clustering of CO emitters around quasars, resulting in a cross-correlation length of $r_{\rm 0,QG}=8.37^{+2.42}_{-2.04}\,h^{-1}\,$cMpc, assuming a fixed slope $γ=1.8$. This contradicts the reported mild overdensities (x1.4) of Lyman alpha emitters (LAEs) in the same fields at scales of $R\lesssim7\,h^{-1}\,$cMpc which is well described by a cross-correlation length 3 times lower than that measured for CO emitters. We discuss some possibilities to explain this discrepancy, including low star formation efficiency, and excess of dust in galaxies around quasars. Finally, we constrain, for the first time, the clustering of CO emitters at $z\sim4$, finding an auto-correlation length of $r_{\rm 0,CO}=3.14\pm1.71\,h^{-1}\,$cMpc (with $γ=1.8$). Our work together with the previous study of LAEs around quasars traces simultaneously the clustering properties of both optical and dusty galaxy populations in quasars fields, stressing the importance of multi-wavelength studies, and highlighting important questions about galaxy properties in high$-z$ dense environments.

preprint2021arXiv

IGM damping wing constraints on reionisation from covariance reconstruction of two $z\gtrsim7$ QSOs

Bright, high redshift ($z>6$) QSOs are powerful probes of the ionisation state of the intervening intergalactic medium (IGM). The detection of Ly$α$ damping wing absorption imprinted in the spectrum of high-z QSOs can provide strong constraints on the epoch of reionisation (EoR). In this work, we perform an independent Ly$α$ damping wing analysis of two known $z>7$ QSOs; DESJ0252-0503 at $z=7.00$ (Wang et al.) and J1007+2115 at $z=7.51$ (Yang et al.). For this, we utilise our existing Bayesian framework which simultaneously accounts for uncertainties in: (i) the intrinsic Ly$α$ emission profile (reconstructed from a covariance matrix of measured emission lines; extended in this work to include NV) and (ii) the distribution of ionised (H\,{\scriptsize II}) regions within the IGM using a $1.6^3$ Gpc$^3$ reionisation simulation. This approach is complementary to that used in the aforementioned works as it focuses solely redward of Ly$α$ ($1218 < λ< 1230$Å) making it more robust to modelling uncertainties while also using a different methodology for (i) and (ii). We find, for a fiducial EoR morphology, $\bar{x}_{\rm HI} = 0.64\substack{+0.19 \\ -0.23}$ (68 per cent) at $z=7$ and $\bar{x}_{\rm HI} = 0.27\substack{+0.21 \\ -0.17}$ at $z=7.51$ consistent within $1σ$ to the previous works above, though both are slightly lower in amplitude. Following the inclusion of NV into our reconstruction pipeline, we perform a reanalysis of ULASJ1120+0641 at $z=7.09$ (Mortlock et al.) and ULASJ1342+0928 at $z=7.54$ (Bañados et al.) finding $\bar{x}_{\rm HI} = 0.44\substack{+0.23 \\ -0.24}$ at $z=7.09$ and $\bar{x}_{\rm HI} = 0.31\substack{+0.18 \\ -0.19}$ at $z=7.54$. Finally, we combine the QSO damping wing constraints for all four $z\gtrsim7$ QSOs to obtain a single, unified constraint of $\bar{x}_{\rm HI} = 0.49\substack{+0.11 \\ -0.11}$ at $z=7.29$.

preprint2021arXiv

Paving the Way for Euclid and JWST via Optimal Selection of High-z Quasars

We introduce a probabilistic approach to select 6<z<8 quasar candidates for spectroscopic follow-up, which is based on density estimation in the high-dimensional space inhabited by the optical and near-infrared photometry. Density distributions are modeled as Gaussian mixtures with principled accounting of errors using the extreme deconvolution (XD) technique, generalizing an approach successfully used to select lower redshift (z<3) quasars. We train the probability density of contaminants on 733,694 7-d flux measurements from the 1076 square degrees overlapping area from the DECaLS (z), VIKING (YJHK), and unWISE (W1W2) imaging surveys, after requiring they dropout of DECaLS g and r, whereas the distribution of high-z quasars are trained on synthetic model photometry. Extensive simulations based on these density distributions and current estimates of the quasar luminosity function indicate that this method achieves a completeness of >75% and an efficiency of >15% for selecting quasars at 6<z<8 with J<21.5. Among the classified sources are 8 known 6<z<7 quasars, of which 2/8 are selected suggesting a completeness ~25%, whereas classifying the 6 known (J<21.5) quasars at z>7 from the entire sky, we select 5/6 or a completeness of ~80%.The failure to select the majority of 6<z<7 quasars arises because our model of quasar SEDs underestimates the scatter in the distribution of fluxes. This new optimal approach to quasar selection paves the way for efficient spectroscopic follow-up of Euclid quasar candidates with ground based telescopes and JWST.

preprint2021arXiv

Probing Early Super-massive Black Hole Growth and Quasar Evolution with Near-infrared Spectroscopy of 37 Reionization-era Quasars at 6.3 < z <= 7.64

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $λ_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects.

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 search for dust and molecular gas in enormous Ly$α$ nebulae at $z\approx 2$

Enormous Ly$α$ nebulae, extending over 300-500\,kpc around quasars, represent the pinnacle of galaxy and cluster formation. Here we present IRAM Plateau de Bure Interferometer observations of the enormous Ly$α$ nebulae `Slug' ($z$=$2.282$) and `Jackpot' ($z$=$2.041$). Our data reveal bright, synchrotron emission associated with the two radio-loud AGN embedded in the targeted nebulae, as well as molecular gas, as traced via the CO(3-2) line, in three galaxies (two sources in the Slug, and one in the Jackpot). All of the CO emission is associated with galaxies detected in their rest-frame UV stellar emission. The total mass in molecular gas of these three galaxies [$\sim (3-5)\times10^{10}$ M$_\odot$] is comparable with the total ionized gas mass responsible for the diffuse nebular emission. Our observations place limits on the molecular gas emission in the nebulae: The molecular gas surface density is $Σ_{\rm H2}<12-25$ M$_\odot$ pc$^{-2}$ for the Slug nebula and $Σ_{\rm H2}<34-68$ M$_\odot$ pc$^{-2}$ for the Jackpot nebula. These are consistent with the expected molecular gas surface densities, as predicted via photoionization models of the rest-frame UV line emission in the nebulae, and via Ly$α$ absorption in the Jackpot nebula. Compared to other radio--loud quasars at $z>1$, and high-redshift radio--loud galaxies, we do not find any strong trends relating the molecular gas reservoirs, the radio power, and the Ly$α$ luminosities of these systems. The significant step in sensitivity required to achieve a detection of the molecular gas from the nebulae, if present, will require a substantial time investment with JVLA, NOEMA, or ALMA.

preprint2020arXiv

A Significantly Neutral Intergalactic Medium Around the Luminous z=7 Quasar J0252-0503

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$α$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $λ_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$α$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

preprint2020arXiv

Detecting and Characterizing Young Quasars I: Systemic Redshifts and Proximity Zones Measurements

In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes (SMBHs) in the early universe, and pose challenges on current black hole formation models to explain the rapid formation of billion solar mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission line luminosities, and star formation rates of their host galaxies. Combined with previous results we estimate the fraction of young objects within the high-redshift quasar population at large to be $5\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

preprint2020arXiv

No Redshift Evolution in the Broad Line Region Metallicity up to z=7.54: Deep NIR Spectroscopy of ULAS J1342+0928

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

preprint2020arXiv

Pōniuā'ena: A Luminous $z=7.5$ Quasar Hosting a 1.5 Billion Solar Mass Black Hole

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("Pōniuā'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$α$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $μ$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

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 Accretion Physics of Supermassive Black Holes at Redshift z~7 with Chandra and Infrared Observations

X-ray emission from quasars has been detected up to redshift $z=7.5$, although only limited to a few objects at $z>6.5$. In this work, we present new Chandra observations of five $z>6.5$ quasars. By combining with archival Chandra observations of six additional $z>6.5$ quasars, we perform a systematic analysis on the X-ray properties of these earliest accreting supermassive black holes (SMBHs). We measure the black hole masses, bolometric luminosities ($L_{\rm bol}$), Eddington ratios ($λ_{\rm Edd}$), emission line properties, and infrared luminosities ($L_{\rm IR}$) of these quasars using infrared and sub-millimeter observations. Correlation analysis indicates that the X-ray bolometric correction (the factor that converts from X-ray luminosity to bolometric luminosity) decreases with increasing $L_{\rm bol}$, and that the UV/optical-to-X-ray ratio, $α_{\rm ox}$, strongly correlates with $L_{\rm 2500}$, and moderately correlates with $λ_{\rm Edd}$ and blueshift of CIV emission lines. These correlations are consistent with those found in lower-$z$ quasars, indicating quasar accretion physics does not evolve with redshift. We also find that $L_{\rm IR}$ does not correlate with $L_{\rm 2-10 keV}$ in these luminous distant quasars, suggesting that the ratio of the SMBH growth rate and their host galaxy growth rate in these early luminous quasars are different from those of local galaxies. A joint spectral analysis of the X-ray detected $z>6.5$ quasars yields an average X-ray photon index of $ Γ=2.32^{+0.31}_{-0.30}$, steeper than that of low-$z$ quasars. By comparing it with the $Γ-λ_{\rm Edd}$ relation, we conclude that the steepening of $Γ$ for quasars at $z>6.5$ is mainly driven by their higher Eddington ratios.

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

The REQUIEM Survey I: A Search for Extended Ly-Alpha Nebular Emission Around 31 z>5.7 Quasars

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

preprint2019arXiv

Time-dependent behaviour of quasar proximity zones at $z \sim 6$

Since the discovery of $z\sim 6$ quasars two decades ago, studies of their Ly$α$-transparent proximity zones have largely focused on their utility as a probe of cosmic reionization. But even when in a highly ionized intergalactic medium, these zones provide a rich laboratory for determining the timescales that govern quasar activity and the concomitant growth of their supermassive black holes. In this work, we use a suite of 1D radiative transfer simulations of quasar proximity zones to explore their time-dependent behaviour for activity timescales from $\sim10^3$ to $10^8$ years. The sizes of the simulated proximity zones, as quantified by the distance at which the smoothed Ly$α$ transmission drops below 10% (denoted $R_p$), are in excellent agreement with observations, with the exception of a handful of particularly small zones that have been attributed to extremely short $\lesssim 10^4$ year lifetimes. We develop a physically motivated semi-analytic model of proximity zones which captures the bulk of their equilibrium and non-equilibrium behaviour, and use this model to investigate how quasar variability on $\lesssim10^5$ year timescales is imprinted on the distribution of observed proximity zone sizes. We show that large variations in the ionizing luminosity of quasars on timescales of $\lesssim10^4$ years are disfavored based on the good agreement between the observed distribution of $R_p$ and our model prediction based on "lightbulb" (i.e. steady constant emission) light curves.

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

Constraining the dynamical importance of hot gas and radiation pressure in quasar outflows using emission line ratios

Quasar feedback models often predict an expanding hot gas bubble which drives a galaxy-scale outflow. In many circumstances this hot gas radiates inefficiently and is therefore difficult to observe directly. We present an indirect method to detect the presence of a hot bubble using hydrostatic photoionization calculations of the cold (~10^4 K) line-emitting gas. We compare our calculations with observations of the broad line region, the inner face of the torus, the narrow line region (NLR), and the extended NLR, and thus constrain the hot gas pressure at distances 0.1 pc -- 10 kpc from the center. We find that emission line ratios observed in the average quasar spectrum are consistent with radiation-pressure-dominated models on all scales. On scales <40 pc a dynamically significant hot gas pressure is ruled out, while on larger scales the hot gas pressure cannot exceed six times the local radiation pressure. In individual quasars, ~25% of quasars exhibit NLR ratios that are inconsistent with radiation-pressure-dominated models, though in these objects the hot gas pressure is also unlikely to exceed the radiation pressure by an order of magnitude or more. The derived upper limits on the hot gas pressure imply that the instantaneous gas pressure force acting on galaxy-scale outflows falls short of the time-averaged force needed to explain the large momentum fluxes \dot{p} >> L_AGN/c inferred for galaxy-scale outflows. This apparent discrepancy can be reconciled if optical quasars previously experienced a buried, fully-obscured phase during which the hot gas bubble was more effectively confined and during which galactic wind acceleration occurred.

preprint2016arXiv

Modeling the Lyman-alpha Forest in Collisionless Simulations

Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present Iteratively Matched Statistics (IMS), a novel method to accurately model the Lyman-alpha forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyman-alpha forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can perfectly reproduce line-of-sight observables, such as the PDF and power spectrum, and accurately reproduce the 3D flux power spectrum (5-20%). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic "mock" skies for Lyman-alpha forest surveys.

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

Characterizing the Pressure Smoothing Scale of the Intergalactic Medium

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Ly-alpha forest is degenerate with the 3-D structure of the IGM on ~100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T~$10^4$ K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected Jeans filtering cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Ly-alpha flux, $F_\mathrm{real}$, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The $F_\mathrm{real}$ power spectrum is accurately described by a simple fitting function with cutoff at $λ_F$, allowing us to rigorously quantify the filtering scale for the first time: we find $λ_F$ = 79 kpc (comoving) at z=3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Ly-alpha forest absorption in close quasar pairs, recently proposed as a method to measure the filtering scale. Our results enable one to quantify the filtering scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude $T_0$ and slope $γ$ of the temperature-density relation $T = T_0Δ^{γ-1}$ should be augmented with a third filtering scale parameter $λ_F$.

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

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

Spatially Resolving the Kinematics of the <100 μas Quasar Broad Line Region using Spectroastrometry

The broad line region (BLR) of luminous active galactic nuclei (AGN) is a prominent observational signature of the accretion flow around supermassive black holes, which can be used to measure their masses (M_BH) over cosmic history. Due to the <100 μas angular size of the BLR, current direct constraints on BLR kinematics are limited to those provided by reverberation mapping studies, which are most efficiently carried out on low-luminosity L and low-redshift z AGN. We analyze the possibility to measure the BLR size and study its kinematic structure using spectroastrometry, whereby one measures the spatial position centroid of emission line photons as a function of velocity. We calculate the expected spectroastrometric signal of a rotation-dominated BLR for various assumptions about the ratio of random to rotational motions, and the radial distribution of the BLR gas. We show that for hyper-luminous quasars at z < 2.5, the size of the low-ionization BLR can already be constrained with existing telescopes and adaptive optics systems, thus providing a novel method to spatially resolve the kinematics of the accretion flow at 10^3 -- 10^4 gravitational radii, and measure M_BH at the high-L end of the AGN family. With a 30m-class telescope, BLR spectroastrometry should be routinely detectable for much fainter quasars out to z ~ 6, and for various emission lines. This will enable kinematic M_BH measurements as a function of luminosity and redshift, providing a compelling science case for next generation telescopes.

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

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

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

Observational Requirements for Lyman-alpha Forest Tomographic Mapping of Large-Scale Structure at z ~ 2

The z > 2 Lyman-alpha (Lya) forest traces the underlying dark-matter distribution on large scales and, given sufficient sightlines, can be used to create 3D maps of large-scale structure. We examine the observational requirements to construct such maps and estimate the signal-to-noise as a function of exposure time and sightline density. Sightline densities at z = 2.25 are n_los = [360, 1200,3300] deg^{-2} at limiting magnitudes of g =[24.0, 24.5,25.0], resulting in transverse sightline separations of d_perp = [3.6, 1.9, 1.2] h^{-1} Mpc, which roughly sets the reconstruction scale. We simulate these reconstructions using mock spectra with realistic noise properties, and find that spectra with S/N = 4 per angstrom can be used to generate maps that clearly trace the underlying dark-matter at overdensities of rho/<rho> ~ 1. For the VLT/VIMOS spectrograph, exposure times t_exp = [4, 6, 10] hrs are sufficient for maps with spatial resolution epsilon_3d = [5.0, 3.2, 2.3] h^{-1} Mpc. Assuming ~ 250 h^{-1} Mpc is probed along the line-of-sight, 1 deg^2 of survey area would cover a comoving volume of ~ 10^6 h^{-3} Mpc^3 at <z>=2.3, enabling efficient mapping of large volumes with 8-10m telescopes. These maps could be used to study galaxy environments, detect proto-clusters, and study the topology of large-scale structure at high-z.

preprint2014arXiv

Protocluster Discovery in Tomographic Ly$α$ Forest Flux Maps

We present a new method of finding protoclusters using tomographic maps of Ly$α$ Forest flux. We review our method of creating tomographic flux maps and discuss our new high performance implementation, which makes large reconstructions computationally feasible. Using a large N-body simulation, we illustrate how protoclusters create large-scale flux decrements, roughly 10 $h^{-1}$Mpc across, and how we can use this signal to find them in flux maps. We test the performance of our protocluster finding method by running it on the ideal, noiseless map and tomographic reconstructions from mock surveys, and comparing to the halo catalog. Using the noiseless map, we find protocluster candidates with about 90% purity, and recover about 75% of the protoclusters that form massive clusters ($> 3 \times 10^{14} \, h^{-1} M_{\odot}$). We construct mock surveys similar to the ongoing COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) survey. While the existing data has an average sightline separation of 2.3 $h^{-1}$Mpc, we test separations of 2 - 6 $h^{-1}$Mpc to see what can be tolerated for our application. Using reconstructed maps from small separation mock surveys, the protocluster candidate purity and completeness are very close what was found in the noiseless case. As the sightline separation increases, the purity and completeness decrease, although they remain much higher than we initially expected. We extended our test cases to mock surveys with an average separation of 15 $h^{-1}$Mpc, meant to reproduce high source density areas of the BOSS survey. We find that even with such a large sightline separation, the method can still be used to find some of the largest protoclusters.

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.

preprint2013arXiv

A New Method to Directly Measure the Jeans Scale of the Intergalactic Medium Using Close Quasar Pairs

Although the baryons in the intergalactic medium (IGM) trace dark matter on Mpc scales, small-scale (~100 kpc) fluctuations are suppressed by pressure support, analogous to the classical Jeans argument. This Jeans filtering scale has fundamental cosmological implications: it provides a thermal record of heat injected by UV photons during reionization events, determines the clumpiness of the IGM, and sets the minimum mass scale for gravitational collapse, a key quantity in galaxy formation. Unfortunately, it is extremely challenging to measure via the standard analysis of purely longitudinal Lyman-alpha forest spectra, because the thermal Doppler broadening of absorption lines is highly degenerate with Jeans pressure smoothing. In this work we show that the Jeans scale can be directly measured by characterizing the coherence of correlated Lyman-alpha absorption in quasar pairs with separations small enough to resolve it. We present a novel technique for this purpose, based on the probability distribution function (PDF) of phase angle differences of homologous longitudinal Fourier modes in close quasar pair spectra. A Bayesian formalism is introduced based on the phase angle PDF, and MCMC techniques are used to characterize the precision of a future Jeans scale measurement, and explore degeneracies with other thermal parameters governing the IGM. A semi-analytical model of the IGM is used to generate a grid of 500 thermal models from a dark matter simulation. Our full parameter study indicates that a realistic sample of only 20 close quasar pair spectra can pinpoint the Jeans scale to ~ 5% precision, independent of the parameters governing the temperature-density relation of the IGM. We show that this new method is insensitive to a battery of systematics such as continuum fitting errors, imprecise knowledge of the noise and spectral resolution, and metal-line absorption.

preprint2013arXiv

Chemical constraints on the contribution of Population III stars to cosmic reionization

Recent studies have highlighted that galaxies at z = 6-8 fall short of producing enough ionizing photons to reionize the IGM, and suggest that Population III stars could resolve this tension, because their harder spectra can produce ~10x more ionizing photons than Population II. But this argument depends critically on the duration of the Population III era, and because Population III stars form from pristine gas, in turn depends on the rate of galactic enrichment. We use a semi-analytic model of galaxy formation which tracks galactic chemical evolution, to gauge the impact of Population III stars on reionization. Population III SNe produce distinct metal abundances, and we argue that the duration of the Population III era can be constrained by precise relative abundance measurements in high-z damped Lyα absorbers (DLAs), which provide a chemical record of past star-formation. We find that a single generation of Population III stars can self-enrich galaxies above the critical metallicity Zcrit=10^-4 Zsun for the Population III-to-II transition, on a very short timescale of ~10^6 yr, owing to the large metal yields and short lifetimes of Population III stars. This subsequently terminates the Population III era, hence they contribute >~ 50% of the ionizing photons only for z >~ 30, and at z=10 contribute <1%. The Population III contribution can be increased by delaying metal mixing into the ISM. However comparing the resulting metal abundance pattern to existing measurements in z <~ 6 DLAs, we show that the fractional contribution of high-mass Population III stars to the ionization rate must be <~ 10% at z = 10. Future abundance measurements of z~7-8 QSOs and GRBs should probe the era when the chemical vestiges of Population III star formation become detectable.

preprint2013arXiv

Chemical enrichment of Damped Lyman Alpha systems as a direct constraint on Population III star formation

Damped Ly-alpha absorbers (DLAs) can be used to measure gas-phase metallicities at large cosmological lookback times with high precision. Relative abundances can still be measured accurately deep into the reionization epoch (z > 6) using transitions redward of Ly-alpha. Here we study chemical evolution of DLAs using a constrained model for evolution of galaxies and IGM to determine the degree to which DLA abundance measurements can probe Population III enrichment. We find that if the critical metallicity of Population III to II transition is < 10^-4 Zsun, the cosmic Population III SFR is zero for z<8. Nevertheless, at high redshift (z ~ 6) Population III chemical signatures are retained in low-mass galaxies (halo mass < 10^9 Msun). This is because photoionization feedback suppresses star formation in these galaxies until relatively low redshift (z ~ 10), and the chemical record of early Population III star formation is retained. We model DLAs as these galaxies by assigning to them a mass-dependent H I absorption cross section and predict distribution of DLA abundance ratios. We find that these distributions are anchored towards abundance ratios set by Population II yields, but exhibit a tail that depends on the Population III IMF for z > 5. Thus, a sample of DLA abundance measurements at high redshift holds the promise to constrain Population III IMF. A sample of just 10 DLAs with relative abundances measured to an accuracy of 0.1 dex is sufficient to constrain the Population III IMF at 4-sigma. These constraints may prove stronger than other probes such as metal-poor stars and individual DLAs. Our results provide a global picture of the cosmic thermal, ionization, and chemical evolution, and can rule out certain Population III scenarios.

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

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

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.

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

HST narrow-band search for extended Ly-alpha emission around two z>6 quasars

We search for extended Ly-alpha emission around two z>6 quasars, SDSS J1030+0524 (z=6.309) and SDSS J1148+5251 (z=6.419) using WFC3 narrow-band filters on board the Hubble Space Telescope. For each quasar, we collected two deep, narrow-band images, one sampling the Ly-alpha line+continuum at the quasar redshifts and one of the continuum emission redwards of the line. After carefully modeling the Point Spread Function, we find no evidence for extended Ly-alpha emission. These observations set 2-sigma limits of L(Ly-alpha, extended) < 3.2 x 10^{44} erg/s for J1030+0524 and L(Ly-alpha, extended) < 2.5 x 10^{44} erg/s for J1148+5251. Given the star formation rates typically inferred from (rest-frame) far-infrared measurements of z~6 quasars, these limits are well below the intrinsic bright Ly-alpha emission expected from the recombination of gas photoionized by the quasars or by the star formation in the host galaxies, and point towards significant Ly-alpha suppression or dust attenuation. However, small extinction values have been observed along the line of sight to the nuclei, thus reddening has to be coupled with other mechanisms for Ly-alpha suppression (e.g., resonance scattering). No Ly-alpha emitting companions are found, down to a 5-sigma sensitivity of ~ 1 x 10^{-17} erg/s/cm^2/arcsec^2 (surface brightness) and ~ 5 x 10^{-17} erg/s/cm^2 (assuming point sources).

preprint2012arXiv

Photometric redshifts and quasar probabilities from a single, data-driven generative model

We describe a technique for simultaneously classifying and estimating the redshift of quasars. It can separate quasars from stars in arbitrary redshift ranges, estimate full posterior distribution functions for the redshift, and naturally incorporate flux uncertainties, missing data, and multi-wavelength photometry. We build models of quasars in flux-redshift space by applying the extreme deconvolution technique to estimate the underlying density. By integrating this density over redshift one can obtain quasar flux-densities in different redshift ranges. This approach allows for efficient, consistent, and fast classification and photometric redshift estimation. This is achieved by combining the speed obtained by choosing simple analytical forms as the basis of our density model with the flexibility of non-parametric models through the use of many simple components with many parameters. We show that this technique is competitive with the best photometric quasar classification techniques---which are limited to fixed, broad redshift ranges and high signal-to-noise ratio data---and with the best photometric redshift techniques when applied to broadband optical data. We demonstrate that the inclusion of UV and NIR data significantly improves photometric quasar--star separation and essentially resolves all of the redshift degeneracies for quasars inherent to the ugriz filter system, even when included data have a low signal-to-noise ratio. For quasars spectroscopically confirmed by the SDSS 84 and 97 percent of the objects with GALEX UV and UKIDSS NIR data have photometric redshifts within 0.1 and 0.3, respectively, of the spectroscopic redshift; this amounts to about a factor of three improvement over ugriz-only photometric redshifts. Our code to calculate quasar probabilities and redshift probability distributions is publicly available.

preprint2012arXiv

The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

preprint2012arXiv

The clustering of intermediate redshift quasars as measured by the Baryon Oscillation Spectroscopic Survey

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

preprint2011arXiv

A Simple Likelihood Method for Quasar Target Selection

We present a new method for quasar target selection using photometric fluxes and a Bayesian probabilistic approach. For our purposes we target quasars using Sloan Digital Sky Survey (SDSS) photometry to a magnitude limit of g=22. The efficiency and completeness of this technique is measured using the Baryon Oscillation Spectroscopic Survey (BOSS) data, taken in 2010. This technique was used for the uniformly selected (CORE) sample of targets in BOSS year one spectroscopy to be realized in the 9th SDSS data release. When targeting at a density of 40 objects per sq-deg (the BOSS quasar targeting density) the efficiency of this technique in recovering z>2.2 quasars is 40%. The completeness compared to all quasars identified in BOSS data is 65%. This paper also describes possible extensions and improvements for this technique

preprint2011arXiv

Combined strong and weak lensing analysis of 28 clusters from the Sloan Giant Arcs Survey

We study the mass distribution of a sample of 28 galaxy clusters using strong and weak lensing observations. The clusters are selected via their strong lensing properties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs is combined with weak lensing measurements from deep Subaru/Suprime-cam images to primarily obtain robust constraints on the concentration parameter and the shape of the mass distribution. We find that the concentration c_vir is a steep function of the mass, c_vir \propto M_vir^-0.59\pm0.12, with the value roughly consistent with the lensing-bias-corrected theoretical expectation for high mass (10^15 h^-1 M_sun) clusters. However, the observationally inferred concentration parameters appear to be much higher at lower masses (10^14 h^-1 M_sun), possibly a consequence of the modification to the inner density profiles provided by baryon cooling. The steep mass-concentration relation is also supported from direct stacking analysis of the tangential shear profiles. In addition, we explore the two-dimensional shape of the projected mass distribution by stacking weak lensing shear maps of individual clusters with prior information on the position angle from strong lens modelling, and find significant evidence for a large mean ellipticity with the best-fit value of e = 0.47 \pm 0.06 for the mass distribution of the stacked sample. We find that the luminous cluster member galaxy distribution traces the overall mass distribution very well, although the distribution of fainter cluster galaxies appears to be more extended than the total mass.

preprint2011arXiv

Gemini/GMOS Spectroscopy of 26 Strong Lensing Selected Galaxy Cluster Cores

We present results from a spectroscopic program targeting 26 strong lensing cluster cores that were visually identified in the Sloan Digital Sky Survey (SDSS) and the Second Red-Sequence Cluster Survey (RCS-2). The 26 galaxy cluster lenses span a redshift range of 0.2 < z < 0.65, and our spectroscopy reveals 69 unique background sources with redshifts as high as z=5.200. We also identify redshifts for 262 cluster member galaxies and measure the velocity dispersions and dynamical masses for 18 clusters where we have redshifts for N \geq 10 cluster member galaxies. We include an accounting for the expected biases in dynamical masses of strong lensing selected clusters as predicted by results from numerical simulations and discuss possible sources of bias in our observations. The median dynamical mass of the 18 clusters with N \geq 10 spectroscopic cluster members is M_{Vir} = 7.84 x 10^14 M_sun h_{0.7}^{-1}, which is somewhat higher than predictions for strong lensing selected clusters in simulations. The disagreement is not significant considering the large uncertainty in our dynamical data, systematic uncertainties in the velocity dispersion calibration, and limitations of the theoretical modeling. Nevertheless our study represents an important first step toward characterizing large samples of clusters that are identified in a systematic way as systems exhibiting dramatic strong lensing features.

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 Redshift Distribution of Giant Arcs in the Sloan Giant Arcs Survey

We measure the redshift distribution of a sample of 28 giant arcs discovered as a part of the Sloan Giant Arcs Survey (SGAS). Gemini/GMOS-North spectroscopy provides precise redshifts for 24 arcs, and "redshift desert" constraints for the remaining four. This is a direct measurement of the redshift distribution of a uniformly selected sample of bright giant arcs, which is an observable that can be used to inform efforts to predict giant arc statistics. Our primary giant arc sample has a median redshift z=1.821 and nearly two thirds of the arcs - 64% - are sources at z \gtrsim 1.4, indicating that the population of background sources that are strongly lensed into bright giant arcs resides primarily at high redshift. We also analyze the distribution of redshifts for 19 secondary strongly lensed background sources that are not visually apparent in SDSS imaging, but were identified in deeper follow-up imaging of the lensing cluster fields. Our redshift sample for the secondary sources is not spectroscopically complete, but combining it with our primary giant arc sample suggests that a large fraction of all background galaxies which are strongly lensed by foreground clusters reside at z \gtrsim 1.4. Kolmogorov-Smirnov (KS) tests indicate that our well-selected, spectroscopically complete primary giant arc redshift sample can be reproduced with a model distribution that is constructed from a combination of results from studies of strong lensing clusters in numerical simulations, and observational constraints on the galaxy luminosity function.

preprint2011arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: Quasar Target Selection for Data Release Nine

The SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), a five-year spectroscopic survey of 10,000 deg^2, achieved first light in late 2009. One of the key goals of BOSS is to measure the signature of baryon acoustic oscillations in the distribution of Ly-alpha absorption from the spectra of a sample of ~150,000 z>2.2 quasars. Along with measuring the angular diameter distance at z\approx2.5, BOSS will provide the first direct measurement of the expansion rate of the Universe at z > 2. One of the biggest challenges in achieving this goal is an efficient target selection algorithm for quasars over 2.2 < z < 3.5, where their colors overlap those of stars. During the first year of the BOSS survey, quasar target selection methods were developed and tested to meet the requirement of delivering at least 15 quasars deg^-2 in this redshift range, out of 40 targets deg^-2. To achieve these surface densities, the magnitude limit of the quasar targets was set at g <= 22.0 or r<=21.85. While detection of the BAO signature in the Ly-alpha absorption in quasar spectra does not require a uniform target selection, many other astrophysical studies do. We therefore defined a uniformly-selected subsample of 20 targets deg^-2, for which the selection efficiency is just over 50%. This "CORE" subsample will be fixed for Years Two through Five of the survey. In this paper we describe the evolution and implementation of the BOSS quasar target selection algorithms during the first two years of BOSS operations. We analyze the spectra obtained during the first year. 11,263 new z>2.2 quasars were spectroscopically confirmed by BOSS. Our current algorithms select an average of 15 z > 2.2 quasars deg^-2 from 40 targets deg^-2 using single-epoch SDSS imaging. Multi-epoch optical data and data at other wavelengths can further improve the efficiency and completeness of BOSS quasar target selection. [Abridged]

preprint2011arXiv

Think Outside the Color Box: Probabilistic Target Selection and the SDSS-XDQSO Quasar Targeting Catalog

We present the SDSS-XDQSO quasar targeting catalog for efficient flux-based quasar target selection down to the faint limit of the Sloan Digital Sky Survey (SDSS) catalog, even at medium redshifts (2.5 <~ z <~ 3) where the stellar contamination is significant. We build models of the distributions of stars and quasars in flux space down to the flux limit by applying the extreme-deconvolution method to estimate the underlying density. We convolve this density with the flux uncertainties when evaluating the probability that an object is a quasar. This approach results in a targeting algorithm that is more principled, more efficient, and faster than other similar methods. We apply the algorithm to derive low-redshift (z < 2.2), medium-redshift (2.2 <= z <= 3.5), and high-redshift (z > 3.5) quasar probabilities for all 160,904,060 point sources with dereddened i-band magnitude between 17.75 and 22.45 mag in the 14,555 deg^2 of imaging from SDSS Data Release 8. The catalog can be used to define a uniformly selected and efficient low- or medium-redshift quasar survey, such as that needed for the SDSS-III's Baryon Oscillation Spectroscopic Survey project. We show that the XDQSO technique performs as well as the current best photometric quasar-selection technique at low redshift, and outperforms all other flux-based methods for selecting the medium-redshift quasars of our primary interest. We make code to reproduce the XDQSO quasar target selection publicly available.

preprint2010arXiv

Selecting Quasars by their Intrinsic Variability

We present a new and simple technique for selecting extensive, complete and pure quasar samples, based on their intrinsic variability. We parametrize the single-band variability by a power-law model for the light-curve structure function, with amplitude A and power-law index gamma. We show that quasars can be efficiently separated from other non-variable and variable sources by the location of the individual sources in the A-gamma plane. We use ~60 epochs of imaging data, taken over ~5 years, from the SDSS stripe 82 (S82) survey, where extensive spectroscopy provides a reference sample of quasars, to demonstrate the power of variability as a quasar classifier in multi-epoch surveys. For UV-excess selected objects, variability performs just as well as the standard SDSS color selection, identifying quasars with a completeness of 90% and a purity of 95%. In the redshift range 2.5<z<3, where color selection is known to be problematic, variability can select quasars with a completeness of 90% and a purity of 96%. This is a factor of 5-10 times more pure than existing color-selection of quasars in this redshift range. Selecting objects from a broad griz color box without u-band information, variability selection in S82 can afford completeness and purity of 92%, despite a factor of 30 more contaminants than quasars in the color-selected feeder sample. This confirms that the fraction of quasars hidden in the 'stellar locus' of color-space is small. To test variability selection in the context of Pan-STARRS 1 (PS1) we created mock PS1 data by down-sampling the S82 data to just 6 epochs over 3 years. Even with this much sparser time sampling, variability is an encouragingly efficient classifier. For instance, a 92% pure and 44% complete quasar candidate sample is attainable from the above $griz$-selected catalog.

preprint2010arXiv

The Sloan Digital Sky Survey Quasar Lens Search. IV. Statistical Lens Sample from the Fifth Data Release

We present the second report of our systematic search for strongly lensed quasars from the data of the Sloan Digital Sky Survey (SDSS). From extensive follow-up observations of 136 candidate objects, we find 36 lenses in the full sample of 77,429 spectroscopically confirmed quasars in the SDSS Data Release 5. We then define a complete sample of 19 lenses, including 11 from our previous search in the SDSS Data Release 3, from the sample of 36,287 quasars with i<19.1 in the redshift range 0.6<z<2.2, where we require the lenses to have image separations of 1"<θ<20" and i-band magnitude differences between the two images smaller than 1.25 mag. Among the 19 lensed quasars, 3 have quadruple-image configurations, while the remaining 16 show double images. This lens sample constrains the cosmological constant to be Ω_Λ=0.84^{+0.06}_{-0.08}(stat.)^{+0.09}_{-0.07}(syst.) assuming a flat universe, which is in good agreement with other cosmological observations. We also report the discoveries of 7 binary quasars with separations ranging from 1.1" to 16.6", which are identified in the course of our lens survey. This study concludes the construction of our statistical lens sample in the full SDSS-I data set.

preprint2010arXiv

Two Lensed Lyman-alpha Emitting Galaxies at z~5

We present observations of two strongly lensed $z\sim5$ Lyman-$α$ Emitting (LAE) galaxies that were discovered in the Sloan Giant Arcs Survey (SGAS). We identify the two sources as SGAS J091541+382655, at $z=5.200$, and SGAS J134331+415455 at $z=4.994$. We measure their AB magnitudes at $(i,z)=(23.34\pm0.09,23.29\pm0.13$) mags and $(i,z)=(23.78\pm0.18,24.24^{+0.18}_{-0.16}$) mags, and the rest-frame equivalent widths of the Lyman-$α$ emission at $25.3\pm4.1$Å~and $135.6\pm20.3$Å~for SGAS J091541+382655 and SGAS J134331+415455, respectively. Each source is strongly lensed by a massive galaxy cluster in the foreground, and the magnifications due to gravitational lensing are recovered from strong lens modeling of the foreground lensing potentials. We use the magnification to calculate the intrinsic, unlensed Lyman-$α$ and UV continuum luminosities for both sources, as well as the implied star formation rates (SFR). We find SGAS J091541+382655 and SGAS J134341+415455 to be galaxies with (L$_{Ly-α}$, L$_{UV})\leq(0.6$L$_{Ly-α}^{*}, 2$L$_{UV}^{*}$) and (L$_{Ly-α}$, L$_{UV})=(0.5$L$_{Ly-α}^{*}, 0.9$L$_{UV}^{*}$), respectively. Comparison of the spectral energy distributions (SEDs) of both sources against stellar population models produces estimates of the mass in young stars in each galaxy: we report an upper limit of M$_{stars} \leq 7.9^{+3.7}_{-2.5} \times 10^{7}$ M$_{\sun} h_{0.7}^{-1}$ for SGAS J091531+382655, and a range of viable masses for SGAS J134331+415455 of $2\times10^{8}$ M$_{\sun} h_{0.7}^{-1} <$ M$_{stars} < 6\times10^{9}$ M$_{\sun} h_{0.7}^{-1}$.

preprint2010arXiv

Two Lensed z~3 Lyman Break Galaxies Discovered in the SDSS Giant Arcs Survey

We report the discovery of two strongly-lensed z ~ 3 Lyman Break Galaxies (LBGs) discovered as u-band dropouts as part of the SDSS Giant Arcs Survey (SGAS). The first, SGAS J122651.3+215220 at z = 2.9233 is lensed by one of several sub-clusters, SDSS J1226+2152, in a complex massive cluster at z = 0.43. Its (g, r, i) magnitudes are (21.14, 20.60, 20.51) which translate to surface brightnesses, mu_{g,r,i}, of (23.78, 23.11, 22.81). The second, SGAS J152745.1+065219, is an LBG at z = 2.7593 lensed by the foreground SDSS J1527+0652 at z = 0.39, with (g, r, z)=(20.90, 20.52, 20.58) and mu_{g,r,z}=(25.15, 24.52, 24.12). Moderate resolution spectroscopy confirms the redshifts suggested by photometric breaks, and shows both absorption and emission features typical of LBGs. Lens mass models derived from combined imaging and spectroscopy reveal that SGAS J122651.3+215220 is a highly magnified source (M ~40), while SGAS J152745.1+065219 is magnified by no more than M ~ 15. Compared to LBG survey results (Steidel et al. 2003), the luminosities and lensing-corrected magnitudes suggest that SGAS J122651.3+215220 is among the faintest 20% of LBGs in that sample. SGAS J152745.1+065219, on the other hand, appears to be more representative of the average LBG, similar to the "Cosmic Eye".

preprint2009arXiv

Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4

The clustering of quasars on small scales yields fundamental constraints on models of quasar evolution and the buildup of supermassive black holes. This paper describes the first systematic survey to discover high redshift binary quasars. Using color-selection and photometric redshift techniques, we searched 8142 deg^2 of SDSS imaging data for binary quasar candidates, and confirmed them with follow-up spectroscopy. Our sample of 27 high redshift binaries (24 of them new discoveries) at redshifts 2.9 < z < 4.3 with proper transverse separations 10 kpc < R_{\perp} < 650 kpc increases the number of such objects known by an order of magnitude. Eight members of this sample are very close pairs with R_{\perp} < 100 kpc, and of these close systems four are at z > 3.5. The completeness and efficiency of our well-defined selection algorithm are quantified using simulated photometry and we find that our sample is ~ 50% complete. Our companion paper uses this knowledge to make the first measurement of the small scale clustering (R < 1 Mpc/h comoving) of high-redshift quasars. High redshift binaries constitute exponentially rare coincidences of two extreme (M >~ 10^9 Msun) supermassive black holes. At z ~ 4 there is about one close binary per 10 Gpc^3, thus these could be the highest sigma peaks, the analogs of superclusters, in the early Universe.

preprint2009arXiv

Binary Quasars at High Redshift II: Sub-Mpc Clustering at z ~ 3-4

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

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

Subaru Weak Lensing Measurements of Four Strong Lensing Clusters: Are Lensing Clusters Over-Concentrated?

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

preprint2008arXiv

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

Using a homogenous sample of 38,208 quasars with a sky coverage of $4000 {\rm deg^2}$ drawn from the SDSS Data Release Five quasar catalog, we study the dependence of quasar clustering on luminosity, virial black hole mass, quasar color, and radio loudness. At $z<2.5$, quasar clustering depends weakly on luminosity and virial black hole mass, with typical uncertainty levels $\sim 10%$ for the measured correlation lengths. These weak dependences are consistent with models in which substantial scatter between quasar luminosity, virial black hole mass and the host dark matter halo mass has diluted any clustering difference, where halo mass is assumed to be the relevant quantity that best correlates with clustering strength. However, the most luminous and most massive quasars are more strongly clustered (at the $\sim 2σ$ level) than the remainder of the sample, which we attribute to the rapid increase of the bias factor at the high-mass end of host halos. We do not observe a strong dependence of clustering strength on quasar colors within our sample. On the other hand, radio-loud quasars are more strongly clustered than are radio-quiet quasars matched in redshift and optical luminosity (or virial black hole mass), consistent with local observations of radio galaxies and radio-loud type 2 AGN. Thus radio-loud quasars reside in more massive and denser environments in the biased halo clustering picture. Using the Sheth et al.(2001) formula for the linear halo bias, the estimated host halo mass for radio-loud quasars is $\sim 10^{13} h^{-1}M_\odot$, compared to $\sim 2\times 10^{12} h^{-1}M_\odot$ for radio-quiet quasar hosts at $z\sim 1.5$.

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.

preprint2006arXiv

QSO Absorption Lines from QSOs

We present the results of a search for metal absorption lines in the spectra of background QSOs whose sightlines pass close to foreground QSOs. We detect MgII(2796,2803) absorption in Sloan Digital Sky Survey (SDSS) spectra of four z>1.5 QSOs whose lines of sight pass within 26-98 kpc of lower redshift (z~0.5-1.5) QSOs. The 100% [4/4 pairs] detection of MgII in the background QSOs is clearly at odds with the incidence of associated (z_abs ~ z_em) systems -- absorbers which exist towards only a few percent of QSOs. Although the quality of our foreground QSO spectra is not as high as the SDSS data, absorption seen towards one of the background QSOs clearly does not show up at the same strength in the spectrum of the corresponding foreground QSO. This implies that the absorbing gas is distributed inhomogeneously around the QSO, presumably as a direct consequence of the anisotropic emission from the central AGN. We discuss possible origins for the MgII lines, including: absorption by gas from the foreground QSO host galaxy; companion galaxies fuelling the QSO through gravitational interactions; and tidal debris left by galaxy mergers or interactions which initiated the QSO activity. No single explanation is entirely satisfactory, and we may well be seeing a mixture of phenomena.

preprint2006arXiv

Quasars Probing Quasars I: Optically Thick Absorbers Near Luminous 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 search 149 moderate resolution background quasar spectra, from Gemini, Keck, the MMT, and the SDSS to survey Lyman Limit Systems (LLSs) and Damped Ly-alpha systems (DLAs) in the vicinity of 1.8 < z < 4.0 luminous foreground quasars. A sample of 27 new quasar-absorber pairs is uncovered with column densities, 17.2 < log (N_HI/cm^2) < 20.9, and transverse (proper) distances of 22 kpc/h < R < 1.7 Mpc/h, from the foreground quasars. If they emit isotropically, the implied ionizing photon fluxes are a factor of ~ 5-8000 times larger than the ambient extragalactic UV background over this range of distances. The observed probability of intercepting an absorber is very high for small separations: six out of eight projected sightlines with transverse separations R < 150 kpc/h have an absorber coincident with the foreground quasar, of which four have log N_HI > 10^19. The covering factor of log N_HI > 10^19 absorbers is thus ~ 50 % (4/8) on these small scales, whereas < 2% would have been expected at random. There are many cosmological applications of these new sightlines: they provide laboratories for studying fluorescent Ly-alpha recombination radiation from LLSs, constrain the environments, emission geometry, and radiative histories of quasars, and shed light on the physical nature of LLSs and DLAs.

preprint2006arXiv

The DEEP2 Galaxy Redshift Survey: Clustering of Quasars and Galaxies at z=1

We present the clustering of DEEP2 galaxies at 0.7<z<1.4 around quasars identified using both the SDSS and DEEP2 surveys. We measure the two-point cross-correlation of a sample of 36 optically-selected, spectroscopically-identified quasars from the SDSS and 16 more found in the DEEP2 survey with the full DEEP2 galaxy sample over scales 0.1 < r_p < 10 Mpc/h. The clustering amplitude is found to be similar to the auto-correlation function of DEEP2 galaxies, with a relative bias of b=0.89 +/-0.24 between quasars and DEEP2 galaxies at z~1. No significant dependence is found on scale, quasar luminosity, or redshift over the ranges we probe here. The clustering amplitude errors are comparable to those from significantly larger quasar samples, such as the 2dF QSO Redshift Survey. This results from the statistical power of cross-correlation techniques, which exploit the fact that galaxies are much more numerous than quasars. We also measure the local environments of quasars using the 3rd-nearest-neighbor surface density of surrounding DEEP2 galaxies. Quasars are found in regions of similar mean overdensity as blue DEEP2 galaxies; they differ in environment from the red DEEP2 galaxy population at 2 sigma significance. Our results imply that quasars do not reside in particularly massive dark matter halos at these redshifts, with a mean dark matter halo mass of M_200 ~ 3 10^12 solar masses in a concordance LCDM cosmology.

preprint2005arXiv

A Survey of z>5.7 Quasars in the Sloan Digital Sky Survey IV: Discovery of Seven Additional Quasars

We present the discovery of seven quasars at z>5.7, selected from ~2000 deg^2 of multicolor imaging data of the Sloan Digital Sky Survey (SDSS). The new quasars have redshifts z from 5.79 to 6.13. Five are selected as part of a complete flux-limited sample in the SDSS Northern Galactic Cap; two have larger photometric errors and are not part of the complete sample. One of the new quasars, SDSS J1335+3533 (z=5.93), exhibits no emission lines; the 3-sigma limit on the rest-frame equivalent width of Ly alpha+NV line is 5 A. It is the highest redshift lineless quasar known, and could be a gravitational lensed galaxy, a BL Lac object or a new type of quasar. Two new z>6 quasars, SDSS 1250+3130 (z=6.13) and SDSS J1137+3549 (z=6.01), show deep Gunn-Peterson absorption gaps in Ly alpha. These gaps are narrower the complete Gunn-Peterson absorption troughs observed among quasars at z>6.2 and do not have complete Ly beta absorption.

preprint2005arXiv

Binary Quasars in the Sloan Digital Sky Survey: Evidence for Excess Clustering on Small Scales

We present a sample of 218 new quasar pairs with proper transverse separations R_prop < 1 Mpc/h over the redshift range 0.5 < z < 3.0, discovered from an extensive follow up campaign to find companions around the Sloan Digital Sky Survey and 2dF Quasar Redshift Survey quasars. This sample includes 26 new binary quasars with separations R_prop < 50 kpc/h (theta < 10 arcseconds), more than doubling the number of such systems known. We define a statistical sample of binaries selected with homogeneous criteria and compute its selection function, taking into account sources of incompleteness. The first measurement of the quasar correlation function on scales 10 kpc/h < R_prop < 400 kpc/h is presented. For R_prop < 40 kpc/h, we detect an order of magnitude excess clustering over the expectation from the large scale R_prop > 3 Mpc/h quasar correlation function, extrapolated down as a power law to the separations probed by our binaries. The excess grows to ~ 30 at R_prop ~ 10 kpc/h, and provides compelling evidence that the quasar autocorrelation function gets progressively steeper on sub-Mpc scales. This small scale excess can likely be attributed to dissipative interaction events which trigger quasar activity in rich environments. Recent small scale measurements of galaxy clustering and quasar-galaxy clustering are reviewed and discussed in relation to our measurement of small scale quasar clustering.

preprint2005arXiv

Discovery of a Fifth Image of the Large Separation Gravitationally Lensed Quasar SDSS J1004+4112

We report the discovery of a fifth image in the large separation lensed quasar system SDSS J1004+4112. A faint point source located 0.2'' from the center of the brightest galaxy in the lensing cluster is detected in images taken with the Advanced Camera for Surveys (ACS) and the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) on the Hubble Space Telescope. The flux ratio between the point source and the brightest lensed component in the ACS image is similar to that in the NICMOS image. The location and brightness of the point source are consistent with lens model predictions for a lensed image. We therefore conclude that the point source is likely to be a fifth image of the source quasar. In addition, the NICMOS image reveals the lensed host galaxy of the source quasar, which can strongly constrain the structure of the lensing critical curves and thereby the mass distribution of the lensing cluster.

preprint2004arXiv

Microlensing of the Broad Emission Line Region in the Quadruple Lens SDSS J1004+4112

We present seven epochs of spectroscopy on the quadruply imaged quasar SDSS J1004+4112, spanning observed-frame time delays from 1 to 322 days. The spectra reveal differences in the emission lines between the lensed images. Specifically, component A showed a strong enhancement in the blue wings of several high-ionization lines relative to component B, which lasted at least 28 days (observed frame) then faded. Since the predicted time delay between A and B is <30 days, our time coverage suggests that the event was not intrinsic to the quasar. We attribute these variations to microlensing of part of the broad emission line region of the quasar, apparently resolving structure in the source plane on a scale of ~10^{16} cm at z=1.734. In addition, we observed smaller differences in the emission line profiles between components A and B that persisted throughout the time span, which may also be due to microlensing or millilensing. Further spectroscopic monitoring of this system holds considerable promise for resolving the structure of the broad emission line region in quasars.

Joseph F. Hennawi

What is connected

Connect this record

See the researcher in context

Building this map preview

79 published item(s)

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

(Nearly) Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at $z\sim 5-7$ Using Dark Pixel Fractions in Ly$α$ and Ly$β$ Forests

Constraining IGM enrichment and metallicity with the C IV forest correlation function

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

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference

A Luminous Quasar at Redshift 7.642

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

ALMA reveals large overdensity and strong clustering of galaxies in quasar environments at $z\sim4$

IGM damping wing constraints on reionisation from covariance reconstruction of two $z\gtrsim7$ QSOs

Paving the Way for Euclid and JWST via Optimal Selection of High-z Quasars

Probing Early Super-massive Black Hole Growth and Quasar Evolution with Near-infrared Spectroscopy of 37 Reionization-era Quasars at 6.3 < z <= 7.64

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

A search for dust and molecular gas in enormous Ly$α$ nebulae at $z\approx 2$

A Significantly Neutral Intergalactic Medium Around the Luminous z=7 Quasar J0252-0503

Detecting and Characterizing Young Quasars I: Systemic Redshifts and Proximity Zones Measurements

No Redshift Evolution in the Broad Line Region Metallicity up to z=7.54: Deep NIR Spectroscopy of ULAS J1342+0928

Pōniuā'ena: A Luminous $z=7.5$ Quasar Hosting a 1.5 Billion Solar Mass Black Hole

PypeIt: The Python Spectroscopic Data Reduction Pipeline

Revealing the Accretion Physics of Supermassive Black Holes at Redshift z~7 with Chandra and Infrared Observations

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$

The REQUIEM Survey I: A Search for Extended Ly-Alpha Nebular Emission Around 31 z>5.7 Quasars

Time-dependent behaviour of quasar proximity zones at $z \sim 6$

A Universal Density Structure for Circum-Galactic Gas

Constraining the dynamical importance of hot gas and radiation pressure in quasar outflows using emission line ratios

Modeling the Lyman-alpha Forest in Collisionless Simulations

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

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

The UV-bright Quasar Survey (UVQS): DR1

Characterizing the Pressure Smoothing Scale of the Intergalactic Medium

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

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

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

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

Spatially Resolving the Kinematics of the <100 μas Quasar Broad Line Region using Spectroastrometry

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

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

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

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

Observational Requirements for Lyman-alpha Forest Tomographic Mapping of Large-Scale Structure at z ~ 2

Protocluster Discovery in Tomographic Ly$α$ Forest Flux Maps

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

A New Method to Directly Measure the Jeans Scale of the Intergalactic Medium Using Close Quasar Pairs

Chemical constraints on the contribution of Population III stars to cosmic reionization

Chemical enrichment of Damped Lyman Alpha systems as a direct constraint on Population III star formation

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

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

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

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

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

HST narrow-band search for extended Ly-alpha emission around two z>6 quasars

Photometric redshifts and quasar probabilities from a single, data-driven generative model

The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9

The clustering of intermediate redshift quasars as measured by the Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey quasar catalog: ninth data release

A Simple Likelihood Method for Quasar Target Selection

Combined strong and weak lensing analysis of 28 clusters from the Sloan Giant Arcs Survey

Gemini/GMOS Spectroscopy of 26 Strong Lensing Selected Galaxy Cluster Cores

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

The Redshift Distribution of Giant Arcs in the Sloan Giant Arcs Survey

The SDSS-III Baryon Oscillation Spectroscopic Survey: Quasar Target Selection for Data Release Nine

Think Outside the Color Box: Probabilistic Target Selection and the SDSS-XDQSO Quasar Targeting Catalog

Selecting Quasars by their Intrinsic Variability

The Sloan Digital Sky Survey Quasar Lens Search. IV. Statistical Lens Sample from the Fifth Data Release

Two Lensed Lyman-alpha Emitting Galaxies at z~5

Two Lensed z~3 Lyman Break Galaxies Discovered in the SDSS Giant Arcs Survey

Binary Quasars at High Redshift I: 24 New Quasar Pairs at z ~ 3-4

Binary Quasars at High Redshift II: Sub-Mpc Clustering at z ~ 3-4

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

Subaru Weak Lensing Measurements of Four Strong Lensing Clusters: Are Lensing Clusters Over-Concentrated?

Quasar Clustering from SDSS DR5: Dependences on Physical Properties

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

QSO Absorption Lines from QSOs

Quasars Probing Quasars I: Optically Thick Absorbers Near Luminous Quasars