The nature of dark matter sets the timeline for the formation of first collapsed halos and thus affects the sources of reionization. Here, we consider two different models of dark matter: cold dark matter (CDM) and thermal warm dark matter (WDM), and study how they impact the epoch of reionization (EoR) and its 21-cm observables. Using a suite of simulations, we find that in the WDM scenarios, the structure formation on small scales gets suppressed resulting in a smaller number of low mass dark matter halos compared to the CDM scenario. Assuming that the efficiency of sources in producing ionizing photons remain the same, this leads to a lower number of total ionizing photons produced at any given cosmic time and thus in a delay in the reionization process. We also find visual differences in the neutral hydrogen (HI) topology and in 21-cm maps in case of the WDM compared to the CDM. However, differences in the 21-cm power spectra, at the same neutral fraction, are found to be small. Thus, we focus on the non-Gaussianity in the EoR 21-cm signal, quantified through its bispectrum. We find that the 21-cm bispectra (driven by the HI topology) are significantly different in WDM models when compared with CDM, even for same mass averaged neutral fractions. This establishes that the 21-cm bispectrum is a unique and promising way to differentiate between different dark matter models, and can be used to constrain the nature of the dark matter in the future EoR observations.

Galaxy surveys are an important probe for superimposed oscillations on the primordial power spectrum of curvature perturbations, which are predicted in several theoretical models of inflation and its alternatives. In order to exploit the full cosmological information in galaxy surveys it is necessary to study the matter power spectrum to fully non-linear scales. We therefore study the non-linear clustering in models with superimposed linear and logarithmic oscillations to the primordial power spectrum by running high-resolution dark-matter-only N-body simulations. We fit a Gaussian envelope for the non-linear damping of superimposed oscillations in the matter power spectrum to the results of the N-body simulations for $k \lesssim 0.6\ h/$Mpc at $0 \leq z \leq 5$ with an accuracy below the percent. We finally use this fitting formula to forecast the capabilities of future galaxy surveys, such as Euclid and Subaru, to probe primordial oscillation down to non-linear scales alone and in combination with the information contained in CMB anisotropies.

We constrain and update the bounds on the life-time of a decaying dark matter model with a warm massive daughter particle using the most recent low-redshift probes. We use Supernovae Type-Ia, Baryon Acoustic Oscillations and the time delay measurements of gravitationally lensed quasars. These data sets are complemented by the early universe priors taken from the Cosmic Microwave background. For the maximum allowed fraction of the relativistic daughter particle, the updated bounds on the life-time are found to be $\tau > 9\, \rm{Gyr}$ and $\tau >11\,\rm{Gyr}$ at $95\%$ C.L., for the two-body and many-body decay scenarios, respectively. We also comment on the recent proposal that the current two-body decaying dark matter model can provide resolution for the $H_0$-tension, by contrasting against the standard $\Lambda$CDM model. We infer that the current dark matter decaying scenario is unlikely to alleviate the $H_0$-tension. We find that the decaying dark matter is able to reduce the trend of the decreasing $H_0$ values with increasing lens redshifts observed in the strong lensing dataset.

We study the impact of relativistic effects in the 3-dimensional cross-correlation between Lyman-$\alpha$ forest and quasars. Apart from the relativistic effects, which are dominated by the Doppler contribution, several systematic effects are also included in our analysis (intervening metals, unidentified high column density systems, transverse proximity effect and effect of the UV fluctuations). We compute the signal-to-noise ratio for the Baryonic Oscillation Spectroscopic Survey (BOSS), the extended Baryonic Oscillation Spectroscopic Survey (eBOSS) and the Dark Energy Spectroscopic Instrument (DESI) surveys, showing that DESI will be able to detect the Doppler contribution in a Large Scale Structure (LSS) survey for the first time, with a S/N $>7$ for $r_{\rm min} > 10$ Mpc$/h$, where r$_{\rm min}$ denotes the minimum comoving separation between sources. We demonstrate that several physical parameters, introduced to provide a full modelling of the cross-correlation function, are affected by the Doppler contribution. By using a Fisher matrix approach, we establish that if the Doppler contribution is neglected in the data analysis, the derived parameters will be shifted by a non-negligible amount for the upcoming surveys.

We present a computational framework for "painting" galaxies on top of the Dark Matter Halo/Sub-Halo hierarchy obtained from N-body simulations. The method we use is based on the sub-halo clustering and abundance matching (SCAM) scheme which requires observations of the 1- and 2-point statistics of the target (observed) population we want to reproduce. This method is particularly tailored for high redshift studies and thereby relies on the observed high-redshift galaxy luminosity functions and correlation properties. The core functionalities are written in c++ and exploit Object Oriented Programming, with a wide use of polymorphism, to achieve flexibility and high computational efficiency. In order to have an easily accessible interface, all the libraries are wrapped in python and provided with an extensive documentation. We validate our results and provide a simple and quantitative application to reionization, with an investigation of physical quantities related to the galaxy population, ionization fraction and bubble size distribution.

There is a persistent $H_0$-tension, now at more than $\gtrsim 4\sigma$ level, between the local distance ladder value and the \emph{Planck} cosmic microwave background measurement, in the context of flat $\Lambda$CDM model. We reconstruct $H(z)$ in a cosmological-model-independent way using three low-redshift distance probes including the latest data from baryon acoustic oscillation, Type Ia supernova and four gravitational lensing Time-Delay observations. We adopt general parametric models of $H(z)$ and assume a Gaussian prior on the sound horizon at drag epoch, $r_{\mathrm s}$, from \emph{Planck} measurement. The reconstructed $H_0$ using Pantheon SN Ia and BAO data are consistent with the \emph{Planck} flat $\Lambda$CDM value. When including the GLTD data, $H_0$ increases mildly, yet remaining discrepant with the local measurement at $\sim 2.5\sigma$ level. Our reconstructions being blind to the dark sectors at low redshift, we reaffirm the earlier claims that the Hubble tension is not likely to be solved by modifying the energy budget of the low-redshift universe. We further forecast the constraining ability of future realistic mock BAO data from DESI and GLTD data from LSST, combining which, we anticipate that the uncertainty of the inferred $H_0$ would be improved by $\sim 38\%$, reaching $\sigma_{H_0} \approx 0.56$ uncertainty level.

We consider an Early Dark Energy (EDE) cosmological model, and perform an analysis which takes into account both background and perturbation effects via the parameters $c^{2}_{\rm eff}$ and $c^{2}_{\rm vis}$, representing effective sound speed and viscosity, respectively. By using the latest available data we derive constraints on the amount of dark energy at early times and the present value of the equation of state. Our focus is on the effect that early dark energy has on the Cosmic Microwave Background (CMB) data, including polarization and lensing, in a generalized parameter space including a varying total neutrino mass, and tensor to scalar ratio, besides the 6 standard parameters of the minimal cosmological model. We find that the inclusion of Baryonic Acoustic Oscillations (BAO) data and CMB lensing significantly improves the constraints on the EDE parameters, while other high redshift data like the Quasar Hubble diagram and the Lyman-$\alpha$ forest BAO have instead a negligible impact. We find $\Omega_{\rm eDE} < 0.00391 $ and $w_{0} < -0.9500 $ at the $95 \%$ C.L. for EDE accounting for its clustering through the inclusion of perturbation dynamics. This limit becomes weaker $\Omega_{\rm eDE} < 0.00337$ if perturbations are neglected. The constraints on the EDE parameters are remarkably stable even when $\Sigma m_{\nu}$, and $r$ parameters are varied, with weak degeneracies between $\Omega_{\rm eDE}$ and $r$ or $\Sigma m_{\nu}$. Perturbation parameters are not constrained with current data sets, and tensions between the CMB derived $H_0$ and $\sigma_8$ values and those measured with local probes are not eased. This work demonstrates the capability of CMB probes to constrain the total amount of EDE well below the percent level.

Cosmological studies of large-scale structure have relied on two-point statistics, not fully exploiting the rich structure of the cosmic web. In this paper we show how to capture some of this cosmic web information by using the minimum spanning tree (MST), for the first time using it to estimate cosmological parameters in simulations. Discrete tracers of dark matter such as galaxies, $N$-body particles or haloes are used as nodes to construct a unique graph, the MST, that traces skeletal structure. We study the dependence of the MST on cosmological parameters using haloes from a suite of COLA simulations with a box size of $250\ h^{-1}{\rm Mpc}$, varying the amplitude of scalar fluctuations $\left(A_{\rm s}\right)$, matter density $\left(\Omega_{\rm m}\right)$, and neutrino mass $\left(\sum m_{\nu}\right)$. The power spectrum $P$ and bispectrum $B$ are measured for wavenumbers between $0.125$ and $0.5$ $h{\rm Mpc}^{-1}$, while a corresponding lower cut of $\sim12.6$ $h^{-1}{\rm Mpc}$ is applied to the MST. The constraints from the individual methods are fairly similar but when combined we see improved $1\sigma$ constraints of $\sim 17\%$ ($\sim 12\%$) on $\Omega_{\rm m}$ and $\sim 12\%$ ($\sim 10\%$) on $A_{\rm s}$ with respect to $P$ ($P+B$) thus showing the MST is providing additional information. The MST can be applied to current and future spectroscopic surveys (BOSS, DESI, Euclid, PSF, WFIRST, and 4MOST) in 3D and photometric surveys (DES and LSST) in tomographic shells to constrain parameters and/or test systematics.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

The two-point correlation function of the galaxy distribution is a key cosmological observable that allows us to constrain the dynamical and geometrical state of our Universe. To measure the correlation function we need to know both the galaxy positions and the expected galaxy density field. The expected field is commonly specified using a Monte-Carlo sampling of the volume covered by the survey and, to minimize additional sampling errors, this random catalog has to be much larger than the data catalog. Correlation function estimators compare data-data pair counts to data-random and random-random pair counts, where random-random pairs usually dominate the computational cost. Future redshift surveys will deliver spectroscopic catalogs of tens of millions of galaxies. Given the large number of random objects required to guarantee sub-percent accuracy, it is of paramount importance to improve the efficiency of the algorithm without degrading its precision. We show both analytically and numerically that splitting the random catalog into a number of subcatalogs of the same size as the data catalog when calculating random-random pairs, and excluding pairs across different subcatalogs provides the optimal error at fixed computational cost. For a random catalog fifty times larger than the data catalog, this reduces the computation time by a factor of more than ten without affecting estimator variance or bias.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

We present a comprehensive analysis of atomic hydrogen (HI) properties using a semi-analytical model of galaxy formation and N-body simulations covering a large cosmological volume at high resolution. We examine the HI mass function and the HI density, characterizing both their redshift evolution and their dependence on hosting halo mass. We analyze the HI content of dark matter haloes in the local Universe and up to redshift $z=5$, discussing the contribution of different galaxy properties. We find that different assembly history plays a crucial role in the scatter of this relation. We propose new fitting functions useful for constructing mock HI maps with HOD techniques. We investigate the HI clustering properties relevant for future $21$~cm Intensity Mapping (IM) experiments, including the HI bias and the shot noise level. The HI bias increases with redshift and it is roughly flat on the largest scales probed. The scale dependency is found at progressively larger scales with increasing redshift, apart from a dip feature at $z=0$. The shot-noise values are consistent with the ones inferred by independent studies, confirming that shot-noise will not be a limiting factor for IM experiments. We detail the contribution from various galaxy properties on the HI power spectrum and their relation to the halo bias. We find that HI poor satellite galaxies play an important role at the scales of the 1-halo term. Finally, we present the $21$~cm signal in redshift space, a fundamental prediction to be tested against data from future radio telescopes such as SKA.

We present constraints on neutrino masses, the primordial fluctuation spectrum from inflation, and other parameters of the $\Lambda$CDM model, using the one-dimensional Ly$\alpha$-forest power spectrum measured by Palanque-Delabrouille et al. (2013) from SDSS-III/BOSS, complemented by Planck 2015 cosmic microwave background (CMB) data and other cosmological probes. This paper improves on the previous analysis by Palanque-Delabrouille et al. (2015) by using a more powerful set of calibrating hydrodynamical simulations that reduces uncertainties associated with resolution and box size, by adopting a more flexible set of nuisance parameters for describing the evolution of the intergalactic medium, by including additional freedom to account for systematic uncertainties, and by using Planck 2015 constraints in place of Planck 2013. Fitting Ly$\alpha$ data alone leads to cosmological parameters in excellent agreement with the values derived independently from CMB data, except for a weak tension on the scalar index $n_s$. Combining BOSS Ly$\alpha$ with Planck CMB constrains the sum of neutrino masses to $\sum m_\nu < 0.12$ eV (95\% C.L.) including all identified systematic uncertainties, tighter than our previous limit (0.15 eV) and more robust. Adding Ly$\alpha$ data to CMB data reduces the uncertainties on the optical depth to reionization $\tau$, through the correlation of $\tau$ with $\sigma_8$. Similarly, correlations between cosmological parameters help in constraining the tensor-to-scalar ratio of primordial fluctuations $r$. The tension on $n_s$ can be accommodated by allowing for a running ${\mathrm d}n_s/{\mathrm d}\ln k$. Allowing running as a free parameter in the fits does not change the limit on $\sum m_\nu$. We discuss possible interpretations of these results in the context of slow-roll inflation.

The renewed interest in the possibility that primordial black holes (PBHs) may constitute a significant part of the dark matter has motivated revisiting old observational constraints, as well as developing new ones. We present new limits on the PBH abundance, from a comprehensive analysis of high-resolution, high-redshift Lyman-$\alpha$ forest data. Poisson fluctuations in the PBH number density induce a small-scale power enhancement which departs from the standard cold dark matter prediction. Using a grid of hydrodynamic simulations exploring different values of astrophysical parameters, {we obtain a marginalized upper limit on the PBH mass of $f_{\rm PBH}M_{\rm PBH} \sim 60~M_{\odot}$ at $2\sigma$, when a Gaussian prior on the reionization redshift is imposed, preventing its posterior distribution to peak on very high values, which are disfavoured by the most recent estimates obtained both through Cosmic Microwave Background and Inter-Galactic Medium observations. Such bound weakens to $f_{\rm PBH}M_{\rm PBH} \sim 170~M_{\odot}$, when a conservative flat prior is instead assumed. Both limits significantly improves previous constraints from the same physical observable.} We also extend our predictions to non-monochromatic PBH mass distributions, ruling out large regions of the parameter space for some of the most viable PBH extended mass functions.

Nonlinear objects like halos and voids exhibit a scale-dependent bias on linear scales in massive neutrino cosmologies. The shape of this scale-dependent bias is a unique signature of the neutrino masses, but the amplitude of the signal is generally small, of the order of $f_\nu$, the contribution of neutrinos to the total matter content ($\lesssim 1\%$). In this paper, we demonstrate for the first time how the strength of this signal can be substantially enhanced by using information about the halo environment at a range of scales. This enhancement is achieved by using certain combinations of the large scale Cold Dark Matter and total matter environments of halos, both of which are measurable from galaxy clustering and weak lensing surveys.

Several interesting Dark Matter (DM) models invoke a dark sector leading to two types of relic particles, possibly interacting with each other: non-relativistic DM, and relativistic Dark Radiation (DR). These models have interesting consequences for cosmological observables, and could in principle solve problems like the small-scale cold DM crisis, Hubble tension, and/or low $\sigma_8$ value. Their cosmological behaviour is captured by the ETHOS parametrisation, which includes a DR-DM scattering rate scaling like a power-law of the temperature, $T^n$. Scenarios with $n=0$, $2$, or $4$ can easily be realised in concrete dark sector set-ups. Here we update constraints on these three scenarios using recent CMB, BAO, and high-resolution Lyman-$\alpha$ data. We introduce a new Lyman-$\alpha$ likelihood that is applicable to a wide range of cosmological models with a suppression of the matter power spectrum on small scales. For $n=2$ and $4$, we find that Lyman-$\alpha$ data strengthen the CMB+BAO bounds on the DM-DR interaction rate by many orders of magnitude. However, models offering a possible solution to the missing satellite problem are still compatible with our new bounds. For $n=0$, high-resolution Lyman-$\alpha$ data bring no stronger constraints on the interaction rate than CMB+BAO data, except for extremely small values of the DR density. Using CMB+BAO data and a theory-motivated prior on the minimal density of DR, we find that the $n=0$ model can reduce the Hubble tension from $4.1\sigma$ to $2.7\sigma$, while simultaneously accommodating smaller values of the $\sigma_8$ and $S_8$ parameters hinted by cosmic shear data.

Massive neutrinos, due to their free streaming, produce a suppression in the matter power spectrum at intermediate and small scales which could be probed by galaxy clustering and/or weak lensing observables. This effect happens at scales that are also influenced by baryon feedback, i.e. galactic winds or Active Galactic Nuclei (AGN) feedback, which in realistic hydrodynamic simulations has also been shown to produce a suppression of power. Leaving aside, for the moment, the complex issue of galaxy bias, we focus here on matter clustering and tomographic weak lensing, we investigate the possible degeneracy between baryon feedback and neutrinos showing that it is not likely to degrade significantly the measurement of neutrino mass in future surveys. To do so, we generate mock data sets and fit them using the Markov Chain Monte Carlo (MCMC) technique and explore degeneracies between feedback parameters and neutrino mass. We model baryon feedback through fitting functions, while massive neutrinos are accounted for, also in the non-linear regime, using Halofit calibrated against accurate N-body neutrino simulations. In the error budget, we include the uncertainty in the modelling of non-linearities. For both matter clustering and weak lensing, we always recover the input neutrino mass within $\sim 0.25\sigma$ confidence level. Finally, we also take into account the intrinsic alignment effect in the weak lensing mock data. Even in this case, we are able to recover the right parameters: in particular, we find a significant degeneracy pattern between $M_\nu$ and the intrinsic alignment parameter $A_\mathrm{IA}$ .

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

In this paper we present an efficient method for constraining both thermal and non-thermal Dark Matter (DM) scenarios with the Lyman-$\alpha$ forest, based on a simple and flexible parametrisation capable to reproduce the small scale clustering signal of a large set of non-cold DM (nCDM) models. We extract new limits on the fundamental DM properties, through an extensive analysis of the high resolution, high redshift data obtained by the MIKE/HIRES spectrographs. By using a large suite of hydrodynamical simulations, we determine constraints on both astrophysical, cosmological, and nCDM parameters by performing a full Monte Carlo Markov Chain (MCMC) analysis. We obtain a marginalised upper limit on the largest possible scale at which a power suppression induced by nearly any nCDM scenario can occur, i.e. $\alpha<0.03~{\rm{Mpc}}/h$ (2$\sigma$ C.L.). We explicitly describe how to test several of the most viable nCDM scenarios without the need to run any specific numerical simulations, due to the novel parametrisation proposed, and due to a new scheme that interpolates between the cosmological models explored. The shape of the linear matter power spectrum for standard thermal warm DM models appear to be in mild tension ($\sim 2\sigma$ C.L.) with the data, compared to non-thermal scenarios. We show that a DM fluid composed by both a warm (thermal) and a cold component is also in tension with the Lyman-$\alpha$ forest, at least for large $\alpha$ values. This is the first study that allows to probe the linear small scale shape of the DM power spectrum for a large set of nCDM models.

Fuzzy Dark Matter (FDM) represents an alternative and intriguing description of the standard Cold Dark Matter (CDM) fluid, able to explain the lack of direct detection of dark matter particles in the GeV sector and to alleviate small scales tensions in the cosmic large-scale structure formation. Cosmological simulations of FDM models in the literature were performed either with very expensive high-resolution grid-based simulations of individual haloes or through N-body simulations encompassing larger cosmic volumes but resorting on significant approximations in the FDM non-linear dynamics to reduce their computational cost. With the use of the new N-body cosmological hydrodynamical code AX-GADGET, we are now able not only to overcome such numerical problems, but also to combine a fully consistent treatment of FDM dynamics with the presence of gas particles and baryonic physical processes, in order to quantify the FDM impact on specific astrophysical observables. In particular, in this paper we perform and analyse several hydrodynamical simulations in order to constrain the FDM mass by quantifying the impact of FDM on Lyman-$\alpha$ forest observations, as obtained for the first time in the literature in a N-body setup without approximating the FDM dynamics. We also study the statistical properties of haloes, exploiting the large available sample, to extract information on how FDM affects the abundance, the shape, and density profiles of dark matter haloes.

We derive a generalized luminosity distance versus redshift relation for a linearly perturbed FLRW (Friedmann-Lemaitre-Robertson-Walker) metric with two scalar mode excitations. We use two equivalent approaches, based on the Jacobi map and the van Vleck determinant respectively. We apply the resultant formula to two simple models - an exact FLRW universe and an approximate FLRW universe perturbed by a single scalar mode sinusoidally varying with time. For both models we derive a cosmographic expansion for d_L in terms of z. We comment on the interpretation of our results and their possible application to more realistic cosmological models.

We characterize the properties of the intergalactic medium (IGM) around a sample of galaxies extracted from state-of-the-art hydrodynamical simulations of structure formation in a cosmological volume of 25 Mpc comoving at $z\sim 2$. The simulations are based on two different subresolution schemes for star formation and supernova feedback: the MUlti-Phase Particle Integrator (MUPPI) scheme and the Effective Model. We develop a quantitative and probabilistic analysis based on the apparent optical depth method of the properties of the absorbers as a function of impact parameter from their nearby galaxies: in such a way we probe different environments from circumgalactic medium to low-density filaments. Absorbers' properties are then compared with a spectroscopic observational data set obtained from high-resolution quasar spectra. Our main focus is on the N$_{\rm{ CIV}}$-N$_{\rm { HI}}$ relation around simulated galaxies: the results obtained with MUPPI and the Effective model are remarkably similar, with small differences only confined to regions at impact parameters $b = [1-3] \times r_{\rm {vir}}$. Using $\mbox{C IV}$ as a tracer of the metallicity, we obtain evidence that the observed metal absorption systems have the highest probability to be confined in a region of 150-400 kpc around galaxies. Near-filament environments have instead metallicities too low to be probed by present-day telescopes, but could be probed by future spectroscopical studies. Finally we compute $\mbox{C IV}$ covering fractions which are in agreement with observational data.

The notion of luminosity distance is most often defined in purely FLRW (Friedmann-Lemaitre-Robertson-Walker) cosmological spacetimes, or small perturbations thereof. However, the abstract notion of luminosity distance is actually much more robust than this, and can be defined non-perturbatively in almost arbitrary spacetimes. Some quite general results are already known, in terms of $dA_\mathrm{observer}/d\Omega_\mathrm{source}$, the cross-sectional area per unit solid angle of a null geodesic spray emitted from some source and subsequently detected by some observer. We shall reformulate these results in terms of a suitably normalized null geodesic affine parameter and the van Vleck determinant, $\Delta_{vV}$. The contribution due to the null geodesic affine parameter is effectively the inverse square law for luminosity, and the van Vleck determinant can be viewed as providing a measure of deviations from the inverse square law. This formulation is closely related to the so-called Jacobi determinant, but the van Vleck determinant has somewhat nicer analytic properties and wider and deeper theoretical base in the general relativity, quantum physics, and quantum field theory communities. In the current article we shall concentrate on non-perturbative results, leaving near-FLRW perturbative investigation for future work.

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 investigate the future perspectives of the detection of the relativistic dipole by cross-correlating the 21 cm emission in Intensity Mapping (IM) and galaxy surveys at low redshift. We model the neutral hydrogen (HI) and the galaxy population by means of the halo model to relate the parameters that affect the dipole signal such as the biases of the two tracers and the Poissonian noise. We investigate the behavior of the signal-to-noise as a function of the galaxy and magnification biases, for two fixed models of the neutral hydrogen. In both cases we found that the signal-to-noise does not grow by increasing the difference between the biases of the two tracers, due to the larger shot-noise yields by highly biased tracers. We also study and provide an optimal luminosity-threshold galaxy catalogue to enhance the signal-to-noise ratio of the relativistic dipole. Interestingly, we show that the maximum magnitude provided by the survey does not lead to the maximum signal-to-noise for detecting relativistic effects and we predict the optimal value for the limiting magnitude. Our work suggests that an optimal analysis could increase the signal-to-noise ratio up to a factor five compared to a standard one.

We investigate the possibility of performing cosmological studies in the redshift range $2.5<z<5$ through suitable extensions of existing and upcoming radio-telescopes like CHIME, HIRAX and FAST. We use the Fisher matrix technique to forecast the bounds that those instruments can place on the growth rate, the BAO distance scale parameters, the sum of the neutrino masses and the number of relativistic degrees of freedom at decoupling, $N_{\rm eff}$. We point out that quantities that depend on the amplitude of the 21cm power spectrum, like $f\sigma_8$, are completely degenerate with $\Omega_{\rm HI}$ and $b_{\rm HI}$, and propose several strategies to independently constraint them through cross-correlations with other probes. Assuming $5\%$ priors on $\Omega_{\rm HI}$ and $b_{\rm HI}$, $k_{\rm max}=0.2~h{\rm Mpc}^{-1}$ and the primary beam wedge, we find that a HIRAX extension can constrain, within bins of $\Delta z=0.1$: 1) the value of $f\sigma_8$ at $\simeq4\%$, 2) the value of $D_A$ and $H$ at $\simeq1\%$. In combination with data from Euclid-like galaxy surveys and CMB S4, the sum of the neutrino masses can be constrained with an error equal to $23$ meV ($1\sigma$), while $N_{\rm eff}$ can be constrained within 0.02 ($1\sigma$). We derive similar constraints for the extensions of the other instruments. We study in detail the dependence of our results on the instrument, amplitude of the HI bias, the foreground wedge coverage, the nonlinear scale used in the analysis, uncertainties in the theoretical modeling and the priors on $b_{\rm HI}$ and $\Omega_{\rm HI}$. We conclude that 21cm intensity mapping surveys operating in this redshift range can provide extremely competitive constraints on key cosmological parameters.

We investigate constraints on scalar dark matter (DM) by analyzing the Lyman-alpha forest, which probes structure formation at medium and small scales, and also by studying its cosmological consequences at high and low redshift. For scalar DM that constitutes more than 30% of the total DM density, we obtain a lower limit m >~ 10^{-21} eV for the mass of scalar DM. This implies an upper limit on the initial field displacement (or the decay constant for an axion-like field) of phi <~ 10^{16} GeV. We also derive limits on the energy scale of cosmic inflation and establish an upper bound on the tensor-to-scalar ratio of r < 10^{-3} in the presence of scalar DM. Furthermore, we show that there is very little room for ultralight scalar DM to solve the "small-scale crisis" of cold DM without spoiling the Lyman-alpha forest results. The constraints presented in this paper can be used for testing generic theories that contain light scalar fields.

Structure formation at small cosmological scales provides an important frontier for dark matter (DM) research. Scenarios with small DM particle masses, large momenta or hidden interactions tend to suppress the gravitational clustering at small scales. The details of this suppression depend on the DM particle nature, allowing for a direct link between DM models and astrophysical observations. However, most of the astrophysical constraints obtained so far refer to a very specific shape of the power suppression, corresponding to thermal warm dark matter (WDM), i.e., candidates with a Fermi-Dirac or Bose-Einstein momentum distribution. In this work we introduce a new analytical fitting formula for the power spectrum, which is simple yet flexible enough to reproduce the clustering signal of large classes of non-thermal DM models, which are not at all adequately described by the oversimplified notion of WDM. We show that the formula is able to fully cover the parameter space of sterile neutrinos (whether resonantly produced or from particle decay), mixed cold and warm models, fuzzy dark matter, as well as other models suggested by effective theory of structure formation (ETHOS). Based on this fitting formula, we perform a large suite of N-body simulations and we extract important nonlinear statistics, such as the matter power spectrum and the halo mass function. Finally, we present first preliminary astrophysical constraints, based on linear theory, from both the number of Milky Way satellites and the Lyman-$\alpha$ forest. This paper is a first step towards a general and comprehensive modeling of small-scale departures from the standard cold DM model.

Gravitational non-linear evolution induces a shift in the position of the baryon acoustic oscillations (BAO) peak together with a damping and broadening of its shape that bias and degrades the accuracy with which the position of the peak can be determined. BAO reconstruction is a technique developed to undo part of the effect of non-linearities. We present and analyse a reconstruction method that consists of displacing pixels instead of galaxies and whose implementation is easier than the standard reconstruction method. We show that this method is equivalent to the standard reconstruction technique in the limit where the number of pixels becomes very large. This method is particularly useful in surveys where individual galaxies are not resolved, as in 21cm intensity mapping observations. We validate this method by reconstructing mock pixelated maps, that we build from the distribution of matter and halos in real- and redshift-space, from a large set of numerical simulations. We find that this method is able to decrease the uncertainty in the BAO peak position by 30-50% over the typical angular resolution scales of 21 cm intensity mapping experiments.

Large-scale extragalactic magnetic fields may induce conversions between very-high-energy photons and axionlike particles (ALPs), thereby shielding the photons from absorption on the extragalactic background light. However, in simplified "cell" models, used so far to represent extragalactic magnetic fields, this mechanism would be strongly suppressed by current astrophysical bounds. Here we consider a recent model of extragalactic magnetic fields obtained from large-scale cosmological simulations. Such simulated magnetic fields would have large enhancement in the filaments of matter. As a result, photon-ALP conversions would produce a significant spectral hardening for cosmic TeV photons. This effect would be probed with the upcoming Cherenkov Telescope Array detector. This possible detection would give a unique chance to perform a tomography of the magnetized cosmic web with ALPs.

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 imprint of non-standard dark energy (DE) and dark matter (DM) models on the 21cm intensity map power spectra from high-redshift neutral hydrogen (HI) gas. To this purpose we use halo catalogs from N-body simulations of dynamical DE models and DM scenarios which are statistically indistinguishable from the standard Cold Dark Matter model with Cosmological Constant (LCDM) using currently available cosmological observations. We limit our analysis to halo catalogs at redshift z = 1 and 2.3 which are common to all simulations. For each catalog we model the HI distribution by using a simple prescription to associate the HI gas mass to N-body halos. We find that the DE models leave a distinct signature on the HI spectra across a wide range of scales, which correlates with differences in the halo mass function and the onset of the non-linear regime of clustering. In the case of the non-standard DM model significant differences of the HI spectra with respect to the LCDM model only arise from the suppressed abundance of low mass halos. These cosmological model dependent features also appear in the 21cm spectra. In particular, we find that future SKA measurements can distinguish the imprints of DE and DM models at high statistical significance.

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 investigate the cross-correlation signal between 21cm intensity mapping maps and the Lyman-alpha forest in the fully non-linear regime using state-of-the-art hydrodynamic simulations. The cross-correlation signal between these fields can provide a coherent and comprehensive picture of the neutral hydrogen (HI) content of our Universe in the post-reionization era, probing both its mass content and volume distribution. We compute the auto-power spectra of both fields together with their cross-power spectrum at z = 2.4 and find that on large scales the fields are completely anti-correlated. This anti-correlation arises because regions with high (low) 21cm emission, such as those with a large (low) concentration of damped Lyman-alpha systems, will show up as regions with low (high) transmitted flux. We find that on scales smaller than k = 0.2 h/Mpc the cross-correlation coefficient departs from -1, at a scale where non-linearities show up. We use the anisotropy of the power spectra in redshift-space to determine the values of the bias and of the redshift-space distortion parameters of both fields; we find that the errors on the value of these parameters could decrease by 30% when adding data from the cross-power spectrum in a conservative analysis. Our results point out that linear theory is capable of reproducing the shape and amplitude of the cross-power up to rather non-linear scales. Finally, we find that the 21cm-Lya cross-power spectrum can be detected by combining data from a BOSS-like survey together with 21cm intensity mapping observations by SKA1-MID with a S/N ratio higher than 3 in the range 0.06< k <1 h/Mpc. We emphasize that while the shape and amplitude of the 21cm auto-power spectrum can be severely affected by residual foreground contamination, cross-power spectra will be less sensitive to that and therefore can be used to identify systematics in the 21cm maps.

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

We investigate the Alcock Paczynski (AP) test applied to the Baryon Acoustic Oscillation (BAO) feature in the galaxy correlation function. By using a general formalism that includes relativistic effects, we quantify the importance of the linear redshift space distortions and gravitational lensing corrections to the galaxy number density fluctuation. We show that redshift space distortions significantly affect the shape of the correlation function, both in radial and transverse directions, causing different values of galaxy bias to induce offsets up to 1% in the AP test. On the other hand, we find that the lensing correction around the BAO scale modifies the amplitude but not the shape of the correlation function and therefore does not introduce any systematic effect. Furthermore, we investigate in details how the AP test is sensitive to redshift binning: a window function in transverse direction suppresses correlations and shifts the peak position toward smaller angular scales. We determine the correction that should be applied in order to account for this effect, when performing the test with data from three future planned galaxy redshift surveys: Euclid, the Dark Energy Spectroscopic Instrument (DESI) and the Square Kilometer Array (SKA).

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

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.

We report the detection of a cross-correlation signal between {\it Fermi} Large Area Telescope diffuse gamma-ray maps and catalogs of clusters. In our analysis, we considered three different catalogs: WHL12, redMaPPer and PlanckSZ. They all show a positive correlation with different amplitudes, related to the average mass of the objects in each catalog, which also sets the catalog bias. The signal detection is confirmed by the results of a stacking analysis. The cross-correlation signal extends to rather large angular scales, around 1 degree, that correspond, at the typical redshift of the clusters in these catalogs, to a few to tens of Mpc, i.e. the typical scale-length of the large scale structures in the Universe. Most likely this signal is contributed by the cumulative emission from AGNs associated to the filamentary structures that converge toward the high peaks of the matter density field in which galaxy clusters reside. In addition, our analysis reveals the presence of a second component, more compact in size and compatible with a point-like emission from within individual clusters. At present, we cannot distinguish between the two most likely interpretations for such a signal, i.e. whether it is produced by AGNs inside clusters or if it is a diffuse gamma-ray emission from the intra-cluster medium. We argue that this latter, intriguing, hypothesis might be tested by applying this technique to a low redshift large mass cluster sample.

In these proceedings we describe the WEAVE-QSO survey, which will observe around 400,000 high redshift quasars starting in 2018. This survey is part of a broader WEAVE survey to be conducted at the 4.2m William Herschel Telescope. We will focus on chiefly on the science goals, but will also briefly summarise the target selection methods anticipated and the expected survey plan. Understanding the apparent acceleration in the expansion of the Universe is one of the key scientific challenges of our time. Many experiments have been proposed to study this expansion, using a variety of techniques. Here we describe a survey that can measure this acceleration and therefore help elucidate the nature of dark energy: a survey of the Lyman-alpha forest (and quasar absorption in general) in spectra towards z>2 quasars (QSOs). Further constraints on neutrino masses and warm dark matter are also anticipated. The same data will also shed light on galaxy formation via study of the properties of inflowing/outflowing gas associated with nearby galaxies and in a cosmic web context. Gas properties are sensitive to density, temperature, UV radiation, metallicity and abundance pattern, and so constraint galaxy formation in a variety of ways. WEAVE-QSO will study absorbers with a dynamic range spanning more than 8 orders of magnitude in column density, their thermal broadening, and a host of elements and ionization species. A core principal of the WEAVE-QSO survey is the targeting of QSOs with near 100% efficiency principally through use of the J-PAS (r < 23.2) and Gaia (r < 20) data.

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $\Delta\Sigma$ and show that several of these effects may be non-negligible given the precision of our measurement. Disentangling cosmological effects from the details of the galaxy-halo connection, the effects of baryons, and massive neutrinos, is the next challenge facing joint lensing and clustering analyses. This is especially true in the context of large galaxy samples from Baryon Acoustic Oscillation surveys with precise measurements but complex selection functions.

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 compare the measured angular cross-correlation between the Fermi-LAT gamma-ray sky and catalogues of extra-galactic objects with the expected signal induced by weakly interacting massive particle (WIMP) dark matter (DM). We include a detailed description of the contribution of astrophysical gamma-ray emitters such as blazars, misaligned AGN and star forming galaxies, and perform a global fit to the measured cross-correlation. Five catalogues are considered: SDSS-DR6 quasars, 2MASS galaxies, NVSS radio galaxies, SDSS-DR8 Luminous Red Galaxies and SDSS-DR8 main galaxy sample. To model the cross-correlation signal we use the halo occupation distribution formalism to estimate the number of galaxies of a given catalogue in DM halos and their spatial correlation properties. We discuss uncertainties in the predicted cross-correlation signal arising from the DM clustering and WIMP microscopic properties, which set the DM gamma-ray emission. The use of different catalogues probing objects at different redshifts reduces significantly, though not completely, the degeneracy among the different gamma-ray components. We find that the presence of a significant WIMP DM signal is allowed by the data but not significantly preferred by the fit, although this is mainly due to a degeneracy with the misaligned AGN component. With modest substructure boost, the sensitivity of this method excludes thermal annihilation cross sections at 95% C.L. for WIMP masses up to few tens of GeV. Constraining the low-redshift properties of astrophysical populations with future data will further improve the sensitivity to DM.

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 quantitatively investigate the possibility of detecting baryonic acoustic oscillations (BAO) using single-dish 21cm intensity mapping observations in the post-reionization era. We show that the telescope beam smears out the isotropic BAO signature and, in the case of the Square Kilometer Array (SKA) instrument, makes it undetectable at redshifts $z\gtrsim1$. We however demonstrate that the BAO peak can still be detected in the radial 21cm power spectrum and describe a method to make this type of measurements. By means of numerical simulations, containing the 21cm cosmological signal as well as the most relevant Galactic and extra-Galactic foregrounds and basic instrumental effect, we quantify the precision with which the radial BAO scale can be measured in the 21cm power spectrum. We systematically investigate the signal-to-noise and the precision of the recovered BAO signal as a function of cosmic variance, instrumental noise, angular resolution and foreground contamination. We find that the expected noise levels of SKA would degrade the final BAO errors by $\sim5\%$ with respect to the cosmic-variance limited case at low redshifts, but that the effect grows up to $\sim65\%$ at $z\sim2-3$. Furthermore, we find that the radial BAO signature is robust against foreground systematics, and that the main effect is an increase of $\sim20\%$ in the final uncertainty on the standard ruler caused by the contribution of foreground residuals as well as the reduction in sky area needed to avoid high-foreground regions. We also find that it should be possible to detect the radial BAO signature with high significance in the full redshift range. We conclude that a 21cm experiment carried out by the SKA should be able to make direct measurements of the expansion rate $H(z)$ with competitive per-cent level precision on redshifts $z\lesssim2.5$.

We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.

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.

Clusters of galaxies have the potential of providing powerful constraints on possible deviations from General Relativity. We use the catalogue of Sunyaev-Zel'dovich sources detected by Planck and consider a correction to the halo mass function for a f(R) class of modified gravity models, which has been recently found to reproduce well results from N-body simulations, to place constraints on the scalaron field amplitude at the present time, $f_{R}^0$. We find that applying this correction to different calibrations of the halo mass function produces upper bounds on $f_{R}^0$ tighter by more than an order of magnitude, ranging from $\log_{10}(-f_{R}^0) < -5.81$ to $\log_{10}(-f_{R}^0) < -4.40$ (95 % confidence level). This sensitivity is due to the different shape of the halo mass function, which is degenerate with the parameters used to calibrate the scaling relations between SZ observables and cluster masses. Any claim of constraints more stringent that the weaker limit above, based on cluster number counts, appear to be premature and must be supported by a careful calibration of the halo mass function and by a robust calibration of the mass scaling relations.

The free-streaming of keV-scale particles impacts structure growth on scales that are probed by the Lyman-alpha forest of distant quasars. Using an unprecedentedly large sample of medium-resolution QSO spectra from the ninth data release of SDSS, along with a state-of-the-art set of hydrodynamical simulations to model the Lyman-alpha forest in the non-linear regime, we issue one of the tightest bounds to date, from Ly-$\alpha$ data alone, on pure dark matter particles : $m_X > 4.09 \: \rm{keV}$ (95% CL) for early decoupled thermal relics such as a hypothetical gravitino, and correspondingly $m_s > 24.4 \: \rm{keV}$ (95% CL) for a non-resonantly produced right-handed neutrino. This limit depends on the value on $n_s$, and Planck measures a higher value of $n_s$ than SDSS-III/BOSS. Our bounds thus change slightly when Ly-$\alpha$ data are combined with CMB data from Planck 2016. The limits shift to $m_X > 2.96 \: \rm{keV}$ (95% CL) and $m_s > 16.0 \: \rm{keV}$ (95% CL). Thanks to SDSS-III data featuring smaller uncertainties and covering a larger redshift range than SDSS-I data, our bounds confirm the most stringent results established by previous works and are further at odds with a purely non-resonantly produced sterile neutrino as dark matter.

We use the first release of the VImos Public Extragalactic Redshift Survey of galaxies (VIPERS) of ~50,000 objects to measure the biasing relation between galaxies and mass in the redshift range z=[0.5,1.1]. We estimate the 1-point distribution function [PDF] of VIPERS galaxies from counts in cells and, assuming a model for the mass PDF, we infer their mean bias relation. The reconstruction of the bias relation from PDFs is performed through a novel method that accounts for Poisson noise, redshift distortions, inhomogeneous sky coverage and other selection effects. With this procedure we constrain galaxy bias and its deviations from linearity down to scales as small as 4 Mpc/h and out to z=1.1. We detect small (~3%) but significant deviations from linear bias. The mean biasing function is close to linear in regions above the mean density. The mean slope of the biasing relation is a proxy to the linear bias parameter. It increases both with luminosity, in agreement with results of previous analyses, and with redshift. However, we detect a strong bias evolution only for z>0.9 in agreement with some, but not all, previous studies. We also detected a significant increase of the bias with the scale, from 4 to 8 Mpc/h, now seen for the first time out to z=1. The amplitude of nonlinearity depends on redshift, luminosity and on scales but no clear trend is detected. Thanks to the large cosmic volume probed by VIPERS we find that the mismatch between the previous estimates of bias at z~1 from zCOSMOS and VVDS-Deep galaxy samples is fully accounted for by cosmic variance. The results of our work confirm the importance of going beyond the over-simplistic linear bias hypothesis showing that non-linearities can be accurately measured through the applications of the appropriate statistical tools to existing datasets like VIPERS.

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

The discovery of neutrino masses through the observation of oscillations boosted the importance of neutrinoless double beta decay ($0\nu\beta\beta$). In this paper, we review the main features of this process, underlining its key role both from the experimental and theoretical point of view. In particular, we contextualize the $0\nu\beta\beta$ in the panorama of lepton-number violating processes, also assessing some possible particle physics mechanisms mediating the process. Since the $0\nu\beta\beta$ existence is correlated with neutrino masses, we also review the state-of-art of the theoretical understanding of neutrino masses. In the final part, the status of current $0\nu\beta\beta$ experiments is presented and the prospects for the future hunt for $0\nu\beta\beta$ are discussed. Also, experimental data coming from cosmological surveys are considered and their impact on $0\nu\beta\beta$ expectations is examined.

We present the calculation of the Lyman-alpha (Lyman-$\alpha$) transmitted flux fluctuations with full relativistic corrections to the first order. Even though several studies exist on relativistic effects in galaxy clustering, this is the first study to extend the formalism to a different tracer of underlying matter at unique redshift range ($z=2-5$). Furthermore, we show a comprehensive application of our calculations of the Quasar-Lyman-$\alpha$ cross-correlation function. Our results indicate that the signal of relativistic effects is sizeable at Baryonic Acoustic Oscillation (BAO) scale mainly due to the large difference in density bias factors of our tracers. We construct an observable, the anti-symmetric part of the cross-correlation function, that is dominated by the relativistic signal and offers a new way to measure the relativistic terms at relatively small scales. The analysis shows that relativistic effects are important when considering cross-correlations between tracers with very different biases, and should be included in the data analysis of the current and future surveys. Moreover, the idea presented in this paper is highly complementary to other techniques and observable trying to isolate the effect of the relativistic corrections and thus test the validity of the theory of gravity beyond the Newtonian regime.

Cosmology is making impressive progress and it is producing stringent bounds on the sum of the neutrino masses {\Sigma}, a parameter of great importance for the current laboratory experiments. In this letter, we exploit the potential relevance of the analysis of Palanque-Delabrouille et al. [JCAP 1502, 045 (2015)] to the neutrinoless double beta decay (0{\nu}{\beta}{\beta}) search. This analysis indicates small values for the lightest neutrino mass, since the authors find {\Sigma} < 84 meV at 1{\sigma} C. L., and provides a 1{\sigma} preference for the normal hierarchy. The allowed values for the Majorana effective mass, probed by 0{\nu}{\beta}{\beta}, turn out to be < 75meV at 3{\sigma}C.L. and lower down to less than 20meV at 1{\sigma}C.L.. If this indication is confirmed, the impact on the 0{\nu}{\beta}{\beta} experiments will be tremendous since the possibility of detecting a signal will be out of the reach of the next generation of experiments.

By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $\Omega_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $\Omega_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

We investigate the signatures left by massive neutrinos on the spatial distribution of neutral hydrogen (HI) in the post-reionization era by running hydrodynamic simulations that include massive neutrinos as additional collisionless particles. We find that halos in massive/massless neutrino cosmologies host a similar amount of neutral hydrogen, although for a fixed halo mass, on average, the HI mass increases with the sum of the neutrino masses. Our results show that HI is more strongly clustered in cosmologies with massive neutrinos, while its abundance, $\Omega_{\rm HI}(z)$, is lower. These effects arise mainly from the impact of massive neutrinos on cosmology: they suppress both the amplitude of the matter power spectrum on small scales and the abundance of dark matter halos. Modelling the HI distribution with hydrodynamic simulations at $z > 3$, and a simple analytic model at $z<3$, we use the Fisher matrix formalism to conservatively forecast the constraints that Phase 1 of the Square Kilometre Array (SKA) will place on the sum of neutrino masses, $M_\nu\equiv \Sigma m_{\nu}$. We find that with 10,000 hours of interferometric observations at $3 \lesssim z \lesssim 6$ from a deep and narrow survey with SKA1-LOW, the sum of the neutrino masses can be measured with an error $\sigma(M_\nu)\lesssim0.3$ eV (95% CL). Similar constraints can be obtained with a wide and deep SKA1-MID survey at $z \lesssim 3$, using the single-dish mode. By combining data from MID, LOW, and Planck, plus priors on cosmological parameters from a Stage IV spectroscopic galaxy survey, the sum of the neutrino masses can be determined with an error $\sigma(M_\nu)\simeq0.06$ eV (95% CL).

Cosmic voids are a promising environment to characterize neutrino-induced effects on the large-scale distribution of matter in the universe. We perform a comprehensive numerical study of the statistical properties of voids, identified both in the matter and galaxy distributions, in massive and massless neutrino cosmologies. The matter density field is obtained by running several independent $N$-body simulations with cold dark matter and neutrino particles, while the galaxy catalogs are modeled by populating the dark matter halos in simulations via a halo occupation distribution (HOD) model to reproduce the clustering properties observed by the Sloan Digital Sky Survey (SDSS) II Data Release 7. We focus on the impact of massive neutrinos on the following void statistical properties: number density, ellipticities, two-point statistics, density and velocity profiles. Considering the matter density field, we find that voids in massive neutrino cosmologies are less evolved than those in the corresponding massless neutrinos case: there is a larger number of small voids and a smaller number of large ones, their profiles are less evacuated, and they present a lower wall at the edge. Moreover, the degeneracy between $\sigma_8$ and $\Omega_{\nu}$ is broken when looking at void properties. In terms of the galaxy density field, we find that differences among cosmologies are difficult to detect because of the small number of galaxy voids in the simulations. Differences are instead present when looking at the matter density and velocity profiles around these voids.

The Lyman alpha forest power spectrum has been measured on large scales by the BOSS survey in SDSS-III at $z\sim 2.3$, has been shown to agree well with linear theory predictions, and has provided the first measurement of Baryon Acoustic Oscillations at this redshift. However, the power at small scales, affected by non-linearities, has not been well examined so far. We present results from a variety of hydrodynamic simulations to predict the redshift space non-linear power spectrum of the Lyman Alpha transmission for several models, testing the dependence on resolution and box size. A new fitting formula is introduced to facilitate the comparison of our simulation results with observations and other simulations. The non-linear power spectrum has a generic shape determined by a transition scale from linear to non-linear anisotropy, and a Jeans scale below which the power drops rapidly. In addition, we predict the two linear bias factors of the Lyman Alpha forest and provide a better physical interpretation of their values and redshift evolution. The dependence of these bias factors and the non-linear power on the amplitude and slope of the primordial fluctuations power spectrum, the temperature-density relation of the intergalactic medium, and the mean Lyman Alpha transmission, as well as the redshift evolution, is investigated and discussed in detail. A preliminary comparison to the observations shows that the predicted redshift distortion parameter is in good agreement with the recent determination of Blomqvist et al., but the density bias factor is lower than observed. We make all our results publicly available in the form of tables of the non-linear power spectrum that is directly obtained from all our simulations, and parameters of our fitting formula.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

We study the impact of neutrino masses on the shape and height of the BAO peak of the matter correlation function, both in real and redshift space. In order to describe the nonlinear evolution of the BAO peak we run N-body simulations and compare them with simple analytic formulae. We show that the evolution with redshift of the correlation function and its dependence on the neutrino masses is well reproduced in a simplified version of the Zel'dovich approximation, in which the mode-coupling contribution to the power spectrum is neglected. While in linear theory the BAO peak decreases for increasing neutrino masses, the effect of nonlinear structure formation goes in the opposite direction, since the peak broadening by large scale flows is less effective. As a result of this combined effect, the peak decreases by $\sim 0.6 \%$ for $ \sum m_\nu = 0.15$ eV and increases by $\sim1.2 \%$ for $ \sum m_\nu = 0.3$ eV, with respect to a massless neutrino cosmology with equal value of the other cosmological parameters. We extend our analysis to redshift space and to halos, and confirm the agreement between simulations and the analytic formulae. We argue that all analytical approaches having the Zel'dovich propagator in their lowest order approximation should give comparable performances, irrespectively to their formulation in Lagrangian or in Eulerian space.

If dark matter (DM) is composed by particles which are non-gravitationally coupled to ordinary matter, their annihilations or decays in cosmic structures can result in detectable radiation. We show that the most powerful technique to detect a particle DM signal outside the Local Group is to study the angular cross-correlation of non-gravitational signals with low-redshift gravitational probes. This method allows to enhance signal-to-noise from the regions of the Universe where the DM-induced emission is preferentially generated. We demonstrate the power of this approach by focusing on GeV-TeV DM and on the recent cross-correlation analysis between the 2MASS galaxy catalogue and the Fermi-LAT gamma-ray maps. We show that this technique is more sensitive than other extragalactic gamma-ray probes, such as the energy spectrum and angular autocorrelation of the extragalactic background, and emission from clusters of galaxies. Intriguingly, we find that the measured cross-correlation can be well fitted by a DM component, with thermal annihilation cross section and mass between 10 and 100 GeV, depending on the small-scale DM properties and gamma-ray production mechanism. This solicits further data collection and dedicated analyses.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.

Building on our previous cross-correlation analysis (Xia et al. 2011) between the isotropic gamma-ray background (IGRB) and different tracers of the large-scale structure of the universe, we update our results using 60-months of data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. We perform a cross-correlation analysis between the IGRB and objects that may trace the astrophysical sources of the IGRB: SDSS-DR6 QSOs, the SDSS-DR8 Main Galaxy Sample, Luminous Red Galaxies (LRGs) in the SDSS catalog, 2MASS galaxies, and radio NVSS galaxies. The benefit of correlating the Fermi-LAT signal with catalogs of objects at various redshifts is to provide tomographic information on the IGRB which is crucial to separate the various contributions and to clarify its origin. We observe a significant (>3.5 sigma) cross-correlation signal on angular scales smaller than 1 deg in the NVSS, 2MASS and QSO cases and, at lower statistical significance (~3.0 sigma), with SDSS galaxies. These results are robust against the choice of the statistical estimator, estimate of errors, map cleaning procedure and instrumental effects. Finally, we test the hypothesis that the IGRB observed by Fermi-LAT originates from the summed contributions of three types of unresolved extragalactic sources: BL Lacs, FSRQs and Star-Forming Galaxies (SFGs). We find that a model in which the IGRB is mainly produced by SFGs ($72^{+23}_{-37}$% with 2 sigma errors), with BL Lacs and FSRQs giving a minor contribution, provides a good fit to the data. We also consider a possible contribution from Misaligned Active Galactic Nuclei, and we find that, depending on the details of the model and its uncertainty, they can also provide a substantial contribution, partly degenerate with the SFG one. (abridged)

We investigate the impact that warm dark matter (WDM) has in terms of 21cm intensity mapping in the post-reionization Universe at z = 3 - 5. We perform hydrodynamic simulations for 5 different models: cold dark matter and WDM with 1,2,3,4 keV (thermal relic) mass and assign the neutral hydrogen a-posteriori using two different methods that both reproduce observations in terms of column density distribution function of neutral hydrogen systems. Contrary to naive expectations, the suppression of power present in the linear and non-linear matter power spectra, results in an increase of power in terms of neutral hydrogen and 21cm power spectra. This is due to the fact that there is a lack of small mass halos in WDM models with respect to cold dark matter: in order to distribute a total amount of neutral hydrogen within the two cosmological models, a larger quantity has to be placed in the most massive halos, that are more biased compared to the cold dark matter cosmology. We quantify this effect and address significance for the telescope SKA1-LOW, including a realistic noise modeling. The results indicate that we will be able to rule out a 4 keV WDM model with 5000 hours of observations at z > 3, with a statistical significance of > 3 sigma, while a smaller mass of 3 keV, comparable to present day constraints, can be ruled out at more than 2 sigma confidence level with 1000 hours of observations at z > 5.

We provide a quantitative analysis of the halo model in the context of massive neutrino cosmologies. We discuss all the ingredients necessary to model the non-linear matter and cold dark matter power spectra and compare with the results of N-body simulations that incorporate massive neutrinos. Our neutrino halo model is able to capture the non-linear behavior of matter clustering with a $\sim 20\%$ accuracy up to very non-linear scales of $k=10~h/$Mpc (which would be affected by baryon physics). The largest discrepancies arise in the range $k=0.5-1~h/$Mpc where the 1-halo and 2-halo terms are comparable and are present also in a massless neutrino cosmology. However, at scales $k<0.2~h/$Mpc our neutrino halo model agrees with the results of N-body simulations at the level of 8\% for total neutrino masses of $<0.3$ eV. We also model the neutrino non-linear density field as a sum of a linear and clustered component and predict the neutrino power spectrum and the cold dark matter-neutrino cross-power spectrum up to $k=1~h/$Mpc with $\sim$ 30\% accuracy. For masses below 0.15 eV the neutrino halo model captures the neutrino induced suppression, casted in terms of matter power ratios between massive and massless scenarios, with a 2\% agreement with the results of N-body/neutrino simulations. Finally, we provide a simple application of the halo model: the computation of the clustering of galaxies, in massless and massive neutrinos cosmologies, using a simple Halo Occupation Distribution scheme and our halo model extension.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

We present constraints on the parameters of the $\Lambda$CDM cosmological model in the presence of massive neutrinos, using the one-dimensional Ly$\alpha$ forest power spectrum obtained with the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS) by Palanque-Delabrouille et al. (2013), complemented by additional cosmological probes. The interpretation of the measured Ly$\alpha$ spectrum is done using a second-order Taylor expansion of the simulated power spectrum. BOSS Ly$\alpha$ data alone provide better bounds than previous Ly$\alpha$ results, but are still poorly constraining, especially for the sum of neutrino masses $\sum m_\nu$, for which we obtain an upper bound of 1.1~eV (95\% CL), including systematics for both data and simulations. Ly$\alpha$ constraints on $\Lambda$CDM parameters and neutrino masses are compatible with CMB bounds from the Planck collaboration. Interestingly, the combination of Ly$\alpha$ with CMB data reduces the uncertainties significantly, due to very different directions of degeneracy in parameter space, leading to the strongest cosmological bound to date on the total neutrino mass, $\sum m_\nu < 0.15$~eV at 95\% CL (with a best-fit in zero). Adding recent BAO results further tightens this constraint to $\sum m_\nu < 0.14$~eV at 95\% CL. This bound is nearly independent of the statistical approach used, and of the different combinations of CMB and BAO data sets considered in this paper in addition to Ly$\alpha$. Given the measured values of the two squared mass differences $\Delta m^2$, this result tends to favor the normal hierarchy scenario against the inverted hierarchy scenario for the masses of the active neutrino species.

The constraints on neutrino masses led to the revision of their cosmological role, since the existence of a cosmological neutrino background is a clear prediction of the standard cosmological model. In this paper, we study the impact of such background on the spatial distribution of both Dark Matter (DM) and galaxies, by coupling $N$-body numerical simulations with semi-analytic models (SAMs) of galaxy formation. Cosmological simulations including massive neutrinos predict a slower evolution of DM perturbations with respect to the $\Lambda$CDM runs with the same initial conditions and a suppression on the matter power spectrum on small and intermediate scales, thus impacting on the predicted properties of galaxy populations. We explicitly show that most of these deviations are driven by the different $\sigma_8$ predicted for cosmologies including a massive neutrino background. We conclude that independent estimates of $\sigma_8$ are needed, in order to unambiguously characterise the effect of this background on the growth of structures. Galaxy properties alone are a weak tracer of deviations with respect to the $\Lambda$CDM run, but their combination with the overall matter distribution at all scales allows to disentangle between different cosmological models. Moreover, these deviations go on opposite direction with respect to competing models like modified gravity, thus weakening any detectable cosmological signal. Given the ubiquitous presence of a neutrino background, these effects have to be taken into account in future missions aimed at constraining the properties of the "Dark" components of the Universe.

We present a comprehensive analysis of high resolution hydrodynamic simulations in terms of Lyman-alpha and Lyman-beta one dimensional flux power spectra ($P_{\alpha\alpha}$ and $P_{\beta\beta}$). In particular, we focus on the behaviour that the flux auto-power spectra and cross-power spectra ($P_{\alpha\beta}$) display when the intergalactic medium (IGM) thermal history is changed in a range of values that bracket a reference model, while cosmological parameters are kept fixed to best fit the cosmic microwave background data. We present empirical fits that describe at the sub-percent level the dependence of the power spectra on the thermal parameters. At the largest scales, the power spectra show a constant bias between each other that is set by the parameters describing the IGM thermal state. The cross-power spectrum has an oscillatory pattern and crosses zero at a scale which depends on $T_0$, the IGM temperature at the mean density, for reasonable values of the power-law index $\gamma$ of the IGM temperature-density relation ($T=T_0(1+\delta)^{\gamma-1}$). By performing a Fisher matrix analysis, we find that the power spectrum $P_{\beta\beta}$ is more sensitive to the thermal history than $P_{\alpha\alpha}$ alone, due to the fact that it probes denser regions than Lyman-alpha. When we combine the power and cross spectra the constraints on $\gamma$ can be improved by a factor $\sim 4$, while the constraints on $T_0$ improve by a factor of $\sim 2$. We address the role of signal-to-noise and resolution by mocking realistic observations and we conclude that the framework presented in this work can significantly improve the knowledge of the IGM thermal state, which will in turn guarantee better constraints on IGM-derived cosmological parameters.

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

Standard cosmological perturbation theory (SPT) for the Large Scale Structure (LSS) of the Universe fails at small scales (UV) due to strong nonlinearities and to multistreaming effects. In Pietroni et al. 2011 a new framework was proposed in which the large scales (IR) are treated perturbatively while the information on the UV, mainly small scale velocity dispersion, is obtained by nonlinear methods like N-body simulations. Here we develop this approach, showing that it is possible to reproduce the fully nonlinear power spectrum (PS) by combining a simple (and fast) 1-loop computation for the IR scales and the measurement of a single, dominant, correlator from N-body simulations for the UV ones. We measure this correlator for a suite of seven different cosmologies, and we show that its inclusion in our perturbation scheme reproduces the fully non-linear PS with percent level accuracy, for wave numbers up to $k\sim 0.4\, h~{\rm Mpc^{-1}}$ down to $z=0$. We then show that, once this correlator has been measured in a given cosmology, there is no need to run a new simulation for a different cosmology in the suite. Indeed, by rescaling this correlator by a proper function computable in SPT, the reconstruction procedure works also for the other cosmologies and for all redshifts, with comparable accuracy. Finally, we clarify the relation of this approach to the Effective Field Theory methods recently proposed in the LSS context.

We measure and quantify properties of galactic outflows and diffuse gas at $z \geq 1$ in cosmological hydrodynamical simulations. Our novel sub-resolution model, MUPPI, implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at $z = 2$ by analysing central galaxies having a stellar mass higher than $10^{9} M_{\odot}$. The outflow velocity and mass outflow rate ($\dot{M}_{\rm out}$) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor ($\eta$) is between $0.2 - 10$. The comparison Effective model generates a constant outflow velocity, and a negative correlation of $\eta$ with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than $80 \%$ over $z = 1 - 5$. High SF activity at $z \sim 2 - 4$ drives strong outflows, causing the positive and steep correlations of velocity and $\dot{M}_{\rm out}$ with SFR. The outflow velocity correlation with SFR becomes flatter at $z = 1$, and $\eta$ displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and Effective models produce significant outflows at $\sim 1 / 10$ of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the Effective model. Our MUPPI model, using local properties of gas in the sub-resolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties, and consistent with observations.

(Abridged) We discuss the reliability of neutrino mass constraints, either active or sterile, from the combination of WMAP 9-year or Planck CMB data with BAO measurements from BOSS DR11, galaxy shear measurements from CFHTLenS, SDSS Ly-$\alpha$ forest constraints and galaxy cluster mass function from Chandra observations. To avoid model dependence of the constraints we perform a full likelihood analysis for all the datasets employed. As for the cluster data analysis we rely on to the most recent calibration of massive neutrino effects in the halo mass function and we explore the impact of the uncertainty in the mass bias and re-calibration of the halo mass function due to baryonic feedback processes on cosmological parameters. We find that none of the low redshift probes alone provide evidence for massive neutrinos in combination with CMB measurements, while a larger than $2\sigma$ detection of non zero neutrino mass, either active or sterile, is achieved combining cluster or shear data with CMB and BAO measurements. The preference for massive neutrino is larger in the sterile neutrino scenario, and for the combination of Planck, BAO, shear and cluster datasets we find that the vanilla $\Lambda$CDM model is rejected at more than $3\sigma$ and a sterile neutrino mass as motivated by accelerator anomaly is within the $2\sigma$ errors. Finally, results from the full data combination reflect the tension between the $\sigma_8$ constraints obtained from cluster and shear data and that inferred from Ly-$\alpha$ forest measurements; in the active neutrino scenario for both CMB datasets employed, the full data combination yields only an upper limits on $\sum m_\nu$, while assuming an extra sterile neutrino we still get preference for non-vanishing mass, $m_s^{\rm eff}=0.26^{+0.22}_{-0.24}$ eV, and dark contribution to the radiation content, $\Delta N_{\rm eff}=0.82\pm0.55$.

We present results of cosmological N-body hydrodynamic chemistry simulations of primordial structure growth and evolution in a scenario with warm dark matter (WDM) having a mass of 3 keV (thermal relic) and compare with a model consisting of standard cold dark matter (CDM). We focus on the high-redshift universe ($z>6$), where the structure formation process should better reflect the primordial (linear) differences in terms of matter power spectrum. We find that early epochs can be exceptional probes of the dark-matter nature. Non-linear WDM power spectra and mass functions are up to 2 dex lower than in CDM and show spreads of factor of a few persisting in the whole first Gyr. Runaway molecular cooling in WDM haloes results severely inhibited because of the damping of power at large $k$ modes and hence cosmic (population III and II-I) star formation rate (SFR) is usually suppressed with respect to CDM predictions. Luminous objects formed in a WDM background are very rare at $z>10$, due to the sparser and retarded evolution of early WDM mini-haloes during the dark ages and their lack can be fitted with a simple analytical formula depending only on magnitude and redshift. Future high-$z$ observations of faint galaxies have the potential to discriminate between CDM and WDM scenarios by means of cosmic stellar mass density (SMD) and specific SFR, as well. When compared to the effects of alternative cosmologies (e.g. non-Gaussian or dark-energy models) or of high-order corrections at large $z$ (e.g. primordial streaming motions or changes in the pristine IMF) the ones caused by WDM are definitely more dramatic.

We model the distribution of neutral hydrogen (HI) in the post-reionization era and investigate its detectability in 21 cm intensity mapping with the future SKA radio telescope. We rely on high resolution hydrodynamical N-body simulations. The HI is assigned a-posteriori to the gas particles following two different approaches: a halo-based method in which HI is assigned only to gas particles residing within dark matter halos; a particle-based method that assigns HI to all gas particles using a prescription based on the physical properties of the particles. The HI statistical properties are then compared to the observational properties of Damped Lyman-$\alpha$ Absorbers (DLAs) and of lower column density systems and reasonable good agreement is found for all the cases. Among the halo-based method, we further consider two different schemes that aim at reproducing the observed properties of DLAs by distributing HI inside halos: one of this results in a much higher bias for DLAs, in agreement with recent observations, which boosts the 21 cm power spectrum by a factor $\sim 4$ with respect to the other recipe. We compute the 21 cm power spectrum from the simulated HI distribution and calculate the expected signal for both SKA1-mid and SKA1-low configurations at $2.4 \leq z \leq 4$. We find that SKA will be able to detect the 21 cm power spectrum, in the non-linear regime, up to $k\sim 1\,h$/Mpc for SKA1-mid and $k\sim 5\,h$/Mpc for SKA1-low with 100 hours of observations. We also investigate the perspective of imaging the HI distribution. Our findings indicate that SKA1-low could detect the most massive HI peaks with a signal to noise ratio (SNR) higher than 5 for an observation time of about 1000 hours at $z=4$, for a synthesized beam width of $2'$. Detection at redshifts $z\geqslant2.4$ with SKA1-mid would instead require a much longer observation time to achieve a comparable SNR level.

The signature left in quasar spectra by the presence of neutral hydrogen in the Universe allows one to constrain the sum of the neutrino masses with improved sensitivity, with respect to laboratory experiments, and may shed a new light on the neutrino mass hierarchy and on the absolute mass scale of neutrinos. Constraints on cosmological parameters and on the dark energy equation of state can also be derived, from a joint parameter estimation procedure. However, this requires a detailed modeling of the line-of-sight power spectrum of the transmitted flux in the Lyman-Alpha (LyA) forest on scales ranging from a few to hundreds of Mpcs, which in turns demands the inclusion and careful treatment of cosmological neutrinos. To this end, we present here a suite of state-of-the-art hydrodynamical simulations with cold dark matter, baryons and massive neutrinos, specifically targeted for modeling the low-density regions of the IGM as probed by the LyA forest at high-redshift. The simulations span volumes ranging from (25 Mpc/h)^3 to (100 Mpc/h)^3, and are made using either 3 X 192^3~21 millions or 3 X 768^3~1.4 billion particles. The resolution of the various runs can be further enhanced, so that we can reach the equivalent of 3 X 3072^3~87 billion particles in a (100 Mpc/h)^3 box size. The chosen cosmological parameters are compatible with the latest Planck (2013) results, although we also explore the effect of slight variations in the main cosmological and astrophysical parameters. We adopt a particle-type implementation of massive neutrinos, and consider three degenerate species having masses M_nu =0.1, 0.2, 0.3, 0.4 and 0.8 eV, respectively. We improve on previous studies in several ways, in particular with updated routines for IGM radiative cooling and heating processes, and initial conditions based on 2LPT rather than the Zeldovich approximation.

We present a suite of cosmological N-body simulations with cold dark matter and baryons aiming at modeling the low-density regions of the IGM as probed by the Lyman-$\alpha$ forests at high redshift. The simulations are designed to match the requirements imposed by the quality of BOSS and eBOSS data. They are made using either 2x768$^3$ or 2x192$^3$ particles, spanning volumes ranging from (25 Mpc.h$^{-1})^3$ for high-resolution simulations to (100 Mpc.h$^{-1})^3$ for large-volume ones. Using a splicing technique, the resolution is further enhanced to reach the equivalent of simulations with 2x3072$^3$= 58 billion particles in a (100 Mpc.h$^{-1}$)^3 box size, i.e. a mean mass per gas particle of 1.2x10$^5$M_sun.h$^{-1}$. We show that the resulting power spectrum is accurate at the 2% level over the full range from a few Mpc to several tens of Mpc. We explore the effect on the one-dimensional transmitted-flux power spectrum of 4 cosmological parameters ($n_s, \sigma_8, \Omega_m, H_0$) and 2 astrophysical parameters ($T_0, \gamma$) related to the heating rate of the IGM. By varying the input parameters around a central model chosen to be in agreement with the latest Planck results, we built a grid of simulations that allows the study of the impact on the flux power spectrum of these six relevant parameters. We improve upon previous studies by not only measuring the effect of each parameter individually, but also probing the impact of the simultaneous variation of each pair of parameters. We thus provide a full second-order expansion, including cross-terms, around our central model. We check the validity of the second-order expansion with independent simulations obtained either with different cosmological parameters or different seeds. Finally, a comparison to the one-dimensional Ly-$\alpha$ forest power spectrum obtained with BOSS shows an excellent agreement.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We measure the large-scale cross-correlation of quasars with the Lyman alpha forest absorption, using over 164,000 quasars from Data Release 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We extend the previous study of roughly 60,000 quasars from Data Release 9 to larger separations, allowing a measurement of the Baryonic Acoustic Oscillation (BAO) scale along the line of sight $c/(H(z=2.36) ~ r_s) = 9.0 \pm 0.3$ and across the line of sight $D_A(z=2.36) / ~ r_s = 10.8 \pm 0.4$, consistent with CMB and other BAO data. Using the best fit value of the sound horizon from Planck data ($r_s=147.49 Mpc$), we can translate these results to a measurement of the Hubble parameter of $H(z=2.36) = 226 \pm 8 km/s / Mpc$ and of the angular diameter distance of $D_A(z=2.36) = 1590 \pm 60 Mpc$. The measured cross-correlation function and an update of the code to fit the BAO scale (baofit) are made publicly available.

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 measure the sum of the neutrino particle masses using the three-dimensional galaxy power spectrum of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) CMASS galaxy sample. Combined with the cosmic microwave background (CMB), supernova (SN) and additional baryonic acoustic oscillation (BAO) data, we find upper 95 percent confidence limits of the neutrino mass $\Sigma m_{\nu}<0.340$ eV within a flat $\Lambda$CDM background, and $\Sigma m_{\nu}<0.821$ eV, assuming a more general background cosmological model. The number of neutrino species is measured to be $N_{\rm eff}=4.308\pm0.794$ and $N_{\rm eff}=4.032^{+0.870}_{-0.894}$ for these two cases respectively. We study and quantify the effect of several factors on the neutrino measurements, including the galaxy power spectrum bias model, the effect of redshift-space distortion, the cutoff scale of the power spectrum, and the choice of additional data. The impact of neutrinos with unknown masses on other cosmological parameter measurements is investigated. The fractional matter density and the Hubble parameter are measured to be $\Omega_M=0.2796\pm0.0097$, $H_0=69.72^{+0.90}_{-0.91}$ km/s/Mpc (flat $\Lambda$CDM) and $\Omega_M=0.2798^{+0.0132}_{-0.0136}$, $H_0=73.78^{+3.16}_{-3.17}$ km/s/Mpc (more general background model). Based on a Chevallier-Polarski-Linder (CPL) parametrisation of the equation-of-state $w$ of dark energy, we find that $w=-1$ is consistent with observations, even allowing for neutrinos. Similarly, the curvature \Omega_K and the running of the spectral index $\alpha_s$ are both consistent with zero. The tensor-to-scaler ratio is constrained down to $r<0.198$ (95 percent CL, flat $\Lambda$ CDM) and $r<0.440$ (95 percent CL, more general background model).

We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $\Lambda$CDM cosmological scenario - namely an $f(R)$ theory of Modified Gravity (MG) and a cosmological background of massive neutrinos - with the aim to investigate their possible observational degeneracies. We focus on three basic statistics of the large-scale matter distribution, more specifically the nonlinear matter power spectrum, the halo mass function, and the halo bias, for which we determine the deviation with respect to the fiducial $\Lambda$CDM cosmology in the context of both separate and combined simulations of $f(R)$ MG and massive neutrinos scenarios. Our results show that while these two extended models separately determine very prominent and potentially detectable features in all the three statistics, when we allow them to be simultaneously at work these features are strongly suppressed, resulting in much weaker deviations from the standard model's predictions. In particular, when an $f(R)$ gravity model with $f_{R0}=-1\times 10^{-4}$ is combined with a total neutrino mass of $\Sigma _{i}m_{\nu _{i}}=0.4$ eV, the resulting matter power spectrum, halo mass function, and bias at z=0 are found to be consistent with the standard model's predictions at the 10%, 20%, and 5% accuracy levels, respectively. Therefore, our results imply an intrinsic theoretical limit to the effective discriminating power of present and future observational data sets with respect to these widely considered extensions of the standard cosmological scenario in the absence of independent measurements of the neutrino masses from laboratory experiments, even though the high-redshift evolution might still allow to partially break the degeneracy [Abridged].

The acronym ESPRESSO stems for Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations; this instrument will be the next VLT high resolution spectrograph. The spectrograph will be installed at the Combined-Coud\'e Laboratory of the VLT and linked to the four 8.2 m Unit Telescopes (UT) through four optical Coud\'e trains. ESPRESSO will combine efficiency and extreme spectroscopic precision. ESPRESSO is foreseen to achieve a gain of two magnitudes with respect to its predecessor HARPS, and to improve the instrumental radial-velocity precision to reach the 10 cm/s level. It can be operated either with a single UT or with up to four UTs, enabling an additional gain in the latter mode. The incoherent combination of four telescopes and the extreme precision requirements called for many innovative design solutions while ensuring the technical heritage of the successful HARPS experience. ESPRESSO will allow to explore new frontiers in most domains of astrophysics that require precision and sensitivity. The main scientific drivers are the search and characterization of rocky exoplanets in the habitable zone of quiet, nearby G to M-dwarfs and the analysis of the variability of fundamental physical constants. The project passed the final design review in May 2013 and entered the manufacturing phase. ESPRESSO will be installed at the Paranal Observatory in 2016 and its operation is planned to start by the end of the same year.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

In the context of the study of the Integrated Sachs Wolfe effect (ISW), we construct a template of the projected density distribution up to $z\simeq 0.7$ by using the Luminous Galaxies (LGs) from the Sloan Digital Sky Survey DR8. We use a photo-z catalogue trained with more than a hundred thousand galaxies from BOSS in the SDSS DR8 imaging area. We consider two different LG samples whose selection matches that of SDSS-III/BOSS: the LOWZ sample ($z\in [0.15,0.5]$) and the CMASS sample ($z\in[0.4,0.7]$). When building the LG density maps we use the information from star density, survey footprint, seeing conditions, sky emission, dust extinction and airmass to explore the impact of these artifacts on the two LG samples. In agreement with previous studies, we find that the CMASS sample is particularly sensitive to Galactic stars, which dominate the contribution to the auto-angular power spectrum below $\ell=7$. Other potential systematics affect mostly the low multipole range ($\ell\in[2,7]$), but leave fluctuations on smaller scales practically unchanged. The resulting power spectra in the multipole range $\ell\in[2,100]$ for the LOWZ, CMASS and LOWZ+CMASS samples are compatible with linear $\Lambda$CDM expectations and constant bias values of $b=1.98 \pm 0.11$, $2.08\pm0.14$ and $1.88\pm 0.11$, respectively, with no traces of non-Gaussianity: $f_{\rm NL}^{\rm local}=59\pm 75$ at 95% confidence level for the full LOWZ+CMASS sample in the range $\ell\in[4,100]$. After cross-correlating WMAP-9yr data with the LOWZ+CMASS LG density field, the ISW signal is detected at the level of 1.62--1.69$\,\sigma$. While this result is in close agreement with predictions from Monte Carlo simulations in the concordance $\Lambda$CDM model, it cannot rule out by itself an Einstein-de Sitter scenario, and has a moderately low signal compared to previous studies conducted on subsets of this LG sample.

We study the role of feedback from supernovae and black holes in the evolution of the star formation rate function (SFRF) of z~4-7 galaxies. We use a new set of cosmological hydrodynamic simulations, ANGUS (AustraliaN GADGET-3 early Universe Simulations), run with a modified and improved version of the parallel TreePM-smoothed particle hydrodynamics code GADGET-3 called P-GADGET3(XXL), that includes a self-consistent implementation of stellar evolution and metal enrichment. In our simulations both Supernova (SN) driven galactic winds and Active Galactic Nuclei (AGN) act simultaneously in a complex interplay. The SFRF is insensitive to feedback prescription at z>5, meaning that it cannot be used to discriminate between feedback models during reionisation. However, the SFRF is sensitive to the details of feedback prescription at lower redshift. By exploring different SN driven wind velocities and regimes for the AGN feedback, we find that the key factor for reproducing the observed SFRFs is a combination of "strong" SN winds and early AGN feedback in low mass galaxies. Conversely, we show that the choice of initial mass function and inclusion of metal cooling have less impact on the evolution of the SFRF. When variable winds are considered, we find that a non-aggressive wind scaling is needed to reproduce the SFRFs at z>4. Otherwise, the amount of objects with low SFRs is greatly suppressed and at the same time winds are not effective enough in the most massive systems.

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.

We use a suite of N-body simulations that incorporate massive neutrinos as an extra-set of particles to investigate their effect on the halo mass function. We show that for cosmologies with massive neutrinos the mass function of dark matter haloes selected using the spherical overdensity (SO) criterion is well reproduced by the fitting formula of Tinker et al. (2008) once the cold dark matter power spectrum is considered instead of the total matter power, as it is usually done. The differences between the two implementations, i.e. using $P_{\rm cdm}(k)$ instead of $P_{\rm m}(k)$, are more pronounced for large values of the neutrino masses and in the high end of the halo mass function: in particular, the number of massive haloes is higher when $P_{\rm cdm}(k)$ is considered rather than $P_{\rm m}(k)$. As a quantitative application of our findings we consider a Planck-like SZ-clusters survey and show that the differences in predicted number counts can be as large as $30\%$ for $\sum m_\nu = 0.4$ eV. Finally, we use the Planck-SZ clusters sample, with an approximate likelihood calculation, to derive Planck-like constraints on cosmological parameters. We find that, in a massive neutrino cosmology, our correction to the halo mass function produces a shift in the $\sigma_8(\Omega_{\rm m}/0.27)^\gamma$ relation which can be quantified as $\Delta \gamma \sim 0.05$ and $\Delta \gamma \sim 0.14$ assuming one ($N_\nu=1$) or three ($N_\nu=3$) degenerate massive neutrino, respectively. The shift results in a lower mean value of $\sigma_8$ with $\Delta \sigma_8 = 0.01$ for $N_\nu=1$ and $\Delta \sigma_8 = 0.02$ for $N_\nu=3$, respectively. Such difference, in a cosmology with massive neutrinos, would increase the tension between cluster abundance and Planck CMB measurements.

Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.

We review the current state of the theory of large scale structure in a warm dark matter (WDM) cosmological model. In particular, we focus on the non-linear modelling of the matter power spectrum and on the mass function of dark matter haloes. We describe the results of N-body simulations with WDM and mention the effects that could be induced by baryonic physics. We also examine the halo model of large scale structure and its recently suggested modifications for a WDM cosmology, which account for the small scale smoothness of the initial matter density field and better fit the results of N-body simulations. Having described the theoretical models, we discuss the current lower limits on the WDM particle mass, m_w, which correspond to upper limits on the WDM temperature under the assumption that the particles are thermal relics. The best such constraints come from the Ly{\alpha} forest and exclude all masses below 3.3 keV at the 2{\sigma} confidence level. We finally review the forecasts for future lensing surveys, which will be of the same order of magnitude as the already existing constraints from the Ly{\alpha} forest data but explore a different redshift regime.

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$\alpha$ forest. The first method is based on a Fourier transform, and the second on a maximum likelihood estimator. The two methods are independent and have different systematic uncertainties. The determination of the noise level in the data spectra was subject to a novel treatment, because of its significant impact on the derived power spectrum. We applied the two methods to 13,821 quasar spectra from SDSS-III/BOSS DR9 selected from a larger sample of over 60,000 spectra on the basis of their high quality, large signal-to-noise ratio, and good spectral resolution. The power spectra measured using either approach are in good agreement over all twelve redshift bins from $<z> = 2.2$ to $<z> = 4.4$, and scales from 0.001 $\rm(km/s)^{-1}$ to $0.02 \rm(km/s)^{-1}$. We determine the methodological and instrumental systematic uncertainties of our measurements. We provide a preliminary cosmological interpretation of our measurements using available hydrodynamical simulations. The improvement in precision over previously published results from SDSS is a factor 2--3 for constraints on relevant cosmological parameters. For a $\Lambda$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $\sigma_8 =0.83\pm0.03$ and $n_s= 0.97\pm0.02$ based on \ion{H}{i} absorption in the range $2.1<z<3.7$.

We present the Data Release 10 Quasar (DR10Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the first 2.5 years of the survey and that are confirmed as quasars via visual inspection of the spectra. The catalog also includes known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 166,583 quasars (74,454 are new discoveries since SDSS-DR9) detected over 6,373 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with $z>2.15$ (117,668) is $\sim$5 times greater than the number of $z>2.15$ quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII, MgII). The catalog identifies 16,461 broad absorption line quasars and gives their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry with typical accuracy of 0.03 mag and information on the optical morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3,600-10,500\AA\ at a spectral resolution in the range 1,300$<$R$<$2,500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 2,376 quasars that have been identified among the galaxy targets of the SDSS-III/BOSS.

We perform a forecast analysis on how well a Euclid-like photometric galaxy cluster survey will constrain the total neutrino mass and effective number of neutrino species. We base our analysis on the Monte Carlo Markov Chains technique by combining information from cluster number counts and cluster power spectrum. We find that combining cluster data with CMB measurements from Planck improves by more than an order of magnitude the constraint on neutrino masses compared to each probe used independently. For the LCDM+m_nu model the 2 sigma upper limit on total neutrino mass shifts from M_nu < 0.35 eV using cluster data alone to M_nu < 0.031 eV when combined with CMB data. When a non-standard model with N_eff number of neutrino species is considered, we estimate N_eff<3.14 (95% CL), while the bounds on neutrino mass are relaxed to M_nu < 0.040 eV. This accuracy would be sufficient for a 2 sigma detection of neutrino mass even in the minimal normal hierarchy scenario. We also consider scenarios with a constant dark energy equation of state and a non-vanishing curvature. When these models are considered the error on M_nu is only slightly affected, while there is a larger impact of the order of ~ 15 % and ~ 20% respectively on the 2 sigma error bar of N_eff with respect to the standard case. We also treat the LCDM+m_nu+N_eff case with free nuisance parameters, which parameterize the uncertainties on the cluster mass determination. In this case, the upper bounds on M_nu are relaxed by a factor larger than two, M_nu < 0.083 eV (95% CL), hence compromising the possibility of detecting the total neutrino mass with good significance. We thus confirm the potential that a large optical/near-IR cluster survey, like that to be carried out by Euclid, could have in constraining neutrino properties [abridged].

We use a large suite of N-body simulations to study departures from universality in halo abundances and clustering in cosmologies with non-vanishing neutrino masses. To this end, we study how the halo mass function and halo bias factors depend on the scaling variable $\sigma^2(M,z)$, the variance of the initial matter fluctuation field, rather than on halo mass $M$ and redshift $z$ themselves. We show that using the variance of the cold dark matter rather than the total mass field, i.e., $\sigma^2_{cdm}(M,z)$ rather than $\sigma^2_{m}(M,z)$, yields more universal results. Analysis of halo bias yields similar conclusions: When large-scale halo bias is defined with respect to the cold dark matter power spectrum, the result is both more universal, and less scale- or $k$-dependent. These results are used extensively in Papers I and III of this series.

By using a suite of large box-size N-body simulations that incorporate massive neutrinos as an extra set of particles, we investigate the impact of neutrino masses on the spatial distribution of dark matter haloes and galaxies. We compute the bias between the spatial distribution of dark matter haloes and the overall matter and cold dark matter distributions using statistical tools such as the power spectrum and the two-point correlation function. Overall we find a scale-dependent bias on large scales for the cosmologies with massive neutrinos. However, our results indicate that the scale-dependence in the bias is reduced if the latter is computed with respect to the cold dark matter distribution only. We find that the value of the bias on large scales is reasonably well reproduced by the Tinker fitting formula once the linear cold dark matter power spectrum is used, instead of the total matter power spectrum. We investigate whether scale-dependent bias really comes from purely neutrino's effect or from nonlinear gravitational collapse of haloes. For this purpose, we address the $\Omega_\nu$-$\sigma_8$ degeneracy and find that such degeneracy is not perfect, implying that neutrinos imprint a slight scale dependence on the large-scale bias. Finally, by using a simple halo occupation distribution (HOD) model, we investigate the impact of massive neutrinos on the distribution of galaxies within dark matter haloes. We use the main galaxy sample in the Sloan Digital Sky Survey II Data Release 7 to investigate if the small-scale galaxy clustering alone can be used to discriminate among different cosmological models with different neutrino masses. Our results suggest that different choices of the HOD parameters can reproduce the observational measurements relatively well, and we quantify the difference between the values of the HOD parameters between massless and massive neutrino cosmologies.

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)

(Abridged) We investigate two modes of coupling the feedback energy from a central AGN to the neighboring gas in galaxy simulations: kinetic - velocity boost, and thermal - heating. We formulate kinetic feedback models for energy-driven wind (EDW) and momentum-driven wind (MDW), using two free parameters: feedback efficiency epsilon_f, and AGN wind velocity v_w. A novel numerical algorithm is implemented in the SPH code GADGET-3, to prevent the expansion of a hole in the gas distribution around the BH. We perform simulations of isolated evolution and merger of disk galaxies, of Milky-Way mass as well as lower and higher masses. We find that in the isolated galaxy BH kinetic feedback generates intermittent bipolar jet-like gas outflows. We infer that current prescriptions for BH subgrid physics in galaxy simulations can grow the BH to observed values even in an isolated disk galaxy. The BH growth is enhanced in a galaxy merger. Comparing the [M_BH - sigma_star] relation obtained in our simulations with observational data, we conclude that it is possible to find parameter sets for a fit in all the models, except for the case with MDW feedback in a galaxy merger, in which the BH is always too massive. The BH thermal feedback implementation of Springel, Di Matteo & Hernquist (2005) within the multiphase star-formation model is found to have negligible impact on gas properties; and the effect claimed in all previous studies is attributed to gas depletion around the BH by the creation of an artificial hole. The BH mass accretion rate in our simulations exhibit heavy fluctuations. The star formation rate is quenched with feedback by removal of gas. The CGM gas at galactocentric distances (20 - 100)/h kpc are found to give the best metallicity observational diagnostic to distinguish between BH models.

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.

We present the results of observations taken with the X-shooter spectrograph devoted to the study of quasars at z ~ 6. This paper focuses on the properties of metals at high redshift traced, in particular, by the C IV doublet absorption systems. Six objects were observed with resolutions ~27 and 34 km/s in the visual, and 37.5 and 53.5 km/s in the near infrared. We detected 102 C IV lines in the range: 4.35 < z < 6.2 of which 27 are above z ~ 5. Thanks to the characteristics of resolution and spectral coverage of X-shooter, we could also detect 25 Si IV doublets associated with the C IV at z>5. The column density distribution function of the C IV line sample is observed to evolve in redshift for z>~ 5.3, with respect to the normalization defined by low redshift (1.5 < z <4) C IV lines. This behaviour is reflected in the redshift evolution of the C IV cosmic mass density, Omega_CIV, of lines with column density in the range 13.4 < log N(C IV) < 15, which is consistent with a drop of a factor of ~2 for z>~ 5.3. Considering only the stronger C IV lines (13.8 < log N(CIV) < 15), Omega_CIV gently rises by a factor of ~10 between z ~ 6.2 and z~ 1.5 with a possible flattening toward z ~ 0. The increase is well fitted by a power law: Omega_CIV = (2+-1)x10^{-8} [(1+z)/4]^{-3.1+-0.1}. An insight into the properties of the C IV absorbers and their evolution with redshift is obtained by comparing the observed column densities of associated C IV, Si IV and C II absorptions with the output of a set of CLOUDY photo-ionization models. As already claimed by cosmological simulations, we find that C IV is a good tracer of the metallicity in the low-density intergalacitc medium gas at z ~ 5-6 while at z ~ 3 it arises in gas with over density delta ~100.

The injection of turbulence in the circum-galactic medium at redshift z = 2 is investigated using the mesh-based hydrodynamic code Enzo and a subgrid-scale (SGS) model for unresolved turbulence. Radiative cooling and heating by a uniform Ultraviolet (UV) background are included in our runs and compared with the effect of turbulence modelling. Mechanisms of gas exchange between galaxies and the surrounding medium, as well as metal enrichment, are not taken into account, and turbulence is here driven solely by structure formation (mergers and shocks). We find that turbulence, both at resolved and SGS scales, impacts mostly the warm-hot intergalactic medium (WHIM), with temperature between 10^5 and 10^7 K, mainly located around collapsed and shock heated structures, and in filaments. Typical values of the ratio of turbulent to thermal pressure is 0.1 in the WHIM, corresponding to a volume-weighted average of the SGS turbulent to thermal Doppler broadening b_t / b_therm = 0.26, on length scales below the grid resolution of 25 kpc/h. In the diffuse intergalactic medium (IGM), defined in a range of baryon overdensity \delta\ between 1 and 50, the importance of turbulence is smaller, but grows as a function of gas density, and the Doppler broadening ratio is fitted by the function b_t / b_therm = 0.023 \delta^{0.58}.

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

Planck has produced detailed all-sky observations over nine frequency bands between 30 and 857 GHz. These observations allow robust reconstruction of the primordial cosmic microwave background (CMB) temperature fluctuations over nearly the full sky, as well as new constraints on Galactic foregrounds. This paper describes the component separation framework adopted by Planck. We test four foreground-cleaned CMB maps derived using qualitatively different component separation algorithms. The quality of our reconstructions is evaluated through detailed simulations and internal comparisons, and shown through various tests to be internally consistent and robust for CMB power spectrum and cosmological parameter estimation up to l = 2000. The parameter constraints on LambdaCDM cosmologies derived from these maps are consistent with those presented in the cross-spectrum based Planck likelihood analysis. We choose two of the CMB maps for specific scientific goals. We also present maps and frequency spectra of the Galactic low-frequency, CO, and thermal dust emission. The component maps are found to provide a faithful representation of the sky, as evaluated by simulations. For the low-frequency component, the spectral index varies widely over the sky, ranging from about beta = -4 to -2. Considering both morphology and prior knowledge of the low frequency components, the index map allows us to associate a steep spectral index (beta < -3.2) with strong anomalous microwave emission, corresponding to a spinning dust spectrum peaking below 20 GHz, a flat index of beta > -2.3 with strong free-free emission, and intermediate values with synchrotron emission.

We quantify the impact of massive neutrinos on the statistics of low density regions in the intergalactic medium (IGM) as probed by the Lyman-alpha forest at redshifts z=2.2--4. Based on mock but realistic quasar (QSO) spectra extracted from hydrodynamic simulations with cold dark matter, baryons and neutrinos, we find that the probability distribution of weak Lyman-alpha absorption features, as sampled by Lyman-alpha flux regions at high transmissivity, is strongly affected by the presence of massive neutrinos. We show that systematic errors affecting the Lyman-alpha forest reduce but do not erase the neutrino signal. Using the Fisher matrix formalism, we conclude that the sum of the neutrino masses can be measured, using the method proposed in this paper, with a precision smaller than 0.4 eV using a catalog of 200 high resolution (S/N~100) QSO spectra. This number reduces to 0.27 eV by making use of reasonable priors in the other parameters that also affect the statistics of the high transitivity regions of the Lyman-alpha forest. The constraints obtained with this method can be combined with independent bounds from the CMB, large scale structures and measurements of the matter power spectrum from the Lyman-alpha forest to produce tighter upper limits on the sum of the masses of the neutrinos.

Cosmography is the tool that makes possible to untie the interpretation of cosmological observations from the definition of any dynamical prior. We review the constraints on the cosmographic parameter obtained using the most thorough data set ensemble available. We focus on some specific topics about the statistically based selection of the most stringent fitting expansion.

We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations. Our set of simulations explore the properties of neutrinos in a reference $\Lambda$CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass $10^{11}-10^{15}$ $h^{-1}$M$_{\odot}$, over a redshift range $z=0-2$. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than $\sim 10^{13.5}h^{-1}$M$_{\odot}$. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.

We report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $\Lambda$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

(Abridged) We explore new observationally-constrained sub-resolution models of galactic outflows and investigate their impact on the circumgalactic medium (CGM) over redshifts z = 2 - 4. We perform cosmological hydrodynamic simulations, including star formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where the outflow velocity has a positive correlation with galactocentric distance (r), and a RVW with additional dependence on halo mass (RVWb). Overall, we find that the outflows expel metal-enriched gas away from galaxies, significantly quench star formation, and enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy centers are most sensitive to the choice of the wind model for halo masses (10^9 - 10^11) M_sun. We infer that the RVWb model is similar to the NW case, except that it substantially enriches the CGM: the carbon metallicity (Z_C) is 10 times higher in RVWb than in NW at r > R_200; and the warm gas of 10^4 - 10^5 K and delta < 100 is enriched to 50 times higher than in NW. We also find that the impact of models CW and RVWa are similar, with the following differences. RVWa causes a greater suppression of star formation rate at z < 5, and has a higher fraction of low-density (delta < 10), warm-hot (10^4 - 10^6 K) gas than in CW. Outflows in CW produce a higher and earlier enrichment of some IGM phases than RVWa. By visual inspection, we note that the RVWa model shows more pronounced bipolar outflows and galactic disks. We present fitting formulae for [Z_C-delta] and [Z_C-r], also for the abundance of CIV as a function of r. We predict observational diagnostics to distinguish between different outflow scenarios: Z_C of the CGM gas at r = (30 - 300) kpc/h comoving, and CIV fraction of the inner gas at r < (4 - 5) kpc/h comoving.

We describe fitting methods developed to analyze fluctuations in the Lyman-{\alpha} forest and measure the parameters of baryon acoustic oscillations (BAO). We apply our methods to BOSS Data Release 9. Our method is based on models of the three-dimensional correlation function in physical coordinate space, and includes the effects of redshift-space distortions, anisotropic non-linear broadening, and broadband distortions. We allow for independent scale factors along and perpendicular to the line of sight to minimize the dependence on our assumed fiducial cosmology and to obtain separate measurements of the BAO angular and relative velocity scales. Our fitting software and the input files needed to reproduce our main BOSS Data Release 9 results are publicly available.

The forest of Lyman-alpha absorption lines seen in the spectra of distant quasars has become an important probe of the distribution of matter in the Universe. We use large, hydrodynamical simulations from the OWLS project to investigate the effect of feedback from galaxy formation on the probability distribution function and the power spectrum of the Lyman-alpha transmitted flux. While metal-line cooling is unimportant, both galactic outflows from massive galaxies driven by active galactic nuclei and winds from low-mass galaxies driven by supernovae have a substantial impact on the flux statistics. At redshift z=2.25, the effects on the flux statistics are of a similar magnitude as the statistical uncertainties of published data sets. The changes in the flux statistics are not due to differences in the temperature-density relation of the photo-ionised gas. Instead, they are caused by changes in the density distribution and in the fraction of hot, collisionally ionised gas. It may be possible to disentangle astrophysical and cosmological effects by taking advantage of the fact that they induce different redshift dependencies. In particular, the magnitude of the feedback effects appears to decrease rapidly with increasing redshift. Analyses of Lyman-alpha forest data from surveys that are currently in process, such as BOSS/SDSS-III and X-Shooter/VLT, must take galactic winds into account.

We present the BOSS Lyman-alpha (Lya) Forest Sample from SDSS Data Release 9, comprising 54,468 quasar spectra with zqso > 2.15 suitable for Lya forest analysis. This data set probes the intergalactic medium with absorption redshifts 2.0 < z_alpha < 5.7 over an area of 3275 square degrees, and encompasses an approximate comoving volume of 20 h^-3 Gpc^3. With each spectrum, we have included several products designed to aid in Lya forest analysis: improved sky masks that flag pixels where data may be unreliable, corrections for known biases in the pipeline estimated noise, masks for the cores of damped Lya systems and corrections for their wings, and estimates of the unabsorbed continua so that the observed flux can be converted to a fractional transmission. The continua are derived using a principal component fit to the quasar spectrum redwards of restframe Lya (lambda > 1216 Ang), extrapolated into the forest region and normalized by a linear function to fit the expected evolution of the Lya forest mean-flux. The estimated continuum errors are ~5% rms. We also discuss possible systematics arising from uncertain spectrophotometry and artifacts in the flux calibration; global corrections for the latter are provided. Our sample provides a convenient starting point for users to analyze clustering in BOSS Lya forest data, and it provides a fiducial data set that can be used to compare results from different analyses of baryon acoustic oscillations in the Lya forest. The full data set is available from the SDSS-III DR9 web site.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

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

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We measure the quasar two-point correlation function over the redshift range 2.2<z<2.8 using data from the Baryon Oscillation Spectroscopic Survey. We use a homogeneous subset of the data consisting of 27,129 quasars with spectroscopic redshifts---by far the largest such sample used for clustering measurements at these redshifts to date. The sample covers 3,600 square degrees, corresponding to a comoving volume of 9.7(Gpc/h)^3 assuming a fiducial LambdaCDM cosmology, and it has a median absolute i-band magnitude of -26, k-corrected to z=2. After accounting for redshift errors we find that the redshift space correlation function is fit well by a power-law of slope -2 and amplitude s_0=(9.7\pm 0.5)Mpc/h over the range 3<s<25Mpc/h. The projected correlation function, which integrates out the effects of peculiar velocities and redshift errors, is fit well by a power-law of slope -1 and r_0=(8.4\pm 0.6)Mpc/h over the range 4<R<16Mpc/h. There is no evidence for strong luminosity or redshift dependence to the clustering amplitude, in part because of the limited dynamic range in our sample. Our results are consistent with, but more precise than, previous measurements at similar redshifts. Our measurement of the quasar clustering amplitude implies a bias factor of b~3.5 for our quasar sample. We compare the data to models to constrain the manner in which quasars occupy dark matter halos at z~2.4 and infer that such quasars inhabit halos with a characteristic mass of <M>~10^{12}Msun/h with a duty cycle for the quasar activity of 1 per cent.

We have investigated the Cosmic Infrared Background (CIB) anisotropies in the framework of the physical evolutionary model for proto-spheroidal galaxies by Granato et al. (2004). After having re-calibrated the cumulative flux function $dS/dz$ at $\lambda \ge 850\,\mu$m using the available determinations of the shot noise amplitude (the original model already correctly reproduces it at shorter wavelengths) the CIB power spectra at wavelengths from $250\,\mu$m to $2\,$mm measured by {\it Planck}, {\it Herschel}, SPT and ACT experiments have been fitted using the halo model with only 2 free parameters, the minimum halo mass and the power-law index of the mean occupation function of satellite galaxies. The best-fit {\it minimum} halo mass is $\log(M_{\rm min}/M_\odot) = 12.24 \pm 0.06$, higher than, but consistent within the errors, with the estimate by Amblard et al. (2011) and close to the estimate by Planck Collaboration (2011). The redshift evolution of the volume emissivity of galaxies yielded by the model is found to be consistent with that inferred from the data. The derived {\it effective} halo mass, $M_{\rm eff} \simeq 5\times 10^{12}\,M_\odot$, of $z\simeq 2$ sub-millimeter galaxies is close to that estimated for the most efficient star-formers at the same redshift. The effective bias factor and the comoving clustering radius at $z\simeq 2$ yielded by the model are substantially lower than those found for a model whereby the star formation is fueled by steady gas accretion, but substantially higher than those found for a merging-driven galaxy evolution with a top-heavy initial mass function.

We investigate the thermal history of the intergalactic medium (IGM) in the redshift interval z=1.7--3.2 by studying the small-scale fluctuations in the Lyman alpha forest transmitted flux. We apply a wavelet filtering technique to eighteen high resolution quasar spectra obtained with the Ultraviolet and Visual Echelle Spectrograph (UVES), and compare these data to synthetic spectra drawn from a suite of hydrodynamical simulations in which the IGM thermal state and cosmological parameters are varied. From the wavelet analysis we obtain estimates of the IGM thermal state that are in good agreement with other recent, independent wavelet-based measurements. We also perform a reanalysis of the same data set using the Lyman alpha forest flux probability distribution function (PDF), which has previously been used to measure the IGM temperature-density relation. This provides an important consistency test for measurements of the IGM thermal state, as it enables a direct comparison of the constraints obtained using these two different methodologies. We find the constraints obtained from wavelets and the flux PDF are formally consistent with each other, although in agreement with previous studies, the flux PDF constraints favour an isothermal or inverted IGM temperature-density relation. We also perform a joint analysis by combining our wavelet and flux PDF measurements, constraining the IGM thermal state at z=2.1 to have a temperature at mean density of T0/[10^3 K]=17.3 +/- 1.9 and a power-law temperature-density relation exponent gamma=1.1 +/- 0.1 (1 sigma). Our results are consistent with previous observations that indicate there may be additional sources of heating in the IGM at z<4.

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 a measurement of the Lyman alpha flux probability distribution function (PDF) measured from a set of eight high resolution quasar spectra with emission redshifts at 3.3 < z < 3.8. We carefully study the effect of metal absorption lines on the shape of the PDF. Metals have a larger impact on the PDF measurements at lower redshift, where there are fewer Lyman alpha absorption lines. This may be explained by an increase in the number of metal lines which are blended with Lyman alpha absorption lines toward higher redshift, but may also be due to the presence of fewer metals in the intergalactic medium with increasing lookback time. We also provide a new measurement of the redshift evolution of the effective optical depth, tau_eff, at 2.8 < z < 3.6, and find no evidence for a deviation from a power law evolution in the log(tau_eff)-log(1+z) plane. The flux PDF measurements are furthermore of interest for studies of the thermal state of the intergalactic medium (IGM) at z ~ 3 . By comparing the PDF to state-of-the-art cosmological hydrodynamical simulations, we place constraints on the temperature of the IGM and compare our results with previous measurements of the PDF at lower redshift. At redshift z=3, our new PDF measurements are consistent with an isothermal temperature-density relation, T=T_0 Delta^{gamma-1}, with a temperature at the mean density of T_0 = 19250 +/- 4800 K and a slope gamma=0.90+/-0.21 (1 sigma uncertainties). In comparison, joint constraints with previous PDF measurements at z<3 favour an inverted (gamma<1) temperature-density relation with T_0=17900 +/- 3500 K and gamma=0.70 +/- 0.12, in broad agreement with previous analyses.

We perform a cosmographic analysis using several cosmological observables such as the luminosity distance moduli, the volume distance, the angular diameter distance and the Hubble parameter. These quantities are determined using different data sets: Supernovae type Ia and Gamma Ray Bursts, the Baryonic Acoustic Oscillations, the cosmic microwave background power spectrum and the Hubble parameter as measured from surveys of galaxies. This data set allows to put constraints on the cosmographic expansion with unprecedented precision. We also present forecasts for the coefficients of the kinematic expansion using future but realistic data sets: constraints on the coefficients of the expansions are likely to improve by a factor ten with the upcoming large scale structure probes. Finally, we derive the set of the cosmographic parameters for several cosmological models (including $\Lambda$CDM) and compare them with our best fit set. While distance measurements are unable to discriminate among these models, we show that the inclusion of the Hubble data set leads to strong constraints on the lowest order coefficients and in particular it is incompatible with $\Lambda$CDM at 3-$\sigma$ confidence level. We discuss the reliability of this determination and suggest further observations which might be of crucial importance for the viability of cosmographic tests in the next future.

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

We study the effect of the non-Gaussianity induced by gravitational evolution upon the statistical properties of absorption in quasar (QSO) spectra. Using the generic hierarchical ansatz and the lognormal approximation we derive the analytical expressions for the one-point PDF as well as for the joint two-point probability distribution (2PDF) of transmitted fluxes in two neighbouring QSOs. These flux PDFs are constructed in 3D as well as in projection (i.e. in 2D). The PDFs are constructed by relating the lower-order moments, i.e. cumulants and cumulant correlators, of the fluxes to the 3D neutral hydrogen distribution which is, in turn, expressed as a function of the underlying dark matter distribution. The lower-order moments are next modelled using a generating function formalism in the context of a {\em minimal tree-model} for the higher-order correlation hierarchy. These different approximations give nearly identical results for the range of redshifts probed, and we also find a very good agreement between our predictions and outputs of hydrodynamical simulations. The formalism developed here for the joint statistics of flux-decrements concerning two lines of sight can be extended to multiple lines of sight, which could be particularly important for the 3D reconstruction of the cosmic web from QSO spectra (e.g. in the BOSS survey). These statistics probe the underlying projected neutral hydrogen field and are thus linked to "hot-spots" of absorption. The results for the PDF and the bias presented here use the same functional forms of scaling functions that have previously been employed for the modelling of other cosmological observation such as the Sunyaev-Zel'dovich effect.

Cosmological neutrinos strongly affect the evolution of the largest structures in the Universe, i.e. galaxies and galaxy clusters. We use large box-size full hydrodynamic simulations to investigate the non-linear effects that massive neutrinos have on the spatial properties of cold dark matter (CDM) haloes. We quantify the difference with respect to the concordance LambdaCDM model of the halo mass function and of the halo two-point correlation function. We model the redshift-space distortions and compute the errors on the linear distortion parameter beta introduced if cosmological neutrinos are assumed to be massless. We find that, if not taken correctly into account and depending on the total neutrino mass, these effects could lead to a potentially fake signature of modified gravity. Future nearly all-sky spectroscopic galaxy surveys will be able to constrain the neutrino mass if it is larger than 0.6 eV, using beta measurements alone and independently of the value of the matter power spectrum normalisation. In combination with other cosmological probes, this will strengthen neutrino mass constraints and help breaking parameter degeneracies.

We investigate the non-linear evolution of the matter power spectrum by using a large set of high-resolution N-body/hydrodynamic simulations. The linear matter power in the initial conditions is consistently modified to accommodate warm dark matter particles which induce a small scale cut-off in the power as compared to standard cold dark matter scenarios. The impact of such thermal relics is addressed at small scales with k > 1 h/Mpc and at z < 5, which are particularly important for the next generation of Lyman-alpha forest, weak lensing and galaxy clustering surveys. We quantify the mass and redshift dependence of the warm dark matter non-linear matter power and we provide a fitting formula which is accurate at the ~2% level below z=3 and for masses m_wdm > 0.5 keV. The role of baryonic physics (cooling, star formation and feedback recipes) on the warm dark matter induced suppression is also quantified. Furthermore, we compare our findings with the halo model and show their impact on the cosmic shear power spectra.

Semi-analytical methods, based on Eulerian perturbation theory, are a promising tool to follow the time evolution of cosmological perturbations at small redshifts and at mildly nonlinear scales. All these schemes are based on two approximations: the existence of a smoothing scale and the single-stream approximation, where velocity dispersion of the dark matter fluid, as well as higher moments of the particle distributions, are neglected. Despite being widely recognized, these two assumptions are, in principle, incompatible, since any finite smoothing scale gives rise to velocity dispersion and higher moments at larger scales. We describe a new approach to perturbation theory, where the Vlasov and fluid equations are derived in presence of a finite coarse-graining scale: this allows a clear separation between long and short distance modes and leads to a hybrid approach where the former are treated perturbatively and the effect of the latter is encoded in external source terms for velocity, velocity dispersion, and all the higher order moments, which can be computed from N-body simulations. We apply the coarse-grained perturbation theory to the computation of the power spectrum and the cross-spectrum between density and velocity dispersion, and compare the results with N-body simulations, finding good agreement.

In this work, starting from 21 months of data from the Fermi-Large Area Telescope, we derive maps of the residual isotropic gamma-ray emission, a relevant fraction of which is expected to be contributed by the extragalactic diffuse gamma-ray background. We compute the angular two-point auto-correlation function of the residual Fermi-LAT maps at energies E>1GeV, E>3GeV and E>30GeV well above the Galactic plane and find no significant correlation signal. This is, indeed, what is expected if the EGB were contributed by BL Lacertae, Flat Spectrum Radio Quasars or star-forming galaxies, since, in this case, the predicted signal is very weak. Then, we search for the Integrated Sachs-Wolfe signature by cross-correlating the Fermi-LAT maps with the WMAP7-Cosmic Microwave Background map. We find a cross-correlation consistent with zero, even though the expected signal is larger than that of the EGB auto-correlation. Finally, in an attempt to constrain the nature of the gamma-ray background we cross-correlate the Fermi-LAT maps with the angular distributions of objects that may contribute to the EGB: QSOs in the SDSS-DR6 catalog, NVSS galaxies, 2MASS galaxies and LRG in the SDSS catalog. The cross-correlation is always consistent with zero, in agreement with theoretical expectations, but we find (with low statistical significance) some interesting features that may indicate that some specific classes of objects contribute to the EGB. A chi2 analysis confirms that the correlation properties of the 21-month data do not provide strong constraints of the EGB origin. However, the results suggest that the situation will significantly improve with the 5- and 10-year Fermi-LAT data. The future EGB analysis will then allow placing significant constraints on the nature of the EGB and might provide in addition a detection of the ISW signal. (Abridged)

In this paper we measure the angular power spectra $C_\ell$ of three high-redshift large-scale structure probes: the radio sources from the NRAO VLA Sky Survey (NVSS), the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs) and the MegaZ-LRG (DR7), the final SDSS II Luminous Red Galaxy (LRG) photometric redshift survey. We perform a global analysis of the constraints on the amplitude of primordial non-Gaussianity from these angular power spectra, as well as from their cross-correlation power spectra with the cosmic microwave background (CMB) temperature map. In particular, we include non-Gaussianity of the type arising from single-field slow roll, multifields, curvaton (local type), and those which effects on the halo clustering can be described by the equilateral template (related to higher-order derivative type non-Gaussianity) and by the enfolded template (related to modified initial state or higher-derivative interactions). When combining all data sets, we obtain limits of $f_{\rm NL}=48\pm20$, $f_{\rm NL}=50\pm265$ and $f_{\rm NL}=183\pm95$ at 68% confidence level for local, equilateral and enfolded templates, respectively. Furthermore, we explore the constraint on the cubic correction $g_{\rm NL}\phi^3$ on the bias of dark matter haloes and obtain a limit of $-1.2\times10^5<g_{\rm NL}<11.3\times10^5$ at 95% confidence level.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

We assess the possibility to detect the warm-hot intergalactic medium (WHIM) in emission and to characterize its physical conditions and spatial distribution through spatially resolved X-ray spectroscopy, in the framework of the recently proposed DIOS, EDGE, Xenia, and ORIGIN missions, all of which are equipped with microcalorimeter-based detectors. For this purpose we analyze a large set of mock emission spectra, extracted from a cosmological hydrodynamical simulation. These mock X-ray spectra are searched for emission features showing both the OVII K alpha triplet and OVIII Ly alpha line, which constitute a typical signature of the warm hot gas. Our analysis shows that 1 Ms long exposures and energy resolution of 2.5 eV will allow us to detect about 400 such features per deg^2 with a significance >5 sigma and reveals that these emission systems are typically associated with density ~100 above the mean. The temperature can be estimated from the line ratio with a precision of ~20%. The combined effect of contamination from other lines, variation in the level of the continuum, and degradation of the energy resolution reduces these estimates. Yet, with an energy resolution of 7 eV and all these effects taken into account, one still expects about 160 detections per deg^2. These line systems are sufficient to trace the spatial distribution of the line-emitting gas, which constitute an additional information, independent from line statistics, to constrain the poorly known cosmic chemical enrichment history and the stellar feedback processes.

Current results from the Lyman alpha forest assume that the primordial power spectrum of density perturbations follows a simple power law form, with running. We present the first analysis of Lyman alpha data to study the effect of relaxing this strong assumption on primordial and astrophysical constraints. We perform a large suite of numerical simulations, using them to calibrate a minimally parametric framework for describing the power spectrum. Combined with cross-validation, a statistical technique which prevents over-fitting of the data, this framework allows us to reconstruct the power spectrum shape without strong prior assumptions. We find no evidence for deviation from scale-invariance; our analysis also shows that current Lyman alpha data do not have sufficient statistical power to robustly probe the shape of the power spectrum at these scales. In contrast, the ongoing Baryon Oscillation Sky Survey (BOSS) will be able to do so with high precision. Furthermore, this near-future data will be able to break degeneracies between the power spectrum shape and astrophysical parameters.

We present an analysis of high-resolution hydrodynamical N-body simulations of coupled dark energy cosmologies which focusses on the statistical properties of the transmitted Lyman-alpha flux in the high-redshift intergalactic medium (IGM). In these models the growth of the diffuse cosmic web differs from the standard LCDM case: the density distribution is skewed towards underdense regions and the matter power spectra are typically larger (in a scale dependent way). These differences are also appreciable in the Lyman-alpha flux and are larger than 5% (10%) at z=2-4 in the flux probability distribution function (pdf) for high transmissivity regions and for values of the coupling parameter \beta = 0.08 (\beta = 0.2). The flux power spectrum is also affected at the ~2% (~ 5-10%) level for \beta = 0.08 (\beta = 0.2) in a redshift dependent way. We infer the behaviour of flux pdf and flux power for a reasonable range of couplings and present constraints using present high and low resolution data sets. We find an upper limit \beta < 0.15 (at 2 sigma confidence level), which is obtained using only IGM data and is competitive with those inferred from other large scale structure probes.

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.

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.

X-shooter, with its characteristics of resolution, spectral coverage and efficiency, provides a unique opportunity to obtain spectra of the highest-redshift quasars (z ~ 6) that will allow us to carry out successful investigations on key cosmological issues, from the details of the re-ionization process, to the evolution of the first galaxies and AGNs. In this paper, we present the spectra of three z ~ 6 quasars: one obtained during the commissioning of X-shooter and two in the context of our ongoing GTO programme. Combining this sample with data in the literature, we update the value of the C IV cosmic mass density in the range 4.5 < z < 5, confirming the constant trend with redshift between 2.5 and 5.

Observations of quasar pairs represent a powerful tool in the study of the intergalactic medium, providing in particular information on the quasar radiative feedback not only along the line of sight, but also in the transverse direction. In this paper, we present the spectra of the four quasar pairs we observed so far with X-shooter, and we discuss two individual spectra affected by the presence of a close quasar in the foreground.

We investigate the properties of triply ionized Carbon (CIV) in the Intergalactic Medium using a set of high-resolution and large box-size cosmological hydrodynamic simulations of a $\Lambda$CDM model. We rely on a modification of the GADGET-2 code that self-consistently follows the metal enrichment mechanism by means of a detailed chemical evolution model. We focus on several numerical implementations of galactic feedback: galactic winds in the energy driven and momentum driven prescriptions and Active Galactic Nuclei (AGN) powered by gas accretion onto massive black holes. We extract mock IGM transmission spectra in neutral hydrogen (HI) and CIV and perform Voigt profile fitting. The results are then compared with high-resolution quasar (QSO) spectra obtained with the UVES spectrograph at the VLT and the HIRES spectrograph at Keck. We find that feedback has little impact on statistics related to the neutral hydrogen, while CIV is more affected by galactic winds and/or AGN feedback. When the same analysis is performed over observed and simulated CIV lines, we find reasonables good agreement between data and simulations over the column density range $N_{\rm CIV}=10^{12.5-15}$ cm$^{-2}$. Also the CIV line-widths distribution appears to be in agreement with the observed values, while the HI Doppler parameters, $b_{\rm HI}$, are in general too large showing that the diffuse cosmic web is heated more than what is inferred by observations. The simulation without feedback fails in reproducing the CIV systems at high column densities at all redshift, while the AGN feedback case agrees with observations only at $z<3$, when this form of feedback is particularly effective. We also present scatter plots in the $b-N$ and in the $N_{\rm CIV}-N_{\rm HI}$ planes, showing that there is rough agreement between observations and simulations only when feedback is taken into account.

We present an analysis of the constraints on the amplitude of primordial non-Gaussianity of local type described by the dimensionless parameter $f_{\rm NL}$. These constraints are set by the auto-correlation functions (ACFs) of two large scale structure probes, the radio sources from NRAO VLA Sky Survey (NVSS) and the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs), as well as by their cross-correlation functions (CCFs) with the cosmic microwave background (CMB) temperature map (Integrated Sachs-Wolfe effect). Several systematic effects that may affect the observational estimates of the ACFs and of the CCFs are investigated and conservatively accounted for. Our approach exploits the large-scale scale-dependence of the non-Gaussian halo bias. The derived constraints on {$f_{\rm NL}$} coming from the NVSS CCF and from the QSO ACF and CCF are weaker than those previously obtained from the NVSS ACF, but still consistent with them. Finally, we obtain the constraints on $f_{\rm NL}=53\pm25$ ($1\,\sigma$) and $f_{\rm NL}=58\pm24$ ($1\,\sigma$) from NVSS data and SDSS DR6 QSO data, respectively.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

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.

The NRAO VLA Sky Survey (NVSS) is the only dataset that allows an accurate determination of the auto-correlation function (ACF) on angular scales of several degrees for Active Galactic Nuclei (AGNs) at typical redshifts $z \simeq 1$. Surprisingly, the ACF is found to be positive on such large scales while, in the framework of the standard hierarchical clustering scenario with Gaussian primordial perturbations it should be negative for a redshift-independent effective halo mass of order of that found for optically-selected quasars. We show that a small primordial non-Gaussianity can add sufficient power on very large scales to account for the observed NVSS ACF. The best-fit value of the parameter $f_{\rm NL}$, quantifying the amplitude of primordial non-Gaussianity of local type is $f_{\rm NL}=62 \pm 27$ ($1\,\sigma$ error bar) and $25<f_{\rm NL}<117$ ($2\,\sigma$ confidence level), corresponding to a detection of non-Gaussianity significant at the $\sim 3\,\sigma$ confidence level. The minimal halo mass of NVSS sources is found to be $M_{\rm min}=10^{12.47\pm0.26}h^{-1}M_{\odot}$ ($1\,\sigma$) strikingly close to that found for optically selected quasars. We discuss caveats and possible physical and systematic effects that can impact on the results.

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.

We constrain the parameters describing the kinematical state of the universe using a cosmographic approach, which is fundamental in that it requires a very minimal set of assumptions (namely to specify a metric) and does not rely on the dynamical equations for gravity. On the data side, we consider the most recent compilations of Supernovae and Gamma Ray Bursts catalogues. This allows to further extend the cosmographic fit up to $z = 6.6$, i.e. up to redshift for which one could start to resolve the low z degeneracy among competing cosmological models. In order to reliably control the cosmographic approach at high redshifts, we adopt the expansion in the improved parameter $y = z/(1+z)$. This series has the great advantage to hold also for $z > 1$ and hence it is the appropriate tool for handling data including non-nearby distance indicators. We find that Gamma Ray Bursts, probing higher redshifts than Supernovae, have constraining power and do require (and statistically allow) a cosmographic expansion at higher order than Supernovae alone. Exploiting the set of data from Union and GRBs catalogues, we show (for the first time in a purely cosmographic approach parametrized by deceleration $q_0$, jerk $j_0$, snap $s_0$) a definitively negative deceleration parameter $q_0$ up to the 3$\sigma$ confidence level. We present also forecasts for realistic data sets that are likely to be obtained in the next few years.

The cosmic evolution of the metal content of the intergalactic medium puts stringent constraints on the properties of galactic outflows and on the nature of UV background. In this paper, we present a new measure of the redshift evolution of the mass density of C IV, Omega_CIV, in the interval 1.5 < z < 4 based on a sample of more than 1500 C IV lines with column densities 10^{12} < N(CIV) < 10^{15} cm^{-2}. This sample more than doubles the absorption redshift path covered in the range z<2.5 by previous samples. The result shows a significant increase of Omega_CIV towards the lower redshifts at variance with the previously pictured constant behaviour.

We analyse the evolution of the properties of the low-redshift Intergalactic Medium (IGM) using high-resolution hydrodynamic simulations that include a detailed chemical evolution model. We focus on the effects that two different forms of energy feedback, strong galactic winds driven by supernova explosion and Active Galactic Nuclei (AGN) powered by gas accretion onto super-massive black holes (BHs), have on the thermo- and chemo-dynamical properties of of the low redshift IGM. We find that feedback associated to winds (W) and BHs leave distinct signatures in both the chemical and thermal history of the baryons, especially at redshift z<3 [..] We present results for the enrichment in terms of mass and metallicity distributions for the WHIM phase, both as a function of density and temperature. Finally, we compute the evolution of the relative abundances between different heavy elements, namely Oxygen, Carbon and Iron. While both C/O and O/Fe evolve differently at high redshifts for different feedback models, their values are similar at z=0 [..]. The sensitivity of WHIM properties on the implemented feedback scheme could be important both for discriminating between different feedback physics and for detecting the WHIM with future far-UV and X-ray telescopes.

We expect a detectable correlation between two seemingly unrelated quantities: the four point function of the cosmic microwave background (CMB) and the amplitude of flux decrements in quasar (QSO) spectra. The amplitude of CMB convergence in a given direction measures the projected surface density of matter. Measurements of QSO flux decrements trace the small-scale distribution of gas along a given line-of-sight. While the cross-correlation between these two measurements is small for a single line-of-sight, upcoming large surveys should enable its detection. This paper presents analytical estimates for the signal to noise (S/N) for measurements of the cross-correlation between the flux decrement and the convergence and for measurements of the cross-correlation between the variance in flux decrement and the convergence. For the ongoing BOSS (SDSS III) and Planck surveys, we estimate an S/N of 30 and 9.6 for these two correlations. For the proposed BigBOSS and ACTPOL surveys, we estimate an S/N of 130 and 50 respectively. Since the cross-correlation between the variance in flux decrement and the convergence is proportional to the fourth power of $\sigma_8$, the amplitude of these cross-correlations can potentially be used to measure the amplitude of $\sigma_8$ at z~2 to 2.5% with BOSS and Planck and even better with future data sets. These measurements have the potential to test alternative theories for dark energy and to constrain the mass of the neutrino. The large potential signal estimated in our analytical calculations motivate tests with non-linear hydrodynamical simulations and analyses of upcoming data sets.

In this paper we rely on the quasar (QSO) catalog of the Sloan Digital Sky Survey Data Release Six (SDSS DR6) of about one million photometrically selected QSOs to compute the Integrated Sachs-Wolfe (ISW) effect at high redshift, aiming at constraining the behavior of the expansion rate and thus the behaviour of dark energy at those epochs. This unique sample significantly extends previous catalogs to higher redshifts while retaining high efficiency in the selection algorithm. We compute the auto-correlation function (ACF) of QSO number density from which we extract the bias and the stellar contamination. We then calculate the cross-correlation function (CCF) between QSO number density and Cosmic Microwave Background (CMB) temperature fluctuations in different subsamples: at high z>1.5 and low z<1.5 redshifts and for two different choices of QSO in a conservative and in a more speculative analysis. We find an overall evidence for a cross-correlation different from zero at the 2.7\sigma level, while this evidence drops to 1.5\sigma at z>1.5. We focus on the capabilities of the ISW to constrain the behaviour of the dark energy component at high redshift both in the \LambdaCDM and Early Dark Energy cosmologies, when the dark energy is substantially unconstrained by observations. At present, the inclusion of the ISW data results in a poor improvement compared to the obtained constraints from other cosmological datasets. We study the capabilities of future high-redshift QSO survey and find that the ISW signal can improve the constraints on the most important cosmological parameters derived from Planck CMB data, including the high redshift dark energy abundance, by a factor \sim 1.5.

We present a theoretical estimate for a new observable: the cross-correlation between the Lyman-alpha flux fluctuations in quasar (QSO) spectra and the convergence of the cosmic microwave background (CMB) as measured along the same line-of-sight. As a first step toward the assessment of its detectability, we estimate the signal-to-noise ratio using linear theory. Although the signal-to-noise is small for a single line-of-sight and peaks at somewhat smaller redshifts than those probed by the Lyman-alpha forest, we estimate a total signal-to-noise of 9 for cross-correlating QSO spectra of SDSS-III with Planck and of 20 for cross-correlating with a future polarization based CMB experiment. The detection of this effect would be a direct measure of the neutral hydrogen-matter cross-correlation and could provide important information on the growth of structures at large scales in a redshift range which is still poorly probed by observations.

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 revisit Lyman-alpha bounds on the dark matter mass in Lambda Warm Dark Matter (Lambda-WDM) models, and derive new bounds in the case of mixed Cold plus Warm models (Lambda-CWDM), using a set up which is a good approximation for several theoretically well-motivated dark matter models. We combine WMAP5 results with two different Lyman-alpha data sets, including observations from the Sloan Digital Sky Survey. We pay a special attention to systematics, test various possible sources of error, and compare the results of different statistical approaches. Expressed in terms of the mass of a non-resonantly produced sterile neutrino, our bounds read m_NRP > 8 keV (frequentist 99.7% confidence limit) or m_NRP > 12.1 keV (Bayesian 95% credible interval) in the pure Lambda-WDM limit. For the mixed model, we obtain limits on the mass as a function of the warm dark matter fraction F_WDM. Within the mass range studied here (5 keV < m_NRP < infinity), we find that any mass value is allowed when F_WDM < 0.6 (frequentist 99.7% confidence limit); similarly, the Bayesian joint probability on (F_WDM, 1/m_NRP) allows any value of the mass at the 95% confidence level, provided that F_WDM < 0.35.

We investigate the properties of Damped Lyman-alpha systems (DLAs) using high-resolution and large box-size cosmological hydrodynamical simulations of a LambdaCDM model. The numerical code used is a modification of GADGET-2 with a self consistent implementation of the metal enrichment mechanism (Tornatore et al. 2007). We explore the numerical convergence of some relevant physical quantities and we vary the parameters describing the properties of galactic winds; the initial stellar mass function; the linear dark matter power spectrum and the metal enrichment pattern of the IGM (Intergalactic Medium) around DLAs. We find that strong galactic winds with speed of about 600 km/s, in an energy-driven wind scenario, are needed in order to match the observed column density distribution function for DLAs and the evolution of the neutral hydrogen content with redshift. The momentum-driven implementation of the galactic wind model, that relates the speed and mass load in the wind to the properties of the dark matter haloes, shows a behaviour which is intermediate between the energy-driven galactic winds of small (100 km/s) and large (600 km/s) velocities. At z=3 the contribution of haloes of masses between 10^9 and 10^10 h^-1 M_sun, for DLAs below 10^20.8 cm^-2, to the column density distribution function, is significant. By interpolating physical quantities along line-of-sights through massive haloes we qualitatively show how different galactic wind models impact on the IGM around DLAs. Furthermore, we analyse statistics related to the velocity widths of SiII associated to DLAs: the metallicity in the wind seems to be rather clumpy and this produces an underestimation of the observed velocity widths. We outline possible solutions to this problem.

Early dark energy models, for which the contribution to the dark energy density at high redshifts is not negligible, influence the growth of cosmic structures and could leave observable signatures that are different from the standard cosmological constant cold dark matter ($\Lambda$CDM) model. In this paper, we present updated constraints on early dark energy using geometrical and dynamical probes. From WMAP five-year data, baryon acoustic oscillations and type Ia supernovae luminosity distances, we obtain an upper limit of the dark energy density at the last scattering surface (lss), $\Omega_{\rm EDE}(z_{\rm lss})<2.3\times10^{-2}$ (95% C.L.). When we include higher redshift observational probes, such as measurements of the linear growth factors, Gamma-Ray Bursts (GRBs) and Lyman-$\alpha$ forest (\lya), this limit improves significantly and becomes $\Omega_{\rm EDE}(z_{\rm lss})<1.4\times10^{-3}$ (95% C.L.). Furthermore, we find that future measurements, based on the Alcock-Paczy\'nski test using the 21cm neutral hydrogen line, on GRBs and on the \lya forest, could constrain the behavior of the dark energy component and distinguish at a high confidence level between early dark energy models and pure $\Lambda$CDM. In this case, the constraints on the amount of early dark energy at the last scattering surface improve by a factor ten, when compared to present constraints. We also discuss the impact on the parameter $\gamma$, the growth rate index, which describes the growth of structures in standard and in modified gravity models.

We assess the possibility to detect and characterize the physical state of the missing baryons at low redshift by analyzing the X-ray absorption spectra of the Gamma Ray Burst [GRB] afterglows, measured by a micro calorimeters-based detector with 3 eV resolution and 1000 cm2 effective area and capable of fast re-pointing, similar to that on board of the recently proposed X-ray satellites EDGE and XENIA. For this purpose we have analyzed mock absorption spectra extracted from different hydrodynamical simulations used to model the properties of the Warm Hot Intergalactic Medium [WHIM]. These models predict the correct abundance of OVI absorption lines observed in UV and satisfy current X-ray constraints. According to these models space missions like EDGE and XENIA should be able to detect about 60 WHIM absorbers per year through the OVII line. About 45 % of these have at least two more detectable lines in addition to OVII that can be used to determine the density and the temperature of the gas. Systematic errors in the estimates of the gas density and temperature can be corrected for in a robust, largely model-independent fashion. The analysis of the GRB absorption spectra collected in three years would also allow to measure the cosmic mass density of the WHIM with about 15 % accuracy, although this estimate depends on the WHIM model. Our results suggest that GRBs represent a valid, if not preferable, alternative to Active Galactic Nuclei to study the WHIM in absorption. The analysis of the absorption spectra nicely complements the study of the WHIM in emission that the spectrometer proposed for EDGE and XENIA would be able to carry out thanks to its high sensitivity and large field of view.

Previous fits of sterile neutrino dark matter models to cosmological data assumed a peculiar production mechanism, which is not representative of the best-motivated particle physics models given current data on neutrino oscillations. These analyses ruled out sterile neutrino masses smaller than 8-10 keV. Here we focus on sterile neutrinos produced resonantly. We show that their cosmological signature can be approximated by that of mixed Cold plus Warm Dark Matter (CWDM). We use recent results on LambdaCWDM models to show that for each mass greater than or equal to 2 keV, there exists at least one model of sterile neutrino accounting for the totality of dark matter, and consistent with Lyman-alpha and other cosmological data. Resonant production occurs in the framework of the nuMSM (the extension of the Standard Model with three right-handed neutrinos). The models we checked to be allowed correspond to parameter values consistent with neutrino oscillation data, baryogenesis and all other dark matter bounds.

We investigate the effects of non-Gaussianity in the primordial density field on the reionization history. We rely on a semi-analytic method to describe the processes acting on the intergalactic medium (IGM), relating the distribution of the ionizing sources to that of dark matter haloes. Extending previous work in the literature, we consider models in which the primordial non-Gaussianity is measured by the dimensionless non-linearity parameter f_NL, using the constraints recently obtained from cosmic microwave background data. We predict the ionized fraction and the optical depth at different cosmological epochs assuming two different kinds of non-Gaussianity, characterized by a scale-independent and a scale-dependent f_NL and comparing the results to those for the standard Gaussian scenario. We find that a positive f_NL enhances the formation of high-mass haloes at early epochs, when reionization begins, and, as a consequence, the IGM ionized fraction can grow by a factor up to 5 with respect to the corresponding Gaussian model. The increase of the filling factor has a small impact on the reionization optical depth and is of order ~ 10 per cent if a scale-dependent non-Gaussianity is assumed. Our predictions for non-Gaussian models are in agreement with the latest WMAP results within the error bars, but a higher precision is required to constrain the scale dependence of non-Gaussianity.

In the study of the process of cosmic structure formation numerical simulations are crucial tools to interface observational data to theoretical models and to investigate issues otherwise unexplored. Enormous advances have been achieved in the last years thanks to the availability of sophisticated codes, now allowing to tackle the problem of cosmic structure formation and subsequent evolution by covering larger and larger dynamical ranges. Moreover, computational cosmology is the ideal interpretative framework for the overwhelming amount of new data from extragalactic surveys and from large sample of individual objects. The Workshop Novicosmo 2008 "The Impact of Simulations in Cosmology and Galaxy Formation' held in SISSA was aimed at providing the state-of-the-art on the latest numerical simulations in Cosmology and in Galaxy Formation. Particular emphasis was given to the implementation of new physical processes in simulation codes, to the comparison between different codes and numerical schemes and how to use best supercomputing facilities of the next generation. Finally, the impact on our knowledge on the Physics of the Universe brought by this new channel of investigation has also been focused. The Workshop was divided in three sections corresponding (roughly) to three main areas of study: Reionization and Intergalactic medium; Dark and Luminous matter in galaxies; Clusters of galaxies and Large scale Structures. This paper will provide i) a short resume' of the scientific results of the Workshop ii) the complete list of the talks and the instructions on how to retrieve the .pdf of the related (powerpoint) presentations iii) a brief presentation of the associated Exhibition "Space Art"

We present results from the first high-resolution hydrodynamical simulations of non-Gaussian cosmological models. We focus on the statistical properties of the transmitted Lyman-alpha flux in the high redshift intergalactic medium. Imprints of non-Gaussianity are present and are larger at high redshifts. Differences larger than 20 % at z>3 in the flux probability distribution function for high transmissivity regions (voids) are expected for values of the non linearity parameter f_NL=\pm 100 when compared to a standard LCDM cosmology with f_NL=0. We investigate also the one-dimensional flux bispectrum: at the largest scales (corresponding to tens of Mpc) we expect deviations in the flux bispectrum up to 20% at z~4 (for f_NL=\pm 100), significantly larger than deviations of ~ 3% in the flux power spectrum. We briefly discuss possible systematic errors that can contaminate the signal. Although challenging, a detection of non-Gaussianities in the interesting regime of scales and redshifts probed by the Lyman-alpha forest, could be possible with future data sets.

We analyse the evolution of the Intergalactic Medium (IGM) by means of an extended set of large box size hydrodynamical simulations which include pre-heating. We focus on the properties of the z~2 Lyman-alpha forest and on the population of clusters and groups of galaxies at z=0. We investigate the distribution of voids in the Lyman-alpha flux and the entropy-temperature relation of galaxy groups, comparing the simulation results to recent data from high-resolution quasar spectra and from X-ray observations. Pre-heating is included through a simple phenomenological prescription, in which at z=4 the entropy of all gas particles, whose overdensity exceeds a threshold value delta_h is increased to a minimum value K_fl. While the entropy level observed in the central regions of galaxy groups requires a fairly strong pre-heating, with K_fl>100 keV cm^2, the void statistics of the Lyman-alpha forest impose that this pre-heating should take place only in relatively high-density regions, in order not to destroy the cold filaments that give rise to the forest. We conclude that any injection of non-gravitational energy in the diffuse baryons should avoid low-density regions at high redshift and/or take place at relatively low redshift.

We revisit the proximity effect produced by QSOs at redshifts 2.1-3.3 applying the FLO approach (Saitta et al. 2008) to a sample of ~6300 Ly-alpha lines fitted in 21 high resolution, high signal-to-noise spectra. This new technique allows to recover the hydrogen density field from the HI column densities of the lines in the Ly-alpha forest, on the basis of simple assumptions on the physical state of the gas. To minimize the systematic uncertainties that could affect the density recovering in the QSO vicinity, we carefully determined the redshifts of the QSOs in our sample and modelled in detail their spectra to compute the corresponding ionising fluxes. The mean density field obtained from the observed spectra shows a significant over-density in the region within 4 proper Mpc from the QSO position, confirming that QSOs are hosted in high density peaks. The absolute value of rho/<rho> for the peak is uncertain by a factor of ~3, depending on the assumed QSO spectral slope and the minimum HI column density detectable in the spectra. We do not confirm the presence of a significant over-density extending to separations of ~15 proper Mpc from the QSO, claimed in previous works at redshifts <z>=2.5 and 3.8. Our best guess for the UV background ionisation rate based on the IGM mean density recovered by FLO is Gamma_UVB ~ 10^{-12} s^{-1}. However, values of Gamma_UVB ~ 3x10^{-12} s^{-1} could be viable if an inverted temperature-density relation with index alpha=-0.5 is adopted.

We use high-resolution hydrodynamical simulations to investigate the spatial correlation between weak (N(HI) <10^{15} cm^{-2}) Ly-alpha absorbers and gas-rich galaxies in the local universe. We confirm that Ly-alpha absorbers are preferentially expected near gas-rich galaxies and that the degree of correlation increases with the column density of the absorber. The real-space galaxy auto-correlation is stronger than the cross-correlation (correlation lengths r_{0,gg}=3.1 \pm 0.1 Mpc/h and r_{0,ag}=1.4 \pm 0.1 Mpc/h, respectively), in contrast with the recent results of Ryan-Weber (2006, RW06), and the auto-correlation of absorbers is very weak. These results are robust to the presence of strong galactic winds in the hydrodynamical simulations. In redshift-space a further mismatch arises since at small separations the distortion pattern of the simulated galaxy-absorber cross-correlation function is different from the one measured by RW06. However, when sampling the intergalactic medium along a limited number of lines-of-sight, as in the real data, uncertainties in the cross correlation estimates are large enough to account for these discrepancies. Our analysis suggests that the statistical significance of difference between the cross-correlation and auto-correlation signal in current datasets is ~ 1-sigma only.

We investigate the properties of one--dimensional flux ``voids'' (connected regions in the flux distribution above the mean flux level) by comparing hydrodynamical simulations of large cosmological volumes with a set of observed high--resolution spectra at z ~ 2. After addressing the effects of box size and resolution, we study how the void distribution changes when the most significant cosmological and astrophysical parameters are varied. We find that the void distribution in the flux is in excellent agreement with predictions of the standard LCDM cosmology, which also fits other flux statistics remarkably well. We then model the relation between flux voids and the corresponding one--dimensional gas density field along the line--of--sight and make a preliminary attempt to connect the one--dimensional properties of the gas density field to the three--dimensional dark matter distribution at the same redshift. This provides a framework that allows statistical interpretations of the void population at high redshift using observed quasar spectra, and eventually it will enable linking the void properties of the high--redshift universe with those at lower redshifts, which are better known.

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.

In this paper we investigate the effects that a dynamic dark energy component dominant in the universe at late epochs has on reionization. We follow the evolution of HII regions with the analytic approach of Furlanetto and Oh (2005) in two different universes for which we assume the Peebles and Ratra (2003) and Brax and Martin (2000) quintessence models and we compare our results to the LCDM scenario. We show that, for a fixed ionization efficiency, at the same cosmological epoch the topology of bubbles is dominated by high-mass objects and the characteristic size of the ionized regions is slightly smaller than in the LCDM model, especially at the latest stages of reionization, due to the higher recombination efficiency. As a consequence, the bubbles' `epoch of overlap' happens earlier than in LCDM. Finally, we show how the different evolution of the HII regions affects the transmission of the high-z QSO spectra, reducing the Lyman flux absorption at small optical depths.

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.

[Abridged] We present FLO (From Lines to Over-densities), a new technique to reconstruct the hydrogen density field for the Lya forest lines observed in high resolution QSO spectra. The method is based on the hypothesis that the Lya lines arise in the low to intermediate density intergalactic gas and that the Jeans length is the typical size of the Lya absorbers. The reliability of FLO is tested against mock spectra obtained from cosmological simulations. The recovering algorithm gives satisfactory results in the range from the mean density to over-densities of ~30 and reproduces correctly the correlation function of the density field and the 1D power spectrum on scales between ~20 and 60 comoving Mpc. A sample of Lya forests from 22 high resolution QSO spectra is analysed, covering the redshift range 1.7<z<3.5. For each line of sight, we fit Voigt profiles to the lines of the Lya forest, providing the largest, homogeneous sample of fitted Lya lines ever studied. The line number density evolution with redshift follows a power-law relation: dn/dz=(166 +/- 4) [(1+z)/3.5]^{(2.8 +/- 0.2)} (1 sigma errors). The two-point correlation function of lines shows a signal up to separations of ~2 comoving Mpc; weak lines (log N(HI)<13.8) also show a significant clustering but on smaller scales (r<1.5 comoving Mpc). We estimate with FLO the hydrogen density field toward the 22 observed lines of sight. The redshift distribution of the average densities computed for each QSO is consistent with the cosmic mean hydrogen density in the analysed redshift range. The two-point correlation function and the 1D power spectrum of the delta field are estimated. The correlation function shows clustering signal up to ~4 comoving Mpc.

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

Several recent studies suggest a correlation between dark matter halo mass and the shape of the density profile. We re-analyze simulations from Ricotti (2003) in which such a correlation was proposed. We use a standard analysis of the halo density profiles and compare the old simulations to new ones performed with Gadget2, including higher resolution runs. We confirm Ricotti's result that, at virialization, the central log slopes alpha, at 5%-10% of the virial radius are correlated with the halo mass and that the halo concentration is a universal constant. Our results do not contradict the majority of published papers: when using a split power law to fit the density profiles, due to the alpha-concentration degeneracy, the fits are consistent with halos having a universal shape with alpha=1 or 1.5 and concentrations that depend on the mass, in agreement with results published elsewhere. Recently, several groups have found no evidence for convergence of the inner halo profile to a constant power law. The choice of a split power law parameterization used in this letter is motivated by the need to compare our results to previous ones and is formally valid because we are not able to resolve regions where the slope of the fitting function reaches its asymptotic constant value. Using a non-parameterized technique, we also show that the density profiles of dwarf galaxies at z ~ 10 have a log slope shallower than 0.5 within 5% of the virial radius.

The incidence of strong MgII systems in gamma-ray burst (GRB) spectra is a few times higher than in quasar (QSO) spectra. We investigate several possible explanations for this effect, including: dust obscuration bias, clustering of the absorbers, different beam sizes of the sources, multiband magnification bias of GRBs, association of the absorbers with the GRB event or the circumburst environment. We find that: i) the incidence rate of MgII systems in QSO spectra could be underestimated by a factor 1.3-2 due to dust obscuration; ii) the equivalent-width distribution of the MgII absorbers along GRBs is consistent with that observed along QSOs thus suggesting that the absorbers are more extended than the beam sizes of the sources; iii) on average, GRB afterglows showing more than one MgII system are a factor of 1.7 brighter than the others, suggesting a lensing origin of the observed discrepancy; iv) gravitational lensing (in different forms, from galaxy lensing to microlensing) can bias high the counts of MgII systems along GRBs if the luminosity functions of the prompt gamma-ray emission and of the optical afterglows have a mean faint-end slope approaching -5/3 -- -2; v) some of the absorbers can be associated with the circumburst environment or produced by supernova remnants unrelated to the GRB event itself but lying in the same star-forming region. With the possible exception of magnification bias, it is unlikely that one of these effects on its own can fully account for the observed counts. However, the combined action of some of them can substantially reduce the statistical significance of the discrepancy.

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

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.

Using analytic arguments and a suite of very high resolution (10^3 Msun per particle) cosmological hydro-dynamical simulations, we argue that high redshift, z ~ 10, M ~ 10^8 Msun halos, form the smallest `baryonic building block' (BBB) for galaxy formation. These halos are just massive enough to efficiently form stars through atomic line cooling and to hold onto their gas in the presence of supernovae winds and reionisation. These combined effects, in particular that of the supernovae feedback, create a sharp transition: over the mass range 3-10x10^7 Msun, the BBBs drop two orders ofmagnitude in stellar mass. Below ~2x10^7 Msun, galaxies will be dark with almost no stars and no gas. Above this scale is the smallest unit of galaxy formation: the BBB. A small fraction (~100) of these gas rich BBBs fall in to a galaxy the size of the Milky Way. Ten percent of these survive to become the observed LG dwarf galaxies at the present epoch. Those in-falling halos on benign orbits which keep them far away from the Milky Way or Andromeda manage to retain their gas and slowly form stars - these become the smallest dwarf irregular galax ies; those on more severe orbits lose their gas faster than they can form stars and become the dwarf spheroidals. The remaining 90% of the BBBs will be accreted. We show that this gives a metallicity and total stellar mass consistent with the Milky Way old stellar halo (abridged).

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.

Various recent studies proved that cosmological models with a significant contribution from cold dark matter isocurvature perturbations are still compatible with most recent data on cosmic microwave background anisotropies and on the shape of the galaxy power spectrum, provided that one allows for a very blue spectrum of primordial entropy fluctuations (n_iso > 2). However, such models predict an excess of matter fluctuations on small scales, typically below 40 Mpc/h. We show that the proper inclusion of high-resolution high signal-to-noise Lyman-alpha forest data excludes most of these models. The upper bound on the isocurvature fraction alpha=f_iso^2/(1+f_iso^2), defined at the pivot scale k_0=0.05/Mpc, is pushed down to alpha<0.4, while n_iso=1.9+-1.0 (95% confidence limits). We also study the bounds on curvaton models characterized by maximal correlation between curvature and isocurvature modes, and a unique spectral tilt for both. We find that f_iso<0.05 (95% c.l.) in that case. For double inflation models with two massive inflatons coupled only gravitationally, the mass ratio should obey R < 3 (95% c.l.).

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

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.

We use LUQAS, a sample of 27 high resolution high signal-to-noise UVES quasar (QSO) spectra (Kim et al. 2004), and the Croft et al. (2002) sample together with a set of high resolution large box size hydro-dynamical simulations run with the code GADGET-II to recover the linear dark matter power spectrum at z>2 and at scales of 1-40 Mpc/h. These scales cannot be probed by any other observable. We address some of the uncertainties in the theoretical modelling of the Lyman-alpha forest structures such as the difference between full hydro-dynamical simulations and a more simplified scheme based on the HPM (Hydro-Particle-Mesh) technique. We combine these data sets with the Cosmic Microwave Background (CMB) data of WMAP in order to get tighter constraints on cosmological parameters. We focus on the recovered values of the power spectrum amplitude, the primordial spectral index and the running of the primordial spectral index. By considering models of slow-roll inflation we give constraints on inflationary parameters. We explore implications for the mass of a warm dark matter particle, including light gravitinos and sterile neutrinos, and the perspectives of constraining the neutrino mass.

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.

We present a simple method for tracing the spatial distribution and predicting the physical properties of the Warm-Hot Intergalactic Medium (WHIM), from the map of galaxy light in the local universe. Under the assumption that biasing is local and monotonic we map the ~ 2 Mpc/h smoothed density field of galaxy light into the mass density field from which we infer the spatial distribution of the WHIM in the local supercluster. Taking into account the scatter in the WHIM density-temperature and density-metallicity relation, extracted from the z=0 outputs of high-resolution and large box size hydro-dynamical cosmological simulations, we are able to quantify the probability of detecting WHIM signatures in the form of absorption features in the X-ray spectra, along arbitrary directions in the sky. To illustrate the usefulness of this semi-analytical method we focus on the WHIM properties in the Virgo Cluster region.

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

In the WMAP era of high precision cosmology an accurate determination of the matter power spectrum from lyman-alpha forest data becomes crucial. When combining the matter power spectrum derived from CMB experiments with that inferred from lyman-alpha absorption an evidence for a running spectral index and a primordial index n < 1 arises (Verde et al. 2003). In this talk, I will describe some results obtained from a sample of 27 high resolution and high signal-to-noise quasar spectra (the LUQAS sample) and I will address possible systematic errors that can affect the estimate of the flux power spectrum.

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.

In quintessence models, the dark energy content of the universe is described by a slowly rolling scalar field whose pressure and energy density obey an equation of state of the form p=w $\rho$; w is in general a function of time such that w<-1/3, in order to drive the observed acceleration of the Universe today. The cosmological constant model LCDM corresponds to the limiting case w=-1. In this paper, we explore the prospects of using the Lyman-alpha forest to constrain w, using semi-analytical techniques to model the intergalactic medium (IGM). A different value of w changes both the growth factor and the Hubble parameter as a function of time. The resulting change in the optical depth distribution affects the optical depth power spectrum, the number of regions of high transmission per unit redshift and the cross-correlation coefficient of spectra of quasar pairs. These can be detected in current data, provided we have independent estimates of the thermal state of the IGM, its ionization parameter and the baryon density.

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.

We have performed hydrodynamical simulations to investigate the effects of galactic winds on the high-redshift (z=3) universe. Strong winds suppress the formation of low-mass galaxies significantly, and the metals carried by them produce C IV absorption lines with properties in reasonable agreement with observations. The winds have little effect on the statistics of the H I absorption lines, because the hot gas bubbles blown by the winds fill only a small fraction of the volume and because they tend to escape into the voids, thereby leaving the filaments that produce these lines intact.

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.

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