BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matteo Viel

Matteo Viel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAgr-qchep-phastro-ph.IMastro-ph.HEhep-thastro-ph

Catalog footprint

What is connected

99works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

99 published item(s)

preprint2025arXiv

Detecting the IGM metal enrichment with the 2-point correlation function of the flux. Application to the UVES deep spectrum

The distribution and the abundance of metals in the intergalactic medium (IGM) have strong implications on galaxy formation and evolution models. The ionic transitions of heavy elements in quasar spectra can probe both the mechanisms and the sources of chemical pollution, but high resolution, high signal-to-noise ratio (S/N) spectra are required, as the IGM absorbers can be too weak for direct detection. In this work, we investigate the IGM metallicity, focusing on the detection of the weak absorption lines. We exploited the cosmological tool of the two-point correlation function (TPCF) and applied it to the transmitted flux in the CIV forest of the ultra-high S/N UVES spectrum of the quasar HE0940-1050 (z~3). We also "deabsorbed" the strongest circum-galactic medium (CGM) systems to reveal the underlying IGM signal. For each test, we generated 1000 mock spectra with shuffled line positions to derive an estimate for the TPCF noise level. The TPCF shows a clear peak at the characteristic velocity separation of the CIV doublet, but when deabsorbing the CGM contribution (i.e. all metal lines and CIV lines associated with logN(HI)>14.0), the peak is no longer significant at 1$σ$, although 7 weak CIV systems remain visible. Even adding up to 135 additional weak mock CIV systems (logN(HI)<14.0) to the spectrum does not produce a significant CIV peak. Eventually, when we create a synthetic spectrum with gaussian distributed noise and same S/N as the complete spectrum, we remove the signal caused by the spectral intrinsic features and thus find a peak compatible with a metallicity of -3.80<[C/H]<-3.50. We conclude that the TPCF method is not sensitive to the presence of the weakest systems in the real spectrum, despite the extremely high S/N and high resolution of the data. However, the results of this statistical technique would possibly change when combining more than one line of sight.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cosmic voids and BAO with relative baryon-CDM perturbations

We study the statistics of various large-scale structure tracers in gravity-only cosmological simulations including baryons and cold dark matter (CDM) initialized with two different transfer functions, and simulated as two distinct fluids. This allows us to study the impact of baryon-CDM relative perturbations on these statistics. In particular, we focus on the statistics of cosmic voids, as well as on the matter and halo real-space 2-point correlation function and baryon acoustic oscillations (BAO) peak. We find that the void size function is affected at the 1-2% level at maximum, and that the impact is more important at higher redshift, while the void density profile and void bias are roughly unaffected. We do not detect a sizeable impact of relative baryon-CDM perturbations on the real-space correlation functions of matter and halos or the BAO peak, which is in line with results from previous works. Our results imply that it would be hard to use voids or real-space correlation functions to constrain baryon-CDM relative perturbations, but also that we might not have to include them in models for the analysis of future cosmological surveys data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational waves $\times$ HI intensity mapping: cosmological and astrophysical applications

Two of the most rapidly growing observables in cosmology and astrophysics are gravitational waves (GW) and the neutral hydrogen (HI) distribution. In this work, we investigate the cross-correlation between resolved gravitational wave detections and HI signal from intensity mapping (IM) experiments. By using a tomographic approach with angular power spectra, including all projection effects, we explore possible applications of the combination of the Einstein Telescope and the SKAO intensity mapping surveys. We focus on three main topics: \textit{(i)} statistical inference of the observed redshift distribution of GWs; \textit{(ii)} constraints on dynamical dark energy models as an example of cosmological studies; \textit{(iii)} determination of the nature of the progenitors of merging binary black holes, distinguishing between primordial and astrophysical origin. Our results show that: \textit{(i)} the GW redshift distribution can be calibrated with good accuracy at low redshifts, without any assumptions on cosmology or astrophysics, potentially providing a way to probe astrophysical and cosmological models; \textit{(ii)} the constrains on the dynamical dark energy parameters are competitive with IM-only experiments, in a complementary way and potentially with less systematics; \textit{(iii)} it will be possible to detect a relatively small abundance of primordial black holes within the gravitational waves from resolved mergers. Our results extend towards $\mathrm{GW \times IM}$ the promising field of multi-tracing cosmology and astrophysics, which has the major advantage of allowing scientific investigations in ways that would not be possible by looking at single observables separately.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Limits on non-canonical heating and turbulence in the intergalactic medium from the low redshift Lyman-alpha forest

We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5σ$) with the shape and amplitude of the observed CDDF at HI column densities $10^{13.3}\rm\,cm^{-2}\leq N_{\rm HI}\leq 10^{14.5}\rm\,cm^{-2}$. However, the Doppler widths of the simulated lines remain too narrow with respect to the COS data. We argue that invoking AGN feedback does not resolve this discrepancy. We also disfavour enhanced photoheating rates as a potential solution, as this requires an unphysically hard UV background spectrum. If instead appealing to a non-canonical source of heating, an additional specific heat injection of $u \lesssim 6.9\rm\,eV\,m_{\rm p}^{-1}$ is required at $z\lesssim 2.5$ for gas that has $N_{\rm HI}\simeq 10^{13.5}\rm\,cm^{-2}$ by $z=0.1$. Alternatively, there may be an unresolved line of sight turbulent velocity component of $v_{\rm turb}\lesssim 8.5\rm\,km\,s^{-1}(N_{\rm HI}/10^{13.5}\rm\,cm^{-2})^{0.21}$ for the coldest gas in the diffuse IGM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effect of inhomogeneous reionisation on the Lyman-$α$ forest power spectrum at redshift $z>4$: implications for thermal parameter recovery

We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k<0.03\rm\,km^{-1}\,s$, the power spectrum is instead enhanced by $10$--$50$ per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Ly$α$ forest power spectrum at $4.1\leq z \leq 5.4$ that can be easily applied within other parameter inference frameworks using similar reionisation models. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Ly$α$ forest data at $z>4$, $0.006 \rm \,km^{-1}\,s\leq k \leq 0.2 \rm\, km^{-1}\,s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $\sim 10$ per cent. However, for $5$ per cent uncertainties, $\sim 1σ$ shifts between the estimated and true parameters occur.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Weighing Cosmic Structures with Clusters of Galaxies and the Intergalactic Medium

We present an analysis aimed at combining cosmological constraints from number counts of galaxy clusters identified through the Sunyaev-Zeldovich effect, obtained with the South Pole Telescope (SPT), and from Lyman-$α$ spectra obtained with the MIKE/HIRES and X-shooter spectrographs. The SPT cluster analysis relies on mass calibration based on weak lensing measurements, while the Lyman-$α$ analysis is built over a suite of hydrodynamical simulations for the extraction of mock spectra. The resulting constraints exhibit a tension ($\sim 3.3σ$) between the low $σ_8$ values preferred by the low-redshift cluster data, $σ_8=0.74 ^{+0.03}_{-0.04}$, and the higher one preferred by the high-redshift Lyman-$α$ data, $σ_8=0.91 ^{+0.03}_{-0.03}$. We present a detailed analysis in order to understand the origin of this tension and, in particular, to establish whether it arises from systematic uncertainties related to the assumptions underlying the analyses of cluster counts and/or Lyman-$α$ forest. We found this tension to be robust with respect to the choice of modeling of the IGM, even when including possible systematics from unaccounted sub-Damped Lyman-$α$ (DLA) and Lyman-limit systems (LLS) in the Lyman-$α$ data. We conclude that to solve this tension from the SPT side would require a large bias on the cluster mass estimate, or from the Lyman-$α$ side large unaccounted errors on the Lyman-$α$ mean fluxes, respectively. Our results have important implications for future analyses based on cluster number counts from future large photometric surveys (e.g. Euclid and LSST) and on larger samples of high-redshift quasar spectra (e.g. DESI and WEAVE surveys). If confirmed at the much higher statistical significance reachable by such surveys, this tension could represent a significant challenge for the standard $Λ$CDM paradigm.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Constraining beyond $Λ$CDM models with 21cm intensity mapping forecast observations combined with latest CMB data

We explore constraints on dark energy and modified gravity with forecast 21cm intensity mapping measurements using the Effective Field Theory approach. We construct a realistic mock data set forecasting a low redshift 21cm signal power spectrum $P_{21}(z,k)$ measurement from the MeerKAT radio-telescope. We compute constraints on cosmological and model parameters through Monte Carlo Markov chain techniques, testing both the constraining power of $P_{21}(k)$ alone and its effect when combined with the latest Planck 2018 CMB data. We complement our analysis by testing the effects of tomography from an ideal mock data set of observations in multiple redshift bins. We conduct our analysis numerically with the codes EFTCAMB/EFTCosmoMC, which we extend by implementing a likelihood module fully integrated with original codes. We find that adding $P_{21}(k)$ to CMB data provides significantly tighter constraints on $Ω_ch^2$ and $H_0$, with a reduction of the error with respect to Planck results at the level of more than $60\%$. For the parameters describing beyond $Λ$CDM theories, we observe a reduction in the error with respect to the Planck constraints at the level of $\lesssim 10\%$. The improvement increases up to $\sim 35\%$ when we constrain the parameters using ideal, tomographic mock observations. We conclude that the power spectrum of the 21cm signal is sensitive to variations of the parameters describing the examined beyond $Λ$CDM models and, thus, $P_{21}(k)$ observations could help to constrain dark energy. The constraining power on such theories is improved significantly by tomography.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$H_0$ Reconstruction with Type Ia Supernovae, Baryon Acoustic Oscillation and Gravitational Lensing Time-Delay

There is a persistent $H_0$-tension, now at more than $\gtrsim 4σ$ level, between the local distance ladder value and the \emph{Planck} cosmic microwave background measurement, in the context of flat $Λ$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 $Λ$CDM value. When including the GLTD data, $H_0$ increases mildly, yet remaining discrepant with the local measurement at $\sim 2.5σ$ 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 $σ_{H_0} \approx 0.56$ uncertainty level.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraints on the Spacetime Dynamics of an Early Dark Energy Component

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-$α$ forest BAO have instead a negligible impact. We find $Ω_{\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 $Ω_{\rm eDE} < 0.00337$ if perturbations are neglected. The constraints on the EDE parameters are remarkably stable even when $Σm_ν$, and $r$ parameters are varied, with weak degeneracies between $Ω_{\rm eDE}$ and $r$ or $Σm_ν$. Perturbation parameters are not constrained with current data sets, and tensions between the CMB derived $H_0$ and $σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of dark matter models on the EoR 21-cm signal bispectrum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Late-time decaying dark matter: constraints and implications for the $H_0$-tension

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 $τ> 9\, \rm{Gyr}$ and $τ>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 $Λ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-linear damping of superimposed primordial oscillations on the matter power spectrum in galaxy surveys

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The impact of relativistic effects on the 3D Quasar-Lyman-$α$ cross-correlation

We study the impact of relativistic effects in the 3-dimensional cross-correlation between Lyman-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The atomic Hydrogen content of the post-reionization Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Weighing neutrinos with the halo environment

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_ν$, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dark matter searches in the gamma-ray extragalactic background via cross-correlations with galaxy catalogues

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lyman-alpha Forests cool Warm Dark Matter

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Relativistic effects in Lyman-alpha forest

We present the calculation of the Lyman-alpha (Lyman-$α$) 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-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Sherwood simulation suite: overview and data comparisons with the Lyman-alpha forest at redshifts 2 < z < 5

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraint on neutrino masses from SDSS-III/BOSS Ly$α$ forest and other cosmological probes

We present constraints on the parameters of the $Λ$CDM cosmological model in the presence of massive neutrinos, using the one-dimensional Ly$α$ 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$α$ spectrum is done using a second-order Taylor expansion of the simulated power spectrum. BOSS Ly$α$ data alone provide better bounds than previous Ly$α$ results, but are still poorly constraining, especially for the sum of neutrino masses $\sum m_ν$, for which we obtain an upper bound of 1.1~eV (95\% CL), including systematics for both data and simulations. Ly$α$ constraints on $Λ$CDM parameters and neutrino masses are compatible with CMB bounds from the Planck collaboration. Interestingly, the combination of Ly$α$ 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_ν< 0.15$~eV at 95\% CL (with a best-fit in zero). Adding recent BAO results further tightens this constraint to $\sum m_ν< 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$α$. Given the measured values of the two squared mass differences $Δm^2$, this result tends to favor the normal hierarchy scenario against the inverted hierarchy scenario for the masses of the active neutrino species.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmological implications of baryon acoustic oscillation (BAO) measurements

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 $Ω_m=0.301 \pm 0.008$ and curvature $Ω_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 $χ^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)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmology with a SKA HI intensity mapping survey

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cross-correlating 21cm intensity maps with Lyman Break Galaxies in the post-reionization era

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large-scale clustering of Lyman-alpha emission intensity from SDSS/BOSS

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measuring redshift-space distortions with future SKA surveys

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutral hydrogen in galaxy clusters: impact of AGN feedback and implications for intensity mapping

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 $Ω_{\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 $Ω_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Particle dark matter searches outside the Local Group

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The contribution of light Majorana neutrinos to neutrinoless double beta decay and cosmology

Cosmology is making impressive progress and it is producing stringent bounds on the sum of the neutrino masses Σ, 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νββ) search. This analysis indicates small values for the lightest neutrino mass, since the authors find Σ < 84 meV at 1σ C. L., and provides a 1σ preference for the normal hierarchy. The allowed values for the Majorana effective mass, probed by 0νββ, turn out to be < 75meV at 3σC.L. and lower down to less than 20meV at 1σC.L.. If this indication is confirmed, the impact on the 0νββ experiments will be tremendous since the possibility of detecting a signal will be out of the reach of the next generation of experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of massive neutrinos on the BAO peak

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_ν= 0.15$ eV and increases by $\sim1.2 \%$ for $ \sum m_ν= 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The halo model in a massive neutrino cosmology

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Non-Linear Power Spectrum of the Lyman Alpha Forest

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tomography of the Fermi-LAT gamma-ray diffuse extragalactic signal via cross-correlations with galaxy catalogs

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)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Voids in massive neutrino cosmologies

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 $σ_8$ and $Ω_ν$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Warm dark matter signatures on the 21cm power spectrum: Intensity mapping forecasts for SKA

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weighing neutrinos with cosmic neutral hydrogen

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, $Ω_{\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_ν\equiv Σm_ν$. 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 $σ(M_ν)\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 $σ(M_ν)\simeq0.06$ eV (95% CL).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A coarse grained perturbation theory for the Large Scale Structure, with cosmology and time independence in the UV

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic Degeneracies I: Joint N-body Simulations of Modified Gravity and Massive Neutrinos

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 $Λ$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 $Λ$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 $Σ_{i}m_{ν_{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].

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galactic outflow and diffuse gas properties at z>=1 using different baryonic feedback models

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 ($η$) is between $0.2 - 10$. The comparison Effective model generates a constant outflow velocity, and a negative correlation of $η$ 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 $η$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function

The Ly$α$ 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, $γ$, and temperature at mean-density, $T_0$, where $T(Δ) = T_0 Δ^{γ-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $β_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$α$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $γ=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 ($γ\leq 1$) are disfavored at a significance of over 4$σ$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Massive Neutrinos and the Non-linear Matter Power Spectrum

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling the neutral hydrogen distribution in the post-reionization Universe: intensity mapping

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quasar-Lyman $α$ Forest Cross-Correlation from BOSS DR11 : Baryon Acoustic Oscillations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Semi-Analytic Galaxy Formation in Massive Neutrino Cosmologies

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 $Λ$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 $σ_8$ predicted for cosmologies including a massive neutrino background. We conclude that independent estimates of $σ_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 $Λ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulated star formation rate functions at z~4-7, and the role of feedback in high-z galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Suite of Hydrodynamical Simulations for the Lyman-Alpha Forest with Massive Neutrinos

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: weighing the neutrino mass using the galaxy power spectrum of the CMASS sample

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 $Σm_ν<0.340$ eV within a flat $Λ$CDM background, and $Σm_ν<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 $Ω_M=0.2796\pm0.0097$, $H_0=69.72^{+0.90}_{-0.91}$ km/s/Mpc (flat $Λ$CDM) and $Ω_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 Ω_K and the running of the spectral index $α_s$ are both consistent with zero. The tensor-to-scaler ratio is constrained down to $r<0.198$ (95 percent CL, flat $Λ$ CDM) and $r<0.440$ (95 percent CL, more general background model).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The First Billion Years of a Warm Dark Matter Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Lyman-beta forest as a cosmic thermometer

We present a comprehensive analysis of high resolution hydrodynamic simulations in terms of Lyman-alpha and Lyman-beta one dimensional flux power spectra ($P_{αα}$ and $P_{ββ}$). In particular, we focus on the behaviour that the flux auto-power spectra and cross-power spectra ($P_{αβ}$) 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 $γ$ of the IGM temperature-density relation ($T=T_0(1+δ)^{γ-1}$). By performing a Fisher matrix analysis, we find that the power spectrum $P_{ββ}$ is more sensitive to the thermal history than $P_{αα}$ alone, due to the fact that it probes denser regions than Lyman-alpha. When we combine the power and cross spectra the constraints on $γ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The SDSS-III Baryonic Oscillation Spectroscopic Survey: Constraints on the Integrated Sachs Wolfe effect

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 $Λ$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$\,σ$. While this result is in close agreement with predictions from Monte Carlo simulations in the concordance $Λ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A consistent determination of the temperature of the intergalactic medium at redshift z=2.4

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$α$ forest. We use a large set of high-resolution hydrodynamical simulations to calibrate the relationship between the temperature-density ($T$--$Δ$) relation in the IGM and the distribution of HI column densities, $N_{\rm HI}$, and velocity widths, $b_{\rm HI}$, of discrete Ly$α$ 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$--$Δ$ relation, $T=T_{0}Δ^{γ-1}$, with a temperature at mean density, $T_{0}=[1.00^{+0.32}_{-0.21}]\times10^{4}\rm\,K$ and slope $(γ-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$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Baryon Acoustic Oscillations in the Ly-α forest of BOSS quasars

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmology with massive neutrinos I: towards a realistic modeling of the relation between matter, haloes and galaxies

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 $Ω_ν$-$σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmology with massive neutrinos III: the halo mass function and an application to galaxy clusters

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_ν= 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 $σ_8(Ω_{\rm m}/0.27)^γ$ relation which can be quantified as $Δγ\sim 0.05$ and $Δγ\sim 0.14$ assuming one ($N_ν=1$) or three ($N_ν=3$) degenerate massive neutrino, respectively. The shift results in a lower mean value of $σ_8$ with $Δσ_8 = 0.01$ for $N_ν=1$ and $Δσ_8 = 0.02$ for $N_ν=3$, respectively. Such difference, in a cosmology with massive neutrinos, would increase the tension between cluster abundance and Planck CMB measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fitting Methods for Baryon Acoustic Oscillations in the Lyman-α Forest Fluctuations in BOSS Data Release 9

We describe fitting methods developed to analyze fluctuations in the Lyman-α 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galactic Winds in Cosmological Simulations of the Circumgalactic Medium

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic or thermal AGN feedback in simulations of isolated and merging disk galaxies calibrated by the M-sigma relation

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lyman-α Forest and Cosmic Weak Lensing in a Warm Dark Matter 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α forest and exclude all masses below 3.3 keV at the 2σ 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α forest data but explore a different redshift regime.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Measurement of Baryon Acoustic Oscillations in the Lyman-alpha Forest Fluctuations in BOSS Data Release 9

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutrino Signatures on the High Transmission Regions of the Lyman-alpha Forest

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-linear evolution of the cosmic neutrino background

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The one-dimensional Ly-alpha forest power spectrum from BOSS

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$α$ 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 $Λ$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $σ_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$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Sloan Digital Sky Survey quasar catalog: tenth data release

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Turbulence driven by structure formation in the circum-galactic medium

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 δ 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 δ^{0.58}.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Clustering of Sloan Digital Sky Survey III Photometric Luminous Galaxies: The Measurement, Systematics and Cosmological Implications

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 Ω_Λ= 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 $Ω_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)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on Massive Neutrinos from the CFHTLS Angular Power Spectrum

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_ν}$ 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_ν}<0.29\,$eV and $N_{\rm{eff}} =4.17^{+1.62}_{-1.26}$ (2$\,σ$ confidence levels). If we omit smaller scales which may be affected by non-linearities, these constraints become $\sum{m_ν}<0.41\,$eV and $N_{\rm{eff}} =3.98^{+2.02}_{-1.20}$ (2$\,σ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmography beyond standard candles and rulers

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 $Λ$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 $Λ$CDM at 3-$σ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutrino masses and cosmological parameters from a Euclid-like survey: Markov Chain Monte Carlo forecasts including theoretical errors

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The impact of feedback from galaxy formation on the Lyman-alpha transmitted flux

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Sloan Digital Sky Survey quasar catalog: ninth data release

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 $Λ$CDM cosmology with H0 = 70 km/s/Mpc, $Ω_{\rm M}$ = 0.3, and $Ω_Λ$ = 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A cross-correlation study of the Fermi-LAT $γ$-ray diffuse extragalactic signal

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)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Coarse-Grained Cosmological Perturbation Theory

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraints on Primordial Non-Gaussianity from Large Scale Structure Probes

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}ϕ^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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effects of Massive Neutrinos on the Large-Scale Structure of the Universe

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Galactic Winds and Extended Lya Emission from the Host Galaxies of High Column Density QSO Absorption Systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Minimally Parametric Power Spectrum Reconstruction from the Lyman-alpha Forest

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Studying the Warm-Hot Intergalactic Medium in Emission

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Non-Linear Matter Power Spectrum in Warm Dark Matter Cosmologies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Statistics of Cosmological Lyman-alpha Absorption

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

X-shooter observations of QSO pairs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A decline and fall in the future of Italian Astronomy?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining Primordial Non-Gaussianity with High-Redshift Probes

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\,σ$) and $f_{\rm NL}=58\pm24$ ($1\,σ$) from NVSS data and SDSS DR6 QSO data, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High-Redshift Cosmography

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$σ$ confidence level. We present also forecasts for realistic data sets that are likely to be obtained in the next few years.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Primordial Non-Gaussianity and the NRAO VLA Sky Survey

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\,σ$ error bar) and $25<f_{\rm NL}<117$ ($2\,σ$ confidence level), corresponding to a detection of non-Gaussianity significant at the $\sim 3\,σ$ 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\,σ$) strikingly close to that found for optically selected quasars. We discuss caveats and possible physical and systematic effects that can impact on the results.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of neutrinos on the matter distribution as probed by the Intergalactic Medium

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_ν = 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_ν = 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_ν < 0.9$ eV (2 sigma C.L.), comparable to constraints obtained from the cosmic microwave background data or other large scale structure probes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The rise of the C IV mass density at z<2.5

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cross-correlations of the Lyman-alpha forest with weak lensing convergence I: Analytical Estimates of S/N and Implications for Neutrino Mass and Dark Energy

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 $σ_8$, the amplitude of these cross-correlations can potentially be used to measure the amplitude of $σ_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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Damped Lyman-alpha systems in high-resolution hydrodynamical simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lyman-alpha constraints on warm and on warm-plus-cold dark matter models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The impact of feedback on the low redshift Intergalactic Medium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Strong MgII systems in quasar and gamma-ray burst spectra

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.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

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

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

0 citations0 reviewsNo code linkedOpen access