We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of eight compact star-forming galaxies at redshifts z=0.02811-0.06540, with low oxygen abundances 12+log(O/H)=7.43-7.82 and extremely high emission-line flux ratios O32=[OIII]5007/[OII]3727~22-39, aiming to study the properties of Ly-alpha emission in such conditions. We find a diversity in Ly-alpha properties. In five galaxies Ly-alpha emission line is strong, with equivalent width (EW) in the range 45-190A. In the remaining galaxies, weak Ly-alpha emission with EW(Ly-alpha)~2-7A is superposed on a broad Ly-alpha absorption line, indicating a high neutral hydrogen column density N(HI)~(1-3)x10^21 cm^-2. We examine the relation between the Ly-alpha escape fraction fesc(Ly-alpha) and the Lyman continuum escape fraction fesc(LyC), using direct measures of the latter in eleven low-redshift LyC leakers, to verify whether fesc(Ly-alpha) can be an indirect measure of escaping LyC radiation. The usefulness of O32, of the Ly-alpha equivalent width EW(Ly-alpha) and of the Ly-alpha peak separation Vsep as indirect indicators of Ly-alpha leakage is also discussed. It is shown that there is no correlation between O32 and fesc(Ly-alpha). We find an increase of fesc(Ly-alpha) with increasing EW(Ly-alpha) for EW(Ly-alpha)<100A, but for higher EW(Ly-alpha)>150A the fesc(Ly-alpha) is nearly constant attaining the value of ~0.25. We find an anticorrelation between fesc(Ly-alpha) and Vsep, though not as tight as the one found earlier between fesc(LyC) and Vsep. This finding makes Vsep a promising indirect indicator of both the Ly-alpha and ionizing radiation leakage.

Sub-damped Lyman alpha systems (subDLAs; HI column densities of 19.0<=logN(HI)<20.3) are rarely included in the cosmic HI census performed at redshifts z>=1.5, yet are expected to contribute significantly to the overall HI mass budget of the Universe. In this paper, we present a blindly selected sample of 155 subDLAs found along 100 quasar sightlines (with a redshift path length X=475) in the XQ-100 survey to investigate the contribution of subDLAs to the HI mass density of the Universe. The impact of X-Shooter's spectral resolution on sub-DLA identification is evaluated, and found to be sufficient for reliably finding absorbers with logN(HI)>=18.9. We compared the implications of searching for subDLAs solely using HI absorption versus the use of additional metal lines to confirm the identification, and found that metal-selection techniques would have missed 75 subDLAs. Using a bootstrap-Monte Carlo simulation, we computed the column density distribution function (f(N,X)) and the cosmological HI mass density of subDLAs and compared with our previous work based on the XQ-100 damped Lyman alpha systems. We do not find any significant redshift evolution in f(N,X) or cosmological HI mass density for subDLAs. However, subDLAs contribute 10-20 per cent of the total cosmological HI mass density measured at redshifts 2<z<5 (agreeing with previous measurements), and thus have a small but significant contribution to the HI budget of the Universe.

We present a survey for metal absorption systems traced by neutral oxygen over $3.2 < z < 6.5$. Our survey uses Keck/ESI and VLT/X-Shooter spectra of 199 QSOs with redshifts up to 6.6. In total we detect 74 O I absorbers, of which 57 are separated from the background QSO by more than 5000 km/s. We use a maximum likelihood approach to fit the distribution of O I $\lambda$1302 equivalent widths in bins of redshift, and from this determine the evolution in number density of absorbers with $W_{1302} > 0.05$ angstroms. We find that the number density does not monotonically increase with decreasing redshift, as would naively be expected from the buildup of metal-enriched circumgalactic gas with time. The number density over $4.9 < z < 5.7$ is a factor of 1.7--4.1 lower (68% confidence) than over $5.7 < z < 6.5$, with a lower value at $z < 5.7$ favored with 99% confidence. This decrease suggests that the fraction of metals in a low-ionization phase is larger at $z \sim 6$ than at lower redshifts. Absorption from highly ionized metals traced by C IV is also weaker in higher-redshift O I systems, supporting this picture. The evolution of O I absorbers implies that metal-enriched circumgalactic gas at $z \sim 6$ is undergoing an ionization transition driven by a strengthening ultraviolet background. This in turn suggests that the reionization of the diffuse intergalactic medium may still be ongoing at or only recently ended by this epoch.

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.

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.

X-ray and mid-infrared emission are signposts of the accretion of matter onto the supermassive black holes that reside at the centres of most galaxies. As a major step towards understanding accreting supermassive black holes and their role in the evolution of galaxies, we will use the 4MOST multi-object spectrograph to provide a highly complete census of active galactic nuclei over a large fraction of the extragalactic sky observed in X-rays by eROSITA that is visible to 4MOST. We will systematically follow up all eROSITA point-like extragalactic X-ray sources (mostly active galactic nuclei), and complement them with a heavily obscured active galactic nuclei selection approach using mid-infrared data from the Wide-field Infrared Survey Explorer (WISE). The X-ray and mid-infrared flux limits of eROSITA and WISE are well matched to the spectroscopic capabilities of a 4-metre-class telescope, allowing us to reach completeness levels of ~80-90% for all X-ray selected active galactic nuclei with fluxes $f_{0.5-2\ {\rm keV}} > 10^{-14}$ erg s$^{-1}$ cm$^{-2}$; this is about a factor of 30 deeper than the ROSAT all-sky survey. With these data we will determine the physical properties (redshift, luminosity, line emission strength, masses, etc.) of up to one million supermassive black holes, constrain their cosmic evolution and clustering properties, and explore the connection between active galactic nuclei and large-scale structure over redshifts $0 \le z \le 6$.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

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

The 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 investigate the NV absorption signal along the line of sight of background quasars, in order to test the robustness of the use of this ion as criterion to select intrinsic (i.e. physically related to the quasar host galaxy) narrow absorption lines (NALs). We build composite spectra from a sample of $\sim$ 1000 CIV absorbers, covering the redshift range 2.55 < z < 4.73, identified in 100 individual sight lines from the XQ-100 Legacy Survey. We detect a statistical significant NV absorption signal only within 5000 km s$^{-1}$ of the systemic redshift, z$\rm_{em}$. This absorption trough is $\sim$ 15$\sigma$ when only CIV systems with N(CIV) > 10$^{14}$ cm$^{-2}$ are included in the composite spectrum. This result confirms that NV offers an excellent statistical tool to identify intrinsic systems. We exploit the stacks of 11 different ions to show that the gas in proximity to a quasar exhibits a considerably different ionization state with respect to gas in the transverse direction and intervening gas at large velocity separations from the continuum source. Indeed, we find a dearth of cool gas, as traced by low-ionization species and in particular by MgII, in the proximity of the quasar. We compare our findings with the predictions given by a range of Cloudy ionization models and find that they can be naturally explained by ionization effects of the quasar.

We have obtained the first complete ultraviolet (UV) spectrum of a strong Lyman continuum(LyC) emitter at low redshift -- the compact, low-metallicity, star-forming galaxy J1154+2443 -- with a Lyman continuum escape fraction of 46% discovered recently. The Space Telescope Imaging Spectrograph spectrum shows strong Lya and CIII] 1909 emission, as well as OIII] 1666. Our observations show that strong LyC emitters can have UV emission lines with a high equivalent width (e.g. EW(CIII])$=11.7 \pm 2.9 \AA$ rest-frame), although their equivalent widths should be reduced due to the loss of ionizing photons. The intrinsic ionizing photon production efficiency of J1154+2443 is high, $\log(\xi_{\rm ion}^0)=25.56$ erg$^{-1}$ Hz, comparable to that of other recently discovered $z \sim 0.3-0.4$ LyC emitters. Combining our measurements and earlier determinations from the literature, we find a trend of increasing $\xi_{\rm ion}^0$ with increasing CIII] 1909 equivalent width, which can be understood by a combination of decreasing stellar population age and metallicity. Simple ionization and density-bounded photoionization models can explain the main observational features including the UV spectrum of J1154+2443.

We present new measurements of the free-streaming of warm dark matter (WDM) from Lyman-$\alpha$ flux-power spectra. We use data from the medium resolution, intermediate redshift XQ-100 sample observed with the X-shooter spectrograph ($z=3 - 4.2$) and the high-resolution, high-redshift sample used in Viel et al. (2013) obtained with the HIRES/MIKE spectrographs ($z=4.2 - 5.4$). Based on further improved modelling of the dependence of the Lyman-$\alpha$ flux-power spectrum on the free-streaming of dark matter, cosmological parameters, as well as the thermal history of the intergalactic medium (IGM) with hydrodynamical simulations, we obtain the following limits, expressed as the equivalent mass of thermal relic WDM particles. The XQ-100 flux power spectrum alone gives a lower limit of 1.4 keV, the re-analysis of the HIRES/MIKE sample gives 4.1 keV while the combined analysis gives our best and significantly strengthened lower limit of 5.3 keV (all 2$\sigma$ C.L.). The further improvement in the joint analysis is partly due to the fact that the two data sets have different degeneracies between astrophysical and cosmological parameters that are broken when the data sets are combined, and more importantly on chosen priors on the thermal evolution. These results all assume that the temperature evolution of the IGM can be modelled as a power law in redshift. Allowing for a non-smooth evolution of the temperature of the IGM with sudden temperature changes of up to 5000K reduces the lower limit for the combined analysis to 3.5 keV. A WDM with smaller thermal relic masses would require, however, a sudden temperature jump of $5000\,K$ or more in the narrow redshift interval $z=4.6-4.8$, in disagreement with observations of the thermal history based on high-resolution resolution Lyman-$\alpha$ forest data and expectations for photo-heating and cooling in the low density IGM at these redshifts.

We present observations with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope of five star-forming galaxies at redshifts z in the range 0.2993-0.4317 and with high emission-line flux ratios O32=[OIII]5007/[OII]3727 ~ 8-27 aiming to detect the Lyman continuum (LyC) emission. We detect LyC emission in all galaxies with the escape fractions fesc(LyC) in a range of 2-72 per cent. A narrow Ly-alpha emission line with two peaks in four galaxies and with three peaks in one object is seen in medium-resolution COS spectra with a velocity separation between the peaks Vsep varying from ~153 km/s to ~345 km/s. We find a general increase of the LyC escape fraction with increasing O32 and decreasing stellar mass M*, but with a large scatter of fesc(LyC). A tight anti-correlation is found between fesc(LyC) and Vsep making Vsep a good parameter for the indirect determination of the LyC escape fraction. We argue that one possible source driving the escape of ionizing radiation is stellar winds and radiation from hot massive stars.

We recently discovered that the luminous radio-quiet QSO LBQS0302-0019 at z=3.286 is likely accompanied by an obscured AGN at 20kpc projected distance, which we dubbed Jil. It represents the tightest candidate obscured/unobscured dual AGN system at z>3. To verify the dual AGN scenario we obtained deep $K_s$ band (rest-frame $V$ band) imaging with the VLT/HAWK-I+GRAAL instrument at 0.4" resolution during science verification in January 2018. Indeed, we detect the individual host galaxies of the QSO and Jil with estimated stellar masses of $\log(M_\star/M_{\odot})=11.4\pm0.5$ and $\log(M_\star/M_{\odot})=10.9\pm0.5$, respectively. Near-IR spectra obtained with VLT-KMOS reveal a clear [OIII] $\lambda$5007 line detection at the location of Jil which does not contribute significantly to the $K_s$ band flux. Both observations therefore corroborate the dual AGN scenario. A comparison to Illustris cosmological simulations suggests a parent halo mass of $\log(M_\mathrm{halo}/M_{\odot})=13.2\pm0.5$ for this interacting galaxy system, corresponding to a very massive dark matter halo at that epoch.

A prediction of the current paradigm of the hierarchical assembly of galaxies is the presence of supermassive dual black holes at separations of a few kpc or less. In this context, we report the detection of a narrow-line emitter within the extended LyA nebula (~120kpc diameter) of the luminous radio-quiet quasi-stellar object (QSO) LBQS0302-0019 at z=3.286. We identify several high-ionization narrow emission lines (HeII, CIV, CIII]) associated with this point-like source, which we have named "Jil", which is only ~20kpc (2.9") away from the QSO in projection. Emission-line diagnostics confirm that the source is likely powered by photoionization of an obscured active galactic nucleus (AGN) three orders of magnitude fainter than the QSO. The system represents the tightest unobscured/obscured dual AGN currently known at z>3, highlighting the power of MUSE to detect these elusive systems.

We report the detection of the Lyman continuum (LyC) radiation of the compact star-forming galaxy (SFG) J1154+2443 observed with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope. This galaxy, at a redshift of z=0.3690, is characterized by a high emission-line flux ratio O32=[OIII]5007/[OII]3727=11.5. The escape fraction of the LyC radiation fesc(LyC) in this galaxy is 46 per cent, the highest value found so far in low-redshift SFGs and one of the highest values found in galaxies at any redshift. The narrow double-peaked Lya emission line is detected in the spectrum of J1154+2443 with a separation between the peaks Vsep of 199 km/s, one of the lowest known for Lya-emitting galaxies, implying a high fesc(Lya). Comparing the extinction-corrected Lya/Hb flux ratio with the case B value we find fesc(Lya) = 98 per cent. Our observations, combined with previous detections in the literature, reveal an increase of O32 with increasing fesc(LyC). We also find a tight anticorrelation between fesc(LyC) and Vsep. The surface brightness profile derived from the COS acquisition image reveals a bright star-forming region in the centre and an exponential disc in the outskirts with a disc scale length alpha=1.09 kpc. J1154+2443, compared to other known low-redshift LyC leakers, is characterized by the lowest metallicity, 12+logO/H=7.65+/-0.01, the lowest stellar mass M*=10^8.20 Msun, a similar star formation rate SFR=18.9 Msun/yr and a high specific SFR of 1.2x10^-7 yr^-1.

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\%)$.

Previous studies have shown that the incidence rate of intervening strong MgII absorbers towards GRBs were a factor of 2 - 4 higher than towards quasars. Exploring the similar sized and uniformly selected legacy data sets XQ-100 and XSGRB, each consisting of 100 quasar and 81 GRB afterglow spectra obtained with a single instrument (VLT/X-shooter), we demonstrate that there is no disagreement in the number density of strong MgII absorbers with rest-frame equivalent widths $W_r^{2796} >$ 1 {\AA} towards GRBs and quasars in the redshift range 0.1 < z < 5. With large and similar sample sizes, and path length coverages of $\Delta$z = 57.8 and 254.4 for GRBs and quasars, respectively, the incidences of intervening absorbers are consistent within 1 sigma uncertainty levels at all redshifts. For absorbers at z < 2.3 the incidence towards GRBs is a factor of 1.5$\pm$0.4 higher than the expected number of strong MgII absorbers in SDSS quasar spectra, while for quasar absorbers observed with X-shooter we find an excess factor of 1.4$\pm$0.2 relative to SDSS quasars. Conversely, the incidence rates agree at all redshifts with reported high spectral resolution quasar data, and no excess is found. The only remaining discrepancy in incidences is between SDSS MgII catalogues and high spectral resolution studies. The rest-frame equivalent width distribution also agrees to within 1 sigma uncertainty levels between the GRB and quasar samples. Intervening strong MgII absorbers towards GRBs are therefore neither unusually frequent, nor unusually strong.

To reionize the early universe, high-energy photons must escape the galaxies that produce them. It has been suggested that stellar feedback drives galactic outflows out of star-forming regions, creating low density channels through which ionizing photons escape into the inter-galactic medium. We compare the galactic outflow properties of confirmed Lyman continuum (LyC) leaking galaxies to a control sample of nearby star-forming galaxies to explore whether the outflows from leakers are extreme as compared to the control sample. We use data from the Cosmic Origins Spectrograph on the Hubble Space Telescope to measure the equivalent widths and velocities of Si II and Si III absorption lines, tracing neutral and ionized galactic outflows. We find that the Si II and Si III equivalent widths of the LyC leakers reside on the low-end of the trend established by the control sample. The leakers' velocities are not statistically different than the control sample, but their absorption line profiles have a different asymmetry: their central velocities are closer to their maximum velocities. The outflow kinematics and equivalent widths are consistent with the scaling relations between outflow properties and host galaxy properties -- most notably metallicity -- defined by the control sample. Additionally, we use the Ly\alpha\ profiles to show that the Si II equivalent width scales with the Ly\alpha\ peak velocity separation. We determine that the low equivalent widths of the leakers are likely driven by low metallicities and low H I column densities, consistent with a density-bounded ionization region, although we cannot rule out significant variations in covering fraction. While we do not find that the LyC leakers have extreme outflow velocities, the low maximum-to-central velocity ratios demonstrate the importance of the acceleration and density profiles for LyC and Ly\alpha\ escape. [abridged]

We analyze new far-ultraviolet spectra of 13 quasars from the z~0.2 COS-Halos survey that cover the HI Lyman limit of 14 circumgalactic medium (CGM) systems. These data yield precise estimates or more constraining limits than previous COS-Halos measurements on the HI column densities NHI. We then apply a Monte-Carlo Markov Chain approach on 32 systems from COS-Halos to estimate the metallicity of the cool (T~10^4K) CGM gas that gives rise to low-ionization state metal lines, under the assumption of photoionization equilibrium with the extragalactic UV background. The principle results are: (1) the CGM of field L* galaxies exhibits a declining HI surface density with impact parameter Rperp (at >99.5%$ confidence), (2) the transmission of ionizing radiation through CGM gas alone is 70+/-7%; (3) the metallicity distribution function of the cool CGM is unimodal with a median of 1/3 Z_Sun and a 95% interval from ~1/50 Z_Sun to over 3x solar. The incidence of metal poor (<1/100 Z_Sun) gas is low, implying any such gas discovered along quasar sightlines is typically unrelated to L* galaxies; (4) we find an unexpected increase in gas metallicity with declining NHI (at >99.9% confidence) and, therefore, also with increasing Rperp. The high metallicity at large radii implies early enrichment; (5) A non-parametric estimate of the cool CGM gas mass is M_CGM_cool = 9.2 +/- 4.3 10^10 Msun, which together with new mass estimates for the hot CGM may resolve the galactic missing baryons problem. Future analyses of halo gas should focus on the underlying astrophysics governing the CGM, rather than processes that simply expel the medium from the halo.

We present the Lyman-$\alpha$ flux power spectrum measurements of the XQ-100 sample of quasar spectra obtained in the context of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high redshift universe with VLT/XSHOOTER". Using $100$ quasar spectra with medium resolution and signal-to-noise ratio we measure the power spectrum over a range of redshifts $z = 3 - 4.2$ and over a range of scales $k = 0.003 - 0.06\,\mathrm{s\,km^{-1}}$. The results agree well with the measurements of the one-dimensional power spectrum found in the literature. The data analysis used in this paper is based on the Fourier transform and has been tested on synthetic data. Systematic and statistical uncertainties of our measurements are estimated, with a total error (statistical and systematic) comparable to the one of the BOSS data in the overlapping range of scales, and smaller by more than $50\%$ for higher redshift bins ($z>3.6$) and small scales ($k > 0.01\,\mathrm{s\,km^{-1}}$). The XQ-100 data set has the unique feature of having signal-to-noise ratios and resolution intermediate between the two data sets that are typically used to perform cosmological studies, i.e. BOSS and high-resolution spectra (e.g. UVES/VLT or HIRES). More importantly, the measured flux power spectra span the high redshift regime which is usually more constraining for structure formation models.

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 XQ-100 survey provides optical and near infrared coverage of 36 blindly selected, intervening damped Lyman alpha systems (DLAs) at 2 < z < 4, simultaneously covering the MgII doublet at 2796A, 2803A, and the Ly-alpha transition. Using the XQ-100 DLA sample, we investigate the completeness of selecting DLA absorbers based on their MgII rest-frame equivalent width (W2796) at these redshifts. Of the 29 DLAs with clean MgII profiles, we find that six (20% of DLAs) have W2796 < 0.6A. The DLA incidence rate of W2796 < 0.6A absorbers is a factor of ~5 higher than what is seen in z~1 samples, indicating a potential evolution in the MgII properties of DLAs with redshift. All of the W2796 < 0.6A DLAs have low metallicities (-2.5 < [M/H] < -1.7), small velocity widths (v90 < 50 km/s), and tend to have relatively low N(HI). We demonstrate that the exclusion of these low W2796 DLAs results in a higher mean N(HI) which in turn leads to a ~7% increase in the cosmological gas density of HI of DLAs at 2 < z < 4; and that this exclusion has a minimal effect on the HI-weighted mean metallicity.

The XQ-100 survey has provided high signal-noise spectra of 100 redshift 3-4.5 quasars with the X-Shooter spectrograph. The metal abundances for 13 elements in the 41 damped Lyman alpha systems (DLAs) identified in the XQ-100 sample are presented, and an investigation into abundances of a variety of DLA classes is conducted. The XQ-100 DLA sample contains five DLAs within 5000 km/s of their host quasar (proximate DLAs; PDLAs) as well as three sightlines which contain two DLAs within 10,000 km/s of each other along the same line-of-sight (multiple DLAs; MDLAs). Combined with previous observations in the literature, we demonstrate that PDLAs with logN(HI)<21.0 show lower [S/H] and [Fe/H] (relative to intervening systems with similar redshift and N(HI)), whilst higher [S/H] and [Si/H] are seen in PDLAs with logN(HI)>21.0. These abundance discrepancies are independent of their line-of-sight velocity separation from the host quasar, and the velocity width of the metal lines (v90). Contrary to previous studies, MDLAs show no difference in [alpha/Fe] relative to single DLAs matched in metallicity and redshift. In addition, we present follow-up UVES data of J0034+1639, a sightline containing three DLAs, including a metal-poor DLA with [Fe/H]=-2.82 (the third lowest [Fe/H] in DLAs identified to date) at z=4.25. Lastly we study the dust-corrected [Zn/Fe], emphasizing that near-IR coverage of X-Shooter provides unprecedented access to MgII, CaII and TiII lines (at redshifts 3-4) to provide additional evidence for subsolar [Zn/Fe] ratio in DLAs.

We have recently reported the discovery of five low redshift Lyman continuum (LyC) emitters (LCEs, hereafter) with absolute escape fractions fesc(LyC) ranging from 6 to 13%, higher than previously found, and which more than doubles the number of low redshift LCEs.We use these observations to test theoretical predictions about a link between the characteristics of the Lyman-alpha (Lya) line from galaxies and the escape of ionising photons. We analyse the Lya spectra of eight LCEs of the local Universe observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope (our five leakers and three galaxies from the litterature), and compare their strengths and shapes to the theoretical criteria and comparison samples of local galaxies: the Lyman Alpha Reference Survey, Lyman Break Analogs, Green Peas, and the high-redshift strong LyC leaker Ion2. Our LCEs are found to be strong Lya emitters, with high equivalent widths, EW(Lya)> 70 {\AA}, and large Lya escape fractions, fesc(Lya) > 20%. The Lya profiles are all double-peaked with a small peak separation, in agreement with our theoretical expectations. They also have no underlying absorption at the Lya position. All these characteristics are very different from the Lya properties of typical star-forming galaxies of the local Universe. A subset of the comparison samples (2-3 Green Pea galaxies) share these extreme values, indicating that they could also be leaking. We also find a strong correlation between the star formation rate surface density and the escape fraction of ionising photons, indicating that the compactness of star-forming regions plays a role in shaping low column density paths in the interstellar medium of LCEs. The Lya properties of LCEs are peculiar: Lya can be used as a reliable tracer of LyC escape from galaxies, in complement to other indirect diagnostics proposed in the literature.

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.

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

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.

We have recently discovered five Lyman continuum leaking galaxies at z~0.3, selected for their compactness, intense star-formation, and high [OIII]/[OII] ratio (Izotov et al. 2016ab). Here we derive their ionizing photon production efficiency, a fundamental quantity for inferring the number of photons available to reionize the Universe, for the first time for galaxies with confirmed strong Lyman continuum escape (fesc~6-13%). We find an ionizing photon production per unit UV luminosity, which is a factor 2-6 times higher than the canonical value when reported to their observed UV luminosity. After correction for extinction this value is close to the canonical value. The properties of our five Lyman continuum leakers are found to be very similar to those of the confirmed z=3.218 leaker Ion2 from de Barros et al. (2016) and very similar to those of typical star-forming galaxies at z>~6. Our results suggest that UV bright galaxies at high-z such as Lyman break galaxies can be Lyman continuum leakers and that their contribution to cosmic reionization may be underestimated.

Following our first detection reported in Izotov et al. (2016), we present the detection of Lyman continuum (LyC) radiation of four other compact star-forming galaxies observed with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope (HST). These galaxies, at redshifts of z~0.3, are characterized by high emission-line flux ratios [OIII]5007/[OII]3727 > 5. The escape fractions of the LyC radiation fesc(LyC) in these galaxies are in the range of ~6%-13%, the highest values found so far in low-redshift star-forming galaxies. Narrow double-peaked Lyalpha emission lines are detected in the spectra of all four galaxies, compatible with predictions for Lyman continuum leakers. We find escape fractions of Lyalpha, fesc(Lyalpha) ~20%-40%, among the highest known for Lyalpha emitters (LAEs). Surface brightness profiles produced from the COS acquisition images reveal bright star-forming regions in the center and exponential discs in the outskirts with disc scale lengths alpha in the range ~0.6-1.4 kpc. Our galaxies are characterized by low metallicity, ~1/8-1/5 solar, low stellar mass ~(0.2 - 4)e9 Msun, high star formation rates SFR~14-36 Msun/yr, and high SFR densities Sigma~2-35 Msun/yr/kpc^2. These properties are comparable to those of high-redshift star-forming galaxies. Finally, our observations, combined with our first detection reported in Izotov et al. (2016), reveal that a selection for compact star-forming galaxies showing high [OIII]5007/[OII]3727 ratios appears to pick up very efficiently sources with escaping Lyman continuum radiation: all five of our selected galaxies are LyC leakers.

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

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.

One of the key questions in observational cosmology is the identification of the sources responsible for ionisation of the Universe after the cosmic Dark Ages, when the baryonic matter was neutral. The currently identified distant galaxies are insufficient to fully reionise the Universe by redshift z~6, but low-mass star-forming galaxies are thought to be responsible for the bulk of the ionising radiation. Since direct observations at high redshift are difficult for a variety of reasons, one solution is to identify local proxies of this galaxy population. However, starburst galaxies at low redshifts are generally opaque to their ionising radiation. This radiation with small escape fractions of 1-3% is directly detected only in three low-redshift galaxies. Here we present far-ultraviolet observations of a nearby low-mass star-forming galaxy, J0925+1403, selected for its compactness and high excitation. The galaxy is leaking ionising radiation, with an escape fraction of ~8%. The total number of photons emitted during the starburst phase is sufficient to ionize intergalactic medium material, which is about 40 times more massive than the stellar mass of the galaxy.

We present a sample of 38 intervening Damped Lyman $\alpha$ (DLA) systems identified towards 100 $z>3.5$ quasars, observed during the XQ-100 survey. The XQ-100 DLA sample is combined with major DLA surveys in the literature. The final combined sample consists of 742 DLAs over a redshift range approximately $1.6 < z_{\rm abs} < 5.0$. We develop a novel technique for computing $\Omega_{\rm HI}^{\rm DLA}$ as a continuous function of redshift, and we thoroughly assess and quantify the sources of error therein, including fitting errors and incomplete sampling of the high column density end of the column density distribution function. There is a statistically significant redshift evolution in $\Omega_{\rm HI}^{\rm DLA}$ ($\geq 3 \sigma$) from $z \sim 2$ to $z \sim$ 5. In order to make a complete assessment of the redshift evolution of $\Omega_{\rm HI}$, we combine our high redshift DLA sample with absorption surveys at intermediate redshift and 21cm emission line surveys of the local universe. Although $\Omega_{\rm HI}^{\rm DLA}$, and hence its redshift evolution, remains uncertain in the intermediate redshift regime ($0.1 < z_{\rm abs} < 1.6$), we find that the combination of high redshift data with 21cm surveys of the local universe all yield a statistically significant evolution in $\Omega_{\rm HI}$ from $z \sim 0$ to $z \sim 5$ ($\geq 3 \sigma$). Despite its statistical significance, the magnitude of the evolution is small: a linear regression fit between $\Omega_{\rm HI}$ and $z$ yields a typical slope of $\sim$0.17$\times 10^{-3}$, corresponding to a factor of $\sim$ 4 decrease in $\Omega_{\rm HI}$ between $z=5$ and $z=0$.

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

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.

KA1858+4850 is a narrow-line Seyfert 1 galaxy at redshift 0.078 and is among the brightest active galaxies monitored by the Kepler mission. We have carried out a reverberation mapping campaign designed to measure the broad-line region size and estimate the mass of the black hole in this galaxy. We obtained 74 epochs of spectroscopic data using the Kast Spectrograph at the Lick 3-m telescope from February to November of 2012, and obtained complementary V-band images from five other ground-based telescopes. We measured the H-beta light curve lag with respect to the V-band continuum light curve using both cross-correlation techniques (CCF) and continuum light curve variability modeling with the JAVELIN method, and found rest-frame lags of lag_CCF = 13.53 (+2.03, -2.32) days and lag_JAVELIN = 13.15 (+1.08, -1.00) days. The H-beta root-mean-square line profile has a width of sigma_line = 770 +/- 49 km/s. Combining these two results and assuming a virial scale factor of f = 5.13, we obtained a virial estimate of M_BH = 8.06 (+1.59, -1.72) x 10^6 M_sun for the mass of the central black hole and an Eddington ratio of L/L_Edd ~ 0.2. We also obtained consistent but slightly shorter emission-line lags with respect to the Kepler light curve. Thanks to the Kepler mission, the light curve of KA1858+4850 has among the highest cadences and signal-to-noise ratios ever measured for an active galactic nucleus; thus, our black hole mass measurement will serve as a reference point for relations between black hole mass and continuum variability characteristics in active galactic nuclei.

We present HST/COS observations of highly ionized absorption lines associated with a radio-loud QSO at $z=1.1319$. The absorption system has multiple velocity components, tracing gas that is largely outflowing from the QSO at velocities of a few 100 km s$^{-1}$. There is an unprecedented range in ionization, with detections of HI, NIII, NIV, NV, OIV, OIV*, OV, OVI, NeVIII, MgX, SV and ArVIII. We estimate the total hydrogen number density from the column density ratio N(OIV*)/N(OIV) to be $\log(n_{\textrm{H}}/\textrm{cm}^3)\sim 3$. Assuming photoionization equilibrium, we derive a distance to the absorbing complex of $2.3<R<6.0$ kpc from the centre of the QSO. A range in ionization parameter, covering $\sim 2$ orders of magnitude, suggest absorption path lengths in the range $10^{-4.5}<l_{\textrm{abs}}<1$ pc. In addition, the absorbing gas only partially covers the background emission from the QSO continuum, which suggests clouds with transverse sizes $l_{\textrm{trans}}<10^{-2.5}$ pc. Widely differing absorption path lengths, combined with covering fractions less than unity across all ions pose a challenge to models involving simple cloud geometries. These issues may be mitigated by the presence of non-equilibrium effects, together with the possibility of multiple gas temperatures. The dynamics and expected lifetimes of the gas clouds suggest that they do not originate from close to the AGN, but are instead formed close to their observed location. Their inferred distance, outflow velocities and gas densities are broadly consistent with scenarios involving gas entrainment or condensations in winds driven by either supernovae, or the supermassive black hole accretion disc. In the case of the latter, the present data most likely does not trace the bulk of the outflow by mass, which could instead manifest itself as an accompanying warm absorber, detectable in X-rays.

Previous studies of the 2.2 < z< 2.7 HeII Lyman-alpha forest measured much larger ionizing background fluctuations than are anticipated theoretically. We re-analyze recent Hubble Space Telescope data from the two HeII sightlines that have been used to make these measurements, HE2347-4342 and HS1700+6416, and find that the vast majority of the absorption is actually consistent with a single HeII photoionization rate. We show that the data constrains the RMS fractional fluctuation level smoothed at 1 Mpc to be < 2 and discuss why other studies had found larger fluctuations. Our measurement is consistent with models in which quasars dominate the z=2.5 metagalactic HeII-ionizing background (but it can accommodate less compelling models), and it suggests that quasars (rather than stars) are the dominant contributor to the HI-ionizing background. We detect a HeII transverse proximity effect that is slightly offset in redshift from a known quasar. Its profile and offset may indicate that the quasar turned on 10 Myr ago.

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

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

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

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

We present the first search for galaxy counterparts of intervening high-z (2<z< 3.6) sub-DLAs and DLAs towards GRBs. Our final sample comprises of five intervening sub-DLAs and DLAs in four GRB fields. To identify candidate galaxy counterparts of the absorbers we use deep optical and near-infrared imaging, and low-, mid- and high-resolution spectroscopy acquired with 6 to 10-m class telescopes, the Hubble and the Spitzer space telescopes. Furthermore, we use the spectroscopic information and spectral-energy-distribution fitting techniques to study them in detail. Our main result is the detection and spectroscopic confirmation of the galaxy counterpart of the intervening DLA at z=3.096 in the field of GRB 070721B (z_GRB=3.6298) as proposed by other authors. We also identify good candidates for the galaxy counterparts of the two strong MgII absorbers at z=0.6915 and 1.4288 towards GRB 050820A (z_GRB=2.615). The properties of the detected DLA galaxy are typical for Lyman-break galaxies (LBGs) at similar redshifts; a young, highly starforming galaxy that shows evidence for a galactic outflow. This supports the hypothesis that a DLA can be the gaseous halo of an LBG. In addition, we report a redshift coincidence of different objects associated with metal lines in the same field, separated by 130-161 kpc. The high detection rate of three correlated structures on a length scale as small as ~150 kpc in two pairs of lines of sight is intriguing. The absorbers in each of these are most likely not part of the same gravitationally bound structure. They more likely represent groups of galaxies.

On 2012 May 17.2 UT, only 1.5 +/- 0.2 d after explosion, we discovered SN 2012cg, a Type Ia supernova (SN Ia) in NGC 4424 (d ~ 15 Mpc). As a result of the newly modified strategy employed by the Lick Observatory SN Search, a sequence of filtered images was obtained starting 161 s after discovery. Utilizing recent models describing the interaction of SN ejecta with a companion star, we rule out a ~1 M_Sun companion for half of all viewing angles and a red-giant companion for nearly all orientations. SN 2012cg reached a B-band maximum of 12.09 +/- 0.02 mag on 2012 June 2.0 and took ~17.3 d from explosion to reach this, typical for SNe Ia. Our pre-maximum brightness photometry shows a narrower-than-average B-band light curve for SN 2012cg, though slightly overluminous at maximum brightness and with normal color evolution (including some of the earliest SN Ia filtered photometry ever obtained). Spectral fits to SN 2012cg reveal ions typically found in SNe Ia at early times, with expansion velocities >14,000 km/s at 2.5 d past explosion. Absorption from C II is detected early, as well as high-velocity components of both Si II 6355 Ang. and Ca II. Our last spectrum (13.5 d past explosion) resembles that of the somewhat peculiar SN Ia 1999aa. This suggests that SN 2012cg will have a slower-than-average declining light curve, which may be surprising given the faster-than-average rising light curve.

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

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

Understanding the origin and diversity of emission processes responsible for Gamma-ray Bursts (GRBs) remains a pressing challenge. While prompt and contemporaneous panchromatic observations have the potential to test predictions of the internal-external shock model, extensive multiband imaging has been conducted for only a few GRBs. We present rich, early-time, multiband datasets for two \swift\ events, GRB 110205A and GRB 110213A. The former shows optical emission since the early stages of the prompt phase, followed by the steep rising in flux up to ~1000s after the burst ($t^{-\alpha}$ with $\alpha=-6.13 \pm 0.75$). We discuss this feature in the context of the reverse-shock scenario and interpret the following single power-law decay as being forward-shock dominated. Polarization measurements, obtained with the RINGO2 instrument mounted on the Liverpool Telescope, also provide hints on the nature of the emitting ejecta. The latter event, instead, displays a very peculiar optical to near-infrared lightcurve, with two achromatic peaks. In this case, while the first peak is probably due to the onset of the afterglow, we interpret the second peak to be produced by newly injected material, signifying a late-time activity of the central engine.

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.

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

We have measured redshifts for 243 z ~3 quasars in nine VLT VIMOS LBG redshift survey areas, each of which is centred on a known bright quasar. Using spectra of these quasars, we measure the cross-correlation between neutral hydrogen gas causing the Lya forest and 1020 Lyman-break galaxies at z ~3. We find an increase in neutral hydrogen absorption within 5 h^-1 Mpc of a galaxy in agreement with the results of Adelberger et al. (2003, 2005). The Lya-LBG cross-correlation can be described by a power-law on scales larger than 3 h^-1 Mpc. When galaxy velocity dispersions are taken into account our results at smaller scales (<2 h^-1 Mpc) are also in good agreement with the results of Adelberger et al. (2005). There is little immediate indication of a region with a transmission spike above the mean IGM value which might indicate the presence of star-formation feedback. To measure the galaxy velocity dispersions, which include both intrinsic LBG velocity dispersion and redshift errors, we have used the LBG-LBG redshift space distortion measurements of Bielby et al. (2010). We find that the redshift-space transmission spike implied in the results of Adelberger et al. (2003) is too narrow to be physical in the presence of the likely LBG velocity dispersion and is likely to be a statistical fluke. Nevertheless, neither our nor previous data can rule out the presence of a narrow, real-space transmission spike, given the evidence of the increased Lya absorption surrounding LBGs which can mask the spike's presence when convolved with a realistic LBG velocity dispersion. Finally, we identify 176 CIV systems in the quasar spectra and find an LBG-CIV correlation strength on scales of 10 h^-1 Mpc consistent with the relation measured at ~Mpc scales.

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

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

(Abridged) We present the results from the largest investigation to date of the proximity effect in the HI Lyalpha forest, using the fifth SDSS data release. The sample consists of 1733 QSOs at redshifts z>2.3 and S/N>10. We adopted the flux statistic to infer the evolution of the HI effective optical depth in the Lyalpha forest between 2<z<4.5, finding very good agreement with results from high-resolution QSO samples. We compared the average opacity close to the quasar emission with its expected behavior in the Lyalpha forest and estimated the signature of the proximity effect towards individual objects at high significance in about 98% of the QSOs. Dividing the whole sample of objects in eight subsets according to their emission redshift, we inferred the proximity effect strength distribution (PESD) on each of them finding in all cases a prominent peak and an extending tail towards values associated to a weak effect. We provide for the first time observational evidence for an evolution in the asymmetry of the PESD with redshift. Adopting the modal values of the PESDs as our best and unbiased estimates of the UV background photoionization rate (Gamma_HI), we determine its evolution within the range 2.3<z<4.6. Our measurements do not show any significant decline towards high redshift and are located at Log(Gamma_HI)=-11.78\pm0.07 in units of s^-1. We decompose the observed photoionization rate into two major contributors: quasar and star-forming galaxies. By modeling the quasar contribution with different luminosity functions we estimated their contribution to Gamma_HI, thus putting a constraint on the residuals. We conclude that independently of the assumed luminosity function, stars are dominating the UVB at z>3.

We used 40 high resolution, high S/N QSO spectra at 2.1<z<4.7 to search for the signature of the proximity effect in the HI Lyalpha forest. Comparing the effective optical depth near each QSO with the expected one, we clearly detect the proximity effect on the combined QSO sample and towards each individual QSO. The observed proximity effect strength distribution (PESD) is asymmetric towards a weak effect. We demonstrate that this is not simply an effect of gravitational clustering around QSOs. Comparing simulated PESDs with observations, we argue that the averaging method to determine the UVB intensity J is heavily biased towards high values because of the PESD asymmetry. Using instead the mode of the PESD provides an unbiased estimate of J. For our sample its modal value is log(J)=-21.51+/-0.15 (in units of ergcm^-2s^-1Hz^-1sr^-1) at z=2.73. We estimated the excess HI absorption attributed to gravitational clustering. On scales of ~3 Mpc, only a minority of QSOs shows overdensities of up to a factor of a few in tau_eff; these are exactly the objects with the weakest proximity effects. After removing them, we redetermined the UVB intensity arriving at log(J)=-21.46+0.14-0.21. This is the most accurate measurement of J to date. We present a new diagnostic based on the shape of the PESD which strongly supports our conclusion that there is no systematic overdensity bias for the proximity effect. This additional diagnostic breaks the otherwise unavoidable degeneracy of the proximity effect between UVB and overdensity. We estimated the redshift evolution of J and found tentative evidence for a mild decrease with increasing redshift. Our results are in excellent agreement with predictions for the evolving UVB intensity, supporting the notion of a substantial contribution of star-forming galaxies.

We present the results of the Quasars near Quasars (QNQ) survey, a CCD-based slitless spectroscopic survey for faint V<22 quasars at 1.7<z<3.6 on 18 26.2'x33.5' fields centred on bright quasars at 2.76<z<4.69. In total 169 quasar candidates with emission lines were selected from the extracted flux-calibrated spectra on the basis of well-defined automatic selection criteria followed by visual inspection and verification. With follow-up spectroscopy of 81 candidates that were likely to reside at z>1.7 we were able to confirm 80 new quasars at 0.580<z<3.586 on 16 of our fields. 64 of the newly discovered quasars are located at z>1.7. The overall high success rate implies that most of the remaining 88 candidates are quasars as well, although the majority of them likely resides at z<1.7 on the basis of the observed line shapes and strengths. Due to the insufficient depth of the input source catalogues needed for extraction of the slitless spectra our survey is not well defined in terms of limiting magnitude for faint 2.5<z<3.6 quasars whose Lyman alpha emission is detectable well beyond V=22, albeit at a continuum S/N<1. While not useful for characterising the evolving space density of quasars, our sample provides many new closely spaced quasar sightlines around intensely studied quasars for further investigations on the three-dimensional distribution of the intergalactic medium.

We exploit a set of high signal-to-noise (~70), low-resolution (R~800) quasar spectra to search for the signature of the so-called proximity effect in the HI Ly alpha forest. Our sample consists of 17 bright quasars in the redshift range 2.7<z<4.1. Analysing the spectra with the flux transmission technique, we detect the proximity effect in the sample at high significance. We use this to estimate the average intensity of the metagalactic UV background, assuming it to be constant over this redshift range. We obtain a value of J = (9+-4)x10^{-22}ergcm^{-2}s^{-1}Hz^{-1}sr^{-1}, in good agreement with previous measurements at similar z. We then apply the same procedure to individual lines of sight, finding a significant deficit in the effective optical depth close to the emission redshift in every single object except one (which by a different line of evidence does nevertheless show a noticeable proximity effect). Thus, we clearly see the proximity effect as a universal phenomenon associated with individual quasars. Using extensive Monte-Carlo simulations to quantify the error budget, we assess the expected statistical scatter in the strength of the proximity effect due to shot noise (cosmic variance). The observed scatter is larger than the predicted one, so that additional sources of scatter are required. We rule out a dispersion of spectral slopes as a significant contributor. Possible effects are long time-scale variability of the quasars and/or gravitational clustering of Ly alpha forest lines. We speculate on the possibility of using the proximity effect as a tool to constrain individual quasar ages, finding that ages between ~10^6 and ~10^8 yrs might produce a characteristic signature in the optical depth profile towards the QSO. We identify one possible candidate for this effect in our sample.

We report the discovery of 14 quasars in the vicinity of HE2347-4342, one of the two quasars whose intergalactic HeII forest has been resolved with FUSE. By analysing the HI and the HeII opacity variations separately, no transverse proximity effect is detected near three foreground quasars of HE2347-4342: QSOJ23503-4328 (z=2.282, $\vartheta=3.59$ arcmin), QSOJ23500-4319 (z=2.302, $\vartheta=8.77$ arcmin) and QSOJ23495-4338 (z=2.690, $\vartheta=16.28$ arcmin). This is primarily due to line contamination and overdensities probably created by large-scale structure. By comparing the HI absorption and the corresponding HeII absorption, we estimated the fluctuating spectral shape of the extragalactic UV radiation field along this line of sight. We find that the UV spectral shape near HE2347-4342 and in the projected vicinity of the three foreground quasars is statistically harder than expected from UV background models dominated by quasars. In addition, we find three highly ionised metal line systems near the quasars. However, they do not yield further constraints on the shape of the ionising field. We conclude that the foreground quasars show a transverse proximity effect that is detectable as a local hardening of the UV radiation field, although the evidence is strongest for QSOJ23495-4338. Thus, the relative spectral hardness traces the proximity effect also in overdense regions prohibiting the traditional detection in the HI forest. Furthermore, we emphasise that softening of quasar radiation by radiative transfer in the intergalactic medium is important to understand the observed spectral shape variations. From the transverse proximity effect of QSOJ23495-4338 we obtain a lower limit on the quasar lifetime of ~25 Myr.

We report on our recent detection of the transverse proximity effect of 4 foreground quasars near Q0302-003 (z=3.285) as a local hardness fluctuation in the spectral shape of the intergalactic UV radiation field, found by correlating the HI Ly\alpha absorption with the corresponding HeII Ly\alpha absorption. We argue that the spectral hardness is a sensitive physical measure to reveal the influence of quasars onto the UV background over scales of several Mpc, and that it breaks the density degeneracy hampering the search for the elusive transverse proximity effect. From the transverse proximity effect on scales of several Mpc we obtain minimum quasar lifetimes of ~10-30Myr.

We report the discovery of the faint (V=21.7) quasar QSO03027-0010 at z=2.808 in the vicinity of Q0302-003, one of the few quasars observed with STIS to study intergalactic HeII absorption. Together with another newly discovered QSO at z=2.29, there are now 6 QSOs known near the line of sight towards Q0302-003, of which 4 are located within the redshift region 2.76<=z<=3.28 covered by the STIS spectrum. We correlated the opacity variations in the HI and HeII Lyman forest spectra with the locations of known quasars. There is no significant proximity effect in the HI Lyman alpha forest for any of the QSOs, except for the well-known line of sight effect for Q0302-003 itself. By comparing the absorption properties in HI and HeII, we estimated the fluctuating hardness of the extragalactic UV radiation field along this line of sight. We find that close to each foreground quasar, the ionizing background is considerably harder than on average. In particular, our newly discovered QSO03027-0010 shows such a hardness increase despite being associated with an overdensity in the HI Lyman forest. We argue that the spectral hardness is a sensitive physical measure to reveal the influence of QSOs onto the UV background even over scales of several Mpc, and that it breaks the density degeneracy hampering the traditional transverse proximity effect analysis. We infer from our sample that there is no need for significantly anisotropic UV radiation from the QSOs. From the transverse proximity effect detected in the sample we obtain minimum quasar lifetimes in the range ~10-30 Myr.