The impact of spatial fluctuations in the ultra-violet background on intergalactic carbon and silicon

Spatial inhomogeneities in the spectral shape of the ultra-violet background (UVB) at the tail-end of HeII reionisation are thought to be the primary cause of the large fluctuations observed in the HeII to HI Ly-a forest optical depth ratio, tau_HeII/tau_HI, at z~2-3. These spectral hardness fluctuations will also influence the ionisation balance of intergalactic metals; we extract realistic quasar absorption spectra from a large hydrodynamical simulation to examine their impact on intergalactic SiIV and CIV absorbers. Using a variety of toy UVB models, we find that while the predicted spatial inhomogeneities in spectral hardness have a significant impact on tau_HeII/tau_HI, the longer mean free path for photons with frequencies above and below the HeII ionisation edge means these fluctuations have less effect on the SiIV and CIV ionisation balance. Furthermore, UVB models which produce the largest fluctuations in specific intensity at the HeII ionisation edge also have the softest ionising spectra, and thus result in photo-ionisation rates which are too low to produce significant fluctuations in the observed tau_SiIV/tau_CIV. Instead, we find spatial variations in the IGM metallicity will dominate any scatter in tau_SiIV/tau_CIV. Our results suggest that observational evidence for homogeneity in the observed tau_SiIV/tau_CIV distribution does not rule out the possibility of significant fluctuations in the UVB spectral shape at z~2-3. On the other hand, the scatter in metallicity inferred from observations of intergalactic CIV and SiIV absorption at z~2-3 using spatially uniform ionisation corrections is likely intrinsic, and therefore provides a valuable constraint on intergalactic metal enrichment scenarios at these redshifts.Comment: 13 pages, 7 figures, accepted to MNRA

The Lyman-alpha forest opacity and the metagalactic hydrogen ionization rate at z~2-4

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.Comment: 13 pages, 7 figures, accepted to MNRAS, minor changes to submitted versio

The power spectrum of the flux distribution in the Lyman-alpha forest of a Large sample of UVES QSO Absorption Spectra (LUQAS)

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.Comment: 18 pages, 9 PS figures, 6 tables. Note that the k-values of the 1D flux power spectrum had been erroneously shifted by half a bin size (in log k) in the previous version. All the other results are unaffected. New tables can be found at http://www.ast.cam.ac.uk/~rtnigm/luqas.ht

Probing 3-D matter distribution at z~2 with QSO multiple lines of sight

We investigate the 3-D matter distribution at z~2 with high resolution (R ~ 40000) spectra of QSO pairs and groups obtained with the UVES spectrograph at ESO VLT. Our sample is unique for the number density of objects and the variety of separations, between 0.5 and 7 proper Mpc. We compute the real space cross-correlation function of the Lyman-alpha forest transmitted fluxes. There is a significant clustering signal up to ~2 proper Mpc, which is still present when absorption lines with high column density (log N > 13.8) are excluded.Comment: Poster paper presented at the IAU Colloquium #199 on "Probing Galaxies through Quasar Absorption Lines" held in Shanghai, China from March 14th to 18th, 200

Detecting X-ray filaments in the low redshift Universe with XEUS and Constellation-X

We propose a possible way to detect baryons at low redshifts from the analysis of X-ray absorption spectra of bright AGN pairs. A simple semi-analytical model to simulate the spectra is presented. We model the diffuse warm-hot intergalactic medium (WHIM) component, responsible for the X-ray absorption, using inputs from high-resolution hydro-dynamical simulations and analytical prescriptions. We show that the number of OVII absorbers per unit redshift with column density larger than $10^{13.5}$ cm$^{-2}$ - corresponding to an equivalent width of $\sim$ 1 km/s - which will be possibly detectable by {\it XEUS}, is \magcir 30 per unit redshift. {\it Constellation-X} will detect $\sim 6$ OVII absorptions per unit redshift with an equivalent width of 10 km/s. Our results show that, in a $\Lambda$CDM Universe, the characteristic size of these absorbers at $z\sim 0.1$ is $\sim 1$ $h^{-1}$ Mpc. The filamentary structure of WHIM can be probed by finding coincident absorption lines in the spectra of background AGN pairs. We estimate that at least 20 AGN pairs at separation \mincir 20 arcmin are needed to detect this filamentary structure at a 3$\sigma$ level. Assuming observations of distant sources using {\it XEUS} for exposure times of 500 ksec, we find that the minimum source flux to probe the filamentary structure is $\sim 2\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, in the 0.1-2.4 keV energy band. Thus, most pairs of these extragalactic X-ray bright sources have already been identified in the {\it ROSAT} All-Sky Survey. Re-observation of these objects by future missions could be a powerful way to search for baryons in the low redshift Universe.Comment: 18 pages, 10 Figures. Two figures added, Sections 2 and 3 expanded. More optimistic results for Constellation-X. Accepted by MNRA

Tomography of the intergalactic medium with Ly-alpha forests in close QSO pairs

We study the three-dimensional distribution of non virialised matter at z~2 using high resolution spectra of QSO pairs and simulated spectra drawn from cosmological hydrodynamical simulations. We have collected the largest sample of QSO pairs ever observed with UVES at the ESO-VLT, with angular separations between ~1 and 14 arcmin. The observed correlation functions of the transmitted flux in the HI Lyman alpha forest along and transverse to the lines of sight are in good agreement implying that the distortions in redshift space due to peculiar velocities are small. The clustering signal is significant up to velocity separations of ~200 km/s, or about 3 h^{-1} comoving Mpc. The regions at lower overdensity (rho/ < 6.5) are still clustered but on smaller scales (Delta v < 100 km/s). The observed and simulated correlation functions are compatible at the 3 sigma level. A better concordance is obtained when only the low overdensity regions are selected for the analysis or when the effective optical depth of the simulated spectra is increased artificially, suggesting a deficiency of strong lines in the simulated spectra. We found that also a lower value of the power-law index of the temperature-density relation for the Lyman alpha forest gas improves the agreement between observed and simulated results. If confirmed, this would be consistent with other observations favouring a late HeII reionization epoch (at z~3). We remark the detection of a significant clustering signal in the cross correlation coefficient at a transverse velocity separation Delta v_{\perp} ~500 km/s whose origin needs further investigation.Comment: Accepted for publication in MNRAS, revised version matching the accepted on

Resolving the high redshift Lyman-alpha forest in smoothed particle hydrodynamics simulations

We use a large set of cosmological smoothed particle hydrodynamics (SPH) simulations to examine the effect of mass resolution and box size on synthetic Lya forest spectra at 2 \leq z \leq 5. The mass resolution requirements for the convergence of the mean Lya flux and flux power spectrum at z=5 are significantly stricter than at lower redshift. This is because transmission in the high redshift Lya forest is primarily due to underdense regions in the intergalactic medium (IGM), and these are less well resolved compared to the moderately overdense regions which dominate the Lya forest opacity at z~2-3. We further find that the gas density distribution in our simulations differs significantly from previous results in the literature at large overdensities (\Delta>10). We conclude that studies of the Lya forest at z=5 using SPH simulations require a gas particle mass of M_gas \leq 2x10^5 M_sol/h, which is >8 times the value required at z=2. A box size of at least 40 Mpc/h is preferable at all redshifts.Comment: 5 pages, 5 figures, 2 tables, accepted by MNRA

How cold is cold dark matter? Small scales constraints from the flux power spectrum of the high-redshift Lyman-alpha forest

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.Comment: Some issues clarified (especially resolution related). Conclusions unchanged. Accepted version by PR

High resolution spectroscopy of the three dimensional cosmic web with close QSO groups

We study the three-dimensional distribution of matter at z~2 using high resolution spectra of QSO pairs and simulated spectra drawn from cosmological hydro-dynamical simulations. We present a sample of 15 QSOs, corresponding to 21 baselines of angular separations evenly distributed between ~1 and 14 arcmin, observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) at the European Southern Observatory-Very Large Telescope (ESO-VLT). The observed correlation functions of the transmitted flux in the HI Lya forest transverse to and along the line of sight are in agreement, implying that the distortions in redshift space due to peculiar velocities are relatively small and - within the relatively large error bars - not significant. The clustering signal is significant up to velocity separations of ~300 km/s, corresponding to about 5 h^{-1} comoving Mpc. Compatibility at the 2 sigma level has been found both for the Auto- and Cross-correlation functions and for the set of the Cross correlation coefficients. The analysis focuses in particular on two QSO groups of the sample. Searching for alignments in the redshift space between Lya absorption lines belonging to different lines of sight, it has been possible to discover the presence of a wide HI structures extending over about ten Mpc in comoving space, and give constraints on the sizes of two cosmic under-dense regions in the intergalactic medium.Comment: Accepted by MNRAS, version matching the published on

Fast, large volume, GPU enabled simulations for the Ly-alpha forest: power spectrum forecasts for baryon acoustic oscillation experiments

High redshift measurements of the baryonic acoustic oscillation scale (BAO) from large Ly-alpha forest surveys represent the next frontier of dark energy studies. As part of this effort, efficient simulations of the BAO signature from the Ly-alpha forest will be required. We construct a model for producing fast, large volume simulations of the Ly-alpha forest for this purpose. Utilising a calibrated semi-analytic approach, we are able to run very large simulations in 1 Gpc^3 volumes which fully resolve the Jeans scale in less than a day on a desktop PC using a GPU enabled version of our code. The Ly-alpha forest spectra extracted from our semi-analytical simulations are in excellent agreement with those obtained from a fully hydrodynamical reference simulation. Furthermore, we find our simulated data are in broad agreement with observational measurements of the flux probability distribution and 1D flux power spectrum. We are able to correctly recover the input BAO scale from the 3D Ly-alpha flux power spectrum measured from our simulated data, and estimate that a BOSS-like 10^4 deg^2 survey with ~15 background sources per square degree and a signal-to-noise of ~5 per pixel should achieve a measurement of the BAO scale to within ~1.4 per cent. We also use our simulations to provide simple power-law expressions for estimating the fractional error on the BAO scale on varying the signal-to-noise and the number density of background sources. The speed and flexibility of our approach is well suited for exploring parameter space and the impact of observational and astrophysical systematics on the recovery of the BAO signature from forthcoming large scale spectroscopic surveys.Comment: 16 pages, 11 figures, accepted to MNRA