Modeling the QSO luminosity and spatial clustering at low redshifts

We investigate the ability of hierarchical models of QSO formation and evolution to match the observed luminosity, number counts and spatial clustering of quasars at redshift z<2. These models assume that the QSO emission is triggered by galaxy mergers, that the mass of the central black hole correlates with halo properties and that quasars shine at their Eddington luminosity except, perhaps, during the very early stages of evolution. We find that models based on simple analytic approximations successfully reproduce the observed B-band QSO luminosity function at all redshifts, provided that some mechanisms is advocated to quench mass accretion within haloes larger than about 1e13 Msun that host bright quasars. These models also match the observed strength of QSO clustering at z~0.8. At larger redshifts, however, they underpredict the QSO biasing which, instead, is correctly reproduced by semi-analytic models in which the halo merger history and associated BHs are followed by Monte Carlo realizations of the merger hierarchy. We show that the disagreement between the luminosity function predicted by semi-analytic models and observations can be ascribed to the use of B-band data, which are a biased tracer of the quasar population, due to obscuration.Comment: 13 pages, 9 figures. Accepted by MNRA

The cosmological co-evolution of supermassive black holes, AGN and galaxies

We model the cosmological co-evolution of galaxies and their central supermassive black holes (BHs) within a semi-analytical framework developed on the outputs of the Millennium Simulation (Croton et al., 2006; De Lucia & Blaizot, 2007). In this work, we analyze the model BH scaling relations, fundamental plane and mass function, and compare them with the most recent observational data. Furthermore, we extend the original code developed by Croton et al. (2006) to follow the evolution of the BH mass accretion and its conversion into radiation, and compare the derived AGN bolometric luminosity function with the observed one. We find, for the most part, a very good agreement between predicted and observed BH properties. Moreover, the model is in good agreement with the observed AGN number density in 0<z<5, provided it is assumed that the cold gas fraction accreted by BHs at high redshifts is larger than at low redshifts (Marulli et al., 2008).Comment: Proceedings of "The Central Kiloparsec: Active Galactic Nuclei and Their Hosts", Ierapetra, Crete, 4-6 June, 2008. To appear in Volume 79 of the Memorie della Societa' Astronomica Italiana. 5 pages, 4 figure

A numerical study of the effects of primordial non-Gaussianities on weak lensing statistics

While usually cosmological initial conditions are assumed to be Gaussian, inflationary theories can predict a certain amount of primordial non-Gaussianity which can have an impact on the statistical properties of the lensing observables. In order to evaluate this effect, we build a large set of realistic maps of different lensing quantities starting from light-cones extracted from large dark-matter only N-body simulations with initial conditions corresponding to different levels of primordial local non-Gaussianity strength $f_{\rm NL}$. Considering various statistical quantities (PDF, power spectrum, shear in aperture, skewness and bispectrum) we find that the effect produced by the presence of primordial non-Gaussianity is relatively small, being of the order of few per cent for values of $|f_{\rm NL}|$ compatible with the present CMB constraints and reaching at most 10-15 per cent for the most extreme cases with $|f_{\rm NL}|=1000$. We also discuss the degeneracy of this effect with the uncertainties due to the power spectrum normalization $\sigma_8$ and matter density parameter $\Omega_{\rm m}$, finding that an error in the determination of $\sigma_8$ ($\Omega_{\rm m}$) of about 3 (10) per cent gives differences comparable with non-Gaussian models having $f_{\rm NL}=\pm 1000$. These results suggest that the possible presence of an amount of primordial non-Gaussianity corresponding to $|f_{\rm NL}|=100$ is not hampering a robust determination of the main cosmological parameters in present and future weak lensing surveys, while a positive detection of deviations from the Gaussian hypothesis is possible only breaking the degeneracy with other cosmological parameters and using data from deep surveys covering a large fraction of the sky.Comment: accepted by MNRA

Searching for galaxy clusters in the Kilo-Degree Survey

In this paper, we present the tools used to search for galaxy clusters in the Kilo Degree Survey (KiDS), and our first results. The cluster detection is based on an implementation of the optimal filtering technique that enables us to identify clusters as over-densities in the distribution of galaxies using their positions on the sky, magnitudes, and photometric redshifts. The contamination and completeness of the cluster catalog are derived using mock catalogs based on the data themselves. The optimal signal to noise threshold for the cluster detection is obtained by randomizing the galaxy positions and selecting the value that produces a contamination of less than 20%. Starting from a subset of clusters detected with high significance at low redshifts, we shift them to higher redshifts to estimate the completeness as a function of redshift: the average completeness is ~ 85%. An estimate of the mass of the clusters is derived using the richness as a proxy. We obtained 1858 candidate clusters with redshift 0 < z_c < 0.7 and mass 13.5 < log(M500/Msun) < 15 in an area of 114 sq. degrees (KiDS ESO-DR2). A comparison with publicly available Sloan Digital Sky Survey (SDSS)-based cluster catalogs shows that we match more than 50% of the clusters (77% in the case of the redMaPPer catalog). We also cross-matched our cluster catalog with the Abell clusters, and clusters found by XMM and in the Planck-SZ survey; however, only a small number of them lie inside the KiDS area currently available.Comment: 13 pages, 15 figures. Accepted for publication on Astronomy & Astrophysic

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 Lambda cold dark matter (Lambda CDM) cosmological scenario - namely an f (R) theory of modified gravity 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 non-linear matter power spectrum, the halo mass function, and the halo bias. 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. In particular, when an f (R) gravity model with f(R0) = -1 x 10(-4) is combined with a total neutrino mass of Sigma(i)m(nu i) = 0.4 eV, the resulting matter power spectrum, halo mass function, and bias at z = 0 are found to be consistent with the standard model's predictions at the less than or similar to 10, less than or similar to 20, and less than or similar to 5 per cent 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

Studying the WHIM with Gamma Ray Bursts

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

Gravitational Clustering from Chi^2 Initial Conditions

We consider gravitational clustering from primoridal non-Gaussian fluctuations provided by a $\chi^2$ model, as motivated by some models of inflation. The emphasis is in signatures that can be used to constrain this type of models from large-scale structure galaxy surveys. Non-Gaussian initial conditions provide additional non-linear couplings otherwise forbidden by symmetry that cause non-linear gravitational corrections to become important at larger scales than in the Gaussian case. In fact, the lack of hierarchical scaling in the initial conditions is partially restored by gravitational evolution at scales $k> 0.1$ h/Mpc. However, the bispectrum shows much larger amplitude and residual scale dependence not present in evolution from Gaussian initial conditions that can be used to test this model against observations. We include the effects of biasing and redshift distortions essential to compare this model with galaxy redshift surveys. We also discuss the effects of primordial non-Gaussianity on the redshift-space power spectrum and show that it changes the shape of the quadrupole to monopole ratio through non-linear corrections to infall velocities.Comment: 20 pages, 7 figure

Statistical and systematic errors in redshift-space distortion measurements from large surveys

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society © 2012 RAS © 2012 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.We investigate the impact of statistical and systematic errors on measurements of linear redshift-space distortions (RSD) in future cosmological surveys by analysing large catalogues of dark matter haloes from the baryonic acoustic oscillation simulations at the Institute for Computational Cosmology. These allow us to estimate the dependence of errors on typical survey properties, as volume, galaxy density and mass (i.e. bias factor) of the adopted tracer. We find that measures of the specific growth rate β = f/b using the Hamilton/Kaiser harmonic expansion of the redshift-space correlation function ξ(r p, π) on scales larger than 3h -1 Mpc are typically underestimated by up to 10 per cent for galaxy-sized haloes. This is significantly larger than the corresponding statistical errors, which amount to a few per cent, indicating the importance of non-linear improvements to the Kaiser model, to obtain accurate measurements of the growth rate. The systematic error shows a diminishing trend with increasing bias value (i.e. mass) of the haloes considered. We compare the amplitude and trends of statistical errors as a function of survey parameters to predictions obtained with the Fisher information matrix technique. This is what is usually adopted to produce RSD forecasts, based on the Feldman-Kaiser-Peacock prescription for the errors on the power spectrum. We show that this produces parameter errors fairly similar to the standard deviations from the halo catalogues, provided it is applied to strictly linear scales in Fourier space (k<0.2 h Mpc -1). Finally, we combine our measurements to define and calibrate an accurate scaling formula for the relative error on β as a function of the same parameters, which closely matches the simulation results in all explored regimes. This provides a handy and plausibly more realistic alternative to the Fisher matrix approach, to quickly and accurately predict statistical errors on RSD expected from future surveysEM is supported by the Spanish MICINNs Juan de la Cierva programme (JCI-2010-08112), by CICYT through the project FPA-2009 09017 and by the Community of Madrid through the project HEPHACOS (S2009/ESP-1473) under grant P-ESP-00346. Financial support of PRIN-INAF 2007, PRIN-MIUR 2008 and ASI Contracts I/023/05/0, I/088/06/0, I/016/07/0, I/064/08/0 and I/039/10/0 is gratefully acknowledged. LG is partly supported by ERC Advanced Grant #291521 ‘DARKLIGHT’

The Tully-Fisher Relation and H_not

The use of the Tully-Fisher (TF) relation for the determination of the Hubble Constant relies on the availability of an adequate template TF relation and of reliable primary distances. Here we use a TF template relation with the best available kinematical zero-point, obtained from a sample of 24 clusters of galaxies extending to cz ~ 9,000 km/s, and the most recent set of Cepheid distances for galaxies fit for TF use. The combination of these two ingredients yields H_not = 69+/-5 km/(s Mpc). The approach is significantly more accurate than the more common application with single cluster (e.g. Virgo, Coma) samples.Comment: 10 pages, including 2 figures and 1 table; uses AAS LaTex. Submitted to ApJ Letter

On the trispectrum as a gaussian test for cosmology

In the standard model for structure formation, bound objects originate from the gravitational collapse of small perturbations arising from quantum fluctuations with random phases. In other scenarios, based on defects, structures are seeded by localized energy density. In principle, it is possible to differentiate between these models on the basis of their statistical properties; only in the former case is the initial density field an almost-perfect random gaussian field. In this paper, we investigate the use of the trispectrum of the galaxy density field, which is the connected four-point function in Fourier space, as a discriminant between gaussian and non-gaussian models. It has the advantage of having only weak non-linear growth. We define a related statistic $\tau$ which, as a test of the gaussian hypothesis, is independent of cosmology, the power spectrum and biasing, in real space, and which is, in principle, a measure of the departure from gaussian statistics. For galaxy redshift surveys, the statistic depends on cosmology and bias only through the potentially observable parameter $\beta$. We compute the expected errors on the estimate of $\tau$, and demonstrate with numerical simulations that it can be a useful discriminant of models, with the important proviso that any bias is linear on large scales. Whether it is the most effective method is uncertain and depends on the nature of the departure from gaussianity.Comment: to appear in ApJ, 28 pages, 5 figure