We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$\alpha$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$\beta$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$\beta$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$\alpha$ and Lyman-$\gamma$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively.

We explore a potential new probe of reionization: the cross-correlation of high-redshift galaxies with maps of the thermal Sunyaev-Zel'dovich (tSZ) effect. We consider two types of high redshift galaxies: Lyman break galaxies (LBGs) and Lyman-$\alpha$ emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift ($z \approx$ 4 to 7) by ongoing and future surveys. We consider a LBG sample like that expected from the The Rubin Observatory Legacy Survey of Space and Time (LSST), and a selection of LAEs modelled after the Subaru SILVERRUSH program, but covering a much larger sky fraction. The tSZ effect is sensitive to a line-of-sight integral of the ionized gas pressure, and can be measured across large patches of the sky using multi-frequency CMB surveys. We consider forecast tSZ maps from CMB Stage 4 and more futuristic observations. Using a suite of hydrodynamical simulations, we show that the high-redshift galaxies are correlated with the tSZ signal from reionization, with cross-power amplitude of order $10^{-15}$ at $\ell \sim 1000$. The cross-spectra between LBGs/LAEs with tSZ maps contain information about the reionization history of the Universe, such as the distribution of bubble sizes, and could in principle be used to directly measure the timing of reionization. The amplitude of the signal is small, however, and its detectability is hindered by low-redshift contributions to tSZ maps and by instrumental noise. If the the low-redshift contribution to the observed tSZ signal is suppressed by masking of massive halos, a combination of overlapping futuristic CMB and galaxy surveys could potentially probe this signal.

We outline the experimental concept and key scientific capabilities of AION (Atom Interferometer Observatory and Network), a proposed UK-based experimental programme using cold strontium atoms to search for ultra-light dark matter, to explore gravitational waves in the mid-frequency range between the peak sensitivities of the LISA and LIGO/Virgo/ KAGRA/INDIGO/Einstein Telescope/Cosmic Explorer experiments, and to probe other frontiers in fundamental physics. AION would complement other planned searches for dark matter, as well as probe mergers involving intermediate mass black holes and explore early universe cosmology. AION would share many technical features with the MAGIS experimental programme in the US, and synergies would flow from operating AION in a network with this experiment, as well as with other atom interferometer experiments such as MIGA, ZAIGA and ELGAR. Operating AION in a network with other gravitational wave detectors such as LIGO, Virgo and LISA would also offer many synergies.

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, but it also remains uncertain whether the low-redshift LyC leakers are true "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission-line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, we find that identifying low-redshift galaxies based on [SII] deficiencies does not seem to produce true analogues. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158\mu m}$ and [OIII]$_{\rm 88\mu m}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149\_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

The search for primordial gravitational waves through the $B$-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary $B$-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper, we revise estimates of the reionization-induced $B$-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-$\alpha$ forest. To derive these $B$-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Ly$\alpha$ data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles $\ell\approx 100$, reionization induces $B$-mode power with $\ell(\ell+1)C_\ell^{BB}/2\pi\approx 4\times 10^{-6}\,\mu$K$^2$. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio $r<10^{-4}$, and can increase by a factor of $\sim 50$ if reionization is sourced by highly clustered sources residing in haloes with mass of $\sim 10^{11}$ M$_\odot$. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments.

We compare a sample of five high-resolution, high S/N Ly$\alpha$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

We model the 21-cm signal and LAE population evolution during the epoch of reionization in order to predict the 21cm-LAE cross-power spectrum. We employ high-dynamic-range simulations of the IGM to create models that are consistent with constraints from the CMB, Lyman-$\alpha$ forest and LAE population statistics. Using these models we consider the evolution of the cross-power spectrum for a selection of realistic reionization histories and predict the sensitivity of current and upcoming surveys to measuring this signal. We find that the imprint of a delayed-end to reionization can be observed by future surveys, and that strong constraints can be placed on the progression of reionization as late as $z=5.7$ using a Subaru-SKA survey. We make predictions for the signal-to-noise ratios achievable by combinations of Subaru/PFS with the MWA, LOFAR, HERA and SKA interferometers for an integration time of 1000 hours. We find that a Subaru-SKA survey could measure the cross-power spectrum for a late reionization at $z=6.6$ with a total signal-to-noise greater than 5, making it possible to constrain both the timing and bubble size at the end of reionization. Furthermore, we find that expanding the current Subaru/PFS survey area and depth by a factor of three would double the total signal-to-noise.

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $\Delta z>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

A long (110 cMpc/$h$) and deep absorption trough in the Ly$\alpha$ forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly$\alpha$ transmission is detected, Ly$\beta$ spikes are present which has led to claims that the gas along this trough must be ionized. Using high resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late ($z \sim 5.2$), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly$\beta$. We have also modelled the Ly$\alpha$ emitter population around the simulated troughs, and show that there is a deficit of Ly$\alpha$ emitters close to the trough as is observed.

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

Infrared and nebular lines provide some of our best probes of the physics regulating the properties of the interstellar medium (ISM) at high-redshift. However, interpreting the physical conditions of high-redshift galaxies directly from emission lines remains complicated due to inhomogeneities in temperature, density, metallicity, ionisation parameter, and spectral hardness. We present a new suite of cosmological, radiation-hydrodynamics simulations, each centred on a massive Lyman-break galaxy that resolves such properties in an inhomogeneous ISM. Many of the simulated systems exhibit transient but well defined gaseous disks that appear as velocity gradients in [CII]~158.6$\mu$m emission. Spatial and spectral offsets between [CII]~158.6$\mu$m and [OIII]~88.33$\mu$m are common, but not ubiquitous, as each line probes a different phase of the ISM. These systems fall on the local [CII]-SFR relation, consistent with newer observations that question previously observed [CII]~158.6$\mu$m deficits. Our galaxies are consistent with the nebular line properties of observed $z\sim2-3$ galaxies and reproduce offsets on the BPT and mass-excitation diagrams compared to local galaxies due to higher star formation rate (SFR), excitation, and specific-SFR, as well as harder spectra from young, metal-poor binaries. We predict that local calibrations between H$\alpha$ and [OII]~3727$\AA$ luminosity and galaxy SFR apply up to $z>10$, as do the local relations between certain strong line diagnostics (R23 and [OIII]~5007$\AA$/H$\beta$) and galaxy metallicity. Our new simulations are well suited to interpret the observations of line emission from current (ALMA and HST) and upcoming facilities (JWST and ngVLA).

Mapping the intergalactic medium (IGM) in Lyman-$\alpha$ emission would yield unprecedented tomographic information on the large-scale distribution of baryons and potentially provide new constraints on the UV background and various feedback processes relevant for galaxy formation. Here, we examine the Lyman-$\alpha$ emission of the moderate-density IGM due to collisional excitations and recombinations in the presence of a UV background in the Sherwood simulation suite. We focus on large-scale-structure filaments in which Lyman-$\alpha$ radiative transfer effects are expected to be moderate. At low density the emission is primarily due to fluorescent re-emission of the ionising UV background due to recombinations, while collisional excitations dominate at higher densities. We discuss prospects of current and future observational facilities to detect this emission and focus on VLT/MUSE for a more detailed sensitivity analysis. We construct mock MUSE observations resembling the MUSE Hubble Deep Field South from our simulations and show that our predictions are consistent with recent analyses of diffuse Lyman-$\alpha$ emission using MUSE at redshifts $3 < z < 6$. We find that it should be most feasible to detect the Lyman-$\alpha$ emission from filaments in the IGM in overdense regions, somewhat surprisingly towards the high-redshift end ($z \gtrsim 4.5$) accessible by MUSE, and if narrowband widths close to the MUSE spectral resolution limit are used. This is due to the higher densities and lower temperatures in the IGM at higher redshift. High-redshift protoclusters therefore appear to be the ideal environment to observe filamentary structures in the IGM in Lyman-$\alpha$ emission.

Understanding the escape of Lyman continuum (LyC) and Lyman alpha (Lya) photons from molecular clouds is one of the keys to constraining the reionization history of the Universe and the evolution of galaxies at high redshift. Using a set of radiation-hydrodynamic simulations with adaptive mesh refinement, we investigate how photons propagate and escape from turbulent clouds with different masses, star formation efficiencies (SFEs), and metallicities, as well as with different models of stellar spectra and supernova feedback. We find that the escape fractions in both LyC and Lya are generally increasing with time if the cloud is efficiently dispersed by radiation and supernova feedback. When the total SFE is low (1% of the cloud mass), 0.1-5% of LyC photons leave the metal-poor cloud, whereas the fractions increase to 20-70% in clouds with a 10% SFE. LyC photons escape more efficiently if gas metallicity is lower, if the upper mass limit in the stellar initial mass function is higher, if binary interactions are allowed in the evolution of stars, or if additional strong radiation pressure, such as Lya pressure, is present. As a result, the number of escaping LyC photons can easily vary by a factor of $\sim4$ on cloud scales. The escape fractions of Lya photons are systemically higher (60-80%) than those of LyC photons despite large optical depths at line centre ($\tau_0\sim10^6-10^9$). Scattering of Lya photons is already significant on cloud scales, leading to double-peaked profiles with peak separations of $v_{\rm sep}\sim400\,{\rm km\,s^{-1}}$ during the initial stage of the cloud evolution, while it becomes narrower than $v_{\rm sep} \le 150 \, {\rm km\,s^{-1}}$ in the LyC bright phase. Comparisons with observations of low-redshift galaxies suggest that Lya photons require further interactions with neutral hydrogen to reproduce their velocity offset for a given LyC escape fraction.

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.

Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.

Observations of the high-redshift Universe using the 21 cm line of neutral hydrogen and complimentary emission lines from the first galaxies promise to open a new door for our understanding of the epoch of reionization. We present predictions for the [C II] 158-micron line and H I 21 cm emission from redshifts z=6--9 using high-dynamic-range cosmological simulations combined with semi-analytical models. We find that the CONCERTO experiment should be able to detect the large scale power spectrum of [C II] emission to redshifts of up to z=8 (signal-to-noise ratio ~ 1 at k = 0.1 h/cMpc with 1500 hr of integration). A Stage II experiment similar to CCAT-p should be able to detect [C II] from even higher redshifts to high significance for similar integration times (signal-to-noise ratio of ~50 at k = 0.2 h/cMpc at z=6--9). We study the possibility of combining such future [C II] measurements with 21 cm measurements using LOFAR and SKA to measure the [C II]-21cm cross power spectra, and find that a Stage II experiment should be able to measure the cross-power spectrum for k < 1 h/cMpc to signal-to-noise ratio of better than 10. We discuss the capability of such measurements to constrain astrophysical parameters relevant to reionization and show that a measurement of the [C II]-21cm cross power spectrum helps break the degeneracy between the mass and brightness of ionizing sources.

High-redshift QSO spectra show large spatial fluctuations in the Ly-alpha opacity of the intergalactic medium on surprisingly large scales at z>~5.5. We present a radiative transfer simulation of cosmic reionization driven by galaxies that reproduces this large scatter and the rapid evolution of the Ly-alpha opacity distribution at 5<z<6. The simulation also reproduces the low Thomson scattering optical depth reported by the latest CMB measurement and is consistent with the observed short near-zones and strong red damping wings in the highest-redshift QSOs. It also matches the rapid disappearance of observed Ly-alpha emission by galaxies at z>~6. Reionization is complete at z=5.3 in our model, and 50% of the volume of the Universe is ionized at z=7. Agreement with the Ly-alpha forest data in such a late reionization model requires a rapid evolution of the ionizing emissivity of galaxies that peaks at z~6.8. The late end of reionization results in a large scatter in the photoionisation rate and the neutral hydrogen fraction at redshifts as low as z<~5.5 with large residual neutral 'islands' that can produce very long Gunn-Peterson troughs resembling those seen in the data.

We model the high redshift (z > 5) Lyman-$\alpha$ emitting (LAE) galaxy population using the empirical rest-frame equivalent width distribution. We calibrate to the observed luminosity function and angular correlation function at z = 5.7 as measured by the SILVERRUSH survey. This allows us to populate the high-dynamic-range Sherwood simulation suite with LAEs, and to calculate the transmission of their Ly $\alpha$ emission through the inter-galactic medium (IGM). We use this simulated population to explore the effect of the IGM on high-redshift observations of LAEs, and make predictions for the narrowband filter redshifts at z = 6.6, 7.0 and 7.3. Comparing our model with existing observations, we find a late reionization is suggested, consistent with the recent low optical depth derived from the cosmic microwave background (CMB) by the Planck collaboration and the opacity fluctuations in the Ly $\alpha$ forest. We also explore the role of the circum-galactic medium (CGM) and the larger volume of gas which is infalling onto the host halo versus the IGM in attenuating the Ly $\alpha$ signal, finding that a significant fraction of the attenuation is due to the CGM and infalling gas, which increases towards the end of reionization, albeit with a large scatter across the mock LAE population.

We use the photon flux and absorption tracer algorithm presented in Katz et al.~2018, to characterise the contribution of haloes of different mass and stars of different age and metallicity to the reionization of the Universe. We employ a suite of cosmological multifrequency radiation hydrodynamics AMR simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \geq 6$. In our simulations, haloes with mass $10^9{\rm M_{\odot}}h^{-1}<M<10^{10}{\rm M_{\odot}}h^{-1}$, stars with metallicity $10^{-3}Z_{\odot}<Z<10^{-1.5}Z_{\odot}$, and stars with age $3\,{\rm Myr} < t < 10 \, {\rm Myr}$ dominate reionization by both mass and volume. We show that the sources that reionize most of the volume of the Universe by $z=6$ are not necessarily the same sources that dominate the meta-galactic UV background at the same redshift. We further show that in our simulations, the contribution of each type of source to reionization is not uniform across different gas phases. The IGM, CGM, filaments, ISM, and rarefied supernova heated gas have all been photoionized by different classes of sources. Collisional ionisation contributes at both the lowest and highest densities. In the early stages of the formation of individual HII bubbles, reionization proceeds with the formation of concentric shells of gas ionised by different classes of sources, leading to large temperature variations as a function of galacto-centric radius. The temperature structure of individual HII bubbles may thus give insight into the star formation history of the galaxies acting as the first ionising sources. Our explorative simulations highlight how the complex nature of reionization can be better understood by using our photon tracer algorithm.

We present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role.

We model the transmission of the Lyman-alpha line through the circum- and intergalactic media around dark matter haloes expected to host Lyman-alpha emitters (LAEs) at z > 5.7, using the high-dynamic-range Sherwood simulations. We find very different CGM environments around more massive haloes (~10^11 M_sun) compared to less massive haloes (~10^9 M_sun) at these redshifts, which can contribute to a different evolution of the Lyman-alpha transmission from LAEs within these haloes. Additionally we confirm that part of the differential evolution could result from bright LAEs being more likely to reside in larger ionized regions. We conclude that a combination of the CGM environment and the IGM ionization structure is likely to be responsible for the differential evolution of the bright and faint ends of the LAE luminosity function at z > 6. More generally, we confirm the suggestion that the self-shielded neutral gas in the outskirts of the host halo can strongly attenuate the Lyman-alpha emission from high redshift galaxies. We find that this has a stronger effect on the more massive haloes hosting brighter LAEs. The faint-end of the LAE luminosity function is thus a more reliable probe of the average ionization state of the IGM. Comparing our model for LAEs with a range of observational data we find that the favoured reionization histories are our previously advocated `Late' and `Very Late' reionization histories, in which reionization finishes rather rapidly at around z ~ 6.

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 a new, on-the-fly photon flux and absorption tracer algorithm designed to directly measure the contribution of different source populations to the metagalactic UV background and to the ionisation fraction of gas in the Universe. We use a suite of multifrequency radiation hydrodynamics simulations that are carefully calibrated to reproduce a realistic reionization history and galaxy properties at $z \ge 6$, to disentangle the contribution of photons emitted by different mass haloes and by stars with different metallicities and ages to the UV background during reionization. While at very early cosmic times low mass, metal poor haloes provide most of the LyC photons, their contribution decreases steadily with time. At $z = 6$ it is the photons emitted by massive systems (${\rm M_{halo}}/{\rm M_\odot} > 10^{10} \, {\rm h ^{-1}}$) and by the metal enriched stars ($10^{-3} < Z/Z_{\rm \odot} < 10^{-1.5}$) that provide the largest contribution to the ionising UV background. We demonstrate that there are large variations in the escape fraction depending on the source, with the escape fraction being highest ($\sim 45-60\%$) for photons emitted by the oldest stars that penetrate into the IGM via low opacity channels carved by the ionising photons and supernova from younger stars. Before HII regions begin to overlap, the photoionisation rate strongly fluctuates between different, isolated HII bubbles, depending on the embedded ionising source, which we suggest may result in spatial variations in the properties of dwarf galaxies.

We present cosmological radiation-hydrodynamic simulations, performed with the code Ramses-RT, of radiatively-driven outflows in a massive quasar host halo at $z = 6$. Our simulations include both single- and multi-scattered radiation pressure on dust from a quasar and are compared against simulations performed with thermal feedback. For radiation pressure-driving, we show that there is a critical quasar luminosity above which a galactic outflow is launched, set by the equilibrium of gravitational and radiation forces. While this critical luminosity is unrealistically high in the single-scattering limit for plausible black hole masses, it is in line with a $\approx 3 \times 10^9 \, \rm M_\odot$ black hole accreting at its Eddington limit, if infrared (IR) multi-scattering radiation pressure is included. The outflows are fast ($v \, \gtrsim \, 1000 \, \rm km \, s^{-1}$) and strongly mass-loaded with peak mass outflow rates $\approx 10^3 - 10^4 \, \rm M_\odot \, yr^{-1}$, but short-lived ($< 10 \, \rm Myr$). Outflowing material is multi-phase, though predominantly composed of cool gas, forming via a thermal instability in the shocked swept-up component. Radiation pressure- and thermally-driven outflows both affect their host galaxies significantly, but in different, complementary ways. Thermally-driven outflows couple more efficiently to diffuse halo gas, generating more powerful, hotter and more volume-filling outflows. IR radiation, through its ability to penetrate dense gas via diffusion, is more efficient at ejecting gas from the bulge. The combination of gas ejection through outflows with internal pressurisation by trapped IR radiation leads to a complete shut down of star formation in the bulge. We hence argue that radiation pressure-driven feedback may be an important ingredient in regulating star formation in compact starbursts, especially during the quasar's `obscured' phase.

We present measurements of the mean and scatter of the IGM Lyman-{\alpha} opacity at 4.9 < z < 6.1 along the lines of sight of 62 quasars at z > 5.7, the largest sample assembled at these redshifts to date by a factor of two. The sample size enables us to sample cosmic variance at these redshifts more robustly than ever before. The spectra used here were obtained by the SDSS, DES-VHS and SHELLQs collaborations, drawn from the ESI and X-Shooter archives, reused from previous studies or observed specifically for this work. We measure the effective optical depth of Lyman-{\alpha} in bins of 10, 30, 50 and 70 cMpc h-1, construct cumulative distribution functions under two treatments of upper limits on flux and explore an empirical analytic fit to residual Lyman-{\alpha} transmission. We verify the consistency of our results with those of previous studies via bootstrap re-sampling and confirm the existence of tails towards high values in the opacity distributions, which may persist down to z = 5.2. Comparing our results with predictions from cosmological simulations, we find further strong evidence against models that include a spatially uniform ionizing background and temperature-density relation. We also compare to IGM models that include either a fluctuating UVB dominated by rare quasars or temperature fluctuations due to patchy reionization. Although both models produce better agreement with the observations, neither fully captures the observed scatter in IGM opacity. Our sample of 62 z > 5.7 quasar spectra opens many avenues for future study of the reionisation epoch.

We present radiation-hydrodynamic simulations of radiatively-driven gas shells launched by bright active galactic nuclei (AGN) in isolated dark matter haloes. Our goals are (1) to investigate the ability of AGN radiation pressure on dust to launch galactic outflows and (2) to constrain the efficiency of infrared (IR) multi-scattering in boosting outflow acceleration. Our simulations are performed with the radiation-hydrodynamic code RAMSES-RT and include both single- and multi-scattered radiation pressure from an AGN, radiative cooling and self-gravity. Since outflowing shells always eventually become transparent to the incident radiation field, outflows that sweep up all intervening gas are likely to remain gravitationally bound to their halo even at high AGN luminosities. The expansion of outflowing shells is well described by simple analytic models as long as the shells are mildly optically thick to IR radiation. In this case, an enhancement in the acceleration of shells through IR multi-scattering occurs as predicted, i.e. a force dP/dt = tau_IR L/c is exerted on the gas. For high optical depths tau_IR > 50, however, momentum transfer between outflowing optically thick gas and IR radiation is rapidly suppressed, even if the radiation is efficiently confined. At high tau_IR, the characteristic flow time becomes shorter than the required trapping time of IR radiation such that the momentum flux dP/dt << tau_IR L/c. We argue that while unlikely to unbind massive galactic gaseous haloes, AGN radiation pressure on dust could play an important role in regulating star formation and black hole accretion in the nuclei of massive compact galaxies at high redshift.

We calibrate here cosmological radiative transfer simulation with ATON/RAMSES with a range of measurements of the Lyman alpha opacity from QSO absorption spectra. We find the Lyman alpha opacity to be very sensitive to the exact timing of hydrogen reionisation. Models reproducing the measured evolution of the mean photoionisation rate and average mean free path reach overlap at z ~ 7 and predict an accelerated evolution of the Lyman alpha opacity at z > 6 consistent with the rapidly evolving luminosity function of Lyman alpha emitters in this redshift range. Similar to "optically thin" simulations our full radiative transfer simulations fail, however, to reproduce the high-opacity tail of the Lyman alpha opacity PDF at z > 5. We argue that this is due to spatial UV fluctuations in the post-overlap phase of reionisation on substantially larger scales than predicted by our source model, where the ionising emissivity is dominated by large numbers of sub-L* galaxies. We further argue that this suggests a significant contribution to the ionising UV background by much rarer bright sources at high redshift.

The thermal state of the post-reionization IGM is sensitive to the timing of reionization and the nature of the ionizing sources. We have modelled here the thermal state of the IGM in cosmological radiative transfer simulations of a realistic, extended, spatially inhomogeneous hydrogen reionization process, carefully calibrated with Ly-alpha forest data. We compare these with cosmological simulations run using a spatially homogeneous ionizing background. The simulations with a realistic growth of ionized regions and a realistic spread in reionization redshifts show, as expected, significant spatial fluctuations in the temperature-density relation (TDR) of the post-reionization IGM. The most recently ionized regions are hottest and exhibit a flatter TDR. In simulations consistent with the average TDR inferred from Ly-alpha forest data, these spatial fluctuations have a moderate but noticeable effect on the statistical properties of the Ly-alpha opacity of the IGM at z ~ 4-6. This should be taken into account in accurate measurements of the thermal properties of the IGM and the free-streaming of dark matter from Ly-alpha forest data in this redshift range. The spatial variations of the TDR predicted by our simulations are, however, smaller by about a factor two than would be necessary to explain the observed large spatial opacity fluctuations on large (> 50 comoving Mpc/h) scales at z > 5.5.

High-signal to noise observations of the Ly$\alpha$ forest transmissivity in the z = 7.085 QSO ULAS J1120+0641 show seven narrow transmission spikes followed by a long 240 cMpc/h trough. Here we use radiative transfer simulations of cosmic reionization previously calibrated to match a wider range of Ly$\alpha$ forest data to show that the occurrence of seven transmission spikes in the narrow redshift range z = 5.85 - 6.1 is very sensitive to the exact timing of reionization. Occurrence of the spikes requires the most under dense regions of the IGM to be already fully ionised. The rapid onset of a long trough at z = 6.12 requires a strong decrease of the photo-ionisation rate at z$\sim$6.1 in this line-of-sight, consistent with the end of percolation at this redshift. The narrow range of reionisation histories that we previously found to be consistent with a wider range of Ly$\alpha$ forest data have a reasonable probability of showing seven spikes and the mock absorption spectra provide an excellent match to the spikes and the trough in the observed spectrum of ULAS J1120+0641. Despite the large overall opacity of Ly$\alpha$ at z > 5.8, larger samples of high signal-to-noise observations of rare transmission spikes should therefore provide important further insights into the exact timing of the percolation of HII bubbles at the tail-end of reionization

We investigate self-shielding of intergalactic hydrogen against ionizing radiation in radiative transfer simulations of cosmic reionization carefully calibrated with Lyman alpha forest data. While self-shielded regions manifest as Lyman-limit systems in the post-reionization Universe, here we focus on their evolution during reionization (redshifts z=6-10). At these redshifts, the spatial distribution of hydrogen-ionizing radiation is highly inhomogeneous, and some regions of the Universe are still neutral. After masking the neutral regions and ionizing sources in the simulation, we find that the hydrogen photoionization rate depends on the local hydrogen density in a manner very similar to that in the post-reionization Universe. The characteristic physical hydrogen density above which self-shielding becomes important at these redshifts is about $\mathrm{n_H \sim 3 \times 10^{-3} cm^{-3}}$, or $\sim$ 20 times the mean hydrogen density, reflecting the fact that during reionization photoionization rates are typically low enough that the filaments in the cosmic web are often self-shielded. The value of the typical self-shielding density decreases by a factor of 3 between redshifts z=3 and 10, and follows the evolution of the average photoionization rate in ionized regions in a simple fashion. We provide a simple parameterization of the photoionization rate as a function of density in self-shielded regions during the epoch of reionization.

Understanding the origin of strong galactic outflows and the suppression of star formation in dwarf galaxies is a key problem in galaxy formation. Using a set of radiation-hydrodynamic simulations of an isolated dwarf galaxy embedded in a $10^{10}\,M_\odot$ halo, we show that the momentum transferred from resonantly scattered Lyman-$\alpha$ (Lya) photons is an important source of stellar feedback which can shape the evolution of galaxies. We find that Lya feedback suppresses star formation by a factor of two in metal-poor galaxies by regulating the dynamics of star-forming clouds before the onset of supernova explosions (SNe). This is possible because each Lya photon resonantly scatters and imparts 10-300 times greater momentum than in the single scattering limit. Consequently, the number of star clusters predicted in the simulations is reduced by a factor of $\sim 5$, compared to the model without the early feedback. More importantly, we find that galactic outflows become weaker in the presence of strong Lya radiation feedback, as star formation and associated SNe become less bursty. We also examine a model in which radiation field is arbitrarily enhanced by a factor of up to 10, and reach the same conclusion. The typical mass loading factors in our metal-poor dwarf system are estimated to be $\sim5-10$ near the mid plane, while it is reduced to $\sim1$ at larger radii. Finally, we find that the escape of ionizing radiation and hence the reionization history of the Universe is unlikely to be strongly affected by Lya feedback.

Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$\alpha$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium.

We present two measurements of the temperature-density relationship (TDR) of the intergalactic medium (IGM) in the redshift range $2.55 < z < 2.95$ using a sample of 13 high-quality quasar spectra and high resolution numerical simulations of the IGM. Our approach is based on fitting the neutral hydrogen column density $N_{HI}$ and the Doppler parameter $b$ of the absorption lines in the \mlya\ forest. The first measurement is obtained using a novel Bayesian scheme which takes into account the statistical correlations between the parameters characterising the lower cut-off of the $b-N_{HI}$ distribution and the power-law parameters $T_0$ and $\gamma$ describing the TDR. This approach yields $T_0/ 10^3\, {\rm K}=15.6 \pm 4.4 $ and $\gamma=1.45 \pm 0.17$ independent of the assumed pressure smoothing of the small scale density field. In order to explore the information contained in the overall $b-N_{HI}$ distribution rather than only the lower cut-off, we obtain a second measurement based on a similar Bayesian analysis of the median Doppler parameter for separate column-density ranges of the absorbers. In this case we obtain $T_0/ 10^3\, {\rm K}=14.6 \pm 3.7$ and $\gamma=1.37 \pm 0.17$ in good agreement with the first measurement. Our Bayesian analysis reveals strong anti-correlations between the inferred $T_0$ and $\gamma$ for both methods as well as an anti-correlation of the inferred $T_0$ and the pressure smoothing length for the second method, suggesting that the measurement accuracy can in the latter case be substantially increased if independent constraints on the smoothing are obtained. Our results are in good agreement with other recent measurements of the thermal state of the IGM probing similar (over-)density ranges.

We present constraints on the masses of extremely light bosons dubbed fuzzy dark matter from Lyman-$\alpha$ forest data. Extremely light bosons with a De Broglie wavelength of $\sim 1$ kpc have been suggested as dark matter candidates that may resolve some of the current small scale problems of the cold dark matter model. For the first time we use hydrodynamical simulations to model the Lyman-$\alpha$ flux power spectrum in these models and compare with the observed flux power spectrum from two different data sets: the XQ-100 and HIRES/MIKE quasar spectra samples. After marginalization over nuisance and physical parameters and with conservative assumptions for the thermal history of the IGM that allow for jumps in the temperature of up to $5000\rm\,K$, XQ-100 provides a lower limit of 7.1$\times 10^{-22}$ eV, HIRES/MIKE returns a stronger limit of 14.3$\times 10^{-22}$ eV, while the combination of both data sets results in a limit of 20 $\times 10^{-22}$ eV (2$\sigma$ C.L.). The limits for the analysis of the combined data sets increases to 37.5$\times 10^{-22}$ eV (2$\sigma$ C.L.) when a smoother thermal history is assumed where the temperature of the IGM evolves as a power-law in redshift. Light boson masses in the range $1-10 \times10^{-22}$ eV are ruled out at high significance by our analysis, casting strong doubts that FDM helps solve the "small scale crisis" of the cold dark matter models.

We study the effect of different feedback prescriptions on the properties of the low redshift ($z\leq1.6$) Ly$\alpha$ forest using a selection of hydrodynamical simulations drawn from the Sherwood simulation suite. The simulations incorporate stellar feedback, AGN feedback and a simplified scheme for efficiently modelling the low column density Ly$\alpha$ forest. We confirm a discrepancy remains between Cosmic Origins Spectrograph (COS) observations of the Ly$\alpha$ forest column density distribution function (CDDF) at $z \simeq 0.1$ for high column density systems ($N_{\rm HI}>10^{14}\rm\,cm^{-2}$), as well as Ly$\alpha$ velocity widths that are too narrow compared to the COS data. Stellar or AGN feedback -- as currently implemented in our simulations -- have only a small effect on the CDDF and velocity width distribution. We conclude that resolving the discrepancy between the COS data and simulations requires an increase in the temperature of overdense gas with $\Delta=4$--$40$, either through additional He$\,\rm \scriptstyle II\ $ photo-heating at $z>2$ or fine-tuned feedback that ejects overdense gas into the IGM at just the right temperature for it to still contribute significantly to the Ly$\alpha$ forest. Alternatively a larger, currently unresolved turbulent component to the line width could resolve the discrepancy.

We present a search for metal absorption line systems at the highest redshifts to date using a deep (30h) VLT/X-Shooter spectrum of the z = 7.084 quasi-stellar object (QSO) ULAS J1120+0641. We detect seven intervening systems at z > 5.5, with the highest-redshift system being a C IV absorber at z = 6.51. We find tentative evidence that the mass density of C IV remains flat or declines with redshift at z < 6, while the number density of C II systems remains relatively flat over 5 < z < 7. These trends are broadly consistent with models of chemical enrichment by star formation-driven winds that include a softening of the ultraviolet background towards higher redshifts. We find a larger number of weak ( W_rest < 0.3A ) Mg II systems over 5.9 < z < 7.0 than predicted by a power-law fit to the number density of stronger systems. This is consistent with trends in the number density of weak Mg II systems at z = 2.5, and suggests that the mechanisms that create these absorbers are already in place at z = 7. Finally, we investigate the associated narrow Si IV, C IV, and N V absorbers located near the QSO redshift, and find that at least one component shows evidence of partial covering of the continuum source.

We present here predictions for the spatial distribution of 21 cm brightness temperature fluctuations from high-dynamic-range simulations for AGN-dominated reionization histories that have been tested against available Lyman-alpha and CMB data. We model AGN by extrapolating the observed M-sigma relation to high redshifts and assign them ionizing emissivities consistent with recent UV luminosity function measurements. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen 21 cm spin temperature is significantly larger than the CMB temperature. Our AGN-dominated reionization histories increase the variance of the 21 cm emission by a factor of up to ten compared to similar reionization histories dominated by faint galaxies, to values close to 100 mK^2 at scales accessible to experiments (k < 1 h/cMpc). This is lower than the sensitivity claimed to have been already reached by ongoing experiments by only a factor of about two or less. When reionization is dominated by AGN, the 21 cm power spectrum is enhanced on all scales due to the enhanced bias of the clustering of the more massive haloes and the peak in the large scale 21 cm power is strongly enhanced and moved to larger scales due to bigger characteristic bubble sizes. AGN dominated reionization should be easily detectable by LOFAR (and later HERA and SKA1) at their design sensitivity, assuming successful foreground subtraction and instrument calibration. Conversely, these could become the first non-trivial reionization scenarios to be ruled out by 21 cm experiments, thereby constraining the contribution of AGN to reionization.

We present cosmological, radiation-hydrodynamics simulations of galaxy formation during the epoch of reionization in an effort towards modelling the interstellar medium (ISM) and interpreting ALMA observations. Simulations with and without stellar radiation are compared at large (Mpc), intermediate (tens of kpc), and small (sub kpc) scales. At large scales, the dense regions around galaxies reionize first before UV photons penetrate the voids; however, considerable amounts of neutral gas remain present within the haloes. The spatial distribution of neutral gas is highly dynamic and is anti-correlated with the presence of stars older than a few Myrs. For our specific feedback implementation, most of the metals remain inside the virial radii of haloes and they are proportionally distributed over the ionized and neutral medium by mass. For our most massive galaxy with ${\rm M_h}\sim10^{11}$M$_{\odot}$, the majority of the CII and OI mass are associated with cold neutral clumps. NII is more diffuse and arises in warmer gas while OIII arises in hotter gas with a higher ionization parameter, produced by photo-heating and supernovae. If smaller pockets of high metallicity gas exist in the ISM, the emission from these ions may be observable by ALMA while the low metallicity of the galaxy may cause these systems to fall below the local [CII]-SFR relation. The presence of dust can cause spatial offsets between UV/Ly$\alpha$ and [CII] emission as suggested by the recent observations of Maiolino et al. [OIII] may be spatially offset from both of these components since it arises from a different part of density-temperature phase-space.

We compare the low redshift (z ~ 0.1) Lyman-alpha forest from hydrodynamical simulations with data from the Cosmic Origin Spectrograph (COS). We find tension between the observed number of lines with b-parameters in the range 25-45 km/s and the predictions from simulations that incorporate either vigorous feedback from active galactic nuclei or that exclude feedback altogether. The gas in these simulations is, respectively, either too hot to contribute to the Lyman-alpha absorption or too cold to produce the required line widths. Matching the observed b-parameter distribution therefore requires feedback processes that thermally or turbulently broaden the absorption features without collisionally (over-)ionising hydrogen. This suggests the Lyman-alpha forest b-parameter distribution is a valulable diagnostic of galactic feedback in the low redshift Universe. We furthermore confirm the low redshift Lyman-alpha forest column density distribution is better reproduced by an ultraviolet background with an HI photo-ionisation rate a factor 1.5-3 higher than predicted by Haardt & Madau (2012).

Lyman-alpha forest data probing the post-reionization Universe shows surprisingly large opacity fluctuations over rather large ($\ge$50 comoving Mpc/h) spatial scales. We model these fluctuations using a hybrid approach utilizing the large volume Millennium simulation to predict the spatial distribution of QSOs combined with smaller scale full hydrodynamical simulation performed with RAMSES and post-processed with the radiative transfer code ATON. We produce realictic mock absorption spectra that account for the contribution of galaxies and QSOs to the ionising UV background. This improved models confirm our earlier findings that a significant ($\ge$50%) contribution of ionising photons from QSOs can explain the large reported opacity fluctuations on large scales. The inferred QSO luminosity function is thereby consistent with recent estimates of the space density of QSOs at this redshift. Our simulations still somewhat struggle, however, to reproduce the very long (110 comoving Mpc/h) high opacity absorption through observed in ULAS J0148+0600, perhaps suggesting an even later end of reionization than assumed in our previously favoured model. Medium-deep/medium area QSO surveys as well as targeted searches for the predicted strong transverse QSO proximity effect whould illuminate the origin of the observed large scale opacity fluctuations. They would allow to substantiate whether UV fluctuations due to QSO are indeed primarily responsible, or whether significant contributions from other recently proposed mechansims such as large scale fluctuations in temperature and mean free path (even in the absence of rare bright sources) are required.

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.

Reionisation in the early Universe is likely driven by dwarf galaxies. Using cosmological radiation-hydrodynamic simulations, we study star formation and the escape of Lyman continuum (LyC) photons from mini-haloes with $M_{\rm halo} \le 10^8\,M_\odot$. Our simulations include a new thermo-turbulent star formation model, non-equilibrium chemistry, and relevant stellar feedback processes (photoionisation by young massive stars, radiation pressure, and mechanical supernova explosions). We find that feedback reduces star formation very efficiently in mini-haloes, resulting in the stellar mass consistent with the slope and normalisation reported in Kimm \& Cen and the empirical stellar mass-to-halo mass relation derived in the local Universe. Because star formation is stochastic and dominated by a few gas clumps, the escape fraction in mini-haloes is generally determined by radiation feedback (heating due to photo-ionisation), rather than supernova explosions. We also find that the photon number-weighted mean escape fraction in mini-haloes is higher ($\sim20$-$40\%$) than that in atomic-cooling haloes, although the instantaneous fraction in individual haloes varies significantly. The escape fraction from Pop III stars is found to be significant ($\ge10\%$) only when the mass is greater than $\sim$100\,\msun. Based on simple analytic calculations, we show that LyC photons from mini-haloes are, despite their high escape fractions, of minor importance for reionisation due to inefficient star formation. We confirm previous claims that stars in atomic-cooling haloes with masses $10^8\,M_\odot\le M_{\rm halo} \le 10^{11}\,M_\odot$ are likely to be the most important source of reionisation.

At low densities the standard ionisation history of the intergalactic medium (IGM) predicts a decreasing temperature of the IGM with decreasing density once hydrogen (and helium) reionisation is complete. Heating the high-redshift, low-density IGM above the temperature expected from photo-heating is difficult, and previous claims of high/rising temperatures in low density regions of the Universe based on the probability density function (PDF) of the opacity in Lyman-$\alpha$ forest data at $2<z<4$ have been met with considerable scepticism, particularly since they appear to be in tension with other constraints on the temperature-density relation (TDR). We utilize here an ultra-high signal-to-noise spectrum of the QSO HE0940-1050 and a novel technique to study the low opacity part of the PDF. We show that there is indeed evidence (at 90% confidence level) that a significant volume fraction of the under-dense regions at $z \sim 3$ has temperatures as high or higher than those at densities comparable to the mean and above. We further demonstrate that this conclusion is nevertheless consistent with measurements of a slope of the TDR in over-dense regions that imply a decreasing temperature with decreasing density, as expected if photo-heating of ionised hydrogen is the dominant heating process. We briefly discuss implications of our findings for the need to invoke either spatial temperature fluctuations, as expected during helium reionization, or additional processes that heat a significant volume fraction of the low-density IGM.

We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high resolution, high signal-to-noise observations of the Lyman-alpha forest over the redshift range 2 < z < 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Lyman-alpha forest is a natural consequence of the web-like distribution of matter arising in LCDM cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Lyman-alpha absorption lines with column densities N_HI > 10^14.5 cm^-2 at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z=2.7 and a flatter temperature-density relation is required at z=2.4, consistent with the expected effects of non-equilibrium ionisation during He-II reionisation.

In this work, we investigate the abundance and distribution of metals in the intergalactic medium (IGM) at $\langle z \rangle \simeq 2.8$ through the analysis of an ultra-high signal-to-noise ratio UVES spectrum of the quasar HE0940-1050. In the CIV forest, our deep spectrum is sensitive at $3\,\sigma$ to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60 per cent of the considered redshift range down to $\simeq11.1$. In our sample, all HI lines with $\log N_{\rm HI} \ge 14.8$ show an associated CIV absorption. In the range $14.0 \le \log N_{\rm HI} <14.8$, 43 per cent of HI lines has an associated CIV absorption. At $\log N_{\rm HI} < 14.0$, the detection rates drop to $<10$ per cent, possibly due to our sensitivity limits and not to an actual variation of the gas abundance properties. In the range $\log N_{\rm HI} \ge 14$, we observe a fraction of HI lines with detected CIV a factor of 2 larger than the fraction of HI lines lying in the circum-galactic medium (CGM) of relatively bright Lyman-break galaxies hosted by dark matter haloes with $\langle M\rangle \sim10^{12}$ M$_{\odot}$. The comparison of our results with the output of a grid of photoionization models and of two cosmological simulations implies that the volume filling factor of the IGM gas enriched to a metallicity $\log Z/Z_{\odot} \ge -3$ should be of the order of $\sim 10-13$ percent. In conclusion, our results favour a scenario in which metals are found also outside the CGM of bright star-forming galaxies, possibly due to pollution by lower mass objects and/or to an early enrichment by the first sources.

We present here 21 cm predictions from high dynamic range simulations for a range of reionization histories that have been tested against available Lyman-alpha and CMB data. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen spin temperature is significantly larger than the CMB temperature. Models consistent with the available Lyman-alpha data and CMB measurement of the Thomson optical depth predict typical values of 10--20 mK^2 for the variance of the 21 cm brightness temperature at redshifts z=7--10 at scales accessible to ongoing and upcoming experiments (k < 1 h/cMpc). This is only a factor of a few below the sensitivity claimed to have been already reached by ongoing experiments in the signal rms value. Our different models for the reionization history make markedly different predictions for the redshift evolution and thus frequency dependence of the 21 cm power spectrum and should be easily discernible by LOFAR (and later HERA and SKA1) at their design sensitivity. Our simulations have sufficient resolution to assess the effect of high-density Lyman limit systems that can self-shield against ionizing radiation and stay 21 cm bright even if the hydrogen in their surroundings is highly ionized. Our simulations predict that including the effect of the self-shielded gas in highly ionized regions reduces the large scale 21 cm power by about 30%.

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

We present models of low- and high-ionization metal-line absorbers (O I, C II, C IV and Mg II) during the end of the reionization epoch, at z ~ 6. Using four cosmological hydrodynamical simulations with different feedback schemes (including the Illustris and Sherwood simulations) and two different choices of hydro-solver, we investigate how the overall incidence rate and equivalent width distribution of metal-line absorbers varies with the galactic wind prescription. We find that the O I and C II absorbers are reasonably insensitive to the feedback scheme. All models, however, struggle to reproduce the observations of C IV and Mg II, which are probing down to lower overdensities than O I and C II at z ~ 6, suggesting that the metals in the simulations are not being transported out into the IGM efficiently enough. The situation is improved but not resolved if we choose a harder (but still reasonable) and/or (locally) increased UV background at z ~ 6.

QSO near-zones are an important probe of the the ionization state of the IGM at z ~ 6-7, at the end of reionization. We present here high-resolution cosmological 3D radiative transfer simulations of QSO environments for a wide range of host halo masses, 10$^{10-12.5}$ M_sun. Our simulated near-zones reproduce both the overall decrease of observed near-zone sizes at 6 < z < 7 and their scatter. The observable near-zone properties in our simulations depend only very weakly on the mass of the host halo. The size of the H II region expanding into the IGM is generally limited by (super-)Lyman Limit systems loosely associated with (low-mass) dark matter haloes. This leads to a strong dependence of near-zone size on direction and drives the large observed scatter. In the simulation centred on our most massive host halo, many sightlines show strong red damping wings even for initial volume averaged neutral hydrogen fractions as low as ~ 10$^{-3}$. For QSO lifetimes long enough to allow growth of the central supermassive black hole while optically bright, we can reproduce the observed near-zone of ULAS J1120+0641 only with an IGM that is initially neutral. Our results suggest that larger samples of z > 7 QSOs will provide important constraints on the evolution of the neutral hydrogen fraction and thus on how late reionization ends.

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

We study how runaway stellar collisions in high-redshift, metal-poor star clusters form very massive stars (VMSs) that can directly collapse to intermediate-mass black holes (IMBHs). We follow the evolution of a pair of neighbouring high-redshift mini-haloes with high-resolution, cosmological hydrodynamical zoom-in simulations using the adaptive mesh refinement code RAMSES combined with the non-equilibrium chemistry package KROME. The first collapsing mini-halo is assumed to enrich the central nuclear star cluster (NSC) of the other to a critical metallicity, sufficient for Population II (Pop. II) star formation at redshift $z\approx27$. Using the spatial configuration of the flattened, asymmetrical gas cloud forming in the core of the metal enriched halo, we set the initial conditions for simulations of an initially non-spherical star cluster with the direct summation code NBODY6 which are compared to about 2000 NBODY6 simulations of spherical star clusters for a wide range of star cluster parameters. The final mass of the VMS that forms depends strongly on the initial mass and initial central density of the NSC. For the initial central densities suggested by our RAMSES simulations, VMSs with mass > 400 M$_{\odot}$ can form in clusters with stellar masses of $\approx10^4$ M$_{\odot}$, and this can increase to well over 1000 M$_{\odot}$ for more massive and denser clusters. The high probability we find for forming a VMS in these mini-haloes at such an early cosmic time makes collisional runaway of Pop. II star clusters a promising channel for producing large numbers of high-redshift IMBHs that may act as the seeds of supermassive black holes.

We discuss the recent Baryon Oscillation Spectroscopic Survey measurement of a rather high bias factor for the host galaxies/haloes of Damped Lyman-alpha Absorbers (DLAs), in the context of our previous modelling of the physical properties of DLAs within the $\Lambda$ cold dark matter paradigm. Joint modelling of the column density distribution, the velocity width distribution of associated low ionization metal absorption, and the bias parameter suggests that DLAs are hosted by galaxies with dark matter halo masses in the range $10 < \log M_v < 12$, with a rather sharp cutoff at the lower mass end, corresponding to virial velocities of 35 km/sec. The observed properties of DLAs appear to suggest efficient (stellar) feedback in haloes with masses/virial velocities below the cutoff and a large retained baryon fraction (> 35 %) in haloes above the cutoff.

We investigate the origin of extragalactic continuum emission and its relation to the stellar population of a recently discovered peculiar z=3.344 Lyman alpha emitter. Based on an analysis of the broad-band colors and morphology we find further support for the idea that the underlying galaxy is being fed by a large-scale (L > 35 kpc) accretion stream. Archival HST images show small scale (~5 kpc) tentacular filaments converging near a hot-spot of star-formation, possibly fueled by gas falling in along the filaments. The spectral energy distribution of the tentacles is broadly compatible with either (1) non-ionizing rest-frame far-UV continuum emission from stars formed in an 60 million-year-old starburst; (2) nebular 2-photon-continuum radiation, arising from collisional excitation cooling, or (3) a recombination spectrum emitted by hydrogen fluorescing in response to ionizing radiation escaping from the galaxy. The latter possibility simultaneously accounts for the presence of asymmetric Lyman alpha emission from the large-scale gaseous filament and the nebular continuum in the smaller-scale tentacles as caused by the escape of ionizing radiation from the galaxy. Possible astrophysical explanations for the nature of the tentacles include: a galactic wind powered by the starburst; in-falling gas during cold accretion, or tails of interstellar medium dragged out of the galaxy by satellite halos that have plunged through the main halo. The possibility of detecting extragalactic 2-photon continuum emission in space-based, broad-band images suggests a tool for studying the gaseous environment of high redshift galaxies at much greater spatial detail than possible with Lyman alpha or other resonance line emission.

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called "rVINE", aims to significantly improve the accuracy of close encounters of massive bodies with supermassive black holes in galaxy-scale numerical simulations. We demonstrate the capabilities of the code by studying two test problems, the sinking of a single massive black hole to the centre of a gas-free galaxy due to dynamical friction and the hardening of a supermassive black hole binary due to close stellar encounters. We show that results obtained with rVINE compare well with NBODY7 for problems with particle numbers that can be simulated with NBODY7. In particular, in both NBODY7 and rVINE we find a clear N-dependence of the binary hardening rate, a low binary eccentricity and moderate eccentricity evolution, as well as the conversion of the galaxy's inner density profile from a cusp to a a core via the ejection of stars at high velocity. The much larger number of particles that can be handled by rVINE will open up exciting opportunities to model stellar dynamics close to SMBHs much more accurately in a realistic galactic context. This will help to remedy the inherent limitations of commonly used tree solvers to follow the correct dynamical evolution of black holes in galaxy scale simulations.

We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibility of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.

Observations of the emission from spatially extended cold gas around bright high-redshift QSOs reveal surprisingly large velocity widths exceeding 2000 km s^(-1), out to projected distances as large as 30 kpc. The high velocity widths have been interpreted as the signature of powerful AGN-driven outflows. Naively, these findings appear in tension with hydrodynamic models in which AGN-driven outflows are energy-driven and thus very hot with typical temperatures T = 10^6-7 K. Using the moving-mesh code Arepo, we perform 'zoom-in' cosmological simulations of a z = 6 QSO and its environment, following black hole growth and feedback via energy-driven outflows. In the simulations, the QSO host galaxy is surrounded by a clumpy circum-galactic medium pre-enriched with metals due to supernovae-driven galactic outflows. As a result, part of the AGN-driven hot outflowing gas can cool radiatively, leading to large amounts (> 10^9 M_sun) of cold gas comoving with the hot bipolar outflow. This results in velocity widths of spatially extended cold gas similar to those observed. We caution, however, that gas inflows, random motions in the deep potential well of the QSO host galaxy and cooling of supernovae-driven winds contribute significantly to the large velocity width of the cold gas in the simulations, complicating the interpretation of observational data.

We employ hydrodynamical simulations using the moving-mesh code AREPO to investigate the role of energy and momentum input from Active Galactic Nuclei (AGN) in driving large-scale galactic outflows. We start by reproducing analytic solutions for both energy- and momentum-driven outflowing shells in simulations of a spherical isolated dark matter potential with gas in hydrostatic equilibrium and with no radiative cooling. We confirm that for this simplified setup, galactic outflows driven by a momentum input rate of order L_Edd/c can establish an M_BH - sigma relation with slope and normalisation similar to that observed. We show that momentum input at a rate of L_Edd/c is however insufficient to drive efficient outflows once cooling and gas inflows as predicted by cosmological simulations at resolved scales are taken into account. We argue that observed large-scale AGN-driven outflows are instead likely to be energy-driven and show that such outflows can reach momentum fluxes exceeding 10 L_Edd/c within the innermost 10 kpc of the galaxy. The outflows are highly anisotropic, with outflow rates and a velocity structure found to be inadequately described by spherical outflow models. We verify that the hot energy-driven outflowing gas is expected to be strongly affected by metal-line cooling, leading to significant amounts (>10^9 M_sun) of entrained cold gas.

We employ cosmological hydrodynamical simulations to investigate models in which the supermassive black holes (BHs) powering luminous z ~ 6 QSOs grow from massive seeds. We simulate at high resolution 18 fields sampling regions with densities ranging from the mean cosmic density all the way to the highest sigma peaks in the Millennium simulation volume. Only in the most massive halos, BHs can grow to masses up to ~ 10^9 Msun by z ~ 6 without invoking super-Eddington accretion. Accretion onto the most massive BHs becomes limited by thermal AGN feedback by z ~ 9-8 with further BH growth proceeding in short Eddington limited bursts. Our modelling suggests that current flux-limited surveys of QSOs at high redshift preferentially detect objects at their peak luminosity and therefore miss a substantial population of QSOs powered by similarly massive BHs but with low accretion rates. To test whether the required host halo masses are consistent with the observed galaxy environment of z ~ 6 QSOs, we produce realistic rest-frame UV images of our simulated galaxies. Without strong stellar feedback, our simulations predict numbers of bright galaxies larger than observed by a factor ten or more. Supernova-driven galactic winds reduce the predicted numbers to a level consistent with observations indicating that stellar feedback was already very efficient at high redshifts. We have further investigated the effect of thermal AGN feedback on the surrounding gas. Our adopted AGN feedback prescription drives mostly energy-driven highly anisotropic outflows with gas speeds of >= 1000 km/s to distances of >= 10 kpc consistent with observations. The spatially extended thermal X-ray emission around bright QSOs powered by these outflows can exceed by large factors the emission expected without AGN feedback and is an important diagnostic of the mechanism whereby AGN feedback energy couples to surrounding gas.

We examine the kinematic structure of Damped Lyman-alpha Systems (DLAs) in a series of cosmological hydrodynamic simulations using the AREPO code. We are able to match the distribution of velocity widths of associated low ionisation metal absorbers substantially better than earlier work. Our simulations produce a population of DLAs dominated by halos with virial velocities around 70 km/s, consistent with a picture of relatively small, faint objects. In addition, we reproduce the observed correlation between velocity width and metallicity and the equivalent width distribution of SiII. Some discrepancies of moderate statistical significance remain; too many of our spectra show absorption concentrated at the edge of the profile and there are slight differences in the exact shape of the velocity width distribution. We show that the improvement over previous work is mostly due to our strong feedback from star formation and our detailed modelling of the metal ionisation state.

We examine the abundance, clustering and metallicity of Damped Lyman-alpha Absorbers (DLAs) in a suite of hydrodynamic cosmological simulations using the moving mesh code AREPO. We incorporate models of supernova and AGN feedback, as well as molecular hydrogen formation. We compare our simulations to the column density distribution function at $z=3$, the total DLA abundance at $z=2-4$, the measured DLA bias at $z=2.3$ and the DLA metallicity distribution at $z=2-4$. Our preferred models produce populations of DLAs in good agreement with most of these observations. The exception is the DLA abundance at $z < 3$, which we show requires stronger feedback in $10^{11-12} \, h^{-1} M_\odot$ mass halos. While the DLA population probes a wide range of halo masses, we find the cross-section is dominated by halos of mass $10^{10} - 10^{11} \, h^{-1} M_\odot$ and virial velocities $50 - 100 \;\mathrm{km/s}$. The simulated DLA population has a linear theory bias of $1.7$, whereas the observations require $2.17 \pm 0.2$. We show that non-linear growth increases the bias in our simulations to $2.3$ at $k=1\; \mathrm{Mpc/}h$, the smallest scale observed. The scale-dependence of the bias is, however, very different in the simulations compared against the observations. We show that, of the observations we consider, the DLA abundance and column density function provide the strongest constraints on the feedback model.

We perform an extensive suite of N-body simulations of the matter power spectrum, incorporating massive neutrinos in the range M = 0.15-0.6 eV, probing the non-linear regime at scales k < 10 hMpc-1 at z < 3. We extend the widely used HALOFIT approximation to account for the effect of massive neutrinos on the power spectrum. In the strongly non-linear regime HALOFIT systematically over-predicts the suppression due to the free-streaming of the neutrinos. The maximal discrepancy occurs at k ~ 1 hMpc-1, and is at the level of 10% of the total suppression. Most published constraints on neutrino masses based on HALOFIT are not affected, as they rely on data probing the matter power spectrum in the linear or mildly non-linear regime. However, predictions for future galaxy, Lyman-alpha forest and weak lensing surveys extending to more non-linear scales will benefit from the improved approximation to the non-linear matter power spectrum we provide. Our approximation reproduces the induced neutrino suppression over the targeted scales and redshifts significantly better. We test its robustness with regard to changing cosmological parameters and a variety of modelling effects.

We have performed high-resolution numerical simulations with the hydrodynamical AMR code Enzo to investigate the formation of massive seed black holes in a sample of six dark matter haloes above the atomic cooling threshold. The aim of this study is to illustrate the effects of varying the maximum refinement level on the final object formed. The virial temperatures of the simulated haloes range from $\rm{T} \sim 10000\ \rm{K} - 16000\ \rm{K}$ and they have virial masses in the range $\rm{M} \sim 2 \times 10^7 \rm{M_{\odot}}$ to $\rm{M} \sim 7 \times 10^7 \rm{M_{\odot}}$ at $z \sim 15$. The outcome of our six fiducial simulations is both generic and robust. A rotationally supported, marginally gravitationally stable, disk forms with an exponential profile. The mass and scale length of this disk depends strongly on the maximum refinement level used. Varying the maximum refinement level by factors between 1 / 64 to 256 times the fiducial level illustrates the care that must be taken in interpreting the results. The lower resolution simulations show tentative evidence that the gas may become rotationally supported out to 20 pc while the highest resolution simulations show only weak evidence of rotational support due to the shorter dynamical times for which the simulation runs. The higher resolution simulations do, however, point to fragmentation at small scales of the order of $\sim 100$ AU. In the highest resolution simulations a central object of a few times $10^2\ \rm{M_{\odot}}$ forms with multiple strongly bound, Jeans unstable, clumps of $\sim 10\ \rm{M_{\odot}}$ and radii of 10 - 20 AU suggesting the formation of dense star clusters in these haloes.

We present new measurements of the thermal state of the intergalactic medium (IGM) at $z\sim2.4$ derived from absorption line profiles in the Ly$\alpha$ forest. We use a large set of high-resolution hydrodynamical simulations to calibrate the relationship between the temperature-density ($T$--$\Delta$) relation in the IGM and the distribution of HI column densities, $N_{\rm HI}$, and velocity widths, $b_{\rm HI}$, of discrete Ly$\alpha$ forest absorbers. This calibration is then applied to the measurement of the lower cut-off of the $b_{\rm HI}$--$N_{\rm HI}$ distribution recently presented by Rudie et al. (2012). We infer a power-law $T$--$\Delta$ relation, $T=T_{0}\Delta^{\gamma-1}$, with a temperature at mean density, $T_{0}=[1.00^{+0.32}_{-0.21}]\times10^{4}\rm\,K$ and slope $(\gamma-1)=0.54\pm0.11$. The slope is fully consistent with that advocated by the analysis of Rudie et al (2012); however, the temperature at mean density is lower by almost a factor of two, primarily due to an adjustment in the relationship between column density and physical density assumed by these authors. These new results are in excellent agreement with the recent temperature measurements of Becker et al. (2011), based on the curvature of the transmitted flux in the Ly$\alpha$ forest. This suggests that the thermal state of the IGM at this redshift is reasonably well characterised over the range of densities probed by these methods.

Low ionization metal absorption due to OI has been identified as an important probe of the physical state of the inter-/circumgalactic medium at the tail-end of reionization. We use here high-resolution hydrodynamic simulations to interpret the incidence rate of OI absorbers at z~6 as observed by Becker et al. (2011). We infer weak OI absorbers (EW > 0.1 A) to have typical HI column densities in the range of sub-DLAs, densities of 80 times the mean baryonic density and metallicities of about 1/500 th solar. This is similar to the metallicity inferred at similar overdensities at z~3, suggesting that the metal enrichment of the circumgalactic medium around low-mass galaxies has already progressed considerably by z~6. The apparently rapid evolution of the incidence rates for OI absorption over the redshift range 5<z<6 mirrors that of self-shielded Lyman-Limit systems at lower redshift and is mainly due to the rapid decrease of the meta-galactic photo-ionization rate at z>5. We predict the incidence rate of OI absorbers to continue to rise rapidly with increasing redshift as the IGM becomes more neutral. If the distribution of metals extends to lower density regions, OI absorbers will allow the metal enrichment of the increasingly neutral filamentary structures of the cosmic web to be probed.

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

We present updated constraints on the free-streaming of warm dark matter (WDM) particles derived from an analysis of the Lya flux power spectrum measured from high-resolution spectra of 25 z > 4 quasars obtained with the Keck High Resolution Echelle Spectrometer (HIRES) and the Magellan Inamori Kyocera Echelle (MIKE) spectrograph. We utilize a new suite of high-resolution hydrodynamical simulations that explore WDM masses of 1, 2 and 4 keV (assuming the WDM consists of thermal relics), along with different physically motivated thermal histories. We carefully address different sources of systematic error that may affect our final results and perform an analysis of the Lya flux power with conservative error estimates. By using a method that samples the multi-dimensional astrophysical and cosmological parameter space, we obtain a lower limit mwdm > 3.3 keV (2sigma) for warm dark matter particles in the form of early decoupled thermal relics. Adding the Sloan Digital Sky Survey (SDSS) Lya flux power spectrum does not improve this limit. Thermal relics of masses 1 keV, 2 keV and 2.5 keV are disfavoured by the data at about the 9sigma, 4sigma and 3sigma C.L., respectively. Our analysis disfavours WDM models where there is a suppression in the linear matter power spectrum at (non-linear) scales corresponding to k=10h/Mpc which deviates more than 10% from a LCDM model. Given this limit, the corresponding "free-streaming mass" below which the mass function may be suppressed is 2x10^8 Msun/h. There is thus very little room for a contribution of the free-streaming of WDM to the solution of what has been termed the small scale crisis of cold dark matter.

A spectroscopically detected Lyman alpha emitting halo at redshift 3.216 in the GOODS-N field is found to reside at the convergence of several Lyman alpha filaments. HST images show that some of the filaments are inhabited by galaxies. Several of the galaxies in the field have pronounced head-tail structures, which are partly aligned with each other. The blue colors of most tails suggest the presence of young stars, with the emission from at least one of the galaxies apparently dominated by high equivalent width Lyman alpha. Faint, more diffuse, and similarly elongated, apparently stellar features, can be seen over an area with a linear extent of at least 90 kpc. The region within several arcseconds of the brightest galaxy exhibits spatially extended emission by HeII, NV and various lower ionization metal lines. The gas-dynamical features present are strongly reminiscent of ram-pressure stripped galaxies, including evidence for recent star formation in the stripped contrails. Spatial gradients in the appearance of several galaxies may represent a stream of galaxies passing from a colder to a hotter intergalactic medium. The stripping of gas from the in-falling galaxies, in conjunction with the occurrence of star formation and stellar feedback in the galactic contrails suggests a mechanism for the metal enrichment of the high redshift intergalactic medium that does not depend on long-range galactic winds, at the same time opening a path for the escape of ionizing radiation from galaxies.

In this third paper in a series on the nature of extended, asymmetric Lyman alpha emitters at z ~ 3 we report the discovery, in an ultra-deep, blind, spectroscopic long-slit survey, of a Lyman alpha emitting halo around a QSO at redshift 3.045. The QSO is a previously known, obscured AGN. The halo appears extended along the direction of the slit and exhibits two faint patches separated by 17 proper kpc in projection from the QSO. Comparison of the 2-dimensional spectrum with archival HST ACS images shows that these patches coincide spatially with emission from a peculiar, dumbbell-shaped, faint galaxy. The assumptions that the Lyman alpha emission patches are originating in the galaxy and that the galaxy is physically related to the QSO are at variance with photometric estimates of the galaxy redshift. We show, however, that a population of very young stars at the redshift of the QSO may fit the existing rest frame broad band UV photometry of the galaxy. If this scenario is correct then the symmetry of the galaxy in continuum and Lyman alpha emission, the extension of the QSO's Lyman alpha emission in its direction, and the likely presence of a young stellar population in close proximity to a (short-lived) AGN suggest that this may be an example of AGN feedback triggering external star formation in high redshift galaxies.

We study the X-ray emission from NGC1277, a galaxy in the core of the Perseus cluster, for which van den Bosch et al. have recently claimed the presence of an UltraMassive Black Hole (UMBH) of mass 1.7 times 10^10 Msun, unless the IMF of the stars in the stellar bulge is extremely bottom heavy. The X-rays originate in a power-law component of luminosity 1.3 times 10^40 erg/s embedded in a 1 keV thermal minicorona which has a half-light radius of about 360 pc, typical of many early-type galaxies in rich clusters of galaxies. If Bondi accretion operated onto the UMBH from the minicorona with a radiative efficiency of 10 per cent, then the object would appear as a quasar with luminosity 10^46 erg/s, a factor of almost 10^6 times higher than observed. The accretion flow must be highly radiatively inefficient, similar to past results on M87 and NGC3115. The UMBH in NGC1277 is definitely not undergoing any significant growth at the present epoch. We note that there are 3 UMBH candidates in the Perseus cluster and that the inferred present mean mass density in UMBH could be 10^5 Msun/Mpc^3, which is 20 to 30 per cent of the estimated mean mass density of all black holes. We speculate on the implied growth of UMBH and their hosts, and discuss the possibiity that extreme AGN feedback could make all UMBH host galaxies have low stellar masses at redshifts around 3. Only those which end up at the centres of groups and clusters later accrete large stellar envelopes and become Brightest Cluster Galaxies. NGC1277 and the other Perseus core UMBH, NGC1270, have not however been able to gather more stars or gas owing to their rapid orbital motion in the cluster core.

A variety of independent observational studies have now reported a significant decline in the fraction of Lyman-break galaxies which exhibit Ly-a emission over the redshift interval z=6-7. In combination with the strong damping wing extending redward of Ly-a in the spectrum of the bright z=7.085 quasar ULAS 1120+0641, this has strengthened suggestions that the hydrogen in the intergalactic medium (IGM) is still substantially neutral at z~7. Current theoretical models imply HI fractions as large as 40-90 per cent may be required to explain these data assuming there is no intrinsic evolution in the Ly-a emitter population. We propose that such large neutral fractions are not necessary. Based on a hydrodynamical simulation which reproduces the absorption spectra of high-redshift (z~6-7) quasars, we demonstrate that the opacity of the intervening IGM redward of rest-frame Ly-a can rise rapidly in average regions of the Universe simply because of the increasing incidence of absorption systems which are optically thick to Lyman continuum photons as the tail-end of reionisation is approached. Our simulations suggest these data do not require a large change in the IGM neutral fraction by several tens of per cent from z=6-7, but may instead be indicative of the rapid decrease in the typical mean free path for ionising photons expected during the final stages of reionisation.

Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z=6). Using hydrodynamical simulations of the formation of a massive M_vir=5 10^11 M_sun halo by z=6 (the most massive progenitor of a cluster of M_vir=2 10^15 M_sun at z=0), we evaluate the impact of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and gas distribution inside this massive halo. We find that SN feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.

We present forecasts for the accuracy of determining the parameters of a minimal cosmological model and the total neutrino mass based on combined mock data for a future Euclid-like galaxy survey and Planck. We consider two different galaxy surveys: a spectroscopic redshift survey and a cosmic shear survey. We make use of the Monte Carlo Markov Chains (MCMC) technique and assume two sets of theoretical errors. The first error is meant to account for uncertainties in the modelling of the effect of neutrinos on the non-linear galaxy power spectrum and we assume this error to be fully correlated in Fourier space. The second error is meant to parametrize the overall residual uncertainties in modelling the non-linear galaxy power spectrum at small scales, and is conservatively assumed to be uncorrelated and to increase with the ratio of a given scale to the scale of non-linearity. It hence increases with wavenumber and decreases with redshift. With these two assumptions for the errors and assuming further conservatively that the uncorrelated error rises above 2% at k = 0.4 h/Mpc and z = 0.5, we find that a future Euclid-like cosmic shear/galaxy survey achieves a 1-sigma error on Mnu close to 32 meV/25 meV, sufficient for detecting the total neutrino mass with good significance. If the residual uncorrelated errors indeed rises rapidly towards smaller scales in the non-linear regime as we have assumed here then the data on non-linear scales does not increase the sensitivity to the total neutrino mass. Assuming instead a ten times smaller theoretical error with the same scale dependence, the error on the total neutrino mass decreases moderately from sigma(Mnu) = 18 meV to 14 meV when mildly non-linear scales with 0.1 h/Mpc < k < 0.6 h/Mpc are included in the analysis of the galaxy survey data.

A deep narrow-band survey for Ly-alpha emission carried out on the VLT-FORS2 has revealed 98 Ly-alpha candidates down to a flux limit of 4.e-18 erg/s/cm^2 in a volume of 5500 comoving Mpc^3 at z=2.4 centered on the hyperluminous quasar HE0109-3518. The properties of the detected sources in terms of their i) equivalent width distribution, ii) luminosity function, and iii) the average luminosity versus projected distance from the quasar, all suggest that a large fraction of these objects have been fluorescently "illuminated" by HE0109-3518. This conclusion is supported by comparison with detailed radiative transfer simulations of the effects of the quasar illumination. 18 objects have a rest-frame Equivalent Width (EW0) larger than 240A, the expected limit for Ly-alpha emission powered by Population II star formation and 12 sources among these do not have any continuum counterpart in a deep V-band imaging of the same field. For these, a stacking analysis indicates EW0>800A, effectively ruling out Ly-alpha powered by internal star formation. These sources are thus the best candidates so far for proto-galactic clouds or "dark" galaxies at high-redshift, whose existence has recently been suggested by several theoretical studies. Assuming they are mostly ionized by the quasar radiation, we estimate that their gas masses would be about 10^9 Msun implying that their star formation efficiencies (SFE) are less than 1.e-11 yr^-1 several times below the SFE of the most gas-rich dwarf galaxies locally, and five hundred times lower than typical massive star-forming galaxies at z~2. We have also discovered extended, filamentary gas, also likely illuminated by the quasar, around some of the brightest continuum-detected sources with EW0>240A. This emission is compatible with the expectations for circum-galactic cold streams but other origins, including tidal stripping, are also possible.

We report the observation of a further asymmetric, extended Lyman alpha emitting halo at z=2.63, from our ultra-deep, long-slit spectroscopic survey of faint high redshift emitters, undertaken with Magellan LDSS3 in the GOODS-S field. The Lya emission, detected over more than 30 kpc, is spatially coincident with a concentration of galaxies visible in deep broad-band imaging. While these faint galaxies without spectroscopic redshifts cannot with certainty be associated with one another or with the Lya emission, there are a number of compelling reasons why they very probably form a Milky Way halo-mass group at the Lya redshift. A filamentary structure, possibly consisting of Lya emission at very high equivalent width, and evidence for disturbed stellar populations, suggest that the properties of the emitting region reflect ongoing galaxy assembly, with recent galaxy mergers and star formation occurring in the group. Hence, the Lya provides unique insights into what is probably a key mode of galaxy formation at high redshifts. The Lya emission may be powered by cooling radiation or spatially extended star formation in the halo, but is unlikely to be fluorescence driven by either an AGN or one of the galaxies. The spatial profile of the emission is conspicuously different from that of typical Lya emitters or Lyman break galaxies, which is consistent with the picture that extended emission of this kind represents a different stage in the galaxy formation process. Faint, extended Lya emitters such as these may be lower-mass analogues of the brighter Lya blobs. Our observations provide further, circumstantial evidence that galaxy mergers may promote the production and / or escape of ionizing radiation, and that the halos of interacting galaxies may be significant sources for ionizing photons during the epoch of reionization.

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

We use the galaxy angular power spectrum at $z\sim0.5-1.2$ from the Canada-France-Hawaii-Telescope Legacy Survey Wide fields (CFHTLS-Wide) to constrain separately the total neutrino mass $\sum{m_\nu}$ and the effective number of neutrino species $N_{\rm{eff}}$. This survey has recently benefited from an accurate calibration of the redshift distribution, allowing new measurements of the (non-linear) matter power spectrum in a unique range of scales and redshifts sensitive to neutrino free streaming. Our analysis makes use of a recent model for the effect of neutrinos on the weakly non-linear matter power spectrum derived from accurate N-body simulations. We show that CFHTLS, combined with WMAP7 and a prior on the Hubble constant provides an upper limit of $\sum{m_\nu}<0.29\,$eV and $N_{\rm{eff}} =4.17^{+1.62}_{-1.26}$ (2$\,\sigma$ confidence levels). If we omit smaller scales which may be affected by non-linearities, these constraints become $\sum{m_\nu}<0.41\,$eV and $N_{\rm{eff}} =3.98^{+2.02}_{-1.20}$ (2$\,\sigma$ confidence levels). Finally we show that the addition of other large scale structures probes can further improve these constraints, demonstrating that high redshift large volumes surveys such as CFHTLS are complementary to other cosmological probes of the neutrino mass.

We present measurements of the intergalactic medium (IGM) temperature within ~5 proper Mpc of seven luminous quasars at z~6. The constraints are obtained from the Doppler widths of Lyman-alpha absorption lines in the quasar near-zones and build upon our previous measurement for the z=6.02 quasar SDSS J0818+1722. The expanded data set, combined with an improved treatment of systematic uncertainties, yields an average temperature at the mean density of log (T_0/K) = 4.21 \pm 0.03 (\pm^0.06_0.07) at 68 (95) per cent confidence for a flat prior distribution over 3.2 < log (T_0/K) < 4.8. In comparison, temperatures measured from the general IGM at z~5 are ~0.3 dex cooler, implying an additional source of heating around these quasars which is not yet present in the general IGM at slightly lower redshift. This heating is most likely due to the recent reionisation of He-II in vicinity of these quasars, which have hard, non-thermal ionising spectra. The elevated temperatures may therefore represent evidence for the earliest stages of He-II reionisation in the most biased regions of the high-redshift Universe. The temperature as a function of distance from the quasars is consistent with being constant, log(T_0/K)~4.2, with no evidence for a line-of-sight thermal proximity effect. However, the limited extent of the quasar near-zones prevents the detection of He-III regions larger than ~5 proper Mpc. Under the assumption the quasars have reionised the He-II in their vicinity, we infer the data are consistent with an average optically bright phase of duration in excess of 10^6.5 yr. These measurements represent the highest-redshift IGM temperature constraints to date, and thus provide a valuable data set for confronting models of H-I reionisation.

The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.

We present observations of a peculiar Lyman alpha-emitting galaxy at redshift 3.344, discovered in a deep, blind spectroscopic survey for faint Lyman alpha emitters with the Magellan II telescope in the Hubble Ultra Deep Field (HUDF). The galaxy exhibits complex Lyman alpha emission, including an extended, asymmetric component that is partially suppressed by damped Lyman alpha absorption, and two spatially elongated, narrow emission features. Archival HST ACS imaging shows evidence for tidal disruption of the stellar component. This V=27 galaxy appears to give us unprecedented insights into two fundamental stages in the formation of structure at high redshift: the inflow of gas into ordinary galaxies, and the escape of ionizing radiation into the intergalactic medium. Neutral hydrogen, falling in partly in form of a narrow filament, appears to emit fluorescent Lyman alpha photons induced by the stellar ionizing flux escaping from the disturbed galaxy. The in-falling material may represent primary cold accretion or an interaction-triggered inflow. The rate of ionizing photons required by the observed Lyman alpha emission is consistent with the rate of photons produced by the observed stellar population, with roughly 50 percent of ionizing photons escaping from the immediate galaxy and encountering the in-falling gas. The observational properties of the galaxy lend support to a picture where galaxy interactions facilitate the escape of both Lyman alpha and ionizing radiation. We argue that galaxies like the present object may be common at high redshift. This galaxy may therefore be a late example of an interacting population of dwarf galaxies contributing significantly to the reionization of the universe.

We study the evolution of gravitationally recoiled supermassive black holes (BHs) in massive gas-rich galaxies by means of high-resolution hydrodynamical simulations. We find that the presence of a massive gaseous disc allows recoiled BHs to return to the centre on a much shorter timescale than for purely stellar discs. Also, BH accretion and feedback can strongly modify the orbit of recoiled BHs and hence their return timescale, besides affecting the distribution of gas and stars in the galactic centre. However, the dynamical interaction of kicked BHs with the surrounding medium is in general complex and can facilitate both a fast return to the centre as well as a significant delay. The Bondi-Hoyle-Lyttleton accretion rates of the recoiling BHs in our simulated galaxies are favourably high for the detection of off-centred AGN if kicked within gas-rich discs -- up to a few per cent of the Eddington accretion rate -- and are highly variable on timescales of a few 10^7 yrs. In major merger simulations of gas-rich galaxies, we find that gravitational recoils increase the scatter in the BH mass -- host galaxy relationships compared to simulations without kicks, with the BH mass being more sensitive to recoil kicks than the bulge mass. A generic result of our numerical models is that the clumpy massive discs suggested by recent high-redshift observations, as well as the remnants of gas-rich mergers, exhibit a gravitational potential that falls steeply in the central regions, due to the dissipative concentration of baryons. As a result, supermassive BHs should only rarely be able to escape from massive galaxies at high redshifts, which is the epoch where the bulk of BH recoils is expected to occur.[Abridged]

We present 3D resonant radiative transfer simulations of the spatial and spectral diffusion of the Lya radiation from a central source in the host galaxies of high column density absorption systems at z ~ 3. The radiative transfer simulations are based on a suite of cosmological galaxy formation simulations that reproduce a wide range of observed properties of damped Lya absorption systems. The Lya emission is predicted to be spatially extended up to several arcsec, and the spectral width of the Lya emission is broadened to several hundred (in some case more than thousand) km/s. The distribution and the dynamical state of the gas in the simulated galaxies is complex, the latter with significant contributions from rotation and both in- and out-flows. The emerging Lya radiation extends to gas with column densities of N_HI ~ 10^{18} cm^{-2} and its spectral shape varies strongly with viewing angle. The strong dependence on the central \hi column density and the HI velocity field suggests that the Lya emission will also vary strongly with time on timescales of a few dynamical times of the central region. Such variations with time should be especially pronounced at times where the host galaxy undergoes a major merger and/or starburst. Depending on the pre-dominance of in- or out-flow along a given sightline and the central column density, the spectra show prominent blue peaks, red peaks or double-peaked profiles. Both spatial distribution and spectral shape are very sensitive to details of the galactic wind implementation. Stronger galactic winds result in more spatially extended Lya emission and - somewhat counterintuitively - a narrower spectral distribution.

We present measurements of the ionising ultraviolet background (UVB) at z ~ 5-6 using the quasar proximity effect. The fifteen quasars in our sample cover the range 4.6 < z_q < 6.4, enabling the first proximity effect measurements of the UVB at z > 5. The metagalactic hydrogen ionisation rate, Gamma_bkg, was determined by modelling the combined ionisation field from the quasar and the UVB in the proximity zone on a pixel-by-pixel basis. The optical depths in the spectra were corrected for the expected effect of the quasar until the mean flux in the proximity region equalled that in the average Ly-alpha forest, and from this we make a measurement of Gamma_bkg. A number of systematic effects were tested using synthetic spectra. Noise in the flux was found to be the largest source of bias at z ~ 5, while uncertainties in the mean transmitted Ly-alpha flux are responsible for the largest bias at z ~ 6. The impacts of large-scale overdensities and Lyman limit systems on Gamma_bkg were also investigated, but found to be small at z > 5. We find a decline in Gamma_bkg with redshift, from log(Gamma_bkg) = -12.15 $\pm$ 0.16 at z ~ 5 to log(Gamma_bkg) = -12.84 $\pm$ 0.18 at z ~ 6 (1 sigma errors). Compared to UVB measurements at lower redshifts, our measurements suggest a drop of a factor of five in the HI photoionisation rate between z ~ 4 and z ~ 6. The decline of Gamma_bkg appears to be gradual, and we find no evidence for a sudden change in the UVB at any redshift that would indicate a rapid change in the attenuation length of ionising photons. Combined with recent measurements of the evolution of the mean free path of ionising photons, our results imply decline in the emissivity of ionising photons by roughly a factor of two from z ~ 5 to 6, albeit with significant uncertainty due to the measurement errors in both Gamma_bkg and the mean free path.

We demonstrate that the Lyman alpha emission in the absorption troughs of a large sample of stacked damped Lyman alpha absorption systems (DLAS) presented by Rahmani et al (2010) is consistent with the spectral profiles and luminosities of a recently detected population of faint Lyman alpha emitters at z ~ 3. This result supports the suggestion that the faint emitters are to be identified with the host galaxies of DLAS at these redshifts.

We present new measurements of the temperature of the intergalactic medium (IGM) derived from the Lyman-alpha forest over 2.0 < z < 4.8. The small-scale structure in the forest of 61 high-resolution QSO spectra is quantified using a new statistic, the curvature, and the conversion to temperature calibrated using a suite of hydrodynamic simulations. At each redshift we focus on obtaining the temperature at an optimal overdensity probed by the Lyman-alpha forest, T(Delta), where the temperature is nearly a one-to-one function of the curvature regardless of the slope of the temperature-density relation. The median 2-sigma statistical uncertainty in these measurements is 8 per cent, though there may be comparable systematic errors due to the unknown amount of Jeans smoothing in the IGM. We use our T(Delta) results to infer the temperature at the mean density, T0. Even for a maximally steep temperature-density relation, T0 must increase from ~8000 K at z ~ 4.4 to >~12000 K at z ~ 2.8. This increase is not consistent with the monotonic decline in T0 expected in the absence of He II reionization. We therefore interpret the observed rise in temperature as evidence of He II reionization beginning at z >~ 4.4. The evolution of T0 is consistent with an end to He II reionization at z ~ 3, as suggested by opacity measurements of the He II Lyman-alpha forest, although the redshift at which T0 peaks will depend somewhat on the evolution of the temperature-density relation. These new temperature measurements suggest that the heat input due to the reionization of He II dominates the thermal balance of the IGM over an extended period with Delta_z >~ 1.

We present a suite of full hydrodynamical cosmological simulations that quantitatively address the impact of neutrinos on the (mildly non-linear) spatial distribution of matter and in particular on the neutral hydrogen distribution in the Intergalactic Medium (IGM), which is responsible for the intervening Lyman-alpha absorption in quasar spectra. The free-streaming of neutrinos results in a (non-linear) scale-dependent suppression of power spectrum of the total matter distribution at scales probed by Lyman-alpha forest data which is larger than the linear theory prediction by about 25% and strongly redshift dependent. By extracting a set of realistic mock quasar spectra, we quantify the effect of neutrinos on the flux probability distribution function and flux power spectrum. The differences in the matter power spectra translate into a ~2.5% (5%) difference in the flux power spectrum for neutrino masses with Sigma m_{\nu} = 0.3 eV (0.6 eV). This rather small effect is difficult to detect from present Lyman-alpha forest data and nearly perfectly degenerate with the overall amplitude of the matter power spectrum as characterised by sigma_8. If the results of the numerical simulations are normalized to have the same sigma_8 in the initial conditions, then neutrinos produce a smaller suppression in the flux power of about 3% (5%) for Sigma m_{\nu} = 0.6$ eV (1.2 eV) when compared to a simulation without neutrinos. We present constraints on neutrino masses using the Sloan Digital Sky Survey flux power spectrum alone and find an upper limit of Sigma m_{\nu} < 0.9$ eV (2 sigma C.L.), comparable to constraints obtained from the cosmic microwave background data or other large scale structure probes.

The thermal state of the intergalactic medium (IGM) provides an indirect probe of both the HI and HeII reionisation epochs. Current constraints on the IGM temperature from the Lya forest are restricted to the redshift range 2<z<4.5, limiting the ability to probe the thermal memory of HI reionisation toward higher redshift. In this work, we present the first direct measurement of the IGM temperature around a z=6 quasar by analysing the Doppler widths of Lya absorption lines in the proximity zone of SDSS J0818+1722. We use a high resolution (R= 40000) Keck/HIRES spectrum in combination with detailed numerical modelling to obtain the temperature at mean density, T_0=23600\pm^5000_6900K (\pm^9200_9300K) at 68 (95) per cent confidence assuming a prior probability 13500K<T_0<38500 K following HI and HeII reionisation. This enables us to place an upper limit on the redshift of HI reionisation, z_H, within 33 comoving Mpc of SDSS J0818+1722. If the quasar reionises the HeII in its vicinity, then in the limit of instantaneous reionisation we infer z_H<9.0 (11.0) at 68 (95) per cent confidence assuming photoheating is the dominant heat source and that HI reionisation is driven by ionising sources with soft spectra, typical of population II stars. If the HI and HeII in the IGM around SDSS J0818+1722 are instead reionised simultaneously by a population of massive metal-free stars, characterised by very hard ionising spectra, we obtain a tighter upper limit of z_H<8.4 (9.4). Initiating reionisation at higher redshifts produces temperatures which are too low with respect to our constraint unless the HI ionising sources or the quasar itself have spectra significantly harder than typically assumed.

We use detailed Lyalpha radiative transfer calculations to further test the claim of Rauch et al. (2008) that they have detected spatially extended faint Lyalpha emission from the elusive host population of Damped Lyalpha Absorption systems (DLAs) in their recent ultra-deep spectroscopic survey. We investigate the spatial and spectral distribution of Lyalpha emission due to star-formation at the centre of DLAs, and its dependence on the spatial and velocity structure of the gas. Our model simultaneously reproduces the observed properties of DLAs and the faint Lyalpha emitters, including the velocity width and column density distribution of DLAs and the large spatial extent of the emission of the faint emitters. Our modelling confirms previous suggestions that DLAs are predominately hosted by Dark Matter (DM) halos in the mass range 10^{9.5}-10^{12} M_sun, and are thus of significantly lower mass than those inferred for L_* Lyman Break Galaxies (LBGs). Our modelling suggests that DM halos hosting DLAs retain up to 20% of the cosmic baryon fraction in the form of neutral hydrogen, and that star formation at the centre of the halos is responsible for the faint Lyalpha emission. The scattering of a significant fraction of the Lyalpha emission to the observed radii, which can be as large as 50 kpc or more, requires the amplitude of the bulk motions of the gas at the centre of the halos to be moderate. The observed space density and size distribution of the emitters together with the incidence rate of DLAs suggests that the Lyalpha emission due to star formation has a duty cycle of ~ 25%.

We employ cosmological hydrodynamical simulations to study the growth of massive black holes (BHs) at high redshifts subject to BH merger recoils from gravitational wave emission. We select the most massive dark matter halo at z=6 from the Millennium simulation, and resimulate its formation at much higher resolution including gas physics and a model for BH seeding, growth and feedback. Assuming that the initial BH seeds are relatively massive, of the order of 10^5 Msun, and that seeding occurs around z~15 in dark matter haloes of mass 10^9-10^10 Msun, we find that it is possible to build up supermassive BHs (SMBHs) by z=6 that assemble most of their mass during extended Eddington-limited accretion periods. The properties of the simulated SMBHs are consistent with observations of z=6 quasars in terms of the estimated BH masses and bolometric luminosities, the amount of star formation occurring within the host halo, and the presence of highly enriched gas in the innermost regions of the host galaxy. After a peak in the BH accretion rate at z=6, the most massive BH has become sufficiently massive for the growth to enter into a much slower phase of feedback-regulated accretion. We explore the full range of expected recoils and radiative efficiencies, and also consider models with spinning BHs. In the most `pessimistic' case where BH spins are initially high, we find that the growth of the SMBHs can be potentially hampered if they grow mostly in isolation and experience only a small number of mergers. Whereas BH kicks can expel a substantial fraction of low mass BHs, they do not significantly affect the build up of the SMBHs. On the contrary, a large number of BH mergers has beneficial consequences for the growth of the SMBHs by considerably reducing their spin. [Abridged]

We use the probability distribution function (PDF) of the lya forest flux at z=2-3, measured from high-resolution UVES/VLT data, and hydrodynamical simulations to obtain constraints on cosmological parameters and the thermal state of the intergalactic medium (IGM) at z 2-3. The observed flux PDF at z=3 alone results in constraints on cosmological parameters in good agreement with those obtained from the WMAP data, albeit with about a factor two larger errors. The observed flux PDF is best fit with simulations with a matter fluctuation amplitude of sigma_8=0.8-0.85 \pm 0.07 and an inverted IGM temperature-density relation (gamma ~ 0.5-0.75), consistent with our previous results obtained using a simpler analysis. These results appear to be robust to uncertainties in the quasar (QSO) continuum placement. We further discuss constraints obtained by a combined analysis of the high-resolution flux PDF and the power spectrum measured from the Sloan Digital Sky Survey (SDSS) lya forest data. The joint analysis confirms the suggestion of an inverted temperature-density relation, but prefers somewhat higher values (sigma_8 ~ 0.9) of the matter fluctuation amplitude than the WMAP data and the best fit to the flux PDF alone. The joint analysis of the flux PDF and power spectrum (as well as an analysis of the power spectrum data alone) prefers rather large values for the temperature of the IGM, perhaps suggesting that we have identified a not yet accounted for systematic error in the SDSS flux power spectrum data or that the standard model describing the thermal state of the IGM at z ~ 2-3 is incomplete.

We use a set of semi-numerical simulations based on Zel'dovich approximation, friends-of-friends algorithm and excursion set formalism to generate reionization maps of high dynamic range with a range of assumptions regarding the distribution and luminosity of ionizing sources and the spatial distribution of sinks for the ionizing radiation. We find that ignoring the inhomogeneous spatial distribution of regions of high gas density where recombinations are important -- as is often done in studies of this kind -- can lead to misleading conclusions regarding the topology of reionization, especially if reionization occurs in the photon-starved regime suggested by Lya forest data. The inhomogeneous spatial distribution of recombinations significantly reduces the mean free path of ionizing photons and the typical size of coherently ionized regions. Reionization proceeds then much more as an outside-in process. Low-density regions far from ionizing sources become ionized before regions of high gas density not hosting sources of ionizing radiation. The spatial distribution of sinks of ionization radiation also significantly affects shape and amplitude the power spectrum of fluctuations of 21cm emission. The slope of the 21cm power spectrum as measured by upcoming 21cm experiments should be able to distinguish to what extent the topology of reionization proceeds outside-in or inside-out while the evolution of the amplitude of the power spectrum with increasing ionized mass fraction should be sensitive to the spatial distribution and the luminosity of ionizing sources.

We investigate the population of z=3 damped Lyman alpha absorbers (DLAs) in a recent series of high resolution galaxy formation simulations. The simulations are of interest because they form at z=0 some of the most realistic disk galaxies to date. No free parameters are available in our study: these have been fixed by physical and z=0 observational constraints, and thus our study provides a genuine consistency test. The precise role of DLAs in galaxy formation remains in debate, but they provide a number of strong constraints on the nature of our simulated bound systems at z=3 because of their coupled information on neutral H I densities, kinematics, metallicity and estimates of star formation activity. Our results, without any parameter-tuning, closely match the observed incidence rate and column density distributions of DLAs. Our simulations are the first to reproduce the distribution of metallicities (with a median of Z_{DLA} = Z_{solar}/20) without invoking observationally unsupported mechanisms such as dust biasing. This is especially encouraging given that these simulations have previously been shown to have a realistic 0<z<2 stellar mass-metallicity relation. Additionally, we see a strong positive correlation between sightline metallicity and low-ion velocity width, the normalization and slope of which comes close to matching recent observational results. However, we somewhat underestimate the number of observed high velocity width systems; the severity of this disagreement is comparable to other recent DLA focused studies. By z=0 the majority of the z=3 neutral gas forming the DLAs has been converted into stars, in agreement with rough physical expectations. [Abridged]

Large dynamic range numerical simulations of atomic cooling driven collapse of gas in pre-galactic DM haloes with T_vir ~ 10000 K show that the gas loses 90% and more of its angular momentum before rotational support sets in. In a fraction of these haloes where the metallicity is low and UV radiation suppresses H_2 cooling, conditions are thus very favourable for the rapid build-up of massive black holes. Depending on the progression of metal enrichment, the continued suppression of H_2 cooling by external and internal UV radiation and the ability to trap the entropy produced by the release of gravitational energy, the gas at the centre of the halo is expected to form a supermassive star, a stellar-mass black hole accreting at super-Eddington accretion rates or a compact star-cluster undergoing collisional run-away of massive stars at its centre. In all three cases a massive black hole of initially modest mass finds itself at the center of a rapid inflow of gas with inflow rates of ~ 1 M_solar\yr. The massive black hole will thus grow quickly to a mass of 10^5 to 10^6 M_solar until further inflow is halted either by consumption of gas by star formation or by the increasing energy and momentum feedback from the growing massive black hole. Conditions for the formation of massive seed black holes in this way are most favourable in haloes with T_vir ~ 15000 K and V_vir ~ 20 km\s with less massive haloes not allowing collapse of gas by atomic cooling and more massive haloes being more prone to fragmentation. This should imprint a characteristic mass on the mass spectrum of an early population of massive black hole seeds in pre-galactic haloes which will later grow into the observed population of supermassive black holes in galactic bulges.

We have used the hydrodynamical AMR code ENZO to investigate the dynamical evolution of the gas at the centre of dark matter haloes with virial velocities of ~ 20 - 30 kms and virial temperatures of ~ 13000-30000 K at z ~ 15 in a cosmological context. The virial temperature of the dark matter haloes is above the threshold where atomic cooling by hydrogen allows the gas to cool and collapse. We neglect cooling by molecular hydrogen and metals, as may be plausible if H_2 cooling is suppressed by a meta-galactic Lyman-Werner background or an internal source of Lyman-Werner photons, and metal enrichment has not progressed very far. The gas in the haloes becomes gravitationally unstable and develops turbulent velocities comparable to the virial velocities of the dark matter haloes. Within a few dynamical times it settles into a nearly isothermal density profile over many decades in radius losing most of its angular momentum in the process. About 0.1 - 1 % of the baryons, at the centre of the dark matter haloes, collapse into a self-gravitating, fat, ellipsoidal, centrifugally supported exponential disc with scale-length of ~ 0.075-0.27 pc and rotation velocities of 25-60 kms. We are able to follow the settling of the gas into centrifugal support and the dynamical evolution of the compact disc in each dark matter halo for a few dynamical times. The dynamical evolution of the gas at the centre of the haloes is complex. In one of the haloes the gas at the centre fragments into a triple system leading to strong tidal perturbations and eventually to the in-fall of a secondary smaller clump into the most massive primary clump. The formation of centrifugally supported self-gravitating massive discs is likely to be an important intermediary stage en route to the formation of a massive black hole seed.

The recently discovered population of ultra-faint extended line emitters can account for the majority of the incidence rate of Damped Lyman Alpha systems (DLAs) at z ~ 3 if the line emission is interpreted as Ly alpha. We show here that a model similar to that proposed by Haehnelt, Steinmetz, & Rauch (2000), which explains the incidence rate and kinematics of DLAs in the context of $\Lambda$CDM models for structure formation, also reproduces the size distribution of the new population of faint Ly alpha emitters for plausible parameters. This lends further support to identification of the emitters with the hitherto elusive population of DLA host galaxies. The observed incidence rate of DLAs together with the observed space density and size distribution of the emitters suggest a duty cycle of ~ 0.2 - 0.4 for the Ly alpha emission from DLA host galaxies. We further show that Ly alpha cooling is expected to contribute little to the Ly alpha emission for the majority of emitters. This leaves centrally concentrated star formation at a rate of a few tenths M_sun/yr, surrounded by extended Ly alpha halos with radii up to 30-50 kpc, as the most plausible explanation for the origin of the emission. Both the luminosity function of Ly alpha emission and the velocity width distribution of low ionization absorption require that galaxies inside Dark Matter (DM) halos with virial velocities < 50 - 70 km/s contribute little to the incidence rate of DLAs at z ~ 3, suggesting that energy and momentum input due to star formation efficiently removes gas from these halos. Galaxies with DM halos with virial velocities of 100 - 150 km/s appear to account for the majority of DLA host galaxies. DLA host galaxies at z ~ 3 should thus become the building blocks of typical present-day galaxies.

The next generation of X-ray telescopes have the potential to detect faint quasars at very high redshift and probe the early growth of massive black holes (BHs). We present modelling of the evolution of the optical and X-ray AGN luminosity function at 2 < z < 6 based on a CDM merger-driven model for the triggering of nuclear activity combined with a variety of fading laws. We extrapolate the merger-driven models to z > 6 for a range of BH growth scenarios. We predict significant numbers of sources at z ~ 6 with fluxes just an order of magnitude below the current detection limits and thus detectable with XEUS and Constellation-X, relatively independently of the fading law chosen. The predicted number of sources at even higher redshift depends sensitively on the early growth history of BHs. For passive evolution models in which BHs grow constantly at their Eddington limit, detectable BHs may be rare beyond z ~ 10 even with Generation-X. However, in the more probable scenario that BH growth at z > 6 can be described by passive evolution with a small duty cycle, or by our merger driven accretion model, then we predict that XEUS and Generation-X will detect significant numbers of black holes out to z ~ 10 and perhaps beyond.

We compare the improved measurement of the Lya forest flux probability distribution at 1.7<z<3.2 presented by Kim et al. (2007) to a large set of hydrodynamical simulations of the Lya forest with different cosmological parameters and thermal histories. The simulations are in good agreement with the observational data if the temperature-density relation for the low density intergalactic medium (IGM), T=T_0 Delta^{gamma-1}, is either close to isothermal or inverted (gamma<1). Our results suggest that the voids in the IGM may be significantly hotter and the thermal state of the low density IGM may be substantially more complex than is usually assumed at these redshifts. We discuss radiative transfer effects which alter the spectral shape of ionising radiation during the epoch of HeII reionisation as a possible physical mechanism for achieving an inverted temperature-density relation at z~3.

CODEX and ESPRESSO are concepts for ultra-stable, high-resolution spectrographs at the E-ELT and VLT, respectively. Both instruments are well motivated by distinct sets of science drivers. However, ESPRESSO will also be a stepping stone towards CODEX both in a scientific as well as in a technical sense. Here we discuss this role of ESPRESSO with respect to one of the most exciting CODEX science cases, i.e. the dynamical determination of the cosmic expansion history.

The redshifts of all cosmologically distant sources are expected to experience a small, systematic drift as a function of time due to the evolution of the Universe's expansion rate. A measurement of this effect would represent a direct and entirely model-independent determination of the expansion history of the Universe over a redshift range that is inaccessible to other methods. Here we investigate the impact of the next generation of Extremely Large Telescopes on the feasibility of detecting and characterising the cosmological redshift drift. We consider the Lyman alpha forest in the redshift range 2 < z < 5 and other absorption lines in the spectra of high redshift QSOs as the most suitable targets for a redshift drift experiment. Assuming photon-noise limited observations and using extensive Monte Carlo simulations we determine the accuracy to which the redshift drift can be measured from the Ly alpha forest as a function of signal-to-noise and redshift. Based on this relation and using the brightness and redshift distributions of known QSOs we find that a 42-m telescope is capable of unambiguously detecting the redshift drift over a period of ~20 yr using 4000 h of observing time. Such an experiment would provide independent evidence for the existence of dark energy without assuming spatial flatness, using any other cosmological constraints or making any other astrophysical assumption.

We present constraints on the mass of warm dark matter (WDM) particles derived from the Lyman-alpha flux power spectrum of 55 high- resolution HIRES spectra at 2.0 < z < 6.4. From the HIRES spectra, we obtain a lower limit of mwdm > 1.2 keV 2 sigma if the WDM consists of early decoupled thermal relics and mwdm > 5.6 keV (2 sigma) for sterile neutrinos. Adding the Sloan Digital Sky Survey Lyman-alpha flux power spectrum, we get mwdm > 4 keV and mwdm > 28 keV (2 sigma) for thermal relics and sterile neutrinos. These results improve previous constraints by a factor two.

CODEX, a high resolution, super-stable spectrograph to be fed by the E-ELT, the most powerful telescope ever conceived, will for the first time provide the possibility of directly measuring the change of the expansion rate of the Universe with time and much more, from the variability of fundamental constants to the search for other earths. A study for the implementation at the VLT of a precursor of CODEX, dubbed ESPRESSO, is presently carried out by a collaboration including ESO, IAC, INAF, IoA Cambridge and Observatoire de Geneve. The present talk is focused on the cosmological aspects of the experiment.

We have performed an extensive Voigt profile analysis of the neutral hydrogen (HI) and metal absorption present in a sample of 18 high resolution, high signal-to-noise QSO spectra observed with VLT/UVES. We use this analysis to separate the metal contribution from the HI absorption and present an improved measurement of the flux probability distribution function (PDF) due to HI absorption alone at <z> = 2.07, 2.52, and 2.94. The flux PDF is sensitive to the continuum fit in the normalised flux range 0.8 < F < 1.0 and to metal absorption at 0.2 < F < 0.8. Our new measurements of the flux PDF due to HI absorption alone are systematically lower at 0.2 < F < 0.8 by up to 30% compared to the widely used measurement of McDonald et al. (2000), based on a significantly smaller sample of Keck/HIRES data. This discrepancy is probably due to a combination of our improved removal of the metal absorption and cosmic variance, since variations in the flux PDF between different lines-of-sight are large. The HI effective optical depth at 1.7 < z < 4 is best fit with a single power law, (0.0023 +-0.0007) (1+z)^(3.65 +- 0.21), in good agreement with previous measurements from comparable data. As also found previously, the effect of noise on the flux distribution is not significant in high resolution, high signal-to-noise data.

We have conducted a long slit search for low surface brightness Lyman-alpha emitters at redshift 2.67 < z < 3.75. A 92 hour long exposure with VLT/FORS2 down to a 1-sigma surface brightness detection limit of 8x10^-20 erg/cm2/s/sqarcsec yielded a sample of 27 single line emitters with fluxes of a few times 10^-18 erg/s/cm2. We present arguments that most objects are indeed Lyman-alpha. The large comoving number density, the large covering factor, dN/dz ~ 0.2-1, and the often extended Lyman-alpha emission suggest that the emitters be identified with the elusive host population of damped Lyman-alpha systems (DLAS) and high column density Lyman limit systems. A small inferred star formation rate, perhaps supplanted by cooling radiation, appears to energetically dominate the Lyman-alpha emission, and is consistent with the low metallicity, low dust content, and theoretically inferred low masses of DLAS, and with the relative lack of success of earlier searches for their optical counterparts. (abridged)

We present constraints on the amplitude and shape of the matter power spectrum and the density of dark matter within the framework of a standard LambdaCDM model. We use a Markov Chain Monte Carlo approach to combine independent measurements of the three dimensional weak gravitational lensing shear field by the COSMOS survey, of low and high resolution Ly-alpha forest flux power spectrum by SDSS and LUQAS, and of Cosmic Microwave Background temperature and polarization anisotropies by WMAP. We note good agreement between the amplitude of the matter power spectrum on intermediate and small scales as inferred from Ly-alpha forest and lensing data. The Ly-alpha forest data helps to break the sigma_8-Omega_m degeneracy characteristic of weak lensing results, yielding sigma_8 = 0.876 +- 0.048 for COSMOS plus Ly-alpha SDSS data. This is somewhat larger than the value preferred by the WMAP year three CMB data. Combining all three data sets significantly tightens the constraints on sigma_8, the spectral index of primordial density fluctuation n_s, a possible running of the spectral index n_run and the matter density Omega_m. Assuming no running, the joint constraints for COSMOS, SDSS and WMAP are sigma_8 = 0.800 +- 0.023, n_s = 0.971 +- 0.011, Omega_m = 0.247 +- 0.016 (1-sigma error bars).

We use a large set of hydrodynamical simulations, combined with measurements of the Lyman alpha opacity of the IGM taken from the literature, to obtain robust estimates for the photoionization rate per hydrogen atom at z=5 and 6. We find the photoionization rate drops by a factor of two and four, respectively, compared to our recent measurements at z = 2 - 4. The number of ionizing photons emitted by known sources at z=5 and 6, based on an extrapolation of source numbers below the detection limit and standard assumptions for the relationship between the ionizing emissivity and observed luminosity density at 1500 Angstroms, are in reasonable agreement with the photoionization rates inferred from the Lyman alpha forest if the escape fraction of ionizing photons from galaxies is large (>= 20 per cent). Claims to the contrary may be attributed to the adoption of an unduly high value for the clumping factor of ionized hydrogen. Using physically motivated assumptions for the mean free path of ionizing photons our measurements of the photoionization rate can be turned into an estimate of the ionizing emissivity. In comoving units the inferred ionizing emissivity is nearly constant over the redshift range 2-6 and corresponds to 1.5-3 photons emitted per hydrogen atom over a time interval corresponding to the age of the Universe at z=6. This strongly suggests that the epoch of reionization was photon-starved and extended. [Abridged]

Absorption spectra of high redshift quasars exhibit an increasingly thick Ly-alpha forest towards z~6. However, the interpretation of these spectra is complicated by the fact that the Ly-alpha optical depth is already large for neutral hydrogen fractions in excess of 10^-4, and also because quasars are expected to reside in dense regions of the IGM. We present a model for the evolution of the ionization state of the IGM which is applicable to the dense, biased regions around high-redshift quasars as well as more typical regions in the IGM, and combine this with numerical radiative transfer simulations. Our model is able to simultaneously reproduce the observed Ly-alpha forest opacity at 4<z<6, the ionizing photon mean-free-path at z~4 and the rapid evolution of highly ionized near-zone sizes around high-redshift quasars at 5.8<z<6.4. We find that within 5 physical Mpc of a high redshift quasar, the evolution of the ionization state of the IGM precedes that in more typical regions by around 0.3 redshift units. More importantly, when combined with the rapid increase in the ionizing photon mean-free-path expected shortly after overlap, this offset results in an ionizing background near the quasar which exceeds the value in the rest of the IGM by a factor of ~2-3. We further find that in the post-overlap phase of reionization the size of the observed quasar near-zones is not directly sensitive to the neutral hydrogen fraction of the IGM. Instead, these sizes probe the level of the background ionization rate and the temperature of the surrounding IGM. The observed rapid evolution of the quasar near-zone sizes at 5.8<z<6.4 can thus be explained by the rapid evolution of the ionizing background, which in our model is caused by the completion of overlap at the end of reionization by 6<z<7.

We examine a large set of synthetic quasar spectra to realistically assess the potential of using the relative sizes of highly ionized near-zones in the Lya and Lyb forest as a probe of the neutral hydrogen content of the intergalactic medium (IGM) at z>6. The scatter in the relative near-zone size distribution, induced by underlying fluctuations in the baryonic density field and the filtering of ionizing radiation, is considerable even for fixed assumptions about the IGM neutral fraction. As a consequence, the current observational data cannot distinguish between an IGM which is significantly neutral or highly ionized just above z=6. Under standard assumptions for quasar ages and ionizing luminosities, a future sample of several tens of high resolution Lya and Lyb near-zone spectra should be capable of distinguishing between a volume weighted neutral hydrogen fraction in the IGM which is greater or less than 10 per cent.

We use state-of-the-art hydrodynamical simulations combined with a 1D radiative transfer code to assess the extent to which the highly ionized regions observed close to z~6 quasars, which we refer to as near-zones, can constrain the ionization state of the surrounding IGM. We find the appearance in Lya absorption of a quasar HII ionization front expanding into a neutral IGM can be very similar to a classical proximity zone, produced by the enhancement in ionizing flux close to a quasar embedded in a highly ionized IGM. The observed sizes of these highly ionized near-zones and their redshift evolution can be reproduced for a wide range of IGM neutral hydrogen fractions for plausible values of the luminosity and lifetime of the quasars. The observed near-zone sizes at the highest observed redshifts are equally consistent with a significantly neutral and a highly ionized surrounding IGM. Stronger constraints on the IGM neutral hydrogen fraction can be obtained by considering the relative size of the near-zones in the Lya and Lyb regions of a quasar spectrum. A large sample of high quality quasar absorption spectra with accurate determinations of near-zone sizes and their redshift evolution in both the Lya and Lyb regions should confirm or exclude the possibility that the Universe is predominantly neutral at the highest observed redshifts. The width of the discrete absorption features in these near-zones will contain important additional information on the ionization state and the previous thermal history of the IGM at these redshifts.

We compare simulations of the Lyman-alpha forest performed with two different hydrodynamical codes, Gadget-2 and Enzo. A comparison of the dark matter power spectrum for simulations run with identical initial conditions show differences of 1-3% at the scales relevant for quantitative studies of the Lyman-alpha forest. This allows a meaningful comparison of the effect of the different implementations of the hydrodynamic part of the two codes. Using the same cooling and heating algorithm in both codes the differences in the temperature and the density probability distribution function are of the order of 10%. These differences can be mainly attributed to a slight mismatch in the resolution. The differences are smaller than or equal to the effects of boxsize and resolution on these statistics. Once resolution effects are taken into account the differences in the flux power spectrum - the statistics most widely used for estimating the matter power spectrum and cosmological parameters from Lyman-alpha forest data - are about 5%. This is again smaller than or equal to the effects of boxsize and resolution. Numerical uncertainties due to a particular implementation of solving the hydrodynamic or gravitational equations appear therefore to contribute only moderately to the error budget in estimates of the flux power spectrum from numerical simulations. We further find that the differences in the flux power spectrum for Enzo simulations run with and without adaptive mesh refinement are also of order 5% or smaller. The latter require 10 times less CPU time making the CPU time requirement similar to that of a version of Gadget-2 that is optimised for Lyman-alpha forest simulations.

We present constraints on the mass of warm dark matter (WDM) particles from a combined analysis of the matter power spectrum inferred from the Sloan Digital Sky Survey \lya flux power spectrum at 2.2<z<4.2, cosmic microwave background data, and the galaxy power spectrum. We obtain a lower limit of m~10 keV (2 sigma) if the WDM consists of sterile neutrinos and m~2 keV (2 sigma) for early decoupled thermal relics. If we combine this bound with the constraint derived from x-ray flux observations in the Coma cluster, we find that the allowed sterile neutrino mass is ~10 keV (in the standard production scenario). Adding constraints based on x-ray fluxes from the Andromeda galaxy, we find that dark matter particles cannot be sterile neutrinos, unless they are produced by a nonstandard mechanism (resonant oscillations, coupling with the inflaton) or get diluted by some large entropy release.

HII regions surrounding supermassive black holes (BHs) in an otherwise still neutral intergalactic medium (IGM) are likely to be the most easily detectable sources by future 21cm experiments like LOFAR. We have made predictions for the size distribution of such HII regions for several physically motivated models for BH growth at high redshift and compared this to the expected LOFAR sensitivity to these sources. The number of potentially detectable HII regions does not only depend on the ionisation state of the intergalactic medium and the decoupling of the spin temperature of the neutral hydrogen from the cosmic microwave background (CMB) temperature, but is also strongly sensitive to the rate of growth of BHs at high redshift. If the supermassive BHs at redshift 6 were built up via continuous Eddington-limited accretion from low mass seed BHs at high redshift, then LOFAR is not expected to detect isolated QSO HII regions at redshifts much larger than 6, and only if the IGM is still significantly neutral. If the high-redshift growth of BHs starts with massive seed BHs and is driven by short-lived accretion events following the merging of BH hosting galaxies then the detection of HII regions surrounding supermassive BHs may extend to redshifts as large as 8-9 but is still very sensitive to the redshift to which the IGM remains significantly neutral. The most optimistic predictions are for a model where the supermassive BHs at z>6 have grown slowly. HII regions around supermassive BHs may then be detected to significantly larger redshifts.

We present the transverse flux correlation function of the Lyman-alpha forest in quasar absorption spectra at z~2.1 from VLT-FORS and VLT-UVES observations of a total of 32 pairs of quasars; 26 pairs with separations in the range 0.6<theta<4 arcmin and 6 pairs with 4<theta<10 arcmin. Correlation is detected at the 3sigma level up to separations of the order of ~4 arcmin (or ~4.4 h^{-1} Mpc comoving at z=2.1 for Omega_{m}=0.3 and Omega_{Lambda}=0.7). We have, furthermore, measured the longitudinal correlation function at a somewhat higher mean redshift (z=2.39) from 20 lines of sight observed with high spectral resolution and high signal-to-noise ratio with VLT-UVES. We compare the observed transverse and longitudinal correlation functions to that obtained from numerical simulations and illustrate the effect of spectral resolution, thermal broadening and peculiar motions. The shape and correlation length of the correlation functions are in good agreement with those expected from absorption by the filamentary and sheet-like structures in the photoionized warm intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 C_{iv} systems detected along the lines of sight toward the pairs of quasars we also investigate the transverse correlation of metals on the same scales. The observed transverse correlation function of intervening C_{iv} absorption systems is consistent with that of a randomly distributed population of absorbers. This is likely due to the small number of pairs with separation less than 2 arcmin. We detect, however, a significant overdensity of systems in the sightlines towards the quartet Q0103-294A&B, Q0102-2931 and Q0102-293 which extends over the redshift range 1.5<z<2.2 and an angular scale larger than 10 arcmin.

A combined analysis of Cosmic Microwave Background (CMB) and Lyman-a forest data allows to constrain the matter power spectrum from small scales of about 1 Mpc/h all the way to the horizon scale. The long lever arm and complementarity provided by such an analysis has previously led to a significant tightening of the constraints on the shape and the amplitude of the power spectrum of primordial density fluctuations. We present here a combined analysis of the WMAP three year results with Lyman-a forest data. The amplitude of the matter power spectrum sigma_8 and the spectral index ns inferred from the joint analysis with high resolution Lyman-a forest data and low resolution Lyman-a forest data as analyzed by Viel & Haehnelt (2006) are consistent with the new WMAP results to within 1 sigma. The joint analysis with the mainly low resolution data as analysed by McDonald et al. (2005) suggests a value of sigma_8 which is ~ 2 sigma higher than that inferred from the WMAP three year data alone. The joint analysis of the three year WMAP and the Lyman-a forest data also does not favour a running of the spectral index. The best fit values for a combined analysis of the three year WMAP data, other CMB data, 2dF and the Lyman-a forest data are (sigma_8, ns) = (0.78\pm 0.03,0.96 \pm 0.01).

We implement the hydro-PM (HPM) technique (Gnedin & Hui 1998) in the hydrodynamical simulation code GADGET-II and quantify the differences between this approximate method and full hydrodynamical simulations of the Lyman-alpha forest in a concordance LCDM model. At redshifts z=3 and z=4, the differences between the gas and dark matter (DM) distributions, as measured by the one-point distribution of density fluctuations, the density power spectrum and the flux power spectrum, systematically decrease with increasing resolution of the HPM simulqation. However, reducing these differences to less than a few percent requires a significantly larger number of grid-cells than particles, with a correspondingly larger demand for memory. Significant differences in the flux decrement distribution remain even for very high resolution hydro-PM simulations, particularly at low redshift. At z=2, the differences between the flux power spectra obtained from HPM simulations and full hydrodynamical simulations are generally large and of the order of 20-30 %, and do not decrease with increasing resolution of the HPM simulation. This is due to the presence of large amounts of shock-heated gas, a situation which is not adequately modelled by the HPM approximation. We confirm the results of Gnedin & Hui (1998) that the statistical properties of the flux distribution are discrepant by > 5-20 % when compared to full hydrodynamical simulations. The discrepancies in the flux power spectrum are strongly scale- and redshift-dependent and extend to large scales. Considerable caution is needed in attempts to use calibrated HPM simulations for quantitative predictions of the flux power spectrum and other statistical properties of the Lyman-alpha forest.

We measure the matter power spectrum from the Lyman-alpha forest of 31 quasar spectra spanning the redshift range of 1.6-3.6. The optical depth, $\tau$, for absorption of the intergalactic medium is obtained from the flux using the inversion method of Nusser & Haehnelt (1999). The optical depth is converted to density by using a simple power law relation, $\tau \propto (1+\delta)^\alpha$. The non-linear 1D power spectrum of the gas density is then inferred with a method that makes simultaneous use of the 1 and 2 point statistics of the flux and compared against theoretical models with a likelihood analysis. A Cold Dark Matter (CDM) model with standard cosmological parameters fits the data well. The power spectrum amplitude is measured to be (assuming a flat Universe), $\sigma_8 = (0.9 \pm 0.09)\times (\Omega_{\rm m}/0.3)^{-0.3}$, with $\alpha$ varying in the range of $1.56-1.8$ with redshift. Enforcing the same cosmological parameters in all four redshift bins, the likelihood analysis suggests some evolution in the density-temperature relation and the thermal smoothing length of the gas. The inferred evolution is consistent with that expected if reionization of \hep occurred at $z\sim 3.2$. A joint analysis with the WMAP results together with a prior on the Hubble constant as suggested by the HST key project data, yields values of $\Omega_{\rm m}$ and $\sigma_8$ that are consistent with the cosmological concordance model. We also perform a further inversion to obtain the linear 3D power spectrum of the matter density fluctuations.

The low density hydrogen and helium in the IGM probed by QSO absorption lines is sensitive to the amplitude and spectral shape of the metagalactic UV background. We use realistic HI and HeII Ly-alpha forest spectra, constructed from state-of-the-art hydrodynamical simulations of a Lambda-CDM Universe, to confirm the reliability of using line profile fitting techniques to infer the ratio of the metagalactic HI and HeII ionization rates. We further show that the large spatial variations and the anti-correlation with HI absorber density observed in the ratio of the measured HeII to HI column densities can be explained in a model where the HI ionization rate is dominated by the combined UV emission from young star forming galaxies and QSOs and the HeII ionization rate is dominated by emission from QSOs only. In such a model the large fluctuations in the column density ratio are due to the small number of QSOs expected to contribute at any given point to the HeII ionization rate. A significant contribution to UV emission at the HeII photoelectric edge from hot gas in galaxies and galaxy groups would decrease the expected fluctuations in the column density ratio. Consequently, this model appears difficult to reconcile with the large increase in HeII opacity fluctuations towards higher redshift. Our results further strengthen previous suggestions that observed HeII Ly-alpha forest spectra at z~2-3.5 probe the tail end of the reionization of HeII by QSOs.

Attempts of Magain et al (2005) to detect the host galaxy of the bright QSO HE0450--2958 have not been successful. We suggest that the supermassive black hole powering the QSO was ejected from the observed ULIRG at the same redshift and at 1.5 arcsec distance. Ejection could have either be caused by recoil due to gravitational wave emission from a coalescing binary of supermassive black holes or the gravitational slingshot of three or more supermassive black holes in the ongoing merger of galaxies which triggered the starburst activity in the ULIRG. We discuss implications for the possible hierarchical build-up of supermassive black holes from intermediate and/or stellar mass black holes, and for the detection of coalescing supermassive binary black holes by LISA.

An excess of sight-lines close to Lyman-break galaxies (LBGs) with little or no absorption in QSO absorption spectra has been reported and has been interpreted as the effect of galactic winds on the Intergalactic Medium. We use here numerical simulations to investigate the flux probability function close to plausible sites of LBGs. We show that the flux distribution near our LBGs in the simulation depends strongly on redshift, and is very sensitive to the averaging procedure. We show that a model without galactic winds and a model with a wind bubble size of 0.5Mpc/h (comoving) are equally consistent with the new determination of the conditional flux distribution by Adelberger et al. (2005). Models with the larger bubble sizes suggested by the previous observations of Adelberger et al. (2003) based on a much smaller sample at higher redshift are not consistent with the new data. We therefore argue that the volume filling factor of galactic winds driven by LBGs may be much smaller than previously thought, and that most of the metals responsible for the metal absorption associated with the low column density Lya forest are unlikely to have been ejected by LBGs.

(abridged) The flux power spectrum of the Lyman-alpha forest in quasar (QSO) absorption spectra is sensitive to a wide range of cosmological and astrophysical parameters and instrumental effects. Modelling the flux power spectrum in this large parameter space to an accuracy comparable to the statistical uncertainty of large samples of QSO spectra is very challenging. We use here a coarse grid of hydrodynamical simulations run with GADGET-2 to obtain a ``best guess'' model around which we calculate a finer grid of flux power spectra using a Taylor expansion of the flux power spectrum to first order. We find that the SDSS flux power spectrum alone is able to constrain a wide range of parameters including the amplitude of the matter power spectrum sigma_8, the matter density Omega_m, the spectral index of primordial density fluctuations n, the effective optical depth tau_eff and its evolution. The thermal history of the Intergalactic Medium (IGM) is, however, poorly constrained and the SDSS data favour either an unplausibly large temperature or an unplausibly steep temperature-density relation. By enforcing a thermal history of the IGM consistent with that inferred from high-resolution QSO spectra, we find the following values for the best fitting model (assuming a flat Universe with a cosmological constant and zero neutrino mass): Omega_ m=0.28 \pm 0.03, n=0.95\pm0.04, \sigma_8=0.91\pm0.07 (1\sigma error bars).We argue that the major uncertainties in this measurement are still systematic rather than statistical.

We derive the transverse flux correlation function in the Lyman-alpha forest at z~2.1 from VLT-FORS observations of a total of 32 pairs of quasars. The shape and correlation length of the transverse correlation function are in good agreement with the paradigm of intergalactic medium predicted in CDM-like models for structures formation. Using a sample of 139 Civ systems detected along the lines of sight toward the pairs of quasars we investigate the transverse correlation of metals on the same scales. We find that the correlation function is consistent with that of a randomly distributed population of Civ systems. However, we detect an important overdensity of systems in front of a quartet.

We study the kinematics of the gaseous cosmic web at high redshift with Lyman alpha forest absorption in multiple QSO sightlines. Using a simple analytic model and a cosmological hydrodynamic simulation we constrain the underlying three-dimensional distribution of velocities from the observed line-of-sight distribution of velocity shear across the plane of the sky. The distribution is found to be in good agreement with the intergalactic medium (IGM) undergoing large scale motions dominated by the Hubble flow. Modeling the Lyman alpha clouds analytically and with a hydrodynamics simulation, the average expansion velocity of the gaseous structures causing the Lyman alpha forest in the lower redshift (z = 2) sample appears about 20 percent lower than the local Hubble expansion velocity. We interpret this as tentative evidence for some clouds undergoing gravitational collapse. However, the distribution of velocities is highly skewed, and the majority of clouds at redshifts from 2 to 3.8 expand typically about 5 - 20 percent faster than the Hubble flow. This behavior is explained if most absorbers in the column density range typically detectable are expanding filaments that stretch and drain into more massive nodes. We find no evidence for the observed distribution of velocity shear being significantly influenced by processes other than Hubble expansion and gravitational instability, like galactic winds. To avoid overly disturbing the IGM, winds may be old and/or limp by the time we observe them in the Lyman alpha forest, or they may occupy only an insignificant volume fraction of the IGM. (abridged)

Understanding the sources responsible for reionizing the Universe is a key goal of observational cosmology. A discrepancy has existed between the metagalactic hydrogen ionization rate, Gamma_HI, predicted by early hydrodynamical simulations of the Lyman-alpha forest if scaled to appropriate assumptions for the IGM temperature, when compared to values predicted from the proximity effect. We present new estimates for Gamma_HI in the redshift range 2<z<4 based on hydrodynamical simulations of the Lyman-alpha forest opacity. Within the current concordance cosmology, and assuming updated QSO emissivity rates, a substantial contribution to the UV background from young star-forming galaxies appears to be required over the entire redshift range. Our results are consistent with lower-end estimates from the proximity effect. It is also found that the errors on the ionization rate are dominated by uncertainties in the thermal state of the intergalactic medium and the r.m.s fluctuation amplitude at the Jeans scale.

The matter power spectrum at comoving scales of (1-40) h^{-1} Mpc is very sensitive to the presence of Warm Dark Matter (WDM) particles with large free streaming lengths. We present constraints on the mass of WDM particles from a combined analysis of the matter power spectrum inferred from the large samples of high resolution high signal-to-noise Lyman-alpha forest data of Kim et al. (2004) and Croft et al. (2002) and the cosmic microwave background data of WMAP. We obtain a lower limit of m_wdm > 550 eV (2 sigma for early decoupled thermal relics and m_wdm > 2.0 keV (2 sigma) for sterile neutrinos. We also investigate the case where in addition to cold dark matter a light thermal gravitino with fixed effective temperature contributes significantly to the matter density. In that case the gravitino density is proportional to its mass, and we find an upper limit m_{3/2} < 16 eV (2 sigma). This translates into a bound on the scale of supersymmetry breaking, Lambda_{susy} < 260 TeV, for models of supersymmetric gauge mediation in which the gravitino is the lightest supersymmetric particle.

Estimates of the metagalactic hydrogen ionization rate from the Lyman-alpha forest opacity in QSO absorption spectra depend on the complex density distribution of neutral hydrogen along the line-of-sight. We use a large suite of high resolution hydrodynamical simulations to investigate in detail the dependence of such estimates on physical and numerical parameters in the context of Lambda-CDM models. Adopting fiducial values for cosmological parameters together with published values of the temperature of the IGM and the effective optical depth, the metagalactic ionization rates which reproduce the Lyman-alpha effective optical depth at z=[2,3,4] are Gamma_HI=[1.3\pm^0.8_0.5, 0.9\pm0.3, 1.0\pm^0.5_0.3] \times 10^-12 s^-1, respectively. The errors include estimates of uncertainties in the relevant physical parameters and the numerical accuracy of the simulations. We find the errors are dominated by the uncertainty in the temperature of the low-density IGM. The estimated metagalactic hydrogen ionization rate for the neutral hydrogen distribution in the current concordance Lambda-CDM model is more than four times the value inferred for that in an Einstein-de Sitter model of the same r.m.s. density fluctuation amplitude sigma_8. The estimated ionization rate is also more than double that expected from updated estimates of the emissivity of observed QSOs alone. A substantial contribution from galaxies appears to be required at all redshifts.

We use the LUQAS sample (Kim et al. 2004), a set of 27 high-resolution and high signal-to-noise QSO absorption spectra at a median redshift of z=2.25, and the data from Croft et al. (2002) at a median redshift of z=2.72, together with a large suite of high-resolution large box-size hydro-dynamical simulations, to estimate the linear dark matter power spectrum on scales 0.003 s/km < k <0.03 s/km. Our re-analysis of the Croft et al. data agrees well with their results if we assume the same mean optical depth and gas temperature-density relation. The inferred linear dark matter power spectrum at z=2.72 also agrees with that inferred from LUQAS at lower redshift if we assume that the increase of the amplitude is due to gravitational growth between these redshifts. We further argue that the smaller mean optical depth measured from high-resolution spectra is more accurate than the larger value obtained from low-resolution spectra by Press et al. (1993) which Croft et al. used. For the smaller optical depth we obtain a ~ 20% higher value for the rms fluctuation amplitude of the matter density. By combining the amplitude of the matter power spectrum inferred from the Lyman-alpha forest with the amplitude on large scales inferred from measurements of the CMB we obtain constraints on the primordial spectral index n and the normalisation sigma_8. For values of the mean optical depth favoured by high-resolution spectra, the inferred linear power spectrum is consistent with a LambdaCDM model with a scale-free (n=1) primordial power spectrum.

The combined analysis of the cosmic microwave background on large scales and Lyman-alpha forest on small scales provides a sufficiently long lever arm to obtain strong constraints on the slope and curvature of the power spectrum of primordial density fluctuations. We present results from the combination of the first year WMAP data and the dark matter power spectrum inferred by Viel et al. (2004) for two different sets of high resolution and high signal-to-noise quasar absorption spectra: the Croft et al. (2002) sample with a median redshift z=2.72 and the LUQAS sample (Kim et al. 2004) with a median redshift z=2.125. The best fit value for the {\it rms} fluctuation amplitude of matter fluctuations is sigma_8 =0.94 +- 0.08 and n=0.99 +- 0.03, if we do not include running of the spectral index. The best fit model with a running spectral index has parameters n=0.959 +- 0.036 and n_run=-0.033 +- 0.025. The data is thus consistent with a scale-free primordial power spectrum with no running of the spectral index. We further include tensor modes and constrain the slow-roll parameters of inflation.

We present a deep XMM-Newton observation of the z=6.30 QSO SDSS J1030+0524, the second most distant quasar currently known. The data contain sufficient counts for spectral analysis, demonstrating the ability of XMM-Newton to measure X-ray spectral shapes of z~6 QSOs with integration times >100ks. The X-ray spectrum is well fit by a power law with index Gamma=2.12 +/- 0.11, an optical-X-ray spectral slope of a_{ox}=-1.80, and no absorption excess to the Galactic value, though our data are also consistent with a power law index in the range 2.02 < Gamma < 2.5 and excess absorption in the range 0 < N_H(cm^-2) < 8x10^22. There is also a possible detection (2 sigma) of FeKa emission. The X-ray properties of this QSO are, overall, similar to those of lower-redshift radio-quiet QSOs. This is consistent with the statement that the X-ray properties of radio-quiet QSOs show no evolution over 0<z<6.3. Combined with previous results, this QSO appears indistinguishable in any way from lower redshift QSOs, indicating that QSOs comparable to those seen locally existed less than one Gyr after the Big Bang.

In a recent paper Pell'o et al. have reported a candidate z=10 galaxy, A1835#1916, which was found in a near-infrared survey of the central regions of the gravitational lensing cluster A1835. If this detection is confirmed and the detection rate turns out to be typical, then the volume averaged ultraviolet emissivity must be rising rapidly with increasing redshift. For a magnification due to gravitational lensing by a factor M>25 estimated by Pell'o et al., the inferred star formation rate density at z=10 would be about one order of magnitude higher than estimates of the star formation rate density at z=6. Objects at z=10 would contribute substantially to the total source counts at 1.6 micronm and the estimated space density of sources may exceed the space density of dark matter haloes in a LCDM model. We therefore argue that if A1835#1916 is indeed at z=10 then either the magnification factor may have been overestimated or the galaxy has a top-heavy initial mass function. Sources with the ultraviolet flux and space density of A1835#1916 may be sufficient to reionise most of the diffuse hydrogen in the Universe at z=10. We further use a correlation between the equivalent width and the redshift of the Ly-alpha emission line with respect to the systemic redshift observed in Lyman break galaxies to obtain constraints on the ionisation state of the surrounding intergalactic medium (IGM) from the Gunn-Peterson absorption. These constraints also argue in favour of the surrounding IGM being partially or fully ionised.

We use the LUQAS sample of 27 high resolution, high signal-to-noise QSO absorption spectra (Kim et al. 2004) and the results from Croft et al. (2002) together with a suite of high-resolution hydro-dynamical simulations run with the GADGET-II code, to infer the linear dark matter power spectrum on scales of \~ 0.3-30 Mpc at z=2.125 and at z=2.72.

The flux power spectra of the Lyman-alpha forest from a sample of 27 QSOs taken with the high resolution echelle spectrograph UVES on VLT are presented. We find a similar fluctuation amplitude at the peak of the ``3D'' flux power spectrum at k ~ 0.03 (km/sec)^(-1) as the study by Croft et al. (2002), in the same redshift range. The amplitude of the flux power spectrum increases with decreasing redshift if corrected for the increase in the mean flux level as expected if the evolution of the flux power spectrum is sensitive to the gravitational growth of matter density fluctuations. This is in agreement with the findings of McDonald et al. (2000) at larger redshift. The logarithmic slope of the "3D" flux power spectrum, P_F(k), at large scales k < 0.03 (km/sec)^(-1), is 1.4 +- 0.3, i.e. 0.3 shallower than that found by Croft et al. (2002) but consistent within the errors.

We study the effect of galactic outflows on the Lya forest. The winds are modelled as fully ionised spherical bubbles centered around galactic haloes. The observed flux probability distribution and flux power spectrum limit the volume filling factor of bubbles to be less than 10%. We have compared the mean flux as a function of distance from haloes with the Adelberger etal (ASSP) measurement. A model with bubbles of constant radius of 1.5Mpc/h surrounding the most massive haloes matches the high transmissivity at separations <0.5Mpc/h, but the increase of the transmissivity at small scales is more gradual than observed. At separations <0.5Mpc/h, our estimate of the cosmic variance error is $\Delta F\sim 0.3$, 30% higher than that of ASSP. The difficulty in matching the rise in the transmissivity at small separations is caused by residual absorption of neutral hydrogen lying physically outside the bubbles but having a redshift position similar to the haloes. The flux level is thus sensitive to the amplitude of the coherent velocity shear near haloes which increases with halo mass in the simulation. A model where LBGs are starbursts in small mass haloes matches the observations with smaller bubble size than one where more massive haloes host the LBGs. A starburst model with a bubble radius of 1Mpc/h and a volume filling factor of 2% is consistent with the ASSP measurements at the 1-1.5$\sigma$ level. If this model is correct the sharp rise of the transmissivity at small separation in the ASSP sample is due to cosmic variance and is expected to become more moderate for a larger sample.

We present a determination of the bispectrum of the flux in the Lyman-alpha forest of QSO absorption spectra obtained from LUQAS which consists of spectra observed with the high resolution Echelle spectrograph UVES. Typical errors on the observed bispectrum as obtained from a jack-knife estimator are ~ 50%. For wavenumbers in the range 0.03 s/km <k < 0.1 s/km the observed bispectrum agrees within the errors with that of the synthetic absorption spectra obtained from numerical hydro-simulations of a LCDM model with and without feedback from star formation. Including galactic feedback changes the bispectrum by less than 10%. At smaller wavenumbers the associated metal absorption lines contribute about 50% to the bispectrum and the observed bispectrum exceeds that of the simulations. At wavenumbers k< 0.03 (s/km) second-order perturbation theory applied to the flux spectrum gives a reasonable (errors smaller than 30%) approximation to the bispectra of observed and simulated absorption spectra. The bispectrum of the observed absorption spectra also agrees, within the errors, with that of a randomized set of absorption spectra where a random shift in wavelength has been added to absorption lines identified with VPFIT. This suggests that for a sample of the size presented here, the errors on the bispectrum are too large to discriminate between models with very different 3D distribution of Lyman-alpha absorption. If it were possible to substantially reduce these errors for larger samples of absorption spectra, the bispectrum might become an important statistical tool for probing the growth of gravitational structure in the Universe at redshift z>2.

Artificial absorption spectra are used to test a variety of instrumental and physical effects on the pixel correlation technique for the detection of weak OVI absorption. At HI optical depths <= 0.3-1, the apparent OVI detections are spurious coincidences due to HI absorption at other redshifts. In this range, the apparent OVI optical depth is independent of HI optical depth. At larger HI optical depths, apparent OVI optical depth and HI optical depth are correlated. Detailed modelling is required in order to interpret the significance of this relation. High-resolution spectra of four QSOs together with a large suit of synthetic spectra are used to show that the detection of OVI in individual spectra is only statistically significant for overdensities >= 5. These overdensities are larger than would be naively inferred from the onset of the correlation and a tight optical depth-density relation. The lower limit for the volume filling factor of regions which are enriched by OVI is 4% at 95% confidence. This is no larger than the observed volume filling factor of the winds from Lyman break galaxies. Previous claims that the observed OVI absorption extends to underdense regions and requires a universal metal enrichment with large volume filling factor, as may be expected from population III star formation at very high redshift, appear not to be warranted.

This is the third paper of a series describing the Asiago-ESO/RASS QSO survey (AERQS), a project aimed at the construction of an all-sky statistically well-defined sample of relatively bright QSOs (B<15) at z<0.3. We present here the clustering analysis of the full spectroscopically identified database (392 AGN). The clustering signal at 0.02<z<0.22 is detected at a 3-4 sigma level and its amplitude is measured to be r_0=8.6\pm 2.0 h^{-1} Mpc (in a LambdaCDM model). The comparison with other classes of objects shows that low-redshift QSOs are clustered in a similar way to Radio Galaxies, EROs and early-type galaxies in general, although with a marginally smaller amplitude. The comparison with recent results from the 2QZ shows that the correlation function of QSOs is constant in redshift or marginally increasing toward low redshift. We discuss this behavior with physically motivated models, deriving interesting constraints on the typical mass of the dark matter halos hosting QSOs, M_DMH= 10^{12.7} h^{-1} M_sun (10^{12.0}-10^{13.5}h^{-1} M_sun at 1 sigma confidence level). Finally, we use the clustering data to infer the physical properties of local AGN, obtaining M_BH=2 10^8 h^{-1} M_sun (10^7-3 10^9 h^{-1} M_sun) for the mass of the active black holes, tau_{AGN}= 8 10^6 yr (2 10^{6}-5 10^{7} yr) for their life-time and eta = 0.14 for their efficiency (always for a LambdaCDM model).

We demonstrate that the lyman alpha forest flux power spectrum of ``randomised'' QSO absorption spectra is comparable in shape and amplitude to the flux power spectrum of the original observed spectra. In the randomised spectra a random shift in wavelength has been added to the observed absorption lines as identified and fitted with VPFIT. At 0.03 s/km <k<0.1 s/km the ``3D'' power spectrum of the randomised flux agrees with that of observed spectra within the errors. At larger scales it is still greater than 50 % of that of the observed spectra. At smaller scales the flux power spectrum is dominated by metal lines. Lines of increasing column density contribute to the ``3D'' flux power spectrum at increasingly larger scales. Lines with 13< log (N_HI/cm^{-2})<15 dominate at the peak of the ``3D'' power spectrum while strong absorbers with log (N_HI/cm^{-2})>15 dominate at large scales, k< 0.03 s/km. We further show that a fraction of up to 20 % of the mean flux decrement is contributed by strong absorbers. Analysis of the flux power spectrum which use numerical simulations with too few strong absorption systems calibrated with the observed mean flux are likely to underestimate the inferred rms fluctuation amplitude and the slope of the initial dark matter matter power spectrum. In a combined analysis with other data which constrains the dark matter power spectrum on larger scales this can result in a spurious detection of a running spectral index.

The hierarchical build-up of galactic bulges should lead to the build-up of present-day supermassive black holes by a mixture of gas accretion and merging of supermassive black holes. The tight relation between black hole mass and stellar velocity dispersion is thereby a strong argument that the supermassive black holes in merging galactic bulges do indeed merge. Otherwise the ejection of supermassive black holes by gravitational slingshot would lead to excessive scatter in this relation. At high redshift the coalescence of massive black hole binaries is likely to be driven by the accretion of gas in the major mergers signposted by optically bright QSO activity. If massive black holes only form efficiently by direct collapse of gas in deep galactic potential wells with v_c > 100 km/s as postulated in the model of Kauffmann & Haehnelt (2000) LISA expects to see event rates from the merging of massive binary black holes of about 0.1-1 yr^{-1} spread over the redshift range 0 < z < 5. If, however, the hierarchical build-up of supermassive black holes extends to pre-galactic structures with significantly shallower potential wells event rates may be as high as 10-100 yr^{-1} and will be dominated by events from redshift z > 5.

The tight correlation between black hole mass and velocity dispersion of galactic bulges is strong evidence that the formation of galaxies and supermassive black holes are closely linked. I review the modeling of the joint formation of galaxies and their central supermassive black holes in the context of the hierarchical structure formation paradigm.

Previous modelling has demonstrated that it is difficult to reproduce the SCUBA source counts within the framework of standard hierarchical structure formation models if the sources are assumed to be the high-redshift counterparts of local ultra-luminous infrared galaxies with dust temperatures in the range 40-60 K. Here, we show that the counts are more easily reproduced in a model in which the bulk of the sub-millimetre emission comes from extended, cool (20-25 K) dust in objects with star formation rates of 50-100 solar masses per year. The low temperatures imply typical sizes of about 1 (S_{850}/1mJy})^{1/2} arcsec, a factor two to three larger than those predicted using starburst-like spectral energy distributions. Low dust temperatures also imply a ratio of optical/UV to 850-micron flux which is 30-100 times smaller, for the same optical depth, than expected for objects with a hot, starburst-like SED. This may help explain the small overlap between SCUBA sources and Lyman-break galaxies.

We study the clustering of matter in the intergalactic medium from the Lyman-alpha forests seen in the spectra of pairs or groups of z ~ 2 quasars observed with FORS2 and UVES at the VLT-UT2 Kueyen ESO telescope. The sample consists of five pairs with separations 0.6, 1, 2.1, 2.6 and 4.4 arcmin and a group of four quasars with separations from 2 up to 10 arcmin. This unprecedented data set allows us to measure the transverse flux correlation function for a range of angular scales. Correlations are clearly detectable at separations smaller than 3 arcmin. The shape and correlation length of the transverse correlation function on these scales is in good agreement with those expected from absorption by the photoionized warm intergalactic medium associated with the filamentary and sheet-like structures predicted in CDM-like models for structure formation. At larger separation no significant correlation is detected. Assuming that the absorbing structures are randomly orientated with respect to the line of sight, the comparison of transverse and longitudinal correlation lengths constrains the cosmological parameters (as a modified version of the Alcock & Paczy\'nski test). The present sample is too small to get significant constraints. Using N-body simulations, we investigate the possibility to constrain Ol from future larger samples of QSO pairs with similar separations. The observation of a sample of 30 pairs at 2, 4.5 and 7.5 arcmin should constrain the value of Ol at +- 15 % (2 sigma level). We further use the observed spectra of the group of four quasars, to search for underdense regions in the intergalactic medium. We find a quasi-spherical structure of reduced absorption with radius 12.5 h^{-1} Mpc which we identify as an underdense region.

We study the number and interaction rates of supermassive black holes in galactic bulges as predicted by hierarchical models of galaxy formation in which the spheroidal components of galaxies are formed by mergers. In bright ellipticals, the number of events that can eject a central supermassive binary black hole is large. Central binaries must therefore merge in less than a Hubble time -- otherwise there will be too much scatter in the M-sigma relation and too many off-center galactic nuclei. We propose that binary black holes are able to merge during the major gas accretion events that trigger QSO activity in galaxies. If this is the case, less than 10 percent of faint ellipticals and 40 percent of bright ellipticals are predicted to harbour binary black holes with near equal masses at their centres. This binary may be ejected away from the centre of the galaxy or even into intergalactic space in up to 20 percent of the most luminous ellipticals. The number of low mass black holes that can interact with the central object is predicted to be a strong function of galaxy luminosity. In most faint ellipticals, no black holes fall into the centre of the galaxy after the last major gas accretion event, but in the most luminous ellipticals, an average of 10 low mass black holes interact with the central supermassive object after this time. It is expected that stars will be ejected from galaxy cores as these low mass-ratio binaries harden. Multiple black holes in galactic bulges thus provide a natural explanation for the strong systematic trends in the observed central density profiles of ellipticals as a function of luminosity.

This is the second paper of a series describing the Asiago-ESO/RASS QSO survey, a project aimed at the construction of an all-sky statistically well-defined sample of very bright QSOs (B_J < 15). Such a survey is required to remove the present uncertainties about the properties of the local QSO population and constitutes an homogeneous database for detailed evolutionary studies of AGN. We present here the complete Southern Sample, which comprises 243 bright (12.60 < B_J < 15.13) QSO candidates at high galactic latitudes (|b_{gal}| > 30^{\circ}). The area covered by the survey is 5660 sq. deg. Spectroscopy for the 137 still unidentified objects has been obtained. The total number of AGN turns out to be 111, 63 of which are new identifications. The properties of the selection are discussed. The completeness and the success rate for this survey at the final stage are 63% and 46%, respectively.

A variety of approximate schemes for modelling the low-density Intergalactic Medium (IGM) in the high-redshift Universe is compared to the results of a large high-resolution hydro-dynamical simulation. These schemes use either an analytical description of the dark matter distribution and the IGM or numerical simulations of the DM distributions combined with different approximate relations between dark matter field and the gas distribution. Schemes based on a filtering of the dark matter distribution with a global Jeans scale result in a rather poor description of the gas distribution. An adaptive filtering which takes into account the density/temperature dependence of the Jeans scale is required. A reasonable description of the gas distribution can be achieved using a fit of the mean relation between the dark matter and gas densities in the hydro-dynamical simulation to relate dark matter and gas distribution. In the hydro-dynamical simulations deviations from this mean relation are correlated with gradients in the dark matter peculiar velocity field indicative of shocks in the gas component. A scheme which takes into account this correlation results in a further improved gas distribution. Such adaptive filtering schemes applied to dark matter simulations will be very well suited for studies of statistical properties of the Lyalpha forest which investigate the IGM and the underlying dark matter distribution and require a large dynamic range and/or an extensive parameter study.

We estimate the rate at which stars are captured by supermassive black holes (BHs) in the centres of bulges and elliptical galaxies assuming that these initially had an isothermal cusp (rho ~ r^{-2} with velocity dispersion sigma_*). If radial orbits can be adequately replenished a BH can grow significantly by capturing stars and particles of the dark halo; the latter could be up to 20--40 percent of the total mass captured. A tight relation M_bh 10^8 (sigma_*/200\kms)^5 between BH hole mass and stellar velocity dispersion arises at the high mass end M_bh >= 10^8 M_sun if these giant BHs grow primarily by swallowing stars. For smaller masses a gradual flattening of the relation and an increase of the scatter are expected due to tidal disruption of stars. At moderate redshifts when bright galaxies may still possess isothermal cusps tidal disruption flares may be as frequent as 10^{-4} -- 10^{-2} yr^{-1} per galaxy, a few percent of the supernovae rate. Efficient randomization of the orbits can occur when the dense cores of infalling satellites sink to the sphere of influence of the BH by dynamical friction. Efficient replenishment of the loss cone of radial orbits can also explain why the supermassive binary BHs expected in hierarchically merging galaxies do generally coalesce, as suggested by the observed relation between BH mass and the inferred mass of stars ejected from an isothermal cusp.

(abridged) We present an effective implementation of analytical calculations of the Lyalpha opacity distribution of the Intergalactic Medium (IGM) along multiple lines of sight (LOS) to distant quasars in a cosmological setting. This method fully accounts for the expected correlations between LOS and the cosmic variance in the large-scale modes of the dark matter distribution. Strong correlations extending up to ~300 kpc (proper) and more are found at redshift z~2->3, in agreement with observations. These correlations are investigated using the cross-correlation coefficient and the cross-power spectrum of the flux distribution along different LOS and by identifying coincident absorption features as fitted with a Voigt profile fitting routine. The cross-correlation coefficient between the LOS can be used to constrain the shape-parameter Gamma of the power spectrum if the temperature and the temperature density relation of the IGM can be determined indepedently. We also propose a new technique to recover the 3D linear dark matter power spectrum by integrating over 1D flux cross-spectra which is complementary to the usual `differentiation' of 1D auto-spectra. The cross-power spectrum suffers much less from errors uncorrelated in different LOS, like those introduced by continuum fitting. Investigations of the flux correlations in adjacent LOS should thus allow to extend studies of the dark matter power spectrum with the Lyalpha forest to significantly larger scales than is possible with flux auto-power spectra. 30 pairs with separation of 1-2 arcmin should be sufficient to determine the 1D cross-spectrum at scales of 60 Mpc/h with an accuracy of about 30% if the error is dominated by cosmic variance.

We study the cross-correlation between quasars and galaxies by embedding models for the formation and evolution of the two populations in cosmological N-body simulations. We adopt the quasar evolution model of Kauffmann & Haehnelt (2000), in which supermassive black holes are formed and fuelled during major mergers. We define the `bias' parameter b_QG as the ratio of the cross-correlation function xi_QG to the galaxy auto-correlation function xi_GG. On scales larger than 1 h^{-1} Mpc, the values of b_QG predicted by our models at low redshift depend very little on galaxy selection. They measure the characteristic mass of the dark matter halos that host quasars and can be used to estimate the typical quasar lifetime. In current redshift surveys, such measurements will constrain the lifetimes of low z quasars more accurately than measurements of the quasar auto-correlation function, because galaxies have much higher space densities than quasars. On scales smaller than 1 h^{-1} Mpc, the main contribution to xi_QG comes from quasar/galaxy pairs in the same dark matter halo. The amplitude of xi_QG depends both on the location of the host galaxy and on the density profile of other galaxies within the halo. As a result, measurements on these scales yield information about the processes responsible for fuelling super- massive black holes. At high redshifts, our models predict that quasars of fixed luminosity are located in less massive halos than at low redshift. They are therefore less biased relative to galaxies of given luminosity or stellar mass. We have used the simulations to calculate the evolution of the quasar auto- correlation function. We find that models with quasar lifetimes in the range 10^6-10^7 years provide a good match to the results of the 2dF QSO survey.

We discuss implications of a stellar-dominated UV background at high redshifts. The composite spectrum of 29 LBGs evaluated by Steidel etal (2000) at <z>=3.4 can be well fit by a stellar population with ongoing star formation, a Salpeter initial mass function, modest or negligible dust reddening, and no intrinsic HI photoelectric absorption. Fading starbursts in which star formation has ceased for 10^7 yr or more cannot reproduce the observed flux shortward of 1 Ryd.The escape fraction of H-ionizing photons must be close to 100 percent for the observed sample of LBGs. The spectrum of ionizing photons produced by a stellar population with ongoing star formation is similar to that of QSOs between 1 and 3 Ryd, becomes softer between 3 and 4 Ryd and drops sharply shortward of 4 Ryd. A galaxy-dominated UV background appears inconsistent with the observed HeII/HI opacity ratio at z=2.4, but might be able to explain the SiIV/CIV abundances measured at z>3 in absorption spectra. A scenario may be emerging where star-forming galaxies reionized intergalactic hydrogen at $z>6$ and dominate the 1 Ryd metagalactic flux at z>3, with quasi-stellar sources taking over at low redshifts. If the large amplitude of the H-ionizing flux estimated by Steidel etal is correct, hydrodynamical simulations of structure formation in the IGM require a baryon density (to explain the observed Lyalpha opacity in absorption spectra) which is similar to or larger than that favoured by recent CMB experiments, and is inconsistent with standard nucleosynthesis values.

For the low density filamentary and sheet-like structures in the warm ($\sim 10^4$ to $\sim 10^6$ K) IGM predicted by numerical simulations the resonant line scattering of X-ray background (XRB) photons by He and H-like ions of heavy elements can exceed the ``local'' thermal emission by a factor of a few or more. Due to the conservative nature of scattering this resonantly scattered radiation can only be identified if a significant fraction of the XRB is resolved and removed. While the combined spectrum of the resolved sources will contain X-ray absorption features, the residual background will contain corresponding emission features with the same intensity. At the relevant densities and temperatures the lines of He and H-like oxygen at 0.57 and 0.65 keV are most promising. These lines (which have a typical width of $\sim$ 1-2 eV) may contain up to 50% of the total 0.5-1 keV emission of the filament. For a nearby ($z \la 0.1$) filament with a Thomson optical depth of $10^{-4}$ XMM should detect about 200 photons in the OVII line during a $10^5$ s exposure if the metallicity of the gas is as large as observed in galaxy clusters. On average up to a few percent of the soft XRB could be resonantly scattered by this phase of the IGM and resonantly scattered photons should account for a significant fraction of the truly diffuse background at low energies. Close to bright X-ray sources like galaxy clusters or AGN the flux of scattered radiation will be further enhanced. Off-line blazars are the most promising illuminating sources. The scattered emission from AGN may also constrain the duration of the active phase of these objects.

Recent work has demonstrated that there is a tight correlation between the mass of a black hole and the velocity dispersion of the bulge of its host galaxy. We show that the model of Kauffmann & Haehnelt, in which bulges and supermassive black holes both form during major mergers, produces a correlation between M_bh and sigma with slope and scatter comparable to the observed relation. In the model, the M_bh - sigma relation is significantly tighter than the correlation between black hole mass and bulge luminosity or the correlation between bulge luminosity and velocity dispersion. There are two reasons for this: i) the gas masses of bulge progenitors depend on the velocity dispersion but not on the formation epoch of the bulge, whereas the stellar masses of the progenitors depend on both; ii) mergers between galaxies move black holes along the observed M_bh - sigma relation, even at late times when the galaxies are gas-poor and black holes grow mainly by merging of pre-existing black holes. We conclude that the small scatter in the observed M_bh - sigma relation is consistent with a picture in which bulges and black holes form over a wide range in redshift.

We investigate the dependence of QSO Lyman-alpha absorption features on the temperature of the absorbing gas and on the amplitude of the underlying dark matter fluctuations. We use high-resolution hydrodynamic simulations in cold dark matter dominated cosmological models. In models with a hotter intergalactic medium (IGM), the increased temperature enhances the pressure gradients between low and high density regions and this changes the spatial distribution and the velocity field of the gas. Combined with more thermal broadening, this leads to significantly wider absorption features in hotter models. Cosmological models with little small scale power also have broader absorption features, because fluctuations on the scale of the Jeans length are still in the linear regime . Consequently, both the amplitude of dark matter fluctuations on small scales and thermal smoothing affect the flux decrement distribution in a similar way. However, the b-parameter distribution of Voigt profile fits, obtained by deblending the absorption features into a sum of thermally broadened lines, is largely independent of the amount of small scale power, but does depend strongly on the IGM temperature. The same is true for the two-point function of the flux and for the flux power spectrum on small scales. These three flux statistics are thus sensitive probes of the temperature of the IGM. We compare them computed for our models and obtained from a HIRES spectrum of the quasar Q1422+231 and conclude that the IGM temperature at z=3.25 is fairly high, T_0>= 15000K.

We discuss constraints on the assembly history of supermassive black holes from the observed remnant black holes in nearby galaxies and from the emission caused by accretion onto these black holes. We also summarize the results of a specific model for the evolution of galaxies and their central black holes which traces their hierachical build-up in CDM-like cosmogonies. The model assumes (i) that black holes, ellipticals and starburts form during major mergers of galaxies (ii) that the gas fraction in galaxies decreases with decreasing redshift (iii) that the optical bright phase of a QSO lasts for about 10^7 years. The model succesfully reproduces the evolution of cold gas as traced by damped damped Lyman alpha systems, the evolution of optically bright QSOs, the remnant black hole mass distribution and the host galaxy luminosities of QSOs.

We combine predictions for several hierarchical cosmogonies with observational evidence on damped Lyman alpha systems to establish a correspondence between the high redshift galaxy population and the properties of damped Lyman alpha systems. We assume that high redshift galaxies and damped Lyman alpha systems are hosted by the same dark matter halos and require consistency between the predicted halo space density, the rate of incidence and the velocity width distribution of damped Lyman alpha systems, and the observed galaxy luminosity function at the bright end. We arrive at the following results: (1) predicted impact parameters between the damped absorption system and the luminous part of the absorbing galaxy are expected to be very small (0.3 - 1arcsec) for most galaxies; (2) luminosities of galaxies causing damped absorption are generally fainter than m_R = 25 and damped Lyman alpha systems are predicted to sample preferentially the outer regions of galaxies at the faint end of the galaxy luminosity function at high redshift. Therefore, DLAS should currently provide the best probe of the progenitors of normal present-day galaxies.

A model of the density distribution in the intergalactic medium, motivated by that found in numerical simulations, is used to demonstrate the effect of a clumpy IGM and discrete sources on the reionization of the universe. In an inhomogeneous universe reionization occurs outside-in, starting in voids and gradually penetrating into overdense regions. Reionization should not be sudden but gradual, with a continuous rise of the photon mean free path over a fair fraction of the Hubble time as the emissivity increases. We show that a hydrogen Gunn-Peterson trough should be present at $z\simeq 6$ unless the emissivity increases with redshift at $z>4$. However, the epoch of overlap of cosmological \hii regions could have occurred at a higher redshift if sources of low luminosity reionized the IGM; the Gunn-Peterson trough at $z\sim 6$ would then appear because even the most underdense voids have a large enough neutral fraction in ionization equilibrium to be optically thick to \lya photons. Cosmological \hii regions near the epoch of overlap can produce gaps of transmitted flux only if luminous quasars contributed to the reionization. Despite the clumpiness of the matter distribution, recombinations are not very important during the reionization of hydrogen because the high density gas is not ionized until a late time. We show that the \heii reionization was most likely delayed relative to the hydrogen reionization, but should be completed by $z\sim 3$, the redshift where observations are available. The reported large optical depth fluctuations of \heii are probably not due to an incomplete \heii reionization, but arise from a combination of density fluctuations and the variations in the intensity of the ionizing background due to luminous QSO's.

We incorporate a simple scheme for the growth of supermassive black holes into semi-analytic models that follow the formation and evolution of galaxies in a cold dark matter dominated Universe. We assume that supermassive black holes are formed and fuelled during major mergers. If two galaxies of comparable mass merge, their central black holes coalesce and a few percent of the gas in the merger remnant is accreted by the new black hole over a timescale of a few times 10^7 years. With these simple assumptions, our model not only fits many aspects of the observed evolution of galaxies, but also reproduces quantitatively the observed relation between bulge luminosity and black hole mass in nearby galaxies, the strong evolution of the quasar population with redshift and the relation between the luminosities of nearby quasars and those of their host galaxies. The strong decline in the number density of quasars from z=2 to z=0 is due to the combination of three effects: i) a decrease in the merging rate, ii) a decrease in the amount of cold gas available to fuel black holes, and iii) an increase in the timescale for gas accretion. In a LCDM cosmology the predicted decline in the total content of cold gas in galaxies is consistent with that inferred from observations of damped Lyman-alpha systems. Our results strongly suggest that the evolution of supermassive black holes, quasars and starbursts is inextricably linked to the hierarchical build-up of galaxies.

The real-space optical depth distribution along the line of sight to the QSO Q1422+231 is recovered from two HIRES spectra. The first two moments of the truncated optical depth distribution are used to constrain the density fluctuation amplitude of the intergalactic medium (IGM). The $rms$ of the IGM density at $z\approx 3.25$ estimated from the first spectrum is $\sigma = \sqrt{\exp{[(1.8 \pm 0.27)^2/\alpha^2]}-1}$, with $1.56 <\alpha <2$ for plausible reionization histories. This corresponds to $0.9 \la \sigma \la 2.1$ with $\sigma(\alpha =1.7)= 1.44\pm 0.3 $. The values obtained from the second spectrum are higher by $\approx 20 %$. If the IGM density traces the dark matter (DM) as suggested by numerical simulations we have measured the fluctuation amplitude of the DM density at an effective Jeans scale of about a hundred to two hundred (comoving) kpc. For CDM-like power spectra the amplitude of dark matter fluctuations on these small scales depends on the cosmological density parameter $\Omega$. For power spectra normalized to reproduce the space density of present-day clusters and with a slope parameter of $\Gamma=0.21$ consistent with the observed galaxy power spectrum, the inferred $\Omega$ can be expressed as: $\Omega= 0.61 (\alpha/1.7)^{1.3}(x_{_{\rm J}} /0.62)^{-0.6}$ for a flat universe, and $\Omega= 0.91(\alpha/1.7)^{1.3} (x_{_{\rm J}}/0.62)^{-0.7}$ for a $\lambda=0$ universe. $x_{_{\rm J}}$ is the effective Jeans scale in (comoving) $\mpc$. Based on a suit of detailed mock spectra the 1-$\sigma$ error is $\approx 25 %$. The estimates increase with increasing $\Gamma$. For the second spectrum we obtain 15% lower values.

The thermal history of the intergalactic medium (IGM) after reionization is to a large extent determined by photoheating. Here we demonstrate that calculations of the photoheating rate which neglect radiative transfer effects substantially underestimate the energy input during and after reionization. The neglect of radiative transfer effects results in temperatures of the IGM which are too low by a factor of two after HeII reionization. We briefly discuss implications for the absorption properties of the IGM and the distribution of baryons in shallow potential wells.

The growth of supermassive black holes by merging and accretion in hierarchical models of galaxy formation is studied by means of Monte Carlo simulations. A tight linear relation between masses of black holes and masses of bulges arises if if the mass accreted by supermassive black holes scales linearly with the mass forming stars and if the redshift evolution of mass accretion tracks closely that of star formation. Differences in redshift evolution between black hole accretion and star formation introduce considerable scatter in this relation. A non-linear relation between black hole accretion and star formation results in a non-linear relation between masses of remnant black holes and masses of bulges. The relation of black hole mass to bulge luminosity obseved in nearby galaxies and its scatter are reproduced reasonably well by models in which black hole accretion and star formation are linearly related but do not track each other in redshift. This suggests that a common mechanism determines the efficiency for black hole accretion and the efficiency for star formation, especially for bright bulges.

Some issues relevant for the formation of supermassive black holes are discused and estimates of the event rates for the emission of gravitational waves by coalescing supermassive black hole binaries are given. The models take into account recent improvements in our knowledge of galaxy and star formation in the high-redshift universe. Estimated event rates range from a few to a hundred per year. Typical events will occur at redshift three or larger in galaxies lying at the (very) faint end of the luminosity function at these redshifts.

A method for the recovery of the real space line-of-sight mass density field from Lyman absorption in QSO spectra is presented. The method makes use of a Lucy-type algorithm for the recovery of the HI density. The matter density is inferred from the HI density assuming that the absorption is due to a photoionized intergalactic medium which traces the mass distribution as suggested by recent numerical simulations. Redshift distortions are corrected iteratively from a simultaneous estimate of the peculiar velocity. The method is tested with mock spectra obtained from N-body simulations. The density field is recovered reasonably well up to densities where the absorption features become strongly saturated. The method is an excellent tool to study the density probability distribution and clustering properties of the mass density in the (mildly) non-linear regime. Combined with redshift surveys along QSO sightlines the method will make it possible to relate the clustering of high-redshift galaxies to the clustering of the underlying mass density. We further show that accurate estimates for \Omega_{bar}h^2)^2 J^{-1} H(z)^{-1} and higher order moments of the density probability function can be obtained despite the missing high density tail of the density distribution if a parametric form for the probability distribution of the mass density is assumed.

We demonstrate that the luminosity function of the recently detected population of star-forming galaxies and the QSO luminosity function at z=3 can be matched with the mass function of dark matter haloes predicted by hierarchical cosmogonies for lifetimes of optically bright QSOs in the range 10^6 to 10^8 yr. We suggest that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form M_bh ~(v_halo)^5 and a typical duration of the optically bright QSO phase of a few times 10^7 yr. The high integrated mass density of black holes inferred from recent black hole mass estimates in nearby galaxies may indicate that the overall efficiency of supermassive black holes for producing blue light is smaller than previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion at far above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background and that the ionizing background might be progressively dominated by stars rather than QSOs at higher redshift.

The Doppler parameter distribution of Lyman-alpha absorption is calculated for a set of different reionization histories. The differences in temperature between different reionization histories are as large as a factor three to four depending on the spectrum of the ionizing sources and the redshift of helium reionization. These temperature differences result in observable differences in the Doppler parameter distribution. Best agreement with the observed Doppler parameter distribution between redshift two and four is found if hydrogen and helium are reionized simultaneously at or before redshift five with a quasar-like spectrum.

We investigate the nature of the physical structures giving rise to damped Lyman alpha absorption systems (DLAS) at high redshift. In particular, we examine the suggestion that rapidly rotating large disks are the only viable explanation for the characteristic observed asymmetric profiles of low ionization absorption lines. We demonstrate using hydrodynamic simulations of galaxy formation in a cosmological context that irregular protogalactic clumps can reproduce the observed velocity width distribution and asymmetries of the absorption profiles equally well. The velocity broadening in the simulated clumps is due to a mixture of rotation, random motions, infall and merging. The observed velocity width correlates with the virial velocity of the dark matter halo of the forming protogalactic clump (v_{wid} ~ 0.6 times v_{vir} for the median values with a large scatter of order a factor two between different lines-of-sight). The typical virial velocity of the halos required to give rise to the DLAS population is about 100 km/s and most standard hierarchical structure formation scenarios can easily account even for the largest observed velocity widths. We conclude that the evidence that DLAS at high redshift are related to large rapidly rotating disks with v_circ >= 200 km/s is not compelling.

In hierarchical cosmogonies at redshift three the matter content of a typical present-day galaxy is dispersed over several individual clumps embedded in sheet-like structures, often aligned along filaments. We have used hydrodynamical simulations to investigate the spatial distribution and absorption properties of metal enriched gas in such regions of ongoing galaxy formation. Strong CII and SiIV absorption occurs for lines-of-sight passing the densest regions close to the center of the protogalactic clumps, while CIV is a good tracer of the prominent filamentary structures and OVI becomes the strongest absorption feature for lines-of-sight through low-density regions far away from fully-collapsed objects. Observed column density ratios can be well reproduced if a mean metallicity [Z/H] = -2.5, relative abundances as found in metal-poor stars, a UV background with intensity J_{-22} = 3 at the Lyman limit, and either a power law spectrum (alpha = -1.5) or the spectral shape proposed by Haardt & Madau (1996) are assumed. The observed scatter in [C/H] is about a magnitude larger than that in the simulations suggesting inhomogeneous metal distribution. The detailed small and large scale kinematics observed in CIV and HI are also well reproduced.

Estimates of the linear extent of heavy-element absorption systems along the line-of-sight to a QSO often assume that the cloud is photoionized and that the temperature takes the equilibrium value where photo-heating balances line cooling. We show that rather small deviations from this photoionization equilibrium temperature caused by additional heating processes will lead to an overestimate of the neutral hydrogen fraction and thus to an underestimate of the thickness of the absorber by about two orders of magnitude. Such temperature deviations are indicated both by observations and numerical simulations. This interpretation reconciles the discrepancy between the rather small extent of heavy-element-absorption systems parallel to the line-of-sight obtained from a standard photoionization analysis and the much larger transverse sizes estimates inferred from the observation of common absorption in the spectra of close quasar pairs.

The recently approved COBRAS/SAMBA cosmic microwave background mission will also have a major impact on measurements of the Sunyaev-Zeldovich effect. The frequency of three of the channels (150 GHz, 217 GHz, 353 GHz) are chosen to optimize measurements of the thermal and kinetic SZ effect mainly caused by hot gas in clusters of galaxies. Estimates of the number of detected clusters are somewhat uncertain due to incomplete knowledge of the gas mass function of clusters and the distribution of hot gas at large radii. Straightforward interpolation of X-ray observations of the 200 X-ray brightest clusters gives a firm lower limit of approximately 3000 detected clusters of which roughly half should be resolved. If the gas distribution of clusters should turn out to be favourable and/or the cluster evolution with redshift is weak, this number could be higher by a factor of up to ten. Using an optimal filtering technique a measurement of the bulk peculiar velocity of a sample of about 200 of the resolved clusters will be possible with a 1-sigma error of 100 - 200 km/s.

In a hierarchical cosmogony present-day galaxies build up by continuous merging of smaller structures. At a redshift of three the matter content of a typical present-day galaxy is split into about ten individual protogalactic clumps. Numerical simulations show that these protogalactic clumps have a typical distance of about 100 kpc, are embedded in a sheet-like structure and are often aligned along filaments. Artificial QSO spectra were generated from hydrodynamical simulations of such regions of ongoing galaxy formation. The metal and hydrogen absorption features in the artificial spectra closely resemble observed systems over a wide range in HI column density. Detailed predictions of the column density as a function of impact parameter to protogalactic clumps are presented for HI,CII,CIV,SiIV,NV and OVI. The expected correlations between column densities of different species and their role in understanding the physical properties of the gas from which galaxies form are briefly discussed. The model is able to explain both high-ionization multi-component heavy-element absorbers and damped Lyman alpha systems as groups of small protogalactic clumps.

We investigate the heavy element QSO absorption systems caused by gas condensations at high redshift which evolve into galaxies with circular velocity of 100 to 200 km/s at the present epoch. Artificial QSO spectra were generated for a variety of lines-of-sight through regions of the universe simulated with a hydrodynamics code. The CIV and HI absorption features in these spectra closely resemble observed CIV and HI absorption systems over a wide range in column density. CIV absorption complexes with multiple-component structure and velocity spreads up to about 600 km/s are found. The broadest systems are caused by lines-of-sight passing through groups of protogalactic clumps with individual velocity dispersions of less than 150 km/s aligned along filamentary structures. The temperature of most of the gas does not take the photoionization equilibrium value. This invalidates density and size estimates derived from thermal equilibrium models. Consequences for metal abundance determinations are briefly discussed. We predict occasional exceptionally large ratios of CIV to HI column density (up to a third) for lines-of-sight passing through compact halos of hot gas with temperature close to 3 10^5 K. Our model may be able to explain both high-ionization multi-component heavy-element absorbers and damped Lyman alpha systems as groups of small protogalactic clumps.

The large sizes of Lyman alpha ``clouds'' inferred from coincident absorption in the spectrum of close quasar pairs suggests that these are transient flattened structures of small overdensity. It is argued that the observed absorbers should be preferentially located in underdense regions of the universe and should typically expand faster than the Hubble flow.

We have investigated the possibility of inferring peculiar velocities for clusters of galaxies from the Doppler shift of scattered cosmic microwave background (CMB) photons. We find that if the core radius of the gas distribution or the beam size of the instrument is larger than 3-7 arcminutes, then the maximum attainable signal-to-noise ratio is determined by confusion with primary fluctuations. For smaller angular scales, ``cosmic confusion'' is less important and instrumental noise and/or foreground emission will be the limiting factor. For a cluster with the optical depth of the Coma cluster and for an optimal filtering technique, typical one-sigma errors span the wide range from 400 to 1600 km/s, depending on the cosmological model, the resolution of the instrument and the core radius of the cluster. The results have important implications for the design of future high-resolution surveys of the CMB. Individual peculiar velocities will be measurable only for a few fast moving clusters at intermediate redshift unless cosmic fluctuations are smaller than most standard cosmological scenarios predict. However, a reliable measurement of bulk velocities of ensembles of X-ray bright clusters will be possible on very large scales (100-500 Mpc/h).

Estimates of Omega_baryon from primordial nucleosynthesis together with standard assumptions about the ionization state of low column density Lyman alpha forest clouds can be used to determine an upper limit for the cloud thickness along the line of sight. This upper limit provides significant constraints on the axis ratio of the absorbers and on the overall fraction of the baryonic matter contained in these objects if the recently measured large values for their transversal size are representative. The absorbers have to be considerably flattened structures and may possibly contain a substantial fraction of the baryonic matter at high redshift.

The role of momentum transfer (``radiation pressure'') due to an internal source of ionizing radiation for the formation of baryonic structures is investigated. Fully-ionized self-gravitating gaseous objects can be radiation-pressure supported on a characteristic length scale $D_{\rm rp} \sim 100 \pc - 3 \kpc$. On smaller scales momentum transfer due to ionizing radiation will be the dominant force and for spherical collapsing objects of mass $\la 10^{10}\Msol$ a complete bounce is possible. A population of massive stars and/or accretion onto a central compact object are natural sources of ionizing radiation in newly-forming baryonic structures. Radiation pressure is therefore likely to play an important role for the dynamical and thermal evolution of the intergalactic medium and the mass-to-light ratio of small galaxies. The effect will be especially important for hierarchical cosmogonies where galactic structures build up by merging of smaller objects. Radiation pressure due to ionizing radiation might furthermore be responsible for substructures of size $D_{\rm rp}$ in the stellar component of large galaxies.

Supermassive black holes are investigated as possible sources for low-frequency bursts of gravity waves. The event rate for `known' supermassive black holes at intermediate and high redshifts, inferred from the quasar luminosity function, is low $\sim 0.1 \yr^{-1}$. A space-based interferometer could therefore only see several events per year from supermassive black holes if an additional population of supermassive black holes existed and emitted gravitational waves efficiently. These might reside in the population of dwarf galaxies or in a transient population of small dark matter haloes that have mostly merged into larger haloes hosting the galaxies seen today. The proposed space-based gravitational-wave interferometer LISA/SAGITTARIUS should detect most gravitational-wave events involving supermassive black holes above $10^{4} \Msol$ out to redshifts of $z \sim 100$.