Source author record

Federico Marulli

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

30 published item(s)

preprint2026arXiv

Kinematic lensing with high-resolution spectroscopic surveys. A unique opportunity for transformative cosmology at high redshifts in the 2040s

We present a science case to perform high-redshift cosmic shear surveys for cosmology with next-generation spectroscopic instruments, such as the proposed MegaMapper and Wide-field Spectroscopic Telescope. We argue that by using the novel technique called 'kinematic lensing' (KL) it will be possible to obtain shear catalogues at redshifts between 2 and 5. We show that the signal-to-noise ratio of KL at such high redshifts is on average twice as much that expected from current weak lensing (WL) surveys such as Euclid or LSST, and several times that of the previous generation of WL surveys like DES and KiDS, even with very conservative assumptions about the fraction of spectroscopically-detected sources for which KL shear estimates will be available. This will allow cosmologists to perform joint galaxy clustering-cosmic shear analyses over unprecedented cosmic volumes and to probe the growth of structures deep in the matter-dominated era and across the onset of dark-energy domination, offering a unique opportunity to unveil the mystery of cosmic acceleration.

preprint2026arXiv

Probing Cosmic Expansion and Early Universe with Einstein Telescope

Over the next two decades, gravitational-wave (GW) observations are expected to evolve from a discovery-driven endeavour into a precision tool for astrophysics, cosmology, and fundamental physics. Current second-generation ground-based detectors have established the existence of compact-binary mergers and enabled GW multi-messenger astronomy, but they remain limited in sensitivity, redshift reach, frequency coverage, and duty cycle. These limitations prevent them from addressing many fundamental open questions in cosmology. By the 2040s, wide-field electromagnetic surveys will have mapped the luminous Universe with unprecedented depth and accuracy. Nevertheless, key problems including the nature of dark matter, the physical origin of cosmic acceleration, the properties of gravity on cosmological scales, and the physical conditions of the earliest moments after the Big Bang will remain only partially constrained by electromagnetic observations alone. Progress on these fronts requires access to physical processes and epochs that do not emit light. Gravitational waves provide a unique and complementary observational channel: they propagate over cosmological distances largely unaffected by intervening matter, probe extreme astrophysical environments, and respond directly to the geometry of spacetime. In this context, next-generation GW observatories such as the Einstein Telescope (ET) will be transformative for European astronomy. Operating at sensitivities and frequencies beyond existing detectors, ET will observe binary black holes and neutron stars out to previously inaccessible redshifts, enable continuous high signal-to-noise monitoring of compact sources, and detect gravitational-wave backgrounds of astrophysical and cosmological origin. Together with space-based detectors, ET will play a central role in advancing our understanding of cosmic evolution and fundamental physics.

preprint2026arXiv

Probing General Relativity on Cosmological Scales in the 2040s

General relativity is exquisitely tested in strong-field regimes, yet its validity on cosmological scales remains largely unexplored. Upcoming wide and deep large-scale structure surveys will access the ultra-large, linear scales where relativistic effects - Doppler terms, gravitational redshift, lensing magnification, and potential evolution - leave significant imprints in the clustering of galaxies. These signatures represent unique probes of spacetime that are inaccessible to standard Newtonian analyses but increasingly important as survey volumes grow. We outline the scientific potential of next-generation facilities, such as those envisioned within ESO's Expanding Horizons programme, to deliver the first robust measurements of relativistic effects in large-scale structure through multi-tracer power spectra and the single-tracer bispectrum of high-redshift Lyman-break galaxies. Detecting these contributions would open a new window on gravity, enabling precision tests of general relativity and its alternatives on cosmological scales in the 2040s.

preprint2022arXiv

AMICO galaxy clusters in KiDS-DR3: Measurement of the halo bias and power spectrum normalization from a stacked weak lensing analysis

Galaxy clusters are biased tracers of the underlying matter density field. At very large radii beyond about 10 Mpc/\textit{h}, the shear profile shows evidence of a second-halo term. This is related to the correlated matter distribution around galaxy clusters and proportional to the so-called halo bias. We present an observational analysis of the halo bias-mass relation based on the AMICO galaxy cluster catalog, comprising around 7000 candidates detected in the third release of the KiDS survey. We split the cluster sample into 14 redshift-richness bins and derive the halo bias and the virial mass in each bin by means of a stacked weak lensing analysis. The observed halo bias-mass relation and the theoretical predictions based on the $Λ$CDM standard cosmological model show an agreement within $2σ$. The mean measurements of bias and mass over the full catalog give $M_{200c} = (4.9 \pm 0.3) \times 10^{13} M_{\odot}/\textit{h}$ and $b_h σ_8^2 = 1.2 \pm 0.1$. With the additional prior of a bias-mass relation from numerical simulations, we constrain the normalization of the power spectrum with a fixed matter density $Ω_m = 0.3$, finding $σ_8 = 0.63 \pm 0.10$.

preprint2021arXiv

The XXL survey: XLVI. Forward cosmological analysis of the C1 cluster sample

We present the forward cosmological analysis of an $XMM$ selected sample of galaxy clusters out to a redshift of unity. Following our previous 2018 study based on the dn/dz quantity alone, we perform an upgraded cosmological analysis of the same XXL C1 cluster catalogue (178 objects), with a detailed account of the systematic errors. We follow the ASpiX methodology: the distribution of the observed X-ray properties of the cluster population is analysed in a 3D observable space (count rate, hardness ratio, redshift) and modelled as a function of cosmology. Compared to more traditional methods, ASpiX allows the inclusion of clusters down to a few tens of photons. We obtain an improvement by a factor of 2 compared to the previous analysis by letting the normalisation of the M-T relation and the evolution of the L-T relation free. Adding constraints from the XXL cluster 2-point correlation function and the BAO from various surveys decreases the uncertainties by 23 and 53 % respectively, and 62% when adding both. Switching to the scaling relations from the Subaru analysis, and letting free more parameters, our final constraints are $σ_8$ = $0.99^{+0.14}_{-0.23}$, $Ω_m$ = 0.296 $\pm$ 0.034 ($S_8 = 0.98^{+0.11}_{-0.21}$) for the XXL sample alone. Finally, we combine XXL ASpiX, the XXL cluster 2-point correlation function and the BAO, with 11 free parameters, allowing for the cosmological dependence of the scaling relations in the fit. We find $σ_8$ = $0.793^{+0.063}_{-0.12}$, $Ω_m$ = 0.364 $\pm$ 0.015 ($S_8 = 0.872^{+0.068}_{-0.12}$), but still compatible with Planck CMB at 2.2$σ$. The results obtained by the ASpiX method are promising; further improvement is expected from the final XXL cosmological analysis involving a cluster sample twice as large. Such a study paves the way for the analysis of the eROSITA and future Athena surveys.

preprint2020arXiv

AMICO galaxy clusters in KiDS-DR3: galaxy population properties and their redshift dependence

A catalogue of galaxy clusters was obtained in an area of 414 sq deg up to a redshift $z\sim0.8$ from the Data Release 3 of the Kilo-Degree Survey (KiDS-DR3), using the Adaptive Matched Identifier of Clustered Objects (AMICO) algorithm. The catalogue and the calibration of the richness-mass relation were presented in two companion papers. Here we describe the selection of the cluster central galaxy and the classification of blue and red cluster members, and analyze the main cluster properties, such as the red/blue fraction, cluster mass, brightness and stellar mass of the central galaxy, and their dependence on redshift and cluster richness. We use the Illustris-TNG simulation, which represents the state-of-the-art cosmological simulation of galaxy formation, as a benchmark for the interpretation of the results. A good agreement with simulations is found at low redshifts ($z \le 0.4$), while at higher redshifts the simulations indicate a lower fraction of blue galaxies than what found in the KiDS-AMICO catalogue: we argue that this may be due to an underestimate of star-forming galaxies in the simulations. The selection of clusters with a larger magnitude difference between the two brightest central galaxies, which may indicate a more relaxed cluster dynamical status, improves the agreement between the observed and simulated cluster mass and stellar mass of the central galaxy. We also find that at a given cluster mass the stellar mass of blue central galaxies is lower than that of the red ones.

preprint2020arXiv

CoMaLit -- VI. Intrinsic scatter in stacked relations. The weak lensing AMICO galaxy clusters in KiDS-DR3

Unbiased and precise mass calibration of galaxy clusters is crucial to fully exploit galaxy clusters as cosmological probes. Stacking of weak lensing signal allows us to measure observable-mass relations down to less massive halos halos without extrapolation. We propose a Bayesian inference method to constrain the intrinsic scatter of the mass proxy in stacked analyses. The scatter of the stacked data is rescaled with respect to the individual scatter based on the number of binned clusters. We apply this method to the galaxy clusters detected with the AMICO (Adaptive Matched Identifier of Clustered Objects) algorithm in the third data release of the Kilo-Degree Survey. The results confirm the optical richness as a low scatter mass proxy. Based on the optical richness and the calibrated weak lensing mass-richness relation, mass of individual objects down to ~10^13 solar masses can be estimated with a precision of ~20 per cent.

preprint2020arXiv

The stellar-to-halo mass relation over the past 12 Gyr

Understanding how galaxy properties are linked to the dark matter halos they reside in, and how they co-evolve is a powerful tool to constrain the processes related to galaxy formation. The stellar-to-halo mass relation (SHMR) and its evolution over the history of the Universe provides insights on galaxy formation models and allows to assign galaxy masses to halos in N-body dark matter simulations. We use a statistical approach to link the observed galaxy stellar mass functions on the COSMOS field to dark matter halo mass functions from the DUSTGRAIN simulation and from a theoretical parametrization from z=0 to z=4. We also propose an empirical model to describe the evolution of the stellar-to-halo mass relation as a function of redshift. We calculate the star-formation efficiency (SFE) of galaxies and compare results with previous works and semi-analytical models.

preprint2020arXiv

Validating the methodology for constraining the linear growth rate from clustering anisotropies

Redshift-space clustering distortions provide one of the most powerful probes to test the gravity theory on the largest cosmological scales. We perform a systematic validation study of the state-of-the-art statistical methods currently used to constrain the linear growth rate from redshift-space distortions in the galaxy two-point correlation function. The numerical pipelines are tested on mock halo catalogues extracted from large N-body simulations of the standard cosmological framework. We consider both the monopole and quadrupole multipole moments of the redshift-space two-point correlation function, as well as the radial and transverse clustering wedges, in the comoving scale range $10<r[$\Mpch$]<55$. Moreover, we investigate the impact of redshift measurement errors on the growth rate and linear bias measurements due to the assumptions in the redshift-space distortion model. Considering both the dispersion model and two widely-used models based on perturbation theory, we find that the linear growth rate is underestimated by about $5-10\%$ at $z<1$, while limiting the analysis at larger scales, $r>30$ \Mpch, the discrepancy is reduced below $5\%$. At higher redshifts, we find instead an overall good agreement between measurements and model predictions. Though this accuracy is good enough for clustering analyses in current redshift surveys, the models have to be further improved not to introduce significant systematics in RSD constraints from next generation galaxy surveys. The effect of redshift errors is degenerate with the one of small-scale random motions, and can be marginalised over in the statistical analysis, not introducing any statistically significant bias in the linear growth constraints, especially at $z\geq1$.

preprint2016arXiv

CosmoBolognaLib: C++ libraries for cosmological calculations

We present the CosmoBolognaLib, a large set of Open Source C++ numerical libraries for cosmological calculations. CosmoBolognaLib is a living project aimed at defining a common numerical environment for cosmological investigations of the large-scale structure of the Universe. In particular, one of the primary focuses of this software is to help in handling astronomical catalogues, both real and simulated, measuring one-point, two-point and three-point statistics in configuration space, and performing cosmological analyses. In this paper, we discuss the main features of this software, providing an overview of all the available C++ classes implemented up to now. Both the CosmoBolognaLib and their associated doxygen documentation can be freely downloaded at https://github.com/federicomarulli/CosmoBolognaLib. We provide also some examples to explain how these libraries can be included in either C++ or Python codes.

preprint2015arXiv

Cosmic voids detection without density measurements

Cosmic voids are effective cosmological probes to discriminate among competing world models. Their identification is generally based on density or geometry criteria that, because of their very nature, are prone to shot noise. We propose two void finders that are based on dynamical criterion to select voids in Lagrangian coordinates and minimise the impact of sparse sampling. The first approach exploits the Zel'dovich approximation to trace back in time the orbits of galaxies located in voids and their surroundings, the second uses the observed galaxy-galaxy correlation function to relax the objects' spatial distribution to homogeneity and isotropy. In both cases voids are defined as regions of the negative velocity divergence, that can be regarded as sinks of the back-in-time streamlines of the mass tracers. To assess the performance of our methods we used a dark matter halo mock catalogue CoDECS, and compared the results with those obtained with the ZOBOV void finder. We find that the void divergence profiles are less scattered than the density ones and, therefore, their stacking constitutes a more accurate cosmological probe. The significance of the divergence signal in the central part of voids obtained from both our finders is 60% higher than for overdensity profiles in the ZOBOV case. The ellipticity of the stacked void measured in the divergence field is closer to unity, as expected, than what is found when using halo positions. Therefore our void finders are complementary to the existing methods, that should contribute to improve the accuracy of void-based cosmological tests.

preprint2015arXiv

Cosmic voids in coupled dark energy cosmologies: the impact of halo bias

In this work we analyse the properties of cosmic voids in standard and coupled dark energy cosmologies. Using large numerical simulations, we investigate the effects produced by the dark energy coupling on three statistics: the filling factor, the size distribution and the stacked profiles of cosmic voids. We find that the bias of the tracers of the density field used to identify the voids strongly influences the properties of the void catalogues, and, consequently, the possibility of using the identified voids as a probe to distinguish coupled dark energy models from the standard $Λ$CDM cosmology. In fact, on one hand coupled dark energy models are characterised by an excess of large voids in the cold dark matter distribution as compared to the reference standard cosmology, due to their higher normalisation of linear perturbations at low redshifts. Specifically, these models present an excess of large voids with $R_{eff}>20, 15, 12$ Mpc h^{-1}, at $z=0, 0.55, 1$, respectively. On the other hand, we do not find any significant difference in the properties of the void detected in the distribution of collapsed dark matter halos. These results imply that the tracer bias has a significant impact on the possibility of using cosmic void catalogues to probe cosmology.

preprint2015arXiv

Cosmology and fundamental physics with the Euclid satellite

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

preprint2015arXiv

Forecasts on neutrino mass constraints from the redshift-space two-point correlation function

We provide constraints on the accuracy with which the neutrino mass fraction, $f_ν$, can be estimated when exploiting measurements of redshift-space distortions, describing in particular how the error on neutrino mass depends on three fundamental parameters of a characteristic galaxy redshift survey: density, halo bias and volume. In doing this, we make use of a series of dark matter halo catalogues extracted from the BASICC simulation. The mock data are analysed via a Markov Chain Monte Carlo likelihood analysis. We find a fitting function that well describes the dependence of the error on bias, density and volume, showing a decrease in the error as the bias and volume increase, and a decrease with density down to an almost constant value for high density values. This fitting formula allows us to produce forecasts on the precision achievable with future surveys on measurements of the neutrino mass fraction. For example, a Euclid-like spectroscopic survey should be able to measure the neutrino mass fraction with an accuracy of $δf_ν \approx 6.7\times10^{-4}$, using redshift-space clustering once all the other cosmological parameters are kept fixed to the $Λ$CDM case.

preprint2015arXiv

Measuring the distance-redshift relation with the baryon acoustic oscillations of galaxy clusters

We analyse the largest spectroscopic samples of galaxy clusters to date, and provide observational constraints on the distance-redshift relation from baryon acoustic oscillations. The cluster samples considered in this work have been extracted from the Sloan Digital Sky Survey at three median redshifts, $z=0.2$, $z=0.3$, and $z=0.5$. The number of objects is $12910$, $42215$, and $11816$, respectively. We detect the peak of baryon acoustic oscillations for all the three samples. The derived distance constraints are: $r_s/D_V(z=0.2)=0.18 \pm 0.01$, $r_s/D_V(z=0.3)=0.124 \pm 0.004$ and $r_s/D_V(z=0.5)=0.080 \pm 0.002$. Combining these measurements, we obtain robust constraints on cosmological parameters. Our results are in agreement with the standard $Λ$ cold dark matter model. Specifically, we constrain the Hubble constant in a $Λ$CDM model, $H_0 = 64_{-9}^{+14} \, \mathrm{km} \, \mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, the density of curvature energy, in the $oΛ$CDM context, $Ω_K = -0.015_{-0.36}^{+0.34}$, and finally the parameter of the dark energy equation of state in the $ow$CDM case, $w = -1.01_{-0.44}^{+0.44}$. This is the first time the distance-redshift relation has been constrained using only the peak of baryon acoustic oscillations of galaxy clusters.

preprint2015arXiv

New constraints on $σ_8$ from a joint analysis of stacked gravitational lensing and clustering of galaxy clusters

The joint analysis of clustering and stacked gravitational lensing of galaxy clusters in large surveys can constrain the formation and evolution of structures and the cosmological parameters. On scales outside a few virial radii, the halo bias, $b$, is linear and the lensing signal is dominated by the correlated distribution of matter around galaxy clusters. We discuss a method to measure the power spectrum amplitude $σ_8$ and $b$ based on a minimal modelling. We considered a sample of $\sim 120000$ clusters photometrically selected from the Sloan Digital Sky Survey in the redshift range $0.1<z<0.6$. The auto-correlation was studied through the two-point function of a subsample of $\sim 70000$ clusters; the matter-halo correlation was derived from the weak lensing signal of the subsample of $\sim 1200$ clusters with Canada-France-Hawaii Lensing Survey data. We obtained a direct measurement of $b$, which increases with mass in agreement with predictions of the $Λ$CDM paradigm. Assuming $Ω_\mathrm{M}=0.3$, we found $σ_8=0.78\pm0.16$. We used the same clusters for measuring both lensing and clustering and the estimate of $σ_8$ did require neither the mass-richness relation, nor the knowledge of the selection function, nor the modelling of $b$. With an additional theoretical prior on the bias, we obtained $σ_8=0.75\pm0.08$.

preprint2014arXiv

Disentangling interacting dark energy cosmologies with the three-point correlation function

We investigate the possibility of constraining coupled dark energy (cDE) cosmologies using the three-point correlation function (3PCF). Making use of the CoDECS N-body simulations, we study the statistical properties of cold dark matter (CDM) haloes for a variety of models, including a fiducial $Λ$CDM scenario and five models in which dark energy (DE) and CDM mutually interact. We measure both the halo 3PCF, $ζ(θ)$, and the reduced 3PCF, $Q(θ)$, at different scales ($2<r\,[$Mpc\h$]<40$) and redshifts ($0\leq z\leq2$). In all cDE models considered in this work, $Q(θ)$ appears flat at small scales (for all redshifts) and at low redshifts (for all scales), while it builds up the characteristic V-shape anisotropy at increasing redshifts and scales. With respect to the $Λ$CDM predictions, cDE models show lower (higher) values of the halo 3PCF for perpendicular (elongated) configurations. The effect is also scale-dependent, with differences between $Λ$CDM and cDE models that increase at large scales. We made use of these measurements to estimate the halo bias, that results in fair agreement with the one computed from the two-point correlation function (2PCF). The main advantage of using both the 2PCF and 3PCF is to break the bias$-σ_{8}$ degeneracy. Moreover, we find that our bias estimates are approximately independent of the assumed strength of DE coupling. This study demonstrates the power of a higher-order clustering analysis in discriminating between alternative cosmological scenarios, for both present and forthcoming galaxy surveys, such as e.g. BOSS and Euclid.

preprint2013arXiv

Characterizing dark interactions with the halo mass accretion history and structural properties

We study the halo mass accretion history (MAH) and its correlation with the internal structural properties in coupled dark energy cosmologies (cDE). To accurately predict all the non-linear effects caused by dark interactions, we use the COupled Dark Energy Cosmological Simulations (CoDECS). We measure the halo concentration at z=0 and the number of substructures above a mass resolution threshold for each halo. Tracing the halo merging history trees back in time, following the mass of the main halo, we develope a MAH model that accurately reproduces the halo growth in term of M_{200} in the ΛCDM Universe; we then compare the MAH in different cosmological scenarios. For cDE models with a weak constant coupling, our MAH model can reproduce the simulation results, within 10% of accuracy, by suitably rescaling the normalization of the linear matter power spectrum at z=0, σ_8. However, this is not the case for more complex scenarios, like the "bouncing" cDE model, for which the numerical analysis shows a rapid growth of haloes at high redshifts, that cannot be reproduced by simply rescaling the value of σ_8. Moreover, at fixed value of σ_8, cold dark matter (CDM) haloes in these cDE scenarios tend to be more concentrated and have a larger amount of substructures with respect to ΛCDM predictions. Finally, we present an accurate model that relates the halo concentration to the time at which it assembles half or 4% of its mass. Combining this with our MAH model, we show how halo concentrations change while varying only σ_8 in a ΛCDM Universe, at fixed halo mass.

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.

preprint2013arXiv

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

We investigate the impact of statistical and systematic errors on measurements of linear redshift-space distortions (RSD) in future cosmological surveys, analyzing large catalogues of dark-matter halos from the BASICC simulation. 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 3/h Mpc are typically under-estimated by up to 10% for galaxy sized halos. This is significantly larger than the corresponding statistical errors, which amount to a few percent, indicating the importance of non-linear improvements to the Kaiser model to obtain accurate measurements of the growth rate. We compare the statistical errors to predictions obtained with the Fisher information matrix, based on the usual FKP 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, but only if applied to strictly linear scales in Fourier space (k<0.2 h/Mpc). Finally, we present an accurate scaling formula describing the relative error on β as a function of the survey 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 RSD statistical errors expected from future surveys.

preprint2012arXiv

Clustering and redshift-space distortions in interacting dark energy cosmologies

We investigate the spatial properties of the large scale structure (LSS) of the Universe in the framework of coupled dark energy (cDE) cosmologies. Using the public halo catalogues from the CoDECS simulations -- the largest set of N-body experiments to date for such cosmological scenarios -- we estimate the clustering and bias functions of cold dark matter (CDM) haloes, both in real- and redshift-space. Moreover, we investigate the effects of the dark energy (DE) coupling on the geometric and dynamic redshift-space distortions, quantifying the difference with respect to the concordance LambdaCDM model. At z~0, the spatial properties of CDM haloes in cDE models appear very similar to the LambdaCDM case, even if the cDE models are normalized at last scattering in order to be consistent with the latest Cosmic Microwave Background (CMB) data. At higher redshifts, we find that the DE coupling produces a significant scale-dependent suppression of the halo clustering and bias function. This effect, that strongly depends on the coupling strength, is not degenerate with sigma8 at scales r<5-10 Mpc/h. Moreover, we find that the coupled DE strongly affects both the linear distortion parameter, beta, and the pairwise peculiar velocity dispersion, sigma12. Although the models considered in this work are found to be all in agreement with presently available observational data, the next generation of galaxy surveys will be able to put strong constraints on the level of coupling between DE and CDM exploiting the shape of redshift-space clustering anisotropies.

preprint2012arXiv

Cosmology with clustering anisotropies: disentangling dynamic and geometric distortions in galaxy redshift surveys

We investigate the impact of different observational effects affecting a precise and accurate measurement of the growth rate of fluctuations from the anisotropy of clustering in galaxy redshift surveys. We focus on redshift measurement errors, on the reconstruction of the underlying real-space clustering and on the apparent degeneracy existing with the geometrical distortions induced by the cosmology-dependent conversion of redshifts into distances. We use a suite of mock catalogues extracted from large N-body simulations, focusing on the analysis of intermediate, mildly non-linear scales and apply the standard linear dispersion model to fit the anisotropy of the observed correlation function. We verify that redshift errors up to ~0.2% have a negligible impact on the precision with which the specific growth rate beta can be measured. Larger redshift errors introduce a positive systematic error, which can be alleviated by adopting a Gaussian distribution function of pairwise velocities. This is, in any case, smaller than the systematic error of up to 10% due to the limitations of the linear dispersion model, which is studied in a separate paper. We then show that 50% of the statistical error budget on beta depends on the deprojection procedure through which the real-space correlation function is obtained. Finally, we demonstrate that the degeneracy with geometric distortions can in fact be circumvented. This is obtained through a modified version of the Alcock-Paczynski test in redshift-space, which successfully recovers the correct cosmology by searching for the solution that optimizes the description of dynamical redshift distortions. For a flat cosmology, we obtain largely independent, robust constraints on beta and OmegaM. In a volume of 2.4(Gpc/h)^3, the correct OmegaM is obtained with ~12% error and negligible bias, once the real-space correlation function is properly reconstructed.

preprint2012arXiv

Exploiting the shift of baryonic acoustic oscillations as a dynamical probe for dark interactions

The baryonic acoustic peak in the correlation function of galaxies and galaxy clusters provides a standard ruler to probe the space-time geometry of the Universe, jointly constraining the angular diameter distance and the Hubble expansion rate. Moreover, non-linear effects can systematically shift the peak position, giving us the opportunity to exploit this clustering feature also as a dynamical probe. We investigate the possibility of detecting interactions in the dark sector through an accurate determination of the baryonic acoustic scale. Making use of the public halo catalogues extracted from the CoDECS simulations -- the largest suite of N-body simulations of interacting dark energy models to date -- we determine the position of the baryonic scale fitting a band-filtered correlation function, specifically designed to amplify the signal at the sound horizon. We analyze the shifts due to non-linear dynamics, redshift-space distortions and Gaussian redshift errors, in the range 0 < z < 2. Since the coupling between dark energy and dark matter affects in a particular way the clustering properties of haloes and, specifically, the amplitude and location of the baryonic acoustic oscillations, the cosmic evolution of the baryonic peak position might provide a direct way to discriminate interacting dark energy models from the standard ΛCDM framework. To maximize the efficiency of the baryonic peak as a dynamic probe, the correlation function has to be measured in redshift-space, where the baryonic acoustic shift due to non-linearities is amplified. The typical redshift errors of spectroscopic galaxy surveys do not significantly impact these results.

preprint2012arXiv

Size Evolution of Spheroids in a Hierarchical Universe

We present basic predictions of an updated version of the Munich semi-analytic hierarchical galaxy formation model that grows bulges via mergers and disk instabilities. Overall, we find that while spheroids below Ms ~ 10^11 Msun grow their sizes via a mixture of disk instability and mergers, galaxies above it mainly evolve via mergers. Including gas dissipation in major mergers, efficiently shrinks galaxies, especially those with final mass Ms < 10^11 Msun that are the most gas-rich, improving the match with different observables. We find that the predicted scatter in sizes at fixed stellar mass is still larger than the observed one by up to <40%. Spheroids are, on average, more compact at higher redshifts at fixed stellar mass, and at fixed redshift and stellar mass larger galaxies tend to be more starforming. More specifically, while for bulge-dominated galaxies the model envisages a nearly mass-independent decrease in sizes, the predicted size evolution for intermediate-mass galaxies is more complex. The z=2 progenitors of massive galaxies with mass around Ms and B/T>0.7 at z=0, are found to be mostly disc-dominated galaxies with a median B/T ~ 0.3, with only ~20% remaining bulge-dominated. The model also predicts that central spheroids living in more massive haloes tend to have larger sizes at fixed stellar mass. Including host halo mass dependence in computing velocity dispersions, allows the model to properly reproduce the correlations with stellar mass. We also discuss the fundamental plane, the correlations with galaxy age, the structural properties of pseudobulges, and the correlations with central black holes.

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.

preprint2011arXiv

Expected properties of the Two-Point Autocorrelation Function of the IGM

Recent analyses of the fluctuations of the soft Diffuse X-ray Background (DXB) have provided indirect detection of a component consistent with the elusive Warm Hot Intergalactic Medium (WHIM). In this work we use theoretical predictions obtained from hydrodynamical simulations to investigate the angular correlation properties of the WHIM in emission and assess the possibility of indirect detection with next-generation X-ray missions. Our results indicate that the angular correlation signal of the WHIM is generally weak but dominates the angular correlation function of the DXB outside virialized regions. Its indirect detection is possible but requires rather long exposure times [0.1-1] Ms, large (~1° x1°) fields of view and accurate subtraction of isotropic fore/background contributions, mostly contributed by Galactic emission. The angular correlation function of the WHIM is positive for θ < 5' and provides limited information on its spatial distribution. A satisfactory characterization of the WHIM in 3D can be obtained through spatially resolved spectroscopy. 1 Ms long exposures with next generation detectors will allow to detect ~400 O VII+O VIII X-ray emission systems that we use to trace the spatial distribution of the WHIM. We predict that these observations will allow to estimate the WHIM correlation function with high statistical significance out to ~10 Mpc h^-1 and characterize its dynamical state through the analysis of redshift-space distortions. The detectable WHIM, which is typically associated with the outskirts of virialized regions rather than the filaments has a non-zero correlation function with slope γ = -1.7 \pm 0.1 and correlation length r0 = 4.0 \pm 0.1 Mpc h^-1 in the range r = [4.5, 12] Mpc h^-1. Redshift space distances can be measured to assess the dynamical properties of the gas, typically infalling onto large virialized structures.

preprint2010arXiv

Further constraining galaxy evolution models through the Size Function of SDSS Early-type galaxies

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

preprint2009arXiv

Modeling the cosmological co-evolution of supermassive black holes and galaxies: II. The clustering of quasars and their dark environment

We use semi-analytic modeling on top of the Millennium simulation to study the joint formation of galaxies and their embedded supermassive black holes. Our goal is to test scenarios in which black hole accretion and quasar activity are triggered by galaxy mergers, and to constrain different models for the lightcurves associated with individual quasar events. In the present work we focus on studying the spatial distribution of simulated quasars. At all luminosities, we find that the simulated quasar two-point correlation function is fit well by a single power-law in the range 0.5 < r < 20 h^{-1} Mpc, but its normalization is a strong function of redshift. When we select only quasars with luminosities within the range typically accessible by today's quasar surveys, their clustering strength depends only weakly on luminosity, in agreement with observations. This holds independently of the assumed lightcurve model, since bright quasars are black holes accreting close to the Eddington limit, and are hosted by dark matter haloes with a narrow mass range of a few 10^12 h^{-1} M_sun. Therefore the clustering of bright quasars cannot be used to disentangle lightcurve models, but such a discrimination would become possible if the observational samples can be pushed to significantly fainter limits. Overall, our clustering results for the simulated quasar population agree rather well with observations, lending support to the conjecture that galaxy mergers could be the main physical process responsible for triggering black hole accretion and quasar activity.

preprint2009arXiv

Sizes and ages of SDSS ellipticals: Comparison with hierarchical galaxy formation models

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

preprint2009arXiv

The spatial distribution of X-ray selected AGN in the Chandra deep fields: a theoretical perspective

We study the spatial distribution of X-ray selected AGN in the framework of hierarchical co-evolution of supermassive black holes and their host galaxies and dark matter haloes. To this end, we have applied the model developed by Croton et al.(2006), De Lucia & Blaizot(2007) and Marulli et al.(2008) to the output of the Millennium Run and obtained hundreds of realizations of past light-cones from which we have extracted realistic mock AGN catalogues that mimic the Chandra deep fields. We find that the model AGN number counts are in fair agreement with observations, except at fluxes <1e-15 erg/cm^2/s. The spatial two-point correlation function predicted by the model is well described by a power-law relation out to 20 Mpc/h, in close agreement with observations. Our model matches the correlation length r_0 of AGN in the Chandra Deep Field North but underestimates it in the Chandra Deep Field South. When fixing the slope to gamma = 1.4, as in Gilli et al. (2005), the statistical significance of the mismatch is 2-2.5 sigma, suggesting that the predicted cosmic variance, which dominates the error budget, may not account for the different correlation length of the AGN in the two fields. While our results are robust to changes in the model prescriptions for the AGN lightcurves, the luminosity dependence of the clustering is sensitive to the different lightcurve models adopted. However, irrespective of the model considered, the luminosity dependence of the AGN clustering in our mock fields seems to be weaker than in the real Chandra fields. The significance of this mismatch needs to be confirmed using larger datasets.