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}$.

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.

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'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$\alpha$ 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 $\mu$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.

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.

We present spatially-resolved spectroscopy from the Keck Cosmic Web Imager (KCWI) of a star-forming galaxy at z=0.6942, which shows emission from the Mg II 2796, 2803 Angstrom doublet in the circumgalactic medium (CGM) extending ~37 kpc at 3-sigma significance in individual spaxels (1-sigma detection limit 4.8 x 10^{-19} erg s^-1 cm^-2 arcsec^-2). After deconvolution with the seeing, we obtain 5-sigma detections extending for ~31 kpc measured in 7-spaxel (1.1 arcsec^2) apertures. Spaxels covering the galaxy stellar regions show clear P-Cygni-like emission/absorption profiles with the blueshifted absorption extending to relative velocities of v = -800 km/s; however, the P-Cygni profiles give way to pure emission at large radii from the central galaxy. We have performed three-dimensional radiative transfer modeling to infer the geometry and velocity and density profiles of the outflowing gas. Our observations are most consistent with an isotropic outflow rather than biconical wind models with half-opening angles phi <= 80 deg. Furthermore, our modeling suggests that a wind velocity profile that decreases with radius is necessary to reproduce the velocity widths and strengths of Mg II line emission profiles at large circumgalactic radii. The extent of the Mg II emission we measure directly is further corroborated by our modeling, where we rule out outflow models with extent <30 kpc.

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$\alpha$ 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$\alpha$ 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.

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$\alpha$ 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 $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ 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.

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.

Motivated by the recent discovery of the near-ubiquity of Ly$\alpha$ emission around $z \gtrsim 3$ QSOs, we performed a systematic study of QSO circumgalactic Ly$\alpha$ 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{\nu_{\rm{LL}}}}=10^{31.1-32.3}$ erg s$^{-1}$ Hz$^{-1}$. We found that 14 out of 16 QSOs are associated with Ly$\alpha$ nebulae with projected linear-sizes larger than 50 physical kpc (pkpc). Among them, four nebulae have enormous Ly$\alpha$ emission with the Ly$\alpha$ surface brightness $SB_{\rm{Ly\alpha}} >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$\alpha$ SB profile (SB$_{\rm{Ly\alpha}}$) of QSOs at $z\approx 2$ can be described by a power law of SB$_{\rm{Ly\alpha},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$.

We use MUSE/VLT to conduct a survey of $z\sim3$ physical quasar pairs at close separation with a fast observation strategy. Our aim is twofold: (i) explore the Ly$\alpha$ glow around the faint-end of the quasar population; (ii) take advantage of the combined illumination of a quasar pair to unveil large-scale intergalactic structures extending between the two quasars. Here, we report the results for a quasar pair ($z=3.020,3.008$; $i=21.84,22.15$), separated by 11.6 arcsec (or 89 projected kpc). MUSE reveals filamentary Ly$\alpha$ structures extending between the two quasars with an average surface brightness of SB$_{\rm Ly\alpha}=1.8\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. Photoionization models of the constraints in the Ly$\alpha$, HeII, and CIV line emissions show that the emitting structures are intergalactic bridges with an extent between $\sim89$ and up to $\sim600$ kpc. Our models rule out the possibility that the structure extends for $\sim 2.9$ Mpc, i.e., the separation inferred from the uncertain systemic redshift difference of the quasars if the difference was only due to the Hubble flow. At the current spatial resolution and surface brightness limit, the average projected width of an individual bridge is about 35 kpc. We also detect a strong absorption in HI, NV, and CIV along the background sight-line at higher $z$, which we interpret as due to at least two components of cool, metal enriched, and relatively ionized CGM or IGM surrounding the quasar pair. Two additional HI absorbers are detected along both quasar sight-lines at $\sim -900$ and $-2800$ km s$^{-1}$ from the system, with the latter having associated CIV absorption only along the foreground quasar sight-line. The absence of galaxies in the MUSE field of view at the redshifts of these two absorbers suggests that they trace large-scale structures or expanding shells in front of the quasar pair.

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

The supermassive black holes (SMBHs) observed at the centers of all massive galaxies are believed to have grown via luminous accretion during quasar phases in the distant past. The fraction of inflowing rest mass energy emitted as light, the radiative efficiency, has been inferred to be 10\%, in agreement with expectations from thin disk accretion models. But the existence of billion solar-mass SMBHs powering quasars at $z > 7$ challenges this picture: provided they respect the Eddington limit, there is not enough time to grow $z>7$ SMBHs from stellar remnant seeds unless the radiative efficiency is below 10\%. Here we show that one can constrain the radiative efficiencies of the most distant quasars known using foreground neutral intergalactic gas as a cosmological-scale ionizing photon counter. From the Ly$\alpha$ absorption profiles of ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$), we determine posterior median radiative efficiencies of 0.08\% and 0.1\%, respectively, and the combination of the two measurements rule out the canonical 10\% efficiency at 99.8\% credibility after marginalizing over the unknown obscured fraction. This low radiative efficiency implies rapid mass accretion for the earliest SMBHs, greatly easing the tension between the age of the Universe and the SMBH masses. However, our measured efficiency may instead reflect nearly complete obscuration by dusty gas in the quasar host galaxies over the vast majority of their SMBH growth. Assuming 10\% efficiency during unobscured phases, we find that the obscured fraction would be $>82\%$ at 95\% credibility, and imply a $25.7^{+49.6}_{-16.5}$ times larger obscured than unobscured luminous quasar population at $z>7$.

We measure the intergalactic medium (IGM) opacity in the Ly$\alpha$ as well as in the Ly$\beta$ forest along $19$ quasar sightlines between $5.5\lesssim z_{\rm abs}\lesssim 6.1$, probing the end stages of the reionization epoch. Owing to its lower oscillator strength the Ly$\beta$ transition is sensitive to different gas temperatures and densities than Ly$\alpha$, providing additional constraints on the ionization and thermal state of the IGM. A comparison of our measurements to different inhomogeneous reionization models, derived from post-processing the Nyx cosmological hydrodynamical simulation to include spatial fluctuations in the ultraviolet background (UVB) or the gas temperature field, as well as to a uniform reionization model with varying thermal states of the IGM, leads to two primary conclusions: First, we find that including the effects of spectral noise is key for a proper data to model comparison. Noise effectively reduces the sensitivity to high opacity regions, and thus even stronger spatial inhomogeneities are required to match the observed scatter in the observations than previously inferred. Second, we find that models which come close to reproducing the distribution of Ly$\alpha$ effective optical depths nevertheless underpredict the Ly$\beta$ opacity at the same spatial locations. The origin of this disagreement is not entirely clear but models with an inversion in the temperature-density relation of the IGM just after reionization is completed match our measurements best, although they still do not fully capture the observations at $z\gtrsim 5.8$.

Determinations of the UV luminosity function of AGN at high redshifts are important for constraining the AGN contribution to reionization and understanding the growth of supermassive black holes. Recent inferences of the luminosity function suffer from inconsistencies arising from inhomogeneous selection and analysis of AGN data. We address this problem by constructing a sample of more than 80,000 colour-selected AGN from redshift z=0 to 7.5. While this sample is composed of multiple data sets with spectroscopic redshifts and completeness estimates, we homogenise these data sets to identical cosmologies, intrinsic AGN spectra, and magnitude systems. Using this sample, we derive the AGN UV luminosity function from redshift z=0 to 7.5. The luminosity function has a double power law form at all redshifts. The break magnitude $M_*$ of the AGN luminosity function shows a steep brightening from $M_*\sim -24$ at z=0.7 to $M_*\sim -29$ at z=6. The faint-end slope $\beta$ significantly steepens from $-1.7$ at $z<2.2$ to $-2.4$ at $z\simeq 6$. In spite of this steepening, the contribution of AGN to the hydrogen photoionization rate at $z\sim 6$ is subdominant (< 3%), although it can be non-negligible (~10%) if these luminosity functions hold down to $M_{1450}=-18$. Under reasonable assumptions, AGN can reionize HeII by redshift z=2.9. At low redshifts (z<0.5), AGN can produce about half of the hydrogen photoionization rate inferred from the statistics of HI absorption lines in the IGM. Our global analysis of the luminosity function also reveals important systematic errors in the data, particularly at z=2.2--3.5, which need to be addressed and incorporated in the AGN selection function in future in order to improve our results. We make various fitting functions, luminosity function analysis codes, and homogenised AGN data publicly available.

In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters and combines a semi-numerical reionization model to solve for the topology of reionization and an approximate model of how reionization heats the IGM, with the massively parallel \texttt{Nyx} hydrodynamics code, specifically designed to solve for the structure of diffuse IGM gas. We have produced several large-scale high resolution cosmological hydrodynamical simulations ($2048^3$, $L_{\rm box} = 40$ Mpc/h) with different instantaneous and inhomogeneous HI reionization models that use this new methodology. We study the IGM thermal properties of these models and find that large scale temperature fluctuations extend well beyond the end of reionization. Analyzing the 1D flux power spectrum of these models, we find up to $\sim 50\%$ differences in the large scale properties (low modes, $k\lesssim0.01$ s/km) of the post-reionization power spectrum due to the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios already with existing Ly$\alpha$ forest measurements. Finally, we explore the differences in the small-scale cutoff of the power spectrum and we find that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model happens at the midpoint of the inhomogeneous model.

We perform a Very Large Telescope FOcal Reducer and low dispersion Spectrograph 2 (VLT/FORS2) narrowband imaging search around 5 star-forming galaxies at redshift z=0.67-0.69 in the Great Observatories Origins Deep Survey South (GOODS-S) field to constrain the radial extent of large-scale outflows traced by resonantly scattered MgII emission. The sample galaxies span star formation rates in the range 4 $M_{\odot}/yr$ < SFR < $40\ M_{\odot}/yr$ and have stellar masses $9.9 \lesssim \log M_{*}/M_{\odot} \lesssim 11.0$, and exhibit outflows traced by MgII absorption with velocities ~150-420 km s$^{-1}$ . These observations are uniquely sensitive, reaching surface brightness limits of 5.81 $\times$ $10^{-19}$ ergs sec $^{-1}$ cm$^{-2}$ arcsec$^2$ per 1 arcsec$^2$ aperture (at 5$\sigma$ significance). We do not detect any extended emission around any of the sample galaxies, thus placing 5$\sigma$ upper limits on the brightness of extended MgII emission of $<6.51 \times 10^{-19}$ ergs sec $^{-1}$ cm$^{-2}$ arcsec$^2$ at projected distances $R_{\perp} > 8-21$ kpc. The imaging also resolves the MgII absorption observed toward each galaxy spatially, revealing approximately constant absorption strengths across the galaxy disks. In concert with radiative transfer models predicting the surface brightness of MgII emission for a variety of simple wind morphologies, our detection limits suggest that either (1) the extent of the MgII-emitting material in the outflows from these galaxies is limited to $\lesssim 20$ kpc; or (2) the outflows are anisotropic and/or dusty.

We characterise ionised gas outflows using a large sample of ~330 high-luminosity (45.5 < log(L_bol/erg s^-1) < 49.0), high-redshift (1.5 < z < 4.0) quasars via their [OIII]4960,5008 emission. The median velocity width of the [OIII] emission line is 1540 kms^-1, increasing with increasing quasar luminosity. Broad, blue-shifted wings are seen in the [OIII] profiles of 42 per cent of the sample. Rest-frame ultraviolet spectra with well-characterised CIV 1550 emission line properties are available for more than 210 quasars, allowing an investigation of the relationship between the Broad Line Region (BLR) and Narrow Line Region (NLR) emission properties. The [OIII] blueshift is correlated with CIV blueshift, even when the dependence of both quantities on quasar luminosity has been taken into account. A strong anti-correlation between the [OIII] equivalent width (EW) and CIV blueshift also exists. Furthermore, [OIII] is very weak, with EW<1A, in ~10 per cent of the sample, a factor of 10 higher compared to quasars at lower luminosities and redshifts. If the [OIII] emission originates in an extended NLR, the observations suggest that quasar-driven winds are capable of influencing the host-galaxy environment out to kilo-parsec scales. The mean kinetic power of the ionised gas outflows is then 10^44.7 erg s^-1, which is ~0.15 per cent of the bolometric luminosity of the quasar. These outflow efficiencies are broadly consistent with those invoked in current active galactic nuclei feedback models.

The latest measurements of CMB electron scattering optical depth reported by Planck significantly reduces the allowed space of HI reionization models, pointing toward a later ending and/or less extended phase transition than previously believed. Reionization impulsively heats the intergalactic medium (IGM) to $\sim10^4$ K, and owing to long cooling and dynamical times in the diffuse gas, comparable to the Hubble time, memory of reionization heating is retained. Therefore, a late ending reionization has significant implications for the structure of the $z\sim5-6$ Lyman-$\alpha$ (ly$\alpha$) forest. Using state-of-the-art hydrodynamical simulations that allow us to vary the timing of reionization and its associated heat injection, we argue that extant thermal signatures from reionization can be detected via the Ly$\alpha$ forest power spectrum at $5< z<6$. This arises because the small-scale cutoff in the power depends not only the the IGM temperature at these epochs, but is also particularly sensitive to the pressure smoothing scale set by the IGM full thermal history. Comparing our different reionization models with existing measurements of the Ly$\alpha$ forest flux power spectrum at $z=5.0-5.4$, we find that models satisfying Planck's $\tau_e$, constraint favor a moderate amount of heat injection consistent with galaxies driving reionization, but disfavoring quasar-driven scenarios. We explore the impact of different reionization histories and heating models on the shape of the power spectrum, and find that they can produce similar effects, but argue that this degeneracy can be broken with high enough quality data. We study the feasibility of measuring the flux power spectrum at $z\simeq 6$ using mock quasar spectra and conclude that a sample of $\sim10$ high-resolution spectra with an attainable signal-to-noise ratio will allow distinguishing between different reionization scenarios.

We report high spatial resolution (~0.076", 410pc) Atacama Large Millimeter/submillimeter Array imaging of the dust continuum and the ionized carbon line [CII] in a luminous quasar host galaxy at z=6.6, 800 million years after the big bang. Based on previous studies, this galaxy hosts a ~1x10^9 M_sun black hole and has a star-formation rate of ~1500 M_sun/yr. The unprecedented high resolution of the observations reveals a complex morphology of gas within 3kpc of the accreting central black hole. The gas has a high velocity dispersion with little ordered motion along the line of sight, as would be expected from gas accretion that has yet to settle in a disk. In addition, we find the presence of [CII] cavities in the gas distribution (with diameters of ~0.5kpc), offset from the central black hole. This unique distribution and kinematics cannot be explained by a simple model. Plausible scenarios are that the gas is located in a truncated or warped disk, or the holes are created by interactions with nearby galaxies or due to energy injection into the gas. In the latter case, the energy required to form the cavities must originate from the central active galactic nucleus, as the required energy far exceeds the energy output expected from supernovae. This energy input into the gas, however, does not inhibit the high rate of star-formation. Both star-formation and black hole activity could have been triggered by interactions with satellite galaxies; our data reveal three additional companions detected in [CII] emission around the quasar.

The strong observed clustering of $z>3.5$ quasars indicates they are hosted by massive ($M_{\rm{halo}}\gtrsim10^{12}\,h^{-1}\,\rm{M_{\odot}}$) dark matter halos. Assuming quasars and galaxies trace the same large-scale structures, this should also manifest as strong clustering of galaxies around quasars. Previous work on high-redshift quasar environments, mostly focused at $z>5$, have failed to find convincing evidence for these overdensities. Here we conduct a survey for Lyman alpha emitters (LAEs) in the environs of 17 quasars at $z\sim4$ probing scales of $R\lesssim7\,h^{-1}\,{\rm{Mpc}}$. We measure an average LAE overdensity around quasars of 1.4 for our full sample, which we quantify by fitting the quasar-LAE cross-correlation function. We find consistency with a power-law shape with correlation length of $r^{QG}_{0}=2.78^{+1.16}_{-1.05}\,h^{-1}\,{\rm{cMpc}}$ for a fixed slope of $\gamma=1.8$. We also measure the LAE auto-correlation length and find $r^{GG}_{0}=9.12^{+1.32}_{-1.31}\,h^{-1}$\,cMpc ($\gamma=1.8$), which is $3.3$ times higher than the value measured in blank fields. Taken together our results clearly indicate that LAEs are significantly clustered around $z\sim4$ quasars. We compare the observed clustering with the expectation from a deterministic bias model, whereby LAEs and quasars probe the same underlying dark matter overdensities, and find that our measurements fall short of the predicted overdensities by a factor of 2.1. We discuss possible explanations for this discrepancy including large-scale quenching or the presence of excess dust in galaxies near quasars. Finally, the large cosmic variance from field-to-field observed in our sample (10/17 fields are actually underdense) cautions one from over-interpreting studies of $z\sim6$ quasar environments based on a single or handful of quasar fields.

We present new measurements of the flux power-spectrum P(k) of the $z<0.5$ HI Lyman-$\alpha$ forest spanning scales k ~ 0.001-0.1 s/km. These results were derived from 65 far ultraviolet quasar spectra (resolution R~18000) observed with the Cosmic Origin Spectrograph (COS) on board the Hubble Space Telescope. The analysis required careful masking of all contaminating, coincident absorption from HI and metal-line transitions of the Galactic interstellar medium and intervening absorbers as well as proper treatment of the complex COS line-spread function. From the P(k) measurements, we estimate the HI photoionization rate ($\Gamma_{\rm HI}$) in the z<0.5 intergalactic medium. Our results confirm most of the previous $\Gamma_{\rm HI}$ estimates. We conclude that previous concerns of a photon underproduction crisis are now resolved by demonstrating that the measured $\Gamma_{\rm HI}$ can be accounted for by ultraviolet emission from quasars alone. In a companion paper, we will present constraints on the thermal state of the $z<0.5$ intergalactic medium from the P(k) measurements presented here.

We present measurements of the large-scale (~40 comoving Mpc) effective optical depth of HeII Ly$\alpha$ absorption, $\tau_{\rm eff}$, at 2.54<z<3.86 toward 16 HeII-transparent quasars observed with the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST), to characterize the ionization state of helium in the intergalactic medium (IGM). We provide the first statistical sample of $\tau_{\rm eff}$ measurements in six signal-to-noise ratio >3 HeII sightlines at z>3.5, and study the redshift evolution and sightline-to-sightline variance of $\tau_{\rm eff}$ in 24 HeII sightlines. We confirm an increase of the median $\tau_{\rm eff}$ from ~2 at z=2.7 to >5 at z>3, and a scatter in $\tau_{\rm eff}$ that increases with redshift. The z>3.5 HeII absorption is predominantly saturated, but isolated narrow ($\Delta v<650$ km/s) transmission spikes indicate patches of reionized helium. We compare our measurements to predictions for a range of UV background models applied to outputs of a large-volume (146 comoving Mpc)$^3$ hydrodynamical simulation by forward-modeling our sample's quality and size. At z>2.74 the variance in $\tau_{\rm eff}$ significantly exceeds expectations for a spatially uniform UV background, but is consistent with a fluctuating radiation field sourced by variations in the quasar number density and the mean free path in the post-reionization IGM. We develop a method to infer the approximate median HeII photoionization rate $\Gamma_{\rm HeII}$ of a fluctuating UV background from the median $\tau_{\rm eff}$, finding a factor ~5 decrease in $\Gamma_{\rm HeII}$ between z~2.6 and z~3.1. At z~3.1 a $\Gamma_{\rm HeII}=\left[9.1^{+1.1}_{-1.2}\,\mathrm{(stat.)}\,^{+2.4}_{-3.4}\,\mathrm{(sys.)}\right]\times10^{-16}\,\mathrm{s}^{-1}$ corresponds to a median HeII fraction of ~2.5%, indicating that our data probe the tail end of HeII reionization.

We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$\alpha$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0 (\rho/\rho_0)^{\gamma-1}$, and previous work has inferred $T_0$ and $\gamma$ by fitting its low-$b$ cutoff. This approach discards the majority of available data, and is susceptible to systematics related to cutoff determination. We present a method that exploits all information encoded in the $b$-$N_{\text{HI}}$ distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, then use principal component decomposition to create an emulator which can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D $(T_0, \gamma)$ measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at $z=2$. Using 200 absorbers, equivalent in pathlength to a single Ly$\alpha$ forest spectrum, we measure $\log T_0 =4.092^{+0.050}_{-0.055}$ and $\gamma=1.49^{+0.073}_{-0.074}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of $\log T_0$ and $\gamma$ with precision on $\log T_0$ ($\gamma$) nearly two (three) times higher for current dataset sizes.

Since the discovery of $z\sim 6$ quasars two decades ago, studies of their Ly$\alpha$-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$\alpha$ 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.

Absorption signatures in the spectra of QSOs are one of our most powerful tools for studying galactic and intergalactic environments at high redshifts. With the discovery of QSOs out to z > 7, QSO absorption lines are now tracing the end stages of reionization on multiple fronts using the hydrogen Lyman-$\alpha$ forest and heavy element absorbers. Next-generation QSO absorption line studies with large optical/IR telescopes will reveal in detail how the first galaxies emerged form the cosmic web, transformed their circum- and inter-galactic environments, and completed the last major phase transition of the Universe. These efforts will complement other upcoming studies of reionization, such as those with JWST, ALMA, and redshifted 21cm experiments.

We determine the thermal evolution of the intergalactic medium (IGM) over $3\, \mathrm{Gyr}$ of cosmic time $1.8<z<5.4$ by comparing measurements of the Ly{\alpha} forest power spectrum to a suite of $\sim70$ hydrodynamical simulations. We conduct Bayesian inference of IGM thermal parameters using an end-to-end forward modeling framework whereby mock spectra generated from our simulation grid are used to build a custom emulator which interpolates the power spectrum between thermal grid points. The temperature at mean density $T_0$ rises steadily from $T_0\sim 6000\, \mathrm{K}$ at $z=5.4$, peaks at $14000\, \mathrm{K}$ for $z\sim 3.4$, and decreases at lower redshift reaching $T_0\sim 7000\, \mathrm{K}$ by $z\sim1.8$. This evolution provides conclusive evidence for photoionization heating resulting from the reionization of He II, as well as the subsequent cooling of the IGM due to the expansion of the Universe after all reionization events are complete. Our results are broadly consistent with previous measurements of thermal evolution based on a variety of approaches, but the sensitivity of the power spectrum, the combination of high precision BOSS measurements of large-scale modes ($k\lesssim 0.02\, \mathrm{s/km}$) with our recent determination of the small-scale power, our large grid of models, and our careful statistical analysis allow us to break the well known degeneracy between the temperature at mean density $T_0$ and the slope of the temperature density relation $\gamma$ that has plagued previous analyses. At the highest redshifts $z\geq5$ we infer lower temperatures than expected from the standard picture of IGM thermal evolution leaving little room for additional smoothing of the Ly{\alpha} forest by free streaming of warm dark matter.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

The intense radiation emitted by luminous quasars dramatically alters the ionization state of their surrounding IGM. This so-called proximity effect extends out to tens of Mpc, and manifests as large coherent regions of enhanced Lyman-$\alpha$ (Ly$\alpha$) forest transmission in absorption spectra of background sightlines. Here we present a novel method based on Ly$\alpha$ forest tomography, which is capable of mapping these quasar `light echoes' in three dimensions. Using a dense grid (10-100) of faint ($m_r\approx24.7\,\mathrm{mag}$) background galaxies as absorption probes, one can measure the ionization state of the IGM in the vicinity of a foreground quasar, yielding detailed information about the quasar's radiative history and emission geometry. An end-to-end analysis - combining cosmological hydrodynamical simulations post-processed with a quasar emission model, realistic estimates of galaxy number densities, and instrument + telescope throughput - is conducted to explore the feasibility of detecting quasar light echoes. We present a new fully Bayesian statistical method that allows one to reconstruct quasar light echoes from thousands of individual low S/N transmission measurements. Armed with this machinery, we undertake an exhaustive parameter study and show that light echoes can be convincingly detected for luminous ($M_{1450} < -27.5\,\mathrm{mag}$ corresponding to $m_{1450} < 18.4\,\mathrm{mag}$ at $z\simeq 3.6$) quasars at redshifts $3<z_\mathrm{QSO}<5$, and that a relative precision better than $20\,\%$ on the quasar age can be achieved for individual objects, for the expected range of ages between 1 Myr and 100 Myr. The observational requirements are relatively modest - moderate resolution ($R\gtrsim750$) multi object spectroscopy at low $\rm{}S/N > 5$ is sufficient, requiring three hour integrations using existing instruments on 8m class telescopes.

The duration of quasar accretion episodes is a key quantity for distinguishing between models for the formation and growth of supermassive black holes, the evolution of quasars, and their potential feedback effects on their host galaxies. However, this critical timescale, often referred to as the quasar lifetime, is still uncertain by orders of magnitude ($t_{\rm Q} \simeq 0.01\,{\rm Myr} - 1\,{\rm Gyr}$). Absorption spectra of quasars exhibiting transmission in the He II Ly$\alpha$ forest provide a unique opportunity to make precise measurements of the quasar lifetime. Indeed, the size of a quasar's He II proximity zone, the region near the quasar itself where its own radiation dramatically alters the ionization state of the surrounding intergalactic medium (IGM), depends sensitively on its lifetime for $t_{\rm Q}\lesssim 30$ Myr, comparable to the expected $e$-folding time-scale for SMBH growth $t_{\rm S}=45$ Myr. In this study we compare the sizes of He II proximity zones in the Hubble Space Telescope (HST) spectra of 6 z~4 quasars to theoretical models generated by post-processing cosmological hydrodynamical simulations with $1$D radiative transfer algorithm. We introduce a Bayesian statistical method to infer the lifetimes of individual quasars which allows us to fully marginalize over the unknown ionization state of the surrounding IGM. We measure lifetimes $0.63^{+0.82}_{-0.40}$ Myr and $5.75^{+4.72}_{-2.74}$ Myr for two objects. For the other 4 quasars large redshift uncertainties undermine our sensitivity allowing us to only place upper or lower limits. However a joint analysis of these four systems yields a measurement of their average lifetime of $\langle t_{\rm Q}\rangle = 1.17^{+1.77}_{-0.84}$ Myr. We discuss our short ~1 Myr inferred lifetimes in the context of other quasar lifetime constraints and the growth of SMBHs.

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

Motivated by the recent discovery of rare Enormous Lyman-Alpha Nebulae (ELAN) around z~2 quasars, we have initiated a long-term observational campaign with the MUSE instrument to directly uncover the astrophysics of the gas around quasars. We present here the first 61 targets of our effort under the acronym QSO MUSEUM (Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission). These quasars have a median redshfit of z=3.17, absolute $i$ magnitude in the range $-29.67\leq M_i(z=2)\leq-27.03$, and different levels of radio-loudness. This sample unveils diverse specimens of Ly$\alpha$ nebulosities extending for tens of kiloparsecs around these quasars (on average out to a maximum projected distance of 80 kpc) above a surface brightness SB$>8.8\times10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ ($2\sigma$). The bulk of the extended Ly$\alpha$ emission is within R<50 kpc, and is characterized by relatively quiescent kinematics, with average velocity dispersions of $\langle \sigma_{\rm Ly\alpha}\rangle < 400$~km~s$^{-1}$. Therefore, the motions within all these Ly$\alpha$ nebulosities have amplitudes consistent with gravitational motions expected in dark matter halos hosting quasars at these redshifts, possibly reflecting the complexity in propagating a fast wind on large scales. Our current data suggest a combination of photoionization and resonant scattering as powering mechanisms of the Ly$\alpha$ emission. We discover the first $z\sim3$ ELAN, which confirms a very low probability ($\sim1\%$) of occurrence of such extreme systems at these cosmic epochs. Finally, we discuss the redshift evolution currently seen in extended Ly$\alpha$ emission around radio-quiet quasars from $z\sim3$ to $z\sim2$, concluding that it is possibly linked to a decrease of cool gas mass within the quasars' CGM from $z\sim3$ to $z\sim2$, and thus to the balance of cool vs hot media.

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

We publicly release a new sample of $34$ medium resolution quasar spectra at $5.77\leq z_{\rm em}\leq6.54$ observed with the Echellette Spectrograph and Imager (ESI) on the Keck telescope. This quasar sample represents an ideal laboratory to study the intergalactic medium (IGM) during the end stages of the epoch of reionization, and constrain the timing and morphology of the phase transition. For a subset of $23$ of our highest signal-to-noise ratio spectra (S/N$>7$, per $10\,{\rm km\,s^{-1}}$ pixel), we present a new measurement of the Lyman-$\alpha$ (Ly$\alpha$) forest opacity spanning the redshift range $4.8\lesssim z\lesssim6.3$. We carefully eliminate spectral regions that could be causing biases in our measurements due to additional transmitted flux in the proximity zone of the quasars, or extra absorption caused by strong intervening absorption systems along the line of sight. We compare the observed evolution of the IGM opacity with redshift to predictions from a hydrodynamical simulation with uniform ultraviolet background (UVB) radiation, as well as two semi-numerical patchy reionization models, one with a fluctuating UVB and another with a fluctuating temperature field. Our measurements show a steep rise in opacity at $z\gtrsim5.0$ and an increased scatter and thus support the picture of a spatially inhomogeneous reionization process, consistent with previous work. However, we measure significantly higher optical depths at $5.3\lesssim z\lesssim5.7$ than previous studies, which reduces the contrast between the highest opacity Gunn-Peterson troughs and the average opacity trend of the IGM, which may relieve some of the previously noted tension between these measurements and reionization models.

Faint star-forming galaxies at $z\sim 2-3$ can be used as alternative background sources to probe the Lyman-$\alpha$ forest in addition to quasars, yielding high sightline densities that enable 3D tomographic reconstruction of the foreground absorption field. Here, we present the first data release from the COSMOS Lyman-Alpha Mapping And Mapping Observations (CLAMATO) Survey, which was conducted with the LRIS spectrograph on the Keck-I telescope. Over an observational footprint of 0.157$\mathrm{deg}^2$ within the COSMOS field, we used 240 galaxies and quasars at $2.17<z<3.00$, with a mean comoving transverse separation of $2.37\,h^{-1}\,\mathrm{Mpc^3}$, as background sources probing the foreground Lyman-$\alpha$ forest absorption at $2.05<z<2.55$. The Lyman-$\alpha$ forest data was then used to create a Wiener-filtered tomographic reconstruction over a comoving volume of $3.15\,\times 10^5\,h^{-3}\,\mathrm{Mpc^3}$ with an effective smoothing scale of $2.5\,h^{-1}\,\mathrm{Mpc}$. In addition to traditional figures, this map is also presented as a virtual-reality YouTube360 video visualization and manipulable interactive figure. We see large overdensities and underdensities that visually agree with the distribution of coeval galaxies from spectroscopic redshift surveys in the same field, including overdensities associated with several recently-discovered galaxy protoclusters in the volume. This data release includes the redshift catalog, reduced spectra, extracted Lyman-$\alpha$ forest pixel data, and tomographic map of the absorption.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

We present the most distant detection of cosmic voids ($z \sim 2.3$) and the first detection of three-dimensional voids in the Lyman-$\alpha$ forest. We used a 3D tomographic map of the absorption with effective comoving spatial resolution of $2.5\,h^{-1}\mathrm{Mpc}$ and volume of $3.15\times 10^5\,h^{-3}\mathrm{Mpc}^3$, which was reconstructed from moderate-resolution Keck-I/LRIS spectra of 240 background Lyman-break galaxies and quasars in a $0.16\,\mathrm{deg}^2$ footprint in the COSMOS field. Voids were detected using a spherical overdensity finder calibrated from hydrodynamical simulations of the intergalactic medium. This allows us to identify voids in the IGM corresponding to voids in the underlying matter density field, yielding a consistent volume fraction of voids in both data (19.5%) and simulations (18.2%). We fit excursion set models to the void radius function and compare the radially-averaged stacked profiles of large voids ($r > 5$ $h^{-1}$ Mpc) to stacked voids in mock observations and the simulated density field. Comparing with 432 coeval galaxies with spectroscopic redshifts in the same volume as the tomographic map, we find that the tomography-identified voids are underdense in galaxies by 5.95$\sigma$ compared to random cells.

The baryonic content and physical properties of the warm and hot ($10^5\lesssim T\lesssim 10^7$ K) phases of the circumgalactic medium (CGM) are poorly constrained, owing to the lack of observables probing the requisite range of temperature, spatial scale, halo mass, and redshift. The radiation from a luminous quasar produces a spatially extended emission halo resulting from Thomson scattering off of free electrons in the CGM, which can be used to measure the electron density profile, and therefore, the amount of warm and hot baryonic matter present. We predict the resulting surface brightness profiles and show that they are easily detectable in a three hour integration with the James Webb Space Telescope (JWST), out to $\sim 100$ physical kpc from the centers of individual hyper-luminous quasars. This electron scattering surface brightness is redshift independent, and the signal-to-noise ratio depends only very weakly on redshift, in principle allowing measurements of the warm and hot CGM into the Epoch of Reionization at $z\sim 6.5$. We consider a litany of potential contaminants, and find that for fainter quasars at $z\lesssim1$, extended stellar halos might be of comparable surface brightness. At $z>2$, JWST mid-IR observations start to probe rest-frame optical/UV wavelengths implying that scattering by dust grains in the CGM becomes significant, although multi-color observations should be able to distinguish these scenarios given that Thomson scattering is achromatic.

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

The quasar lifetime $t_{\rm\,Q}$ is one of the most fundamental quantities for understanding quasar evolution and the growth of supermassive black holes (SMBHs), but remains uncertain by several orders of magnitude. In a recent study we uncovered a population of very young quasars ($t_{\rm Q}\lesssim10^4-10^5$ yr), based on the sizes of their proximity zones, which are regions of enhanced Ly$\alpha$ forest transmission near the quasar resulting from its own ionizing radiation. The presence of such young objects poses significant challenges to models of SMBH formation, which already struggle to explain the existence of SMBHs at such early cosmic epochs. We conduct the first comprehensive spectroscopic study of the youngest quasar known, $\rm SDSS\,J1335+3533$ at $z=5.9012$, whose lifetime is $t_{\rm Q}<10^4$ yr ($95\%$ confidence). A careful search of our deep optical and near-infrared spectra for HI and metal absorption lines allows us to convincingly exclude that its small proximity zone results from an associated absorption system rather than a short lifetime. We use the MgII emission line to measure its black hole mass $M_{\rm BH}=(4.09\pm0.58)\times 10^9 M_{\odot}$, implying an Eddington ratio of $0.30\pm0.04$ -- comparable to other co-eval quasars. We similarly find that the relationship between its black hole mass and dynamical mass are consistent with other $z\sim6$ quasars. The only possible anomaly associated with youth are its weak emission lines, but larger samples are needed to shed light on a potential causal connection. We discuss the implications of short lifetimes for various SMBH growth scenarios, and argue that future observations of young quasars with JWST could distinguish between them.

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

Large OVI columns are observed around star-forming, low-redshift ~L* galaxies, with a dependence on impact parameter indicating that most O^5+ particles reside beyond half the halo virial radius (>~100 kpc). In order to constrain the nature of the gas traced by OVI, we analyze additional observables of the outer halo, namely HI to OVI column ratios of 1-10, an absence of low-ion absorption, a mean differential extinction of E(B-V)~10^-3, and a linear relation between OVI column and velocity width. We contrast these observations with two physical scenarios: (1) OVI traces high-pressure (~30 cm^-3 K) collisionally-ionized gas cooling from a virially-shocked phase, and (2) OVI traces low-pressure (<~1 cm^-3 K) gas beyond the accretion shock, where the gas is in ionization and thermal equilibrium with the UV background. We demonstrate that the high-pressure scenario requires multiple gas phases to explain the observations, and a large deposition of energy at >~100 kpc to offset the energy radiated by the cooling gas. In contrast, the low-pressure scenario can explain all considered observations with a single gas phase in thermal equilibrium, provided that the baryon overdensity is comparable to the dark-matter overdensity, and that the gas is enriched to >~Z_sun/3 with an ISM-like dust-to-metal ratio. The low-pressure scenario implies that OVI traces a cool flow with mass flow rate of ~5 M_sun yr^-1, comparable to the star formation rate of the central galaxies. The OVI line widths are consistent with the velocity shear expected within this flow. The low-pressure scenario predicts a bimodality in absorption line ratios at ~100 kpc, due to the pressure jump across the accretion shock.

The rare close projection of two quasars on the sky provides the opportunity to study the host galaxy environment of a foreground quasar in absorption against the continuum emission of a background quasar. For over a decade the "Quasars probing quasars" series has utilized this technique to further the understanding of galaxy formation and evolution in the presence of a quasar at z>2, resolving scales as small as a galactic disc and from bound gas in the circumgalactic medium to the diffuse environs of intergalactic space. Presented here, is the public release of the quasar pair spectral database utilized in these studies. In addition to projected pairs at z>2, the database also includes quasar pair members at z<2, gravitational lens candidates and quasars closely separated in redshift that are useful for small-scale clustering studies. In total the database catalogs 5627 distinct objects, with 4083 lying within 5' of at least one other source. A spectral library contains 3582 optical and near-infrared spectra for 3028 of the cataloged sources. As well as reporting on 54 newly discovered quasar pairs, we outline the key contributions made by this series over the last ten years, summarize the imaging and spectroscopic data used for target selection, discuss the target selection methodologies, describe the database content and explore some avenues for future work. Full documentation for the spectral database, including download instructions, are supplied at the following URL http://specdb.readthedocs.io/en/latest/

We examine the kinematics of the gas in the environments of galaxies hosting quasars at $z\sim2$. We employ 148 projected quasar pairs to study the circumgalactic gas of the foreground quasars in absorption. The sample selects foreground quasars with precise redshift measurements, using emission-lines with precision $\lesssim300\,{\rm km\,s^{-1}}$ and average offsets from the systemic redshift $\lesssim|100\,{\rm km\,s^{-1}}|$. We stack the background quasar spectra at the foreground quasar's systemic redshift to study the mean absorption in \ion{C}{2}, \ion{C}{4}, and \ion{Mg}{2}. We find that the mean absorptions exhibit large velocity widths $\sigma_v\approx300\,{\rm km\,s^{-1}}$. Further, the mean absorptions appear to be asymmetric about the systemic redshifts. The mean absorption centroids exhibit small redshift relative to the systemic $\delta v\approx+200\,{\rm km\,s^{-1}}$, with large intrinsic scatter in the centroid velocities of the individual absorption systems. We find the observed widths are consistent with gas in gravitational motion and Hubble flow. However, while the observation of large widths alone does not require galactic-scale outflows, the observed offsets suggest that the gas is on average outflowing from the galaxy. The observed offsets also suggest that the ionizing radiation from the foreground quasars is anisotropic and/or intermittent.

Models and observations suggest that luminous quasar activity is triggered by mergers, so it should preferentially occur in the most massive primordial dark matter haloes, where the frequency of mergers is expected to be the highest. Since the importance of galaxy mergers increases with redshift, we identify the high-redshift Universe as the ideal laboratory for studying dual AGN. Here we present the X-ray properties of two systems of dual quasars at z=3.0-3.3 selected from the SDSS-DR6 at separations of 6-8 arcsec (43-65kpc) and observed by Chandra for 65ks each. Both members of each pair are detected with good photon statistics to allow us to constrain the column density, spectral slope and intrinsic X-ray luminosity. We also include a recently discovered dual quasar at z=5 (separation of 21 arcsec, 136kpc) for which XMM-Newton archival data allow us to detect the two components separately. Using optical spectra we derived bolometric luminosities, BH masses and Eddington ratios that were compared to those of luminous SDSS quasars in the same redshift ranges. We find that the brighter component of both pairs at z=3.0-3.3 has high luminosities compared to the distribution of SDSS quasars at similar redshift, with J1622A having an order magnitude higher luminosity than the median. This source lies at the luminous end of the z~3.3 quasar luminosity function. While we cannot conclusively state that the unusually high luminosities of our sources are related to their having a close companion, for J1622A there is only a 3% probability that it is by chance.

Recently, Carnall et al. discovered two bright high redshift quasars using the combination of the VST ATLAS and WISE surveys. The technique involved using the 3-D colour plane i-z:z-W1:W1-W2 with the WISE W1 (3.4 micron) and W2 (4.5 micron) bands taking the place of the usual NIR J band to help decrease stellar dwarf contamination. Here we report on our continued search for 5.7<z<6.4 quasars over an ~2x larger area of ~3577 sq. deg. of the Southern Hemisphere. We have found two further z>6 quasars, VST-ATLAS J158.6938-14.4211 at z=6.07 and J332.8017-32.1036 at z=6.32 with magnitudes of z_AB=19.4 and 19.7 mag respectively. J158.6938-14.4211 was confirmed by Keck LRIS observations and J332.8017-32.1036 was confirmed by ESO NTT EFOSC-2 observations. Here we present VLT X-shooter Visible and NIR spectra for the four ATLAS quasars. We have further independently rediscovered two z>5.7 quasars previously found by the VIKING/KiDS and PanSTARRS surveys. This means that in ATLAS we have now discovered a total of six quasars in our target 5.7<z<6.4 redshift range. Making approximate corrections for incompleteness, we find that our quasar space density agrees with the SDSS results of Jiang et al. at M_1450A~-27mag. Preliminary virial mass estimates based on the CIV and MIII emission lines give black hole masses in the range M_BH~1-6x10e9 M_solar for the four ATLAS quasars.

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral Universe, and enables measurements of the hydrogen neutral fraction $x_{\rm HI}(z)$ at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Ly$\alpha$. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift Universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the Universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

Measuring the proximity effect and the damping wing of intergalactic neutral hydrogen in quasar spectra during the epoch of reionization requires an estimate of the intrinsic continuum at rest-frame wavelengths $\lambda_{\rm rest}\sim1200$-$1260$ {\AA}. In contrast to previous works which used composite spectra with matched spectral properties or explored correlations between parameters of broad emission lines, we opted for a non-parametric predictive approach based on principal component analysis (PCA) to predict the intrinsic spectrum from the spectral properties at redder (i.e. unabsorbed) wavelengths. We decomposed a sample of $12764$ spectra of $z\sim2$-$2.5$ quasars from SDSS/BOSS into 10 red-side ($1280$ {\AA} $<\lambda_{\rm rest}<2900$ {\AA}) and 6 blue-side ($1180$ {\AA} $<\lambda_{\rm rest}<1280$ {\AA}) PCA basis spectra, and constructed a projection matrix to predict the blue-side coefficients from a fit to the red-side spectrum. We found that our method predicts the blue-side continuum with $\sim6$-$12\%$ precision and $\lesssim1\%$ bias by testing on the full training set sample. We then computed predictions for the blue-side continua of the two quasars currently known at $z>7$: ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$). Both of these quasars are known to exhibit extreme emission line properties, so we individually calibrated the precision of the continuum predictions from similar quasars in the training set. We find that both $z>7$ quasars, and in particular ULAS J1342+0928, show signs of damping wing-like absorption at wavelengths redward of Ly$\alpha$.

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

The HeII transverse proximity effect - enhanced HeII Ly{\alpha} transmission in a background sightline caused by the ionizing radiation of a foreground quasar - offers a unique opportunity to probe the emission properties of quasars, in particular the emission geometry (obscuration, beaming) and the quasar lifetime. Building on the foreground quasar survey published in Schmidt+2017, we present a detailed model of the HeII transverse proximity effect, specifically designed to include light travel time effects, finite quasar ages, and quasar obscuration. We post-process outputs from a cosmological hydrodynamical simulation with a fluctuating HeII UV background model, plus the added effect of the radiation from a single bright foreground quasar. We vary the age $t_\mathrm{age}$ and obscured sky fractions $\Omega_\mathrm{obsc}$ of the foreground quasar, and explore the resulting effect on the HeII transverse proximity effect signal. Fluctuations in IGM density and the UV background, as well as the unknown orientation of the foreground quasar, result in a large variance of the HeII Ly{\alpha} transmission along the background sightline. We develop a fully Bayesian statistical formalism to compare far UV HeII Ly{\alpha} transmission spectra of the background quasars to our models, and extract joint constraints on $t_\mathrm{age}$ and $\Omega_\mathrm{obsc}$ for the six Schmidt+2017 foreground quasars with the highest implied HeII photoionization rates. Our analysis suggests a bimodal distribution of quasar emission properties, whereby one foreground quasar, associated with a strong HeII transmission spike, is relatively old $(22\,\mathrm{Myr})$ and unobscured $\Omega_\mathrm{obsc}<35\%$, whereas three others are either younger than $(10\,\mathrm{Myr})$ or highly obscured $(\Omega_\mathrm{obsc}>70\%)$.

Quasars are galaxies hosting accreting supermassive black holes; due to their brightness, they are unique probes of the early universe. To date, only few quasars have been reported at $z > 6.5$ ($<$800 Myr after the Big Bang). In this work, we present six additional $z \gtrsim 6.5$ quasars discovered using the Pan-STARRS1 survey. We use a sample of 15 $z \gtrsim 6.5$ quasars to perform a homogeneous and comprehensive analysis of this highest-redshift quasar population. We report four main results: (1) the majority of $z\gtrsim$6.5 quasars show large blueshifts of the broad CIV 1549\AA$\,$emission line compared to the systemic redshift of the quasars, with a median value $\sim$3$\times$ higher than a quasar sample at $z\sim$1; (2) we estimate the quasars' black hole masses (M$\rm_{BH}\sim$0.3$-$5 $\times$ 10$^{9}$ M$_{\odot}$) via modeling of the MgII 2798\AA$\,$emission line and rest-frame UV continuum; we find that quasars at high redshift accrete their material (with $\langle (L_{\mathrm{bol}}/L_{\mathrm{Edd}}) \rangle = 0.39$) at a rate comparable to a luminosity-matched sample at lower$-$redshift, albeit with significant scatter ($0.4$ dex); (3) we recover no evolution of the FeII/MgII abundance ratio with cosmic time; (4) we derive near zone sizes; together with measurements for $z\sim6$ quasars from recent work, we confirm a shallow evolution of the decreasing quasar near zone sizes with redshift. Finally, we present new millimeter observations of the [CII] 158 $\mu$m emission line and underlying dust continuum from NOEMA for four quasars, and provide new accurate redshifts and [CII]/infrared luminosities estimates. The analysis presented here shows the large range of properties of the most distant quasars.

In an effort to search for Ly$\alpha$ emission from circum- and intergalactic gas on scales of hundreds of kpc around $z\sim3$ quasars, and thus characterise the physical properties of the gas in emission, we have initiated an extensive fast-survey with the Multi Unit Spectroscopic Explorer (MUSE): Quasar Snapshot Observations with MUse: Search for Extended Ultraviolet eMission (QSO MUSEUM). In this work, we report the discovery of an enormous Ly$\alpha$ nebula (ELAN) around the quasar SDSS~J102009.99+104002.7 at $z=3.164$, which we followed-up with deeper MUSE observations. This ELAN spans $\sim297$ projected kpc, has an average Ly$\alpha$ surface brightness ${\rm SB}_{\rm Ly\alpha}\sim 6.04\times10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ (within the $2\sigma$ isophote), and is associated with an additional four, previously unknown embedded sources: two Ly$\alpha$ emitters and two faint active galactic nuclei (one Type-1 and one Type-2 quasar). By mapping at high significance the line-of-sight velocity in the entirety of the observed structure, we unveiled a large-scale coherent rotation-like pattern spanning $\sim300$ km s$^{-1}$ with a velocity dispersion of $<270$ km s$^{-1}$, which we interpret as a signature of the inspiraling accretion of substructures within the quasar's host halo. Future multiwavelength data will complement our MUSE observations, and are definitely needed to fully characterise such a complex system. None the less, our observations reveal the potential of new sensitive integral-field spectrographs to characterise the dynamical state of diffuse gas on large scales in the young Universe, and thereby witness the assembly of galaxies.

We utilize the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to search for extended Lyman-Alpha emission around the z~6.6 QSO J0305-3150. After carefully subtracting the point-spread-function, we reach a nominal 5-sigma surface brightness limit of SB = 1.9x10$^{-18}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, collapsing 5 wavelength slices centered at the expected location of the redshifted Lyman-Alpha emission (i.e. at 9256 Ang.). Current data suggest the presence (5-sigma, accounting for systematics) of a Lyman-Alpha nebula that extends for 9 kpc around the QSO. This emission is displaced and redshifted by 155 km/s with respect to the location of the QSO host galaxy traced by the [CII] emission line. The total luminosity is L = 3.0x10$^{42}$ erg/s. Our analysis suggests that this emission is unlikely to rise from optically thick clouds illuminated by the ionizing radiation of the QSO. It is more plausible that the Lyman-Alpha emission is due to fluorescence of the highly ionized optically thin gas. This scenario implies a high hydrogen volume density of n$_H$ ~ 6 cm$^{-3}$. In addition, we detect a Lyman-Alpha emitter (LAE) in the immediate vicinity of the QSO: i.e., with a projected separation of 12.5 kpc and a line-of-sight velocity difference of 560 km/s. The luminosity of the LAE is L = 2.1x10$^{42}$ erg/s and its inferred star-formation-rate is SFR ~ 1.3 M$_\odot$/yr. The probability of finding such a close LAE is one order of magnitude above the expectations based on the QSO-galaxy cross-correlation function. This discovery is in agreement with a scenario where dissipative interactions favour the rapid build-up of super-massive black holes at early Cosmic times.

Galaxy formation depends critically on the physical state of gas in the circumgalactic medium (CGM) and its interface with the intergalactic medium (IGM), determined by the complex interplay between inflows from the IGM and outflows from supernovae or AGN feedback. The average Lyman-alpha (Ly-a) absorption profile around galactic halos represents a powerful tool to probe their gaseous environments. We compare predictions from Illustris and Nyx hydrodynamical simulations with the observed absorption around foreground quasars, damped Ly-a systems, and Lyman-break galaxies. We show how large-scale BOSS and small-scale quasar pair measurements can be combined to precisely constrain the absorption profile over three decades in transverse distance 20kpc$\lesssim b\lesssim$20Mpc. Far from galaxies $\gtrsim2$Mpc, the simulations converge to the same profile and provide a reasonable match to the observations. This asymptotic agreement arises because the $\Lambda$CDM model successfully describes the ambient IGM, and represents a critical advantage of studying the mean absorption profile. However, significant differences between the simulations, and between simulations and observations are present on scales 20kpc$\lesssim b\lesssim$2Mpc, illustrating the challenges of accurately modeling and resolving galaxy formation physics. It is noteworthy that these differences are observed as far out as $\sim2$Mpc, indicating that the `sphere-of-influence' of galaxies could extend to approximately $\sim20$ times the halo virial radius ($\sim100$kpc). Current observations are very precise on these scales and can thus strongly discriminate between different galaxy formation models. We demonstrate that the Ly-a absorption profile is primarily sensitive to the underlying temperature-density relationship of diffuse gas around galaxies, and argue that it thus provides a fundamental test of galaxy formation models.

We present a new method to quantify the value of the escape fraction of ionizing photons, and the existence of ultra-faint galaxies clustered around brighter objects during the epoch of cosmic reionization, using the diffuse Ly$\alpha$, continuum and H$\alpha$ emission observed around galaxies at $z\sim6$. We model the surface brightness profiles of the diffuse halos considering the fluorescent emission powered by ionizing photons escaping from the central galaxies, and the nebular emission from satellite star-forming sources, by extending the formalisms developed in Mas-Ribas & Dijkstra (2016) and Mas-Ribas et al. (2017). The comparison between our predicted profiles and Ly$\alpha$ observations at $z=5.7$ and $z=6.6$ favors a low ionizing escape fraction, $f_{\rm esc}^{\rm ion}\sim5\%$, for galaxies in the range $-19\gtrsim M_{\rm UV} \gtrsim -21.5$. However, uncertainties and possible systematics in the observations do not allow for firm conclusions. We predict H$\alpha$ and rest-frame visible continuum observations with JWST, and show that JWST will be able to detect extended (a few tens of kpc) fluorescent H$\alpha$ emission powered by ionizing photons escaping from a bright, $L\gtrsim 5L^*$, galaxy. Such observations can differentiate fluorescent emission from nebular emission by satellite sources. We discuss how observations and stacking of several objects may provide unique constraints on the escape fraction for faint galaxies and/or the abundance of ultra-faint radiation sources.

We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma-1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq7\%$ and $\simeq10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma-1$ of the IGM temperature-density relation with a precision of $\pm8.6\%$ at $z=3$, $\pm6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.

We study quasar proximity zones in the redshift range $5.77 \leq z \leq 6.54$ by homogeneously analyzing $34$ medium resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift, and argued that this provides evidence for a rapidly evolving intergalactic medium (IGM) neutral fraction during reionization, we measure a much shallower trend $\propto(1+z)^{-1.44}$. We compare our measured proximity zone sizes to predictions from hydrodynamical simulations post-processed with one-dimensional radiative transfer, and find good agreement between observations and theory irrespective of the ionization state of the ambient IGM. This insensitivity to IGM ionization state has been previously noted, and results from the fact that the definition of proximity zone size as the first drop of the smoothed quasar spectrum below the $10\%$ flux transmission level probes locations where the ionizing radiation from the quasar is an order of magnitude larger than the expected ultraviolet ionizing background that sets the neutral fraction of the IGM. Our analysis also uncovered three objects with exceptionally small proximity zones (two have $R_p < 1$proper Mpc), which constitute outliers from the observed distribution and are challenging to explain with our radiative transfer simulations. We consider various explanations for their origin, such as strong absorption line systems associated with the quasar or patchy reionization, but find that the most compelling scenario is that these quasars have been shining for $\lesssim 10^5$yr.

Lyman alpha halos are observed ubiquitously around star-forming galaxies at high redshift, but their origin is still a matter of debate. We demonstrate that the emission from faint unresolved satellite sources, $M_{\rm UV} \gtrsim -17$, clustered around the central galaxies may play a major role in generating spatially extended Ly$\alpha$, continuum (${\rm UV + VIS}$) and H$\alpha$ halos. We apply the analytic formalism developed in Mas-Ribas & Dijkstra (2016) to model the halos around Lyman Alpha Emitters (LAEs) at $z=3.1$, for several different satellite clustering prescriptions. In general, our UV and Ly$\alpha$ surface brightness profiles match the observations well at $20\lesssim r \lesssim 40$ physical kpc from the centers of LAEs. We discuss how our profiles depend on various model assumptions and how these can be tested and constrained with future H$\alpha$ observations by the James Webb Space Telescope (JWST). Our analysis shows how spatially extended halos constrain (i) the presence of otherwise undetectable satellite sources, (ii) the integrated, volumetric production rates of Ly$\alpha$ and LyC photons, and (iii) their population-averaged escape fractions. These quantities are all directly relevant for understanding galaxy formation and evolution and, for high enough redshifts, cosmic reionization.

The distribution of diffuse gas in the intergalactic medium (IGM) imprints a series of hydrogen absorption lines on the spectra of distant background quasars known as the Lyman-$\alpha$ forest. Cosmological hydrodynamical simulations predict that IGM density fluctuations are suppressed below a characteristic scale where thermal pressure balances gravity. We measured this pressure-smoothing scale by quantifying absorption correlations in a sample of close quasar pairs. We compared our measurements to hydrodynamical simulations, where pressure smoothing is determined by the integrated thermal history of the IGM. Our findings are consistent with standard models for photoionization heating by the ultraviolet radiation backgrounds that reionized the universe.

The amplitude of the ionizing background that pervades the intergalactic medium (IGM) at the end of the epoch of reionization provides a valuable constraint on the emissivity of the sources which reionized the Universe. While measurements of the ionizing background at lower redshifts rely on a simulation-calibrated mapping between the photoionization rate and the mean transmission of the Ly$\alpha$ forest, at $z\gtrsim6$ the IGM becomes increasingly opaque, and transmission arises solely in narrow spikes separated by saturated Gunn-Peterson troughs. In this regime, the traditional approach of measuring the average transmission over large $\sim 50$ Mpc$/h$ regions is less sensitive and sub-optimal. Additionally, the five times smaller oscillator strength of the Ly$\beta$ transition implies the Ly$\beta$ forest is considerably more transparent at $z\gtrsim6$, even in the presence of contamination by foreground $z\sim 5$ Ly$\alpha$ forest absorption. In this work we present a novel statistical approach to analyze the joint distribution of transmission spikes in the co-spatial $z\sim 6$ Ly$\alpha$ and Ly$\beta$ forests. Our method relies on Approximate Bayesian Computation (ABC), which circumvents the necessity of computing the intractable likelihood function describing the highly correlated Ly$\alpha$ and Ly$\beta$ transmission. We apply ABC to mock data generated from a large-volume hydrodynamical simulation combined with a state-of-the-art model of ionizing background fluctuations in the post-reionization IGM, and show that it is sensitive to higher IGM neutral hydrogen fractions than previous techniques. As a proof of concept, we apply this methodology to a real spectrum of a $z=6.54$ quasar and measure the ionizing background from $5.4\leq z \leq 6.4$ along this sightline with $\sim0.2$ dex statistical uncertainties.

Despite decades of effort, the timing and duration of He II reionization, as well as its morphology and the properties of the quasars believed to drive it, are still not well constrained. In this paper we present a new method to study both He II reionization and quasars via the thermal proximity effect -- the photoelectric heating of the intergalactic medium around quasars when their hard radiation doubly ionizes helium. We post-process a SPH simulation with 1D radiative calculations, and study how the thermal proximity effect depends on the amount of singly ionized helium, $x_{\rm HeII,0}$, which prevailed in the IGM before the quasar turned on, and the characteristic lifetime $t_{\rm Q}$ for which quasars shine. We find that the amplitude of the temperature boost in the quasar environment depends on $x_{\rm HeII,0}$, with a characteristic value of $\Delta T \simeq 10^4\,{\rm K}$ for an initially singly ionized IGM ($x_{\rm HeII,0} = 1.0$), whereas the size of the thermal proximity zone is sensitive to quasar lifetime $t_{\rm Q}$, with typical sizes of ~100 cMpc for luminous quasars shining for $t_{\rm Q}=10^8$ yr. This temperature boost is manifest as a change in the thermal broadening of H I absorption lines near the quasar. We introduce a new method based on measuring the Ly$\alpha$ forest power spectrum as a function of distance from the quasar, and conduct Bayesian MCMC analysis to demonstrate that the thermal proximity effect should be easily detectable. For a mock dataset of 50 quasars at z~4, we predict that one can measure $x_{\rm HeII,0}$ to a precision $\approx 0.04$, and $t_{\rm Q}$ to a precision of $\approx 0.1$ dex. By applying our formalism to existing high-resolution Ly$\alpha$ forest spectra of quasars at $3.1 \lesssim z \lesssim 5.0$, one should be able to detect the thermal proximity effect, and reconstruct the full reionization history of He II.

The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. We show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier z ~ 15 than they should. This problem arises because at z > 6, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyman-alpha forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. We make our new UVB photoionization and photoheating rates publicly available for use in future simulations.

We present the most precise estimate to date of the clustering of quasars on very small scales, based on a sample of 47 binary quasars with magnitudes of $g<20.85$ and proper transverse separations of $\sim 25\,h^{-1}$\,kpc. Our sample of binary quasars, which is about 6 times larger than any previous spectroscopically confirmed sample on these scales, is targeted using a Kernel Density Estimation technique (KDE) applied to Sloan Digital Sky Survey (SDSS) imaging over most of the SDSS area. Our sample is "complete" in that all of the KDE target pairs with $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc in our area of interest have been spectroscopically confirmed from a combination of previous surveys and our own long-slit observational campaign. We catalogue 230 candidate quasar pairs with angular separations of $<8\arcsec$, from which our binary quasars were identified. We determine the projected correlation function of quasars ($\bar W_{\rm p}$) in four bins of proper transverse scale over the range $17.0 \lesssim R \lesssim 36.2\,h^{-1}$\,kpc. The implied small-scale quasar clustering amplitude from the projected correlation function, integrated across our entire redshift range, is $A=24.1\pm3.6$ at $\sim 26.6 ~h^{-1}$\,kpc. Our sample is the first spectroscopically confirmed sample of quasar pairs that is sufficiently large to study how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ kpc. We find that empirical descriptions of how quasar clustering evolves with redshift at $\sim 25 ~h^{-1}$ Mpc also adequately describe the evolution of quasar clustering at $\sim 25 ~h^{-1}$ kpc.

The HeII transverse proximity effect -- enhanced HeII Ly$\alpha$~transmission in a background sightline caused by the ionizing radiation of a foreground quasar -- offers a unique opportunity to probe the morphology of quasar-driven HeII reionization. We conduct a comprehensive spectroscopic survey to find $z\sim3$ quasars in the foreground of 22 background quasar sightlines with HST/COS HeII Ly$\alpha$~transmission spectra. With our two-tiered survey strategy, consisting of a deep pencil-beam survey and a shallow wide-field survey, we discover 131 new quasars, which we complement with known SDSS/BOSS quasars in our fields. Using a restricted sample of 66 foreground quasars with inferred HeII photoionization rates greater than the expected UV background at these redshifts ($\Gamma_\mathrm{QSO}^\mathrm{HeII} > 5 \times 10^{-16}\,\mathrm{s}^{-1}$) we perform the first statistical analysis of the HeII transverse proximity effect. Our results show qualitative evidence for a large object-to-object variance: among the four foreground quasars with the highest $\Gamma_\mathrm{QSO}^\mathrm{HeII}$ only one (previously known) quasar is associated with a significant HeII transmission spike. We perform a stacking analysis to average down these fluctuations, and detect an excess in the average HeII transmission near the foreground quasars at $3\sigma$ significance. This statistical evidence for the transverse proximity effect is corroborated by a clear dependence of the signal strength on $\Gamma_\mathrm{QSO}^\mathrm{HeII}$. Our detection places a purely geometrical lower limit on the quasar lifetime of $t_\mathrm{Q} > 25\,\mathrm{Myr}$. Improved modeling would additionally constrain quasar obscuration and the mean free path of HeII-ionizing photons.

The CIV broad emission line is visible in optical spectra to redshifts exceeding z~5. CIV has long been known to exhibit significant displacements to the blue and these `blueshifts' almost certainly signal the presence of strong outflows. As a consequence, single-epoch virial black hole (BH) mass estimates derived from CIV velocity-widths are known to be systematically biased compared to masses from the hydrogen Balmer lines. Using a large sample of 230 high-luminosity (log $L_{\rm Bol}$ = 45.5-48 erg/s), redshift 1.5<z<4.0 quasars with both CIV and Balmer line spectra, we have quantified the bias in CIV BH masses as a function of the CIV blueshift. CIV BH masses are shown to be a factor of five larger than the corresponding Balmer-line masses at CIV blueshifts of 3000 km/s and are over-estimated by almost an order of magnitude at the most extreme blueshifts, >5000 km/s. Using the monotonically increasing relationship between the CIV blueshift and the mass ratio BH(CIV)/BH(H$\alpha$) we derive an empirical correction to all CIV BH-masses. The scatter between the corrected CIV masses and the Balmer masses is 0.24 dex at low CIV blueshifts (~0 km/s) and just 0.10 dex at high blueshifts (~3000 km/s), compared to 0.40 dex before the correction. The correction depends only on the CIV line properties - i.e. full-width at half maximum and blueshift - and can therefore be applied to all quasars where CIV emission line properties have been measured, enabling the derivation of un-biased virial BH mass estimates for the majority of high-luminosity, high-redshift, spectroscopically confirmed quasars in the literature.

In the standard picture of structure formation, the first massive galaxies are expected to form at the highest peaks of the density field, which constitute the cores of massive proto-clusters. Luminous quasars (QSOs) at z~4 are the most strongly clustered population known, and should thus reside in massive dark matter halos surrounded by large overdensities of galaxies, implying a strong QSO-galaxy cross-correlation function. We observed six z~4 QSO fields with VLT/FORS exploiting a novel set of narrow band filters custom designed to select Lyman Break Galaxies (LBGs) in a thin redshift slice of Delta_z~0.3, mitigating the projection effects that have limited the sensitivity of previous searches for galaxies around z>~4 QSOs. We find that LBGs are strongly clustered around QSOs, and present the first measurement of the QSO-LBG cross-correlation function at z~4, on scales of 0.1<~R<~9 Mpc/h (comoving). Assuming a power law form for the cross-correlation function xi=(r/r0_QG)^gamma, we measure r0_QG=8.83^{+1.39}_{-1.51} Mpc/h for a fixed slope of gamma=2.0. This result is in agreement with the expected cross-correlation length deduced from measurements of the QSO and LBG auto-correlation function, and assuming a linear bias model. We also measure a strong auto-correlation of LBGs in our QSO fields finding r0_GG=21.59^{+1.72}_{-1.69} Mpc/h for a fixed slope of gamma=1.5, which is ~4 times larger than the LBG auto-correlation length in random fields, providing further evidence that QSOs reside in overdensities of LBGs. Our results qualitatively support a picture where luminous QSOs inhabit exceptionally massive (M_halo>10^12 M_sun) dark matter halos at z~4.

We present the first results from an on-going survey to characterize the circumgalactic medium (CGM) of the massive high-redshift galaxies detected as submillimeter galaxies (SMGs). We constructed a parent sample of 163 SMG-QSO pairs with separations less than $\sim$36 arcsec by cross-matching far-infrared-selected galaxies from Herschel with spectroscopically confirmed QSOs. The Herschel sources were selected to match the properties of SMGs. We determined the sub-arcsecond positions of six Herschel sources with the Very Large Array and obtained secure redshift identification for three of those with near-infrared spectroscopy. The QSO sightlines probe transverse proper distances of 112, 157, and 198 kpc at foreground redshifts of 2.043, 2.515, and 2.184, respectively, which are comparable to the virial radius of the $\sim10^{13}$ Msun halos expected to host SMGs. High-quality absorption-line spectroscopy of the QSOs reveals systematically strong HI Lyman-alpha absorption around all three SMGs, with rest-frame equivalent widths of $\sim2-3$ \AA. However, none of the three absorbers exhibits compelling evidence for optically thick HI gas or metal absorption, in contrast to the dominance of strong neutral absorbers in the CGM of luminous $z \sim 2$ QSOs. The low covering factor of optically thick HI gas around SMGs tentatively indicates that SMGs may not have as prominent cool gas reservoirs in their halos as the co-eval QSOs and that they may inhabit less massive halos than previously thought.

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.

We describe the execution and data reduction of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high-redshift universe with VLT/XSHOOTER" (hereafter `XQ-100'). XQ-100 has produced and made publicly available an homogeneous and high-quality sample of echelle spectra of 100 QSOs at redshifts z~3.5-4.5 observed with full spectral coverage from 315 to 2500 nm at a resolving power ranging from R~4000 to 7000, depending on wavelength. The median signal-to-noise ratios are 33, 25 and 43, as measured at rest-frame wavelengths 1700, 3000 and 3600 Angstrom, respectively. This paper provides future users of XQ-100 data with the basic statistics of the survey, along with details of target selection, data acquisition and data reduction. The paper accompanies the public release of all data products, including 100 reduced spectra. XQ-100 is the largest spectroscopic survey to date of high-redshift QSOs with simultaneous rest-frame UV/optical coverage, and as such enables a wide range of extragalactic research, from cosmology and galaxy evolution to AGN astrophysics.

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 < \rho/\rho_mean < 5x10^5 (\rho_mean is the cosmic mean), while the physical size of individual clouds scales as ~\rho^-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.

MICADO will equip the E-ELT with a first light capability for diffraction limited imaging at near-infrared wavelengths. The instrument's observing modes focus on various flavours of imaging, including astrometric, high contrast, and time resolved. There is also a single object spectroscopic mode optimised for wavelength coverage at moderately high resolution. This contribution provides an overview of the key functionality of the instrument, outlining the scientific rationale for its observing modes. The interface between MICADO and the adaptive optics system MAORY that feeds it is summarised. The design of the instrument is discussed, focussing on the optics and mechanisms inside the cryostat, together with a brief overview of the other key sub-systems.

We present AGNfitter, a publicly available open-source algorithm implementing a fully Bayesian Markov Chain Monte Carlo method to fit the spectral energy distributions (SEDs) of active galactic nuclei (AGN) from the sub-mm to the UV, allowing one to robustly disentangle the physical processes responsible for their emission. AGNfitter makes use of a large library of theoretical, empirical, and semi-empirical models to characterize both the nuclear and host galaxy emission simultaneously. The model consists of four physical emission components: an accretion disk, a torus of AGN heated dust, stellar populations, and cold dust in star forming regions. AGNfitter determines the posterior distributions of numerous parameters that govern the physics of AGN with a fully Bayesian treatment of errors and parameter degeneracies, allowing one to infer integrated luminosities, dust attenuation parameters, stellar masses, and star formation rates. We tested AGNfitter's performace on real data by fitting the SEDs of a sample of 714 X-ray selected AGN from the XMM-COSMOS survey, spectroscopically classified as Type1 (unobscured) and Type2 (obscured) AGN by their optical-UV emission lines. We find that two independent model parameters, namely the reddening of the accretion disk and the column density of the dusty torus, are good proxies for AGN obscuration, allowing us to develop a strategy for classifying AGN as Type1 or Type2, based solely on an SED-fitting analysis. Our classification scheme is in excellent agreement with the spectroscopic classification, giving a completeness fraction of $\sim 86\%$ and $\sim 70\%$, and an efficiency of $\sim 80\%$ and $\sim 77\%$, for Type1 and Type2 AGNs, respectively.

The lifetime of quasars is fundamental for understanding the growth of supermassive black holes, and is an important ingredient in models of the reionization of the intergalactic medium. However, despite various attempts to determine quasar lifetimes, current estimates from a variety of methods are uncertain by orders of magnitude. This work combines cosmological hydrodynamical simulations and 1D radiative transfer to investigate the structure and evolution of the He II Ly$\alpha$ proximity zones around quasars at $z \simeq 3-4$. We show that the time evolution in the proximity zone can be described by a simple analytical model for the approach of the He II fraction $x_{\rm HeII}\left( t \right)$ to ionization equilibrium, and use this picture to illustrate how the transmission profile depends on the quasar lifetime, quasar UV luminosity, and the ionization state of helium in the ambient IGM (i.e. the average He II fraction, or equivalently the metagalactic He II ionizing background). A significant degeneracy exists between the lifetime and the average He II fraction, however the latter can be determined from measurements of the He II Ly$\alpha$ optical depth far from quasars, allowing the lifetime to be measured. We advocate stacking existing He II quasar spectra at $z\sim 3$, and show that the shape of this average proximity zone profile is sensitive to lifetimes as long as $\sim 30$ Myr. At higher redshift $z\sim 4$ where the He II fraction is poorly constrained, degeneracies will make it challenging to determine these parameters independently. Our analytical model for He II proximity zones should also provide a useful description of the properties of H I proximity zones around quasars at $z \simeq 6-7$.

We measure the effective optical depth of HeII Ly\alpha\ absorption \tau$_\mathsf{eff,HeII}$ at 2.3<z<3.5 in 17 UV-transmitting quasars observed with UV spectrographs on the Hubble Space Telescope (HST). The median \tau$_\mathsf{eff,HeII}$ values increase gradually from 1.95 at z=2.7 to 5.17 at z=3.4, but with a strong sightline-to-sightline variance. Many $\simeq$35 comoving Mpc regions of the z>3 intergalactic medium (IGM) remain transmissive (\tau$_\mathsf{eff,HeII}$<4), and the gradual trend with redshift appears consistent with density evolution of a fully reionized IGM. These modest optical depths imply average HeII fractions of x$_\mathsf{HeII}$<0.01 and HeII ionizing photon mean free paths of $\simeq$50 comoving Mpc at z$\simeq$3.4, thus requiring that a substantial volume of the helium in the Universe was already doubly ionized at early times; this stands in conflict with current models of HeII reionization driven by luminous quasars. Along 10 sightlines we measure the coeval HI Ly\alpha\ effective optical depths, allowing us to study the density dependence of \tau$_\mathsf{eff,HeII}$ at z$\sim$3. We establish that the dependence of \tau$_\mathsf{eff,HeII}$ on increasing \tau$_\mathsf{eff,HI}$ is significantly shallower than expected from simple models of an IGM reionized in HeII. This requires higher HeII photoionization rates in overdense regions or underdense regions being not in photoionization equilibrium. Moreover, there are very large fluctuations in \tau$_\mathsf{eff,HeII}$ at all $\tau_\mathsf{eff,HI}$, which greatly exceed the expectations from these simple models. These data present a distinct challenge to scenarios of HeII reionization - an IGM where HeII appears to be predominantly ionized at z$\simeq$3.4, and with a radiation field strength that may be correlated with the density field, but exhibits large fluctuations at all densities.

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$\alpha$ emission from circum- and inter-galactic gas on scales of hundreds of kpc from the central quasar. We do not detect bright giant Ly$\alpha$ 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$\alpha$ 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$\alpha$ 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\alpha}\sim10^{-19}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at $\sim200$ kpc, with a $2.3\sigma$ detection of Ly$\alpha$ emission at SB$_{\rm Ly\alpha}=(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$\alpha$ 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$\alpha$ nebulae around quasars recently discovered by our group.

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.

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.

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.

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 $\lambda_F$, allowing us to rigorously quantify the filtering scale for the first time: we find $\lambda_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 $\gamma$ of the temperature-density relation $T = T_0\Delta^{\gamma-1}$ should be augmented with a third filtering scale parameter $\lambda_F$.

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$\alpha$ forest absorption at $2.2<z<2.5$ with an effective smoothing scale of $\sigma_{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 $\Delta z\sim 0.006$, which does not appear to have a large associated IGM overdensity. These results demonstrate the potential of Ly$\alpha$ 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.

We report the serendipitous detection of a 0.2 L$^*$, Lyman-$\alpha$ 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.

We study how X-rays from stellar binary systems and the hot intracluster medium (ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel implementation of gas cooling in the moving-mesh hydrodynamics code \textsc{arepo}. X-rays from stellar binaries do not affect cooling at all as their emission spectrum is too hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple well with gas in the temperature range $10^4 - 10^6$ K. Idealised simulations show that the hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters because of the high virial temperature ($> 10^7$K), making the interaction between the gas and incident photons very ineffective. Satellite galaxies in cluster environments, on the other hand, experience a high radiation flux due to the emission from the host halo. Low mass satellites ($< 10^{12}\rm{M_\odot}$) in particular have virial temperatures that are exactly in the regime where the effect of the radiation field is maximal. Idealised simulations of satellite galaxies including only the effect of host halo radiation (no ram pressure stripping or tidal effects) fields show a drastic reduction in the amount of cool gas formed ($\sim 40\%$) on a short timescale of about $0.5$ Gyrs. A galaxy merger simulation including all the other environmental quenching mechanisms, shows about $20\%$ reduction in the stellar mass of the satellite and about $\sim 30\%$ reduction in star formation rate after $1$ Gyr due to the host hot halo radiation field. These results indicate that the hot halo radiation fields potentially play an important role in quenching galaxies in cluster environments.

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.

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$\sigma$ 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.

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

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 {\mu}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.

The Ly$\alpha$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $\langle z \rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $\gamma$, and temperature at mean-density, $T_0$, where $T(\Delta) = T_0 \Delta^{\gamma-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $\beta_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$\alpha$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $\gamma=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($\gamma \leq 1$) are disfavored at a significance of over 4$\sigma$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

We present a new method of finding protoclusters using tomographic maps of Ly$\alpha$ 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.

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.

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.

We present the first observations of foreground Lyman-$\alpha$ 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$\alpha$ 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$\alpha$ 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$.

We conduct a deep narrow-band imaging survey of 13 Ly$\alpha$ 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$\alpha$ 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$\sigma$) for the HeII$\lambda$1640 and CIV$\lambda$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$\alpha$ and CIV/Ly$\alpha$ 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$\alpha$ emission is powered by the ionizing luminosity of an obscured AGN. Our models can accommodate HeII/Ly$\alpha$ and CIV/Ly$\alpha$ 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$\sigma$) 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$\alpha$ and CIV/Ly$\alpha$ in a different scenario, where Ly$\alpha$ 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.

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.

We present AGNfitter: a Markov Chain Monte Carlo algorithm developed to fit the spectral energy distributions (SEDs) of active galactic nuclei (AGN) with different physical models of AGN components. This code is well suited to determine in a robust way multiple parameters and their uncertainties, which quantify the physical processes responsible for the panchromatic nature of active galaxies and quasars. We describe the technicalities of the code and test its capabilities in the context of X-ray selected obscured AGN using multiwavelength data from the XMM-COSMOS survey.

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.

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.

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

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{\alpha} 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.

The fraction of AGN luminosity obscured by dust and re-emitted in the mid-IR is critical for understanding AGN evolution, unification, and parsec-scale AGN physics. For unobscured (Type-1) AGN, where we have a direct view of the accretion disk, the dust covering factor can be measured by computing the ratio of re-processed mid-IR emission to intrinsic nuclear bolometric luminosity. We use this technique to estimate the obscured AGN fraction as a function of luminosity and redshift for 513 Type-1 AGN from the XMM-COSMOS survey. The re-processed and intrinsic luminosities are computed by fitting the 18-band COSMOS photometry with a custom SED-fitting code, which jointly models emission from: hot-dust in the AGN torus, the accretion disk, and the host-galaxy. We find a relatively shallow decrease of the luminosity ratio as a function of Lbol, which we interpret as a corresponding decrease in the obscured fraction. In the context of the receding torus model, where dust sublimation reduces the covering factor of more luminous AGN, our measurements require a torus height which increases with luminosity as h ~ Lbol^{0.3-0.4}. Our obscured fraction-luminosity relation agrees with determinations from SDSS censuses of Type-1 and Type-2 quasars, and favors a torus optically thin to mid-IR radiation. We find a much weaker dependence of obscured fraction on 2-10 keV luminosity than previous determinations from X-ray surveys, and argue that X-ray surveys miss a significant population of highly obscured Compton-thick AGN. Our analysis shows no clear evidence for evolution of obscured fraction with redshift.

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.

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.

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.

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.

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]

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.

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.

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 $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 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.

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

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

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.

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.

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.

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

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.

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]

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.

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.

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.

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.

We present observations of two strongly lensed $z\sim5$ Lyman-$\alpha$ 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-$\alpha$ emission at $25.3\pm4.1$\AA~and $135.6\pm20.3$\AA~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-$\alpha$ 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-\alpha}$, L$_{UV})\leq(0.6$L$_{Ly-\alpha}^{*}, 2$L$_{UV}^{*}$) and (L$_{Ly-\alpha}$, L$_{UV})=(0.5$L$_{Ly-\alpha}^{*}, 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}$.

We present tomography of the circum-galactic metal distribution at redshift 1.7 to 4.5 derived from echellete spectroscopy of binary quasars. We find CIV systems at similar redshifts in paired sightlines more often than expected for sightline-independent redshifts. As the separation of the sightlines increases from 36 kpc to 907 kpc, the amplitude of this clustering decreases. At the largest separations, the CIV systems cluster similar to Lyman-break galaxies (Adelberger et al. 2005a). The CIV systems are significantly less correlated than these galaxies, however, at separations less than R_1 ~ 0.42 +/- 0.15 h-1 comoving Mpc. Measured in real space, i.e., transverse to the sightlines, this length scale is significantly smaller than the break scale estimated from the line-of-sight correlation function in redshift space (Scannapieco et al. 2006a). Using a simple model, we interpret the new real-space measurement as an indication of the typical physical size of enriched regions. We adopt this size for enriched regions and fit the redshift-space distortion in the line-of-sight correlation function. The fitted velocity kick is consistent with the peculiar velocity of galaxies as determined by the underlying mass distribution and places an upper limit on the outflow (or inflow) speed of metals. The implied time scale for dispersing metals is larger than the typical stellar ages of Lyman-break galaxies (Shapley et al. 2001), and we argue that enrichment by galaxies at z > 4.3 played a greater role in dispersing metals. To further constrain the growth of enriched regions, we discuss empirical constraints on the evolution of the CIV correlation function with cosmic time. This study demonstrates the potential of tomography for measuring the metal enrichment history of the circum-galactic medium.

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"<\theta<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 \Omega_\Lambda=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.

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

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

We analyze the size evolution of HII regions around 27 quasars between z=5.7 to 6.4 ('quasar near-zones' or NZ). We include more sources than previous studies, and we use more accurate redshifts for the host galaxies, with 8 CO molecular line redshifts and 9 MgII redshifts. We confirm the trend for an increase in NZ size with decreasing redshift, with the luminosity normalized proper size evolving as: R_{NZ,corrected} = (7.4 \pm 0.3) - (8.0 \pm 1.1) \times (z-6) Mpc. While derivation of the absolute neutral fraction remains difficult with this technique, the evolution of the NZ sizes suggests a decrease in the neutral fraction of intergalactic hydrogen by a factor ~ 9.4 from z=6.4 to 5.7, in its simplest interpretation. Alternatively, recent numerical simulations suggest that this rapid increase in near-zone size from z=6.4 to 5.7 is due to the rapid increase in the background photo-ionization rate at the end of the percolation or overlap phase, when the average mean free path of ionizing photons increases dramatically. In either case, the results are consistent with the idea that z ~ 6 to 7 corresponds to the tail end of cosmic reionization. The scatter in the normalized NZ sizes is larger than expected simply from measurement errors, and likely reflects intrinsic differences in the quasars or their environments. We find that the near-zone sizes increase with quasar UV luminosity, as expected for photo-ionization dominated by quasar radiation.

We report the discovery of a multiply lensed Ly Alpha (Lya) emitter at z = 3.90 behind the massive galaxy cluster WARPS J1415.1+3612 at z = 1.026. Images taken by the Hubble Space Telescope(HST) using ACS reveal a complex lensing system that produces a prominent, highly magnified arc and a triplet of smaller arcs grouped tightly around a spectroscopically confirmed cluster member. Spectroscopic observations using FOCAS on Subaru confirm strong Lya emission in the source galaxy and provide redshifts for more than 21 cluster members, from which we obtain a velocity dispersion of 807+/-185 km/s. Assuming a singular isothermal sphere profile, the mass within the Einstein ring (7.13+/-0.38") corresponds to a central velocity dispersion of 686+15-19 km/s for the cluster, consistent with the value estimated from cluster member redshifts. Our mass profile estimate from combining strong lensing and dynamical analyses is in good agreement with both X-ray and weak lensing results.

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.

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.

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.

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.

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\sigma$ 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$.

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.

We report on the discovery of a bright Lyman alpha blob associated with the z=3 quasar SDSSJ124020.91+145535.6 which is also coincident with strong damped Lyman alpha absorption from a foreground galaxy (a so-called proximate damped Lyman alpha system; PDLA). The one dimensional spectrum acquired by the Sloan Digital Sky Survey (SDSS) shows a broad Lyman alpha emission line with a FWHM ~ 500 km/s and a luminosity of L_{Lya} = 3.9e43 erg/s superposed on the trough of the PDLA. Mechanisms for powering this large Lyman alpha luminosity are discussed. We argue against emission from HII regions in the PDLA galaxy since this requires an excessive star-formation rate ~ 500 Msun/yr and would correspond to the largest Lyman alpha luminosity ever measured from a damped Lyman alpha system or starburst galaxy. We use a Monte Carlo radiative transfer simulation to investigate the possibility that the line emission is fluorescent recombination radiation from the PDLA galaxy powered by the ionizing flux of the quasar, but find that the predicted Lyman alpha flux is several orders of magnitude lower than observed. We conclude that the Lyman alpha emission is not associated with the PDLA galaxy at all, but instead is intrinsic to the quasar's host and similar to the extended Lyman alpha 'fuzz' which is detected around many AGN. PDLAs are natural coronagraphs that block their background quasar at Lyman alpha, and we discuss how systems similar to SDSSJ124020.91+145535.6 might be used to image the neutral hydrogen in the PDLA galaxy in silhouette against the screen of extended Lyman alpha emission from the background quasar.

We present the discovery of four gravitationally lensed quasars selected from the spectroscopic quasar catalog of the Sloan Digital Sky Survey. We describe imaging and spectroscopic follow-up observations that support the lensing interpretation of the following four quasars: SDSS J0832+0404 (image separation \theta=1.98", source redshift z_s=1.115, lens redshift z_l=0.659); SDSS J1216+3529 (\theta=1.49", z_s=2.012); SDSS J1322+1052 (\theta=2.00", z_s=1.716); and SDSS J1524+4409 (\theta=1.67", z_s=1.210, z_l=0.320). Each system has two lensed images. We find that the fainter image component of SDSS J0832+0404 is significantly redder than the brighter component, perhaps because of differential reddening by the lensing galaxy. The lens potential of SDSS J1216+3529 might be complicated by the presence of a secondary galaxy near the main lensing galaxy.

We report the first results of our systematic search for strongly lensed quasars using the spectroscopically confirmed quasars in the Sloan Digital Sky Survey (SDSS). Among 46,420 quasars from the SDSS Data Release 3 (~4188 deg^2), we select a subsample of 22,683 quasars that are located at redshifts between 0.6 and 2.2 and are brighter than the Galactic extinction corrected i-band magnitude of 19.1. We identify 220 lens candidates from the quasar subsample, for which we conduct extensive and systematic follow-up observations in optical and near-infrared wavebands, in order to construct a complete lensed quasar sample at image separations between 1'' and 20'' and flux ratios of faint to bright lensed images larger than 10^{-0.5}. We construct a statistical sample of 11 lensed quasars. Ten of these are galaxy-scale lenses with small image separations (~1''-2'') and one is a large separation (15'') system which is produced by a massive cluster of galaxies, representing the first statistical sample of lensed quasars including both galaxy- and cluster-scale lenses. The Data Release 3 spectroscopic quasars contain an additional 11 lensed quasars outside the statistical sample.

We present cosmological results from the statistics of lensed quasars in the Sloan Digital Sky Survey (SDSS) Quasar Lens Search. By taking proper account of the selection function, we compute the expected number of quasars lensed by early-type galaxies and their image separation distribution assuming a flat universe, which is then compared with 7 lenses found in the SDSS Data Release 3 to derive constraints on dark energy under strictly controlled criteria. For a cosmological constant model (w=-1) we obtain \Omega_\Lambda=0.74^{+0.11}_{-0.15}(stat.)^{+0.13}_{-0.06}(syst.). Allowing w to be a free parameter we find \Omega_M=0.26^{+0.07}_{-0.06}(stat.)^{+0.03}_{-0.05}(syst.) and w=-1.1\pm0.6(stat.)^{+0.3}_{-0.5}(syst.) when combined with the constraint from the measurement of baryon acoustic oscillations in the SDSS luminous red galaxy sample. Our results are in good agreement with earlier lensing constraints obtained using radio lenses, and provide additional confirmation of the presence of dark energy consistent with a cosmological constant, derived independently of type Ia supernovae.

We report coincident damped Lyman alpha (DLA) and sub-DLA absorption at z = 2.66 and z = 2.94 towards the z ~ 3 13.8 arcsecond separation binary quasar SDSS 1116+4118 AB. At the redshifts of the absorbers, this angular separation corresponds to a proper transverse separation of ~ 110 kpc. A third absorber, a sub-DLA at z = 2.47, is detected towards SDSS 1116+4118 B, but no corresponding high column density absorber is present towards SDSS 1116+4118 A. We use high resolution galaxy simulations and a clustering analysis to interpret the coincident absorption and its implications for galaxy structure at z ~ 3. We conclude that the common absorption in the two lines of sight is unlikely to arise from a single galaxy, or a galaxy plus satellite system, and is more feasibly explained by a group of two or more galaxies with separations ~ 100 kpc. The impact of these findings on single line of sight observations is also discussed; we show that abundances of DLAs may be affected by up to a few tenths of a dex by line of sight DLA blending. From a Keck ESI spectrum of the two quasars, we measure metal column densities for all five absorbers and determine abundances for the three absorbers with log N(HI) > 20. For the two highest N(HI) absorbers, we determine high levels of metal enrichment, corresponding to 1/3 and 1/5 solar. These metallicities are amongst the highest measured for DLAs at any redshift and are consistent with values measured in Lyman break galaxies at 2 < z < 3. For the DLA at z = 2.94 we also infer an approximately solar ratio of alpha-to-Fe peak elements from [S/Zn] = +0.05, and measure an upper limit for the molecular fraction in this particular line of sight of log f(H_2)< -5.5.

Using the Sloan Digital Sky Survey, Data Release 5, we survey proximate damped Lya systems (PDLAs): absorption line systems with HI column density N(HI)> 2x10^20 cm^{-2} at velocity separation dv < 3000 km/s from their background quasar. These absorbers are physically associated with their background quasars, and their statistics allow us to study quasar environments out to z~5. However, the large ionizing flux emitted by a quasar can ionize the neutral gas in a nearby galaxy possibly giving rise to a ``proximity effect'', analogous to the similar effect observed in the Lya forest. From a sample of 111 PDLAs, we measure the HI frequency distribution f(N), incidence, and gas mass density of the PDLAs near luminous quasars over the redshift interval z= 2.2 to 5. The incidence and mass density of PDLAs at z~3 is approximately twice that of intervening DLAs, but at z<2.5 and z>3.5 the f(N) distribution is consistent with the intervening population. We interpret the observed enhancement of PDLAs around quasars in terms of quasar-galaxy clustering, and compare the strength of the clustering signal to the expectation from independent measures of the respective clustering strengths of DLAs and quasars, as well as a complementary analysis of the clustering of absorbers around quasars in the transverse direction. We find that there are a factor of 5-10 fewer PDLAs around quasars than expected and interpret this result as evidence for the hypothesis that the ionizing flux from the quasars photoevaporates HI in nearby DLA galaxies, thus reducing their cross-section for DLA absorption. This constitues the first detection of a ``proximity effect'' for DLAs.

(Abridged) We study the two-point correlation function of a uniformly selected sample of 4,426 luminous optical quasars with redshift $2.9 \le z\le 5.4$ selected over 4041 deg$^2$ from the Fifth Data Release of the Sloan Digital Sky Survey. For a real-space correlation function of the form $\xi(r)=(r/r_0)^{-\gamma}$, the fitted parameters in comoving coordinates are $r_0 = 15.2 \pm 2.7 h^{-1}$ Mpc and $\gamma = 2.0 \pm 0.3$, over a scale range $4\le r_p\le 150 h^{-1}$ Mpc. Thus high-redshift quasars are appreciably more strongly clustered than their $z \approx 1.5$ counterparts, which have a comoving clustering length $r_0 \approx 6.5 h^{-1}$ Mpc. Dividing our sample into two redshift bins: $2.9\le z\le 3.5$ and $z\ge 3.5$, and assuming a power-law index $\gamma=2.0$, we find a correlation length of $r_0 = 16.9 \pm 1.7 h^{-1}$ Mpc for the former, and $r_0 = 24.3 \pm 2.4 h^{-1}$ Mpc for the latter. Following Martini & Weinberg, we relate the clustering strength and quasar number density to the quasar lifetimes and duty cycle. Using the Sheth & Tormen halo mass function, the quasar lifetime is estimated to lie in the range $4\sim 50$ Myr for quasars with $2.9\le z\le 3.5$; and $30\sim 600$ Myr for quasars with $z\ge 3.5$. The corresponding duty cycles are $0.004\sim 0.05$ for the lower redshift bin and $0.03\sim 0.6$ for the higher redshift bin. The minimum mass of halos in which these quasars reside is $2-3\times 10^{12}\ h^{-1}M_\odot$ for quasars with $2.9\le z\le 3.5$ and $4-6\times 10^{12}\ h^{-1}M_\odot$ for quasars with $z\ge 3.5$.

We report the discovery of a cluster-scale lensed quasar, SDSS J1029+2623, selected from the Sloan Digital Sky Survey. The lens system exhibits two lensed images of a quasar at z_s=2.197. The image separation of 22.5" makes it the largest separation lensed quasar discovered to date. The similarity of the optical spectra and the radio loudnesses of the two components support the lensing hypothesis. Images of the field show a cluster of galaxies at z_l~0.55 that is responsible for the large image separation. The lensed images and the cluster light center are not collinear, which implies that the lensing cluster has a complex structure.

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm.

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.

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 used a sample of 17 Lyman limit systems with column density N_HI > 10^19 cm^-2 selected from 149 projected quasar pair sightlines, to investigate the clustering pattern of optically thick absorbers around luminous quasars at z ~ 2.5. Specifically, we measured the quasar-absorber correlation function in the transverse direction, and found a comoving correlation length of r_0=9.2_{+1.5}_{-1.7} Mpc/h (comoving) assuming a power law correlation function, with gamma=1.6. Applying this transverse clustering strength to the line-of-sight, would predict that ~ 15-50% of all quasars should show a N_HI > 10^19 cm^-2 absorber within a velocity window of v < 3000 km/s. This overpredicts the number of absorbers along the line-of-sight by a large factor, providing compelling evidence that the clustering pattern of optically thick absorbers around quasars is highly anisotropic. The most plausible explanationfor the anisotropy is that the transverse direction is less likely to be illuminated by ionizing photons than the line-of-sight, and that absorbers along the line-of-sight are being photoevaporated. A simple model for the photoevaporation of absorbers subject to the ionizing flux of a quasar is presented, and it is shown that absorbers with volume densities n_H < 0.1 cm^-3 will be photoevaporated if they lie within ~ 1 Mpc (proper) of a luminous quasar. Using this simple model, we illustrate how comparisons of the transverse and line-of-sight clustering around quasars can ultimately be used to constrain the distribution of gas in optically thick absorption line systems.

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.

We present an algorithm for selecting an uniform sample of gravitationally lensed quasar candidates from low-redshift (0.6<z<2.2) quasars brighter than i=19.1 that have been spectroscopically identified in the SDSS. Our algorithm uses morphological and color selections that are intended to identify small- and large-separation lenses, respectively. Our selection algorithm only relies on parameters that the SDSS standard image processing pipeline generates, allowing easy and fast selection of lens candidates. The algorithm has been tested against simulated SDSS images, which adopt distributions of field and quasar parameters taken from the real SDSS data as input. Furthermore, we take differential reddening into account. We find that our selection algorithm is almost complete down to separations of 1'' and flux ratios of 10^-0.5. The algorithm selects both double and quadruple lenses. At a separation of 2'', doubles and quads are selected with similar completeness, and above (below) 2'' the selection of quads is better (worse) than for doubles. Our morphological selection identifies a non-negligible fraction of single quasars: To remove these we fit images of candidates with a model of two point sources and reject those with unusually small image separations and/or large magnitude differences between the two point sources. We estimate the efficiency of our selection algorithm to be at least 8% at image separations smaller than 2'', comparable to that of radio surveys. The efficiency declines as the image separation increases, because of larger contamination from stars. We also present the magnification factor of lensed images as a function of the image separation, which is needed for accurate computation of magnification bias.

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.

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.

We explore the use of the tangential component of weak lensing shear to characterize the ellipticity of clusters of galaxies. We introduce an ellipticity estimator, and quantify its properties for isolated clusters from LCDM N-body simulations. We compare the N-body results to results from smooth analytic models. The expected distribution of the estimator for mock observations is presented, and we show how this distribution is impacted by contaminants such as noise, line of sight projections, and misalignment of the central galaxy used to determine the orientation of the triaxial halo. We examine the radial profile of the estimator and discuss tradeoffs in the observational strategy to determine cluster shape.

We present a detailed investigation into which properties of CDM halos make them effective strong gravitational lenses. Strong lensing cross sections of 878 clusters from an N-body simulation are measured by ray tracing through 13,594 unique projections. We measure concentrations, axis ratios, orientations, and the amount of substructure of each cluster, and compare the lensing weighted distribution of each quantity to that of the cluster population as a whole. The concentrations of lensing clusters are on average 34% larger than the typical cluster in the Universe. Despite this bias, the anomalously high concentrations (c >14) recently measured by several groups, appear to be inconsistent with the concentration distribution in our simulations, which predict < 2% of lensing clusters should have concentrations this high. No correlation is found between lensing cross section and the amount of substructure. We introduce several types of simplified dark matter halos, and use them to isolate which properties of CDM clusters make them effective lenses. Projections of halo substructure onto small radii and the large scale mass distribution of clusters do not significantly influence cross sections. The abundance of giant arcs is primarily determined by the mass distribution within an average overdensity of ~ 10,000. A multiple lens plane ray tracing algorithm is used to show that projections of large scale structure increase the giant arc abundance by a modest amount <7%. We revisit the question of whether there is an excess of giant arcs behind high redshift clusters in the RCS survey and find that the number of high redshift (z > 0.6) lenses is in good agreement with LCDM, although our simulations predict more low redshift (z < 0.6) lenses than were observed. (abridged)

We compute the expected number of quasars multiply imaged by cluster size dark halos for current wide field quasar surveys by carrying out a large ensemble of ray tracing simulations through clusters from a cosmological N-body simulation of the LCDM cosmology. Our calculation predicts ~ 4 quasar lenses with splittings theta > 10" in the SDSS spectroscopic quasar sample, consistent with the recent discovery of the wide separation lens SDSSJ1004+4112 which has theta=14.6". The SDSS faint photometric quasar survey will contain ~12 multiply imaged quasars with splittings theta > 10". Of these, ~ 2 will be lenses with separations theta > 30", and ~ 2 will be at high redshift (z ~ 4).

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.

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.

We study the potential of weak lensing surveys to detect clusters of galaxies, using a fast Particle Mesh cosmological N-body simulation algorithm specifically tailored to investigate the statistics of these mass-selected clusters. In particular, we explore the degree to which the radial positions of galaxy clusters can be determined tomographically, by using photometric redshifts of background source galaxies. We quantify the errors in the tomographic redshifts, and study their dependence on mass, redshift, detection significance, and filtering scheme. A Tomographic Matched Filtering (TMF) scheme, which combines tomography and matched filtering, is introduced, which optimally detects clusters of galaxies in weak lensing surveys. The TMF exploits the extra information provided by photometric redshifts of background source galaxies, neglected in previous studies, to optimally weight the sources. The efficacy and reliability of the TMF is investigated using a large ensemble of mock observations from our simulations and detailed comparisons are made to other filters. Using photometric redshift information with the TMF enhances the number of clusters detected with S/N > 4.5 by as much as 76%, and it increases the dynamic range of weak lensing searches for clusters, detecting more high redshift clusters and extending the mass sensitivity down to the scale of large groups. Furthermore, we find that coarsely binning source galaxies in as few as three redshift bins is sufficient to realize the gains of the TMF. Thus, the tomographic filtering techniques presented here can be applied to current ground based weak lensing data in as few as three bands. Cosmological applications of mass selected cluster samples are also discussed. (abridged)

We report the discovery of two doubly-imaged quasars, SDSS J100128.61+502756.9 and SDSS J120629.65+433217.6, at redshifts of 1.838 and 1.789 and with image separations of 2.86'' and 2.90'', respectively. The objects were selected as lens candidates from the Sloan Digital Sky Survey (SDSS). Based on the identical nature of the spectra of the two quasars in each pair and the identification of the lens galaxies, we conclude that the objects are gravitational lenses. The lenses are complicated; in both systems there are several galaxies in the fields very close to the quasars, in addition to the lens galaxies themselves. The lens modeling implies that these nearby galaxies contribute significantly to the lens potentials. On larger scales, we have detected an enhancement in the galaxy density near SDSS J100128.61+502756.9. The number of lenses with image separation of ~3'' in the SDSS already exceeds the prediction of simple theoretical models based on the standard Lambda-dominated cosmology and observed velocity function of galaxies.

Giant arcs in strong lensing galaxy clusters can provide a purely geometric determination of cosmological parameters, such as the dark energy density and equation of state. We investigate sources of noise in cosmography with giant arcs, focusing in particular on errors induced by density fluctuations along the line-of-sight, and errors caused by modeling uncertainties. We estimate parameter errors in two independent ways, first by developing a Fisher matrix formalism for strong lensing parameters, and next by directly ray-tracing through N-body simulations using a multi-plane lensing code. We show that for reasonable power spectra, density fluctuations from large-scale structure produce > 100% errors in cosmological parameters derived from any single sightline, precluding the use of individual clusters or golden lenses to derive accurate cosmological constraints. Modeling uncertainties similarly can lead to large errors, and we show that the use of parametrized mass models in fitting strong lensing clusters can significantly bias the inferred cosmological parameters. We lastly speculate on means by which these errors may be corrected.

We present the discovery of five new quasars at z>5.7, selected from the multicolor imaging data of the Sloan Digital Sky Survey (SDSS). Three of them, at redshifts 5.93, 6.07, and 6.22, were selected from ~1700 deg^2 of new SDSS Main Survey imaging in the Northern Galactic Cap. An additional quasar, at redshift 5.85, was discovered by coadding the data obtained in the Fall Equatorial Stripe in the SDSS Southern Survey Region. The fifth object, at redshift 5.80, is selected from a non-standard SDSS scan in the Southern Galactic Cap outside the Main Survey area. The spectrum of SDSS J162331.81+311200.5 (z=6.22) shows a complete Gunn-Peterson trough at z_abs > 5.95, similar to the troughs detected in other three z>6.2 quasars known. We present a composite spectrum of the z>5.7 quasars discovered in the SDSS to date. The average emission line and continuum properties of z~6 quasars exhibit no significant evolution compared to those at low redshift. Using a complete sample of nine z>5.7 quasars, we find that the density of quasars with M_1450 < -26.7 at z~6 is (6+/-2) x 10^-10 per Mpc^3 consistent with our previous estimates. The luminosity distribution of the sample is fit with a power law luminosity function Psi(L) ~ L^(-3.2+/-0.7), somewhat steeper than but consistent with our previous estimates.

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.

We report on preliminary results of our X-ray survey of the most massive clusters currently identified from the Deep Lens Survey (DLS). The DLS cluster sample is selected based on weak lensing shear, which makes it possible for the first time to study clusters in a baryon-independent way. In this article we present X-ray properties of a subset of the shear-selected cluster sample.

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

We study the expected properties and statistics of giant arcs produced by galaxy clusters in a LambdaCDM universe and investigate how the characteristics of CDM clusters determine the properties of the arcs they generate. Due to the triaxiality and substructure of CDM halos, the giant arc cross section for individual clusters varies by more than an order of magnitude as a function of viewing angle. In addition, the shallow density cusps and triaxiality of CDM clusters cause systematic alignments of giant arcs which should be testable with larger samples from forthcoming lensing surveys. We compute the predicted statistics of giant arcs for the LambdaCDM model and compare to results from previous surveys. The predicted arc statistics are in excellent agreement with the numbers of giant arcs observed around low redshift (0.2 < z < 0.6) clusters from the EMSS sample, however there are hints of a possible excess of arcs observed around high redshift z > 0.6 clusters. This excess, if real, appears to be due to the presence of highly massive or concentrated clusters at high redshifts.

Secondary anisotropies of the CMB have the potential to reveal intricate details about the history of our universe between the present and recombination epochs. However, because the CMB we observe is the projected sum of a multitude of effects, the interpretation of small scale anisotropies by future high resolution experiments will be marred by uncertainty and speculation without the handles provided by other observations. In this paper we show that cross correlating the CMB with an overlapping weak lensing survey will isolate the elusive kinetic Sunyaev-Zel'dovich Effect from secondary anisotropies generated at higher redshifts. We show that if upcoming high angular resolution CMB experiments, like PLANCK/ACT/SPT, cover the same area of sky as current and future weak lensing surveys, like CFTHLS/SNAP/LSST, the cross correlation of cosmic shear with the kSZ effect will be detected with high signal to noise ratio, increasing the potential science accessible to both sets of surveys. For example, if ACT and a CFHTLS like survey were to overlap this cross-correlation would be detected with a total signal to noise ratio greater than 220, reaching 1.8 per individual multipole around l \sim 5000. Furthermore, this cross-correlation probes the three point coupling between the underlying dark matter and the "momentum" of the ionized baryons in the densest regions of the universe at intermediate redshifts. Similar to the tSZ power spectrum, its strength is extremely sensitive to the power spectrum normalization parameter, \sigma_8, scaling roughly as \sigma_8^7. It provides an effective mechanism to isolate any component of anisotropy due to patchy reionization and rule out primordial small scale anisotropy.

The Deep Lens Survey (DLS) is a deep BVRz' imaging survey of seven 2x2 degree fields, with all data to be made public. The primary scientific driver is weak gravitational lensing, but the survey is also designed to enable a wide array of other astrophysical investigations. A unique feature of this survey is the search for transient phenomena. We subtract multiple exposures of a field, detect differences, classify, and release transients on the Web within about an hour of observation. Here we summarize the scientific goals of the DLS, field and filter selection, observing techniques and current status, data reduction, data products and release, and transient detections. Finally, we discuss some lessons which might apply to future large surveys such as LSST.

The number of mass clusters and their distribution in redshift are very sensitive to the density of matter Omega_m and the equation of state of dark energy w. Using weak lens gravitational tomography one can detect clusters of dark matter, weigh them, image their projected mass distribution, and determine their 3-D location. The degeneracy curve in the Omega_m - w plane is nearly orthogonal to that from CMB or SN measurements. Thus, a combination of CMB data with weak lens tomography of clusters can yield precision measurements of Omega_m and w, independently of the SN observations. The Large Synoptic Survey Telescope (LSST) will repeatedly survey 30,000 square degrees of the sky in multiple wavelengths. LSST will create a 3-D tomographic assay of mass overdensities back to half the age of the universe by measuring the shear and color-redshift of billions of high redshift galaxies. By simultaneously measuring several functions of cosmic shear and mass cluster abundance, LSST will provide a number of independent constraints on the dark energy density and the equation of state.

We consider the consequences of SIDM for a velocity dependent cross section per unit mass. Accretion of SIDM onto seed black holes can produce supermassive black holes that are too large for certain combinations of parameters,which is used to obtain a new constraint on the dark matter interaction. Constraints due to other considerations are presented and previous ones are generalized. The black hole constraint is extremely sensitive to the slope \alpha, of the inner density profile of dark halos. For the most probable value of \alpha=1.3, there exists a narrow range in parameter space, such that all constraints are satisfied. However, the adiabatic compression of the dark halo by baryons as they cool and contract in normal galaxies yields a steeper cusp, \alpha=1.7. This gives a tighter constraint, which would exclude SIDM as a possible solution to the purported problems with CDM in the absence of other dynamical processes. Nevertheless, SIDM with parameters consistent with this stronger constraint, can explain the ubiquity of supermassive black holes in the centers of galaxies. A ``best fit'' model is presented which reproduces the supermassive black hole masses and their observed correlations with the velocity dispersion of the host bulges. Specifically, the fourth power dependence of black hole mass on velocity dispersion is a direct consequence of the power spectrum having an index of n=-2. Although the dark matter collision rates for this model are too small to directly remedy problems with CDM, mergers between dark halos harboring supermassive black holes at high redshift could ameliorate the cuspy halo problem. This scenario also explains the lack of comparable supermassive black holes in bulgeless galaxies like M33.

The gravitational-wave and accretion driven evolution of neutron stars in low mass X-ray binaries and similar systems is analyzed, while the amplitude of the radiating perturbation (here assumed to be an r-mode) remains small. If most of the star is superfluid, with (temperature independent) mutual friction dominating the ordinary (temperature dependent) shear viscosity, the amplitude of the mode and the angular velocity of the star oscillate about their equilibrium values with a period of at least a few hundred years. The resulting oscillation of the neutron star temperature is also computed. For temperature dependent viscosity, the general conditions for the equilibrium to be stable are found.