The Lyman-alpha forest and WMAP year three

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).Comment: 5 pages, 4 figs, 2 tables. MNRAS letters in pres

Cosmological and astrophysical parameters from the SDSS flux power spectrum and hydrodynamical simulations of the Lyman-alpha forest

(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.Comment: 16 pages, 7 figures, 3 tables. Minor changes to match the accepted version. MNRAS, in pres

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

Can sterile neutrinos be ruled out as warm dark matter candidates?

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.Comment: 4 pages, 3 figures, matches published 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

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

Inferring the dark matter power spectrum from the Lyman-alpha forest in high-resolution QSO absorption spectra

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.Comment: 13 pages, 9 figures, 6 tables. Very minor changes to match the accepted version. MNRAS in pres

Constraints on the Primordial Power Spectrum from High Resolution Lyman-alpha Forest Spectra and WMAP

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.Comment: 6 pages, 4 Figures, 1 Table. One figure added and second order extension of the standard slow-roll approximation included. Main results unchanged. MNRAS in pres

An improved cosmological bound on the thermal axion mass

Relic thermal axions could play the role of an extra hot dark matter component in cosmological structure formation theories. By combining the most recent observational data we improve previous cosmological bounds on the axion mass m_a in the so-called hadronic axion window. We obtain a limit on the axion mass m_a < 0.42eV at the 95% c.l. (m_a < 0.72eV at the 99% c.l.). A novel aspect of the analysis presented here is the inclusion of massive neutrinos and how they may affect the bound on the axion mass. If neutrino masses belong to an inverted hierarchy scheme, for example, the above constraint is improved to m_a < 0.38eV at the 95% c.l. (m_a < 0.67eV at the 99% c.l.). Future data from experiments as CAST will provide a direct test of the cosmological bound.Comment: 5 Pages, 3 Figure

Testing the accuracy of the Hydro-PM approximation in numerical simulations of the Lyman-alpha forest

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.Comment: 13 pages, 10 figures, 1 table. MNRAS. Discussion clarified, one more figur