Source author record

Lauro Moscardini

Lauro Moscardini 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

66works
5topics
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

66 published item(s)

preprint2025arXiv

KiDS-Legacy: Cosmological constraints from cosmic shear with the complete Kilo-Degree Survey

We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-infrared, along with an extra $23$ square degrees of KiDS-like calibration observations of deep spectroscopic surveys. Improvements in our redshift distribution estimation methodology, combined with our enhanced calibration data and multi-band image simulations, allow us to extend our lensed sample out to a photometric redshift of $z_{\rm B}\leq2.0$. Compared to previous KiDS analyses, the increased survey area and redshift depth results in a $\sim32\%$ improvement in constraining power in terms of $Σ_8\equivσ_8\left(Ω_{\rm m}/0.3\right)^α= 0.821^{+0.014}_{-0.016}$, where $α= 0.58$ has been optimised to match the revised degeneracy direction of $σ_8$ and $Ω_{\rm m}$. We adopt a new physically motivated intrinsic alignment model that depends jointly on the galaxy sample's halo mass and spectral type distributions, and that is informed by previous direct alignment measurements. We also marginalise over our uncertainty on the impact of baryon feedback on the non-linear matter power spectrum. Comparing to previous KiDS analyses, we conclude that the increase seen in $S_8$ primarily results from our improved redshift distribution estimation and calibration, as well as new survey area and improved image reduction. Our companion paper Stölzner et al. (submitted) presents a full suite of internal and external consistency tests, finding the KiDS-Legacy data set to be the most internally robust sample produced by KiDS to date.

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

preprint2022arXiv

Galaxies in the central regions of simulated galaxy clusters

In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm.

preprint2021arXiv

AMICO galaxy clusters in KiDS-DR3: the impact of estimator statistics on the luminosity-mass scaling relation

As modern-day precision cosmology aims for statistical uncertainties of the percent level or lower, it becomes increasingly important to reconsider estimator assumptions at each step of the process, and their consequences on the statistical variability of the scientific results. We compare $L^1$ regression statistics to the weighted mean, the canonical $L^2$ method based on Gaussian assumptions, for inference of the weak gravitational shear signal from a catalog of background ellipticity measurements around a sample of clusters, in many recent analyses a standard step in the process. We use the shape measurements of background sources around 6925 AMICO clusters detected in the KiDS 3rd data release. We investigate the robustness of our results and the dependence of uncertainties on the signal-to-noise ratios of the background source detections. Using a halo model approach, we derive lensing masses from the estimated excess surface density profiles. The highly significant shear signal allows us to study the scaling relation between the $r$-band cluster luminosity $L_{200}$, and the derived lensing mass $M_{200}$. We show the results of the scaling relations derived in 13 bins in $L_{200}$, with a tightly constrained power law slope of $\sim 1.24\pm 0.08$. We observe a small, but significant relative bias of a few percent in the recovered excess surface density profiles between the two regression methods, which translates to a $1σ$ difference in $M_{200}$. The efficiency of $L^1$ is at least that of the weighted mean, relatively increasing with higher signal-to-noise shape measurements. Our results indicate the relevance of optimizing the estimator for infering the gravitational shear from a distribution of background ellipticities. The interpretation of measured relative biases can be gauged by deeper observations, while increased computation times remain feasible.

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

Testing the Reliability of Fast Methods for Weak Lensing Simulations: WL-MOKA on PINOCCHIO

The generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes.

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

Characterising Strong Lensing Galaxy Clusters using the Millennium-XXL and MOKA simulations

In this paper we investigate the strong lensing statistics in galaxy clusters. We extract dark matter haloes from the Millennium-XXL simulation, compute their Einstein radius distribution, and find a very good agreement with Monte Carlo predictions produced with the MOKA code. The distribution of the Einstein radii is well described by a log-normal distribution, with a considerable fraction of the largest systems boosted by different projection effects. We discuss the importance of substructures and triaxiality in shaping the size of the critical lines for cluster size haloes. We then model and interpret the different deviations, accounting for the presence of a Bright Central Galaxy (BCG) and two different stellar mass density profiles. We present scaling relations between weak lensing quantities and the size of the Einstein radii. Finally we discuss how sensible is the distribution of the Einstein radii on the cosmological parameters Ω_M-σ_8 finding that cosmologies with higher Ω_M and σ_8 possess a large sample of strong lensing clusters. The Einstein radius distribution may help distinguish Planck13 and WMAP7 cosmology at 3σ.

preprint2016arXiv

Comparison of weak lensing by NFW and Einasto halos and systematic errors

Recent N-body simulations have shown that Einasto radial profiles provide the most accurate description of dark matter halos. Predictions based on the traditional NFW functional form may fail to describe the structural properties of cosmic objects at the percent level required by precision cosmology. We computed the systematic errors expected for weak lensing analyses of clusters of galaxies if one wrongly models the lens density profile. Even though the NFW fits of observed tangential shear profiles can be excellent, viral masses and concentrations of very massive halos (>~ 10^{15}M_Sun/h) can be over- and underestimated by ~10 per cent, respectively. Misfitting effects also steepen the observed mass-concentration relation, as observed in multi-wavelength observations of galaxy groups and clusters. Based on shear analyses, Einasto and NFW halos can be set apart either with deep observations of exceptionally massive structures (>~ 2\times10^{15}M_Sun/h) or by stacking the shear profiles of thousands of group-sized lenses (>~ 10^{14}M_Sun/h).

preprint2016arXiv

Cosmology through arc statistics I: sensitivity to $Ω_m$ and $σ_8$

The next generation of large sky photometric surveys will finally be able to use arc statistics as a cosmological probe. Here we present the first of a series of papers on this topic. In particular, we study how arc counts are sensitive to the variation of two cosmological parameters: the (total) matter density parameter, $Ω_m$, and the normalisation of the primordial power spectrum, expressed in terms of $σ_8$. Both these parameters influence the abundances of collapsed structures and their internal structure. We compute the expected number of gravitational arcs with various length-to-width ratios in mock light cones, by varying these cosmological parameters in the ranges $0.1\leqΩ_m\leq0.5$ and $0.6\leqσ_8\leq 1$. We find that the arc counts dependence on $Ω_m$ and $σ_8$ is similar, but not identical, to that of the halo counts. We investigate how the precision of the constraints on the cosmological parameters based on arc counts depends on the survey area. We find that the constraining power of arc statistics degrades critically only for surveys covering an area smaller than $10\%$ of the whole sky. Finally, we consider the case in which the search for arcs is done only in frames where galaxy clusters have been previously identified. Adopting the selection function for galaxy clusters expected to be detected from photometric data in future wide surveys, we find that less than $10\%$ of the arcs will be missed, with only a small degradation of the corresponding cosmological constraints.

preprint2015arXiv

CoMaLit - II. The scaling relation between mass and Sunyaev-Zel'dovich signal for Planck selected galaxy clusters

We discuss the scaling relation between mass and integrated Compton parameter of a sample of galaxy clusters from the all-sky {\it Planck} Sunyaev-Zel'dovich catalogue. Masses were measured with either weak lensing, caustics techniques, or assuming hydrostatic equilibrium. The retrieved $Y_{500}$-$M_{500}$ relation does not strongly depend on the calibration sample. We found a slope of 1.4-1.9, in agreement with self-similar predictions, with an intrinsic scatter of $20\pm10$ per cent. The absolute calibration of the relation can not be ascertained due to systematic differences of $\sim$20-40 per cent in mass estimates reported by distinct groups. Due to the scatter, the slope of the conditional scaling relation, to be used in cosmological studies of number counts, is shallower, $\sim$1.1-1.6. The regression methods employed account for intrinsic scatter in the mass measurements too. We found that Planck mass estimates suffer from a mass dependent bias.

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

Disentangling dark sector models using weak lensing statistics

We perform multi-plane ray-tracing using the GLAMER gravitational lensing code within high-resolution light-cones extracted from the CoDECS simulations: a suite of cosmological runs featuring a coupling between Dark Energy and Cold Dark Matter. We show that the presence of the coupling is evident not only in the redshift evolution of the normalisation of the convergence power spectrum, but also in differences in non-linear structure formation with respect to ΛCDM. Using a tomographic approach under the assumption of a ΛCDM cosmology, we demonstrate that weak lensing measurements would result in a σ8 value that changes with the source redshift if the true underlying cosmology is a coupled Dark Energy one. This provides a generic null test for these types of models. We also find that different models of coupled Dark Energy can show either an enhanced or a suppressed correlation between convergence maps with differing source redshifts as compared to ΛCDM. This would provide a direct way to discriminate between different possible realisations of the coupled Dark Energy scenario. Finally, we discuss the impact of the coupling on several lensing observables for different source redshifts and angular scales with realistic source redshift distributions for current ground-based and future space-based lensing surveys.

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

Neglecting Primordial non-Gaussianity Threatens Future Cosmological Experiment Accuracy

Future galaxy redshift surveys aim at probing the clustering of the cosmic large-scale structure with unprecedented accuracy, thus complementing cosmic microwave background experiments in the quest to deliver the most precise and accurate picture ever of our Universe. Analyses of such measurements are usually performed within the context of the so-called vanilla LCDM model - the six-parameter phenomenological model which, for instance, emerges from best fits against the recent data obtained by the Planck satellite. Here, we show that such an approach is prone to subtle systematics when the Gaussianity of primordial fluctuations is concerned. In particular, we demonstrate that, if we neglect even a tiny amount of primordial non-Gaussianity - fully consistent with current limits - we shall introduce spurious biases in the reconstruction of cosmological parameters. This is a serious issue that must be properly accounted for in view of accurate (as well as precise) cosmology.

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

preprint2015arXiv

Searching for galaxy clusters in the VST-KiDS Survey

We present the methods and first results of the search for galaxy clusters in the Kilo Degree Survey (KiDS). The adopted algorithm and the criterium for selecting the member galaxies are illustrated. Here we report the preliminary results obtained over a small area (7 sq. degrees), and the comparison of our cluster candidates with those found in the RedMapper and SZ Planck catalogues; the analysis to a larger area (148 sq. degrees) is currently in progress. By the KiDS cluster search, we expect to increase the completeness of the clusters catalogue to z = 0.6-0.7 compared to RedMapper.

preprint2015arXiv

The mass-concentration relation in lensing clusters: the role of statistical biases and selection effects

The relation between mass and concentration of galaxy clusters traces their formation and evolution. Massive lensing clusters were observed to be over-concentrated and following a steep scaling in tension with predictions from the concordance $Λ$CDM paradigm. We critically revise the relation in the CLASH, the SGAS, the LOCUSS, and the high-redshift samples of weak lensing clusters. Measurements of mass and concentration are anti-correlated, which can bias the observed relation towards steeper values. We corrected for this bias and compared the measured relation to theoretical predictions accounting for halo triaxiality, adiabatic contraction of the halo, presence of a dominant BCG and, mostly, selection effects in the observed sample. The normalisation, the slope and the scatter of the expected relation are strongly sample-dependent. For the considered samples, the predicted slope is much steeper than that of the underlying relation characterising dark-matter only clusters. We found that the correction for statistical and selection biases in observed relations mostly solve the tension with the $Λ$CDM model.

preprint2015arXiv

The universality of the virial halo mass function and models for non-universality of other halo definitions

The abundance of galaxy clusters can constrain both the geometry and growth of structure in our Universe. However, this probe could be significantly complicated by recent claims of nonuniversality -- non-trivial dependences with respect to the cosmological model and redshift. In this work we analyse the dependance of the mass function on the way haloes are identified and establish if this can cause departures from universality. In order to explore this dependance, we use a set of different N-body cosmological simulations (Le SBARBINE simulations), with the latest cosmological parameters from the Planck collaboration; this first suite of simulations is followed by a lower resolution set, carried out with different cosmological parameters. We identify dark matter haloes using a Spherical Overdensity algorithm with varying overdensity thresholds (virial, 2000rho_c, 1000rho_c, 500rho_c, 200rho_c and 200rho_b) at all redshifts. We notice that, when expressed in term of the rescaled variable nu, the mass functionfor virial haloes is a nearly universal as a function of redshift and cosmology, while this is clearly not the case for the other overdensities we considered. We provide fitting functions for the halo mass function parameters as a function of overdensity, that allow to predict, to within a few percent accuracy, the halo mass function for a wide range of halo definitions, redshifts and cosmological models. We then show how the departures from universality associated with other halo definitions can be derived by combining the universality of the virial definition with the expected shape of the density profile of halos.

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

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.

preprint2014arXiv

Lensing dispersion of supernova flux: a probe of nonlinear structure growth

The scatter in the apparent magnitude of type Ia supernovae induced by stochastic gravitational lensing is highly dependent on the nonlinear growth of cosmological structure. In this paper, we show that such a dependence can potentially be employed to gain significant information about the mass clustering at small scales. While the mass clustering ultimately hinges on cosmology, here we demonstrate that, upon obtaining more precise observational measurements through future cosmological surveys, the lensing dispersion can very effectively be used to gain information on the poorly understood astrophysical aspects of structure formation, such as the clumpiness of dark matter halos and the importance of gas physics and star formation into shaping the large-scale structure. In order to illustrate this point we verify that even the tentative current measurements of the lensing dispersion performed on the Supernova Legacy Survey sample favor a scenario where virialized structures are somewhat less compact than predicted by $n-$body cosmological simulations. Moreover, we are also able to put lower limits on the slope of the concentration-mass relation. By artificially reducing the statistical observational error we argue that with forthcoming data the stochastic lensing dispersion will allow one to importantly improve constraints on the baryonic physics at work during the assembly of cosmological structure.

preprint2014arXiv

Magnification bias as a novel probe for primordial magnetic fields

In this paper we investigate magnetic fields generated in the early Universe. These fields are important candidates at explaining the origin of astrophysical magnetism observed in galaxies and galaxy clusters, whose genesis is still by and large unclear. Compared to the standard inflationary power spectrum, intermediate to small scales would experience further substantial matter clustering, were a cosmological magnetic field present prior to recombination. As a consequence, the bias and redshift distribution of galaxies would also be modified. Hitherto, primordial magnetic fields (PMFs) have been tested and constrained with a number of cosmological observables, e.g. the cosmic microwave background radiation, galaxy clustering and, more recently, weak gravitational lensing. Here, we explore the constraining potential of the density fluctuation bias induced by gravitational lensing magnification onto the galaxy-galaxy angular power spectrum. Such an effect is known as magnification bias. Compared to the usual galaxy clustering approach, magnification bias helps in lifting the pathological degeneracy present amongst power spectrum normalisation and galaxy bias. This is because magnification bias cross-correlates galaxy number density fluctuations of nearby objects with weak lensing distortions of high-redshift sources. Thus, it takes advantage of the gravitational deflection of light, which is insensitive to galaxy bias but powerful in constraining the density fluctuation amplitude. To scrutinise the potentiality of this method, we adopt a deep and wide-field spectroscopic galaxy survey. We show that magnification bias does contain important information on primordial magnetism, which will be useful in combination with galaxy clustering and shear. We find we shall be able to rule out at 95.4% CL amplitudes of PMFs larger than 0.0005 nG for values of the PMF power spectral index ~0.

preprint2014arXiv

Mass and Concentration estimates from Weak and Strong Gravitational Lensing: a Systematic Study

We study how well halo properties of galaxy clusters, like mass and concentration, are recovered using lensing data. In order to generate a large sample of systems at different redshifts we use the code MOKA. We measure halo mass and concentration using weak lensing data alone (WL), fitting to an NFW profile the reduced tangential shear profile, or by combining weak and strong lensing data, by adding information about the size of the Einstein radius (WL+SL). For different redshifts, we measure the mass and the concentration biases and find that these are mainly caused by the random orientation of the halo ellipsoid with respect to the line-of-sight. Since our simulations account for the presence of a bright central galaxy, we perform mass and concentration measurements using a generalized NFW profile which allows for a free inner slope. This reduces both the mass and the concentration biases. We discuss how the mass function and the concentration mass relation change when using WL and WL+SL estimates. We investigate how selection effects impact the measured concentration-mass relation showing that strong lens clusters may have a concentration 20-30% higher than the average, at fixed mass, considering also the particular case of strong lensing selected samples of relaxed clusters. Finally, we notice that selecting a sample of relaxed galaxy clusters, as is done in some cluster surveys, explain the concentration-mass relation biases.

preprint2014arXiv

Raytracing simulations of coupled dark energy models

Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth-force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak-lensing maps using raytracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance $Λ$CDM model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on CMB lensing by Carbone et al., 2013. The most significant differences from the $Λ$CDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40% compared to the $Λ$CDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a $Λ$CDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.

preprint2014arXiv

The effect of massive neutrinos on the Sunyaev-Zeldovich and X-ray observables of galaxy clusters

Massive neutrinos are expected to influence the formation of the large-scale structure of the Universe, depending on the value of their total mass, $Σm_ν$. In particular Planck data indicate that a non-zero $Σm_ν$ may help to reconcile CMB data with Sunyaev-Zel'dovich (SZ) cluster surveys. In order to study the impact of neutrinos on the SZ and X-ray cluster properties we run a set of six very large cosmological simulations (8$h^{-3}$ Gpc$^3$ comoving volume) that include a massive neutrino particle component: we consider the values of $Σm_ν$ = (0, 0.17, 0.34) eV in two cosmological scenarios to test possible degeneracies. Using the halo catalogues extracted from their outputs we produce 50 mock light-cones and, assuming suitable scaling relations, we determine how massive neutrinos affect SZ and X-ray cluster counts, the $y$-parameter and its power spectrum. We provide forecasts for the South Pole Telescope (SPT) and eROSITA cluster surveys, showing that the number of expected detections is reduced by 40 per cent when assuming $Σm_ν$ =0.34 eV with respect to a model with massless neutrinos. However the degeneracy with $σ_8$ and $Ω_m$ is strong, in particular for X-ray data, requiring the use of additional probes to break it. The $y$-parameter properties are also highly influenced by the neutrino mass fraction, $f_ν$, with $<\,y\,>\propto(1-f_ν)^{20}$, considering the cluster component only, and the normalization of the SZ power spectrum is proportional to $(1-f_ν)^{25-30}$. Comparing our findings with SPT and Atacama Cosmology Telescope measurements at $\ell$ = 3000 indicates that, when Planck cosmological parameters are assumed, a value of $Σm_ν\simeq0.34$ eV is required to fit with the data.

preprint2013arXiv

A comparison of structure formation in minimally and non-minimally coupled quintessence models

We study structure formation in non-minimally coupled dark energy models, where there is a coupling in the Lagrangian between a quintessence scalar field and gravity via the Ricci scalar. We consider models with a range of different non-minimal coupling strengths and compare these to minimally coupled quintessence models with time-dependent dark energy densities. The equations of state of the latter are tuned to either reproduce the equation of state of the non-minimally coupled models or their background history. Thereby they provide a reference to study the unique imprints of coupling on structure formation. We show that the coupling between gravity and the scalar field, which effectively results in a time-varying gravitational constant G, is not negligible and its effect can be distinguished from a minimally coupled model. We extend previous work on this subject by showing that major differences appear in the determination of the mass function at high masses, where we observe differences of the order of 40% at z=0. Our new results concern effects on the non-linear matter power spectrum and on the lensing signal (differences of ~10% for both quantities), where we find that non-minimally coupled models could be distinguished from minimally coupled ones.

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

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

Accuracy of photometric redshifts for future weak lensing surveys from space

Photometric redshifts are a key tool to extract as much information as possible from planned cosmic shear experiments. In this work we aim to test the performances that can be achieved with observations in the near-infrared from space and in the optical from the ground. This is done by performing realistic simulations of multi-band observations of a patch of the sky, and submitting these mock images to software usually applied to real images to extract the photometry and then a redshift estimate for each galaxy. In this way we mimic the most relevant sources of uncertainty present in real data analysis, including blending and light pollution between galaxies. As an example we adopt the infrared setup of the ESA-proposed Euclid mission, while we simulate different observations in the optical, modifying filters, exposure times and seeing values. Finally, we consider directly some future ground-based experiments, such as LSST, Pan-Starrs and DES. The results highlight the importance of u-band observations, especially to discriminate between low (z < 0.5) and high (z ~ 3) redshifts, and the need for good observing sites, with seeing FWHM < 1. arcsec. The former of these indications clearly favours the LSST experiment as a counterpart for space observations, while for the other experiments we need to exclude at least 15 % of the galaxies to reach a precision in the photo-zs equal to $<\frac{σ_z}{1+z}>$ < 0.05.

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 in 2D: the concentration-mass relation for galaxy clusters

The aim of this work is to perform a systematic study of the measures of the mass and concentration estimated by fitting the convergence profile of a large sample of mock galaxy cluster size lenses, created with the publicly available code MOKA. We found that the main contribution to the bias in mass and in concentration is due to the halo triaxiality and second to the presence of substructures within the host halo virial radius. We show that knowing the cluster elongation along the line of sight helps in correcting the mass bias, but still keeps a small negative bias for the concentration. If these mass and concentration biases will characterize the galaxy cluster sample of a wide field survey it will be difficult to well recover within one sigma the cosmological parameters that mainly influence the c - M relation, using as reference a 3D c - M relation measured in cosmological N-body simulation. In this work we propose how to correct the c - M relation for projection effects and for adiabatic contraction and suggest to use these as reference for real observed data. Correcting mass and concentration estimates, as we propose, gives a measurement of the cosmological parameter within 1 - σ confidence contours.

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

Counts of high-redshift GRBs as probe of primordial non-Gaussianities

We propose to use high-redshift long $γ$-ray bursts (GRBs) as cosmological tools to constrain the amount of primordial non-Gaussianity in the density field. By using numerical, N-body, hydrodynamic, chemistry simulations of different cosmological volumes with various Gaussian and non-Gaussian models, we self-consistently relate the cosmic star formation rate density to the corresponding GRB rate. Assuming that GRBs are fair tracers of cosmic star formation, we find that positive local non-Gaussianities, described in terms of the non-linear parameter, \fnl, might boost significantly the GRB rate at high redshift, $z \gg 6$. Deviations with respect to the Gaussian case account for a few orders of magnitude if \fnl$\sim 1000$, one order of magnitude for \fnl$\sim 100$, and a factor of $\sim 2$ for \fnl$\sim 50$. These differences are found only at large redshift, while at later times the rates tend to converge. Furthermore, a comparison between our predictions and the observed GRB data at $z > 6$ allows to exclude large negative \fnl, consistently with previous works. Future detections of any long GRB at extremely high redshift ($z\sim 15-20$) could favor non-Gaussian scenarios with positive \fnl. More stringent constraints require much larger high-$z$ GRB complete samples, currently not available in literature. By distinguishing the contributions to the GRB rate from the metal-poor population III regime, and the metal-enriched population II-I regime, we conclude that the latter is a more solid tracer of the underlying matter distribution, while the former is strongly dominated by feedback mechanisms from the first, massive, short-lived stars, rather than by possible non-Gaussian fluctuations. This holds quite independently of the assumed population III initial mass function.

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

Hydrodynamical simulations of galaxy clusters in dark energy cosmologies: I. general properties

We investigate the influence of dark energy on structure formation, within five different cosmological models, namely a concordance $Λ$CDM model, two models with dynamical dark energy, viewed as a quintessence scalar field (using a RP and a SUGRA potential form) and two extended quintessence models (EQp and EQn) where the quintessence scalar field interacts non-minimally with gravity (scalar-tensor theories). We adopted for all models the normalization of the matter power spectrum $σ_{8}$ to match the CMB data. In the models with dynamical dark energy and quintessence, we describe the equation of state with $w_0\approx-0.9$, still within the range allowed by observations. For each model, we have performed hydrodynamical simulations in a cosmological box of $(300 \ {\rm{Mpc}} \ h^{-1})^{3}$ including baryons and allowing for cooling and star formation. The contemporary presence of evolving dark energy and baryon physics allows us to investigate the interplay between the different background cosmology and the evolution of the luminous matter. Since cluster baryon fraction can be used to constrain other cosmological parameters such as $Ω_{m}$, we also analyse how dark energy influences the baryon content of galaxy clusters. We find that, in models with dynamical dark energy, the evolving cosmological background leads to different star formation rates and different formation histories of galaxy clusters, but the baryon physics is not affected in a relevant way. We investigate several proxies of the cluster mass function based on X-ray observables like temperature, luminosity, $M_{gas}$, and $Y_{gas}$. We conclude that the X-ray temperature and $M_{gas}$ functions are better diagnostic to disentangle the growth of structures among different dark energy models. [Abridged]

preprint2012arXiv

Hydrodynamical simulations of galaxy clusters in dark energy cosmologies: II. c-M relation

We use hydrodynamical simulations of different dark energy cosmologies to investigate the concentration-mass (c-M) relation in galaxy clusters. In particular, we consider a reference ΛCDM model, two quintessence models with inverse power-law potentials (RP and SUGRA), and two extended quintessence models, one with positive and one with negative coupling with gravity (EQp and EQn respectively). All the models are normalized in order to match CMB data from WMAP3. We fit both the dark matter only and the total mass profile with a NFW profile, and recover the concentration of each halo from the fit using different definition. We consider both the complete catalog of clusters and groups and subsamples of objects at different level of relaxation. We find that the definition itself of the concentration can lead to differences up to 20% in its value and that these differences are smaller when more relaxed objects are considered. The c-M relation of our reference ΛCDM model is in good agreement with the results in literature, and relaxed objects have a higher normalization and a shallower slope with respect to the complete sample. The inclusion of baryon physics is found to influence more high-mass systems than low-mass ones, due to a higher concentration of baryons in the inner regions of massive halos. For the different dark energy models, we find that for ΛCDM, RP and SUGRA the normalization of the c-M relation is linked to the growth factor, with models having a higher value of σ8 D+ having also a higher normalization. This simple scheme is no longer valid for EQp and EQn because in these models a time dependent effective gravitational interaction, whose redshift evolution depends on the sign of the coupling, is present. This leads to a decrease (increase) of the expected normalization in the EQp (EQn) model. [Abridged]

preprint2012arXiv

Inclusive Constraints on Unified Dark Matter Models from Future Large-Scale Surveys

In the very last years, cosmological models where the properties of the dark components of the Universe - dark matter and dark energy - are accounted for by a single "dark fluid" have drawn increasing attention and interest. Amongst many proposals, Unified Dark Matter (UDM) cosmologies are promising candidates as effective theories. In these models, a scalar field with a non-canonical kinetic term in its Lagrangian mimics both the accelerated expansion of the Universe at late times and the clustering properties of the large-scale structure of the cosmos. However, UDM models also present peculiar behaviours, the most interesting one being the fact that the perturbations in the dark-matter component of the scalar field do have a non-negligible speed of sound. This gives rise to an effective Jeans scale for the Newtonian potential, below which the dark fluid does not cluster any more. This implies a growth of structures fairly different from that of the concordance LCDM model. In this paper, we demonstrate that forthcoming large-scale surveys will be able to discriminate between viable UDM models and LCDM to a good degree of accuracy. To this purpose, the planned Euclid satellite will be a powerful tool, since it will provide very accurate data on galaxy clustering and the weak lensing effect of cosmic shear. Finally, we also exploit the constraining power of the ongoing CMB Planck experiment. Although our approach is the most conservative, with the inclusion of only well-understood, linear dynamics, in the end we also show what could be done if some amount of non-linear information were included.

preprint2012arXiv

Measuring the neutrino mass from future wide galaxy cluster catalogues

We present forecast errors on a wide range of cosmological parameters obtained from a photometric cluster catalogue of a future wide-field Euclid-like survey. We focus in particular on the total neutrino mass as constrained by a combination of the galaxy cluster number counts and correlation function. For the latter we consider only the shape information and the Baryon Acoustic Oscillations (BAO), while marginalising over the spectral amplitude and the redshift space distortions. In addition to the cosmological parameters of the standard LCDM+nu model we also consider a non-vanishing curvature, and two parameters describing a redshift evolution for the dark energy equation of state. For completeness, we also marginalise over a set of "nuisance" parameters, representing the uncertainties on the cluster mass determination. We find that combining cluster counts with power spectrum information greatly improves the constraining power of each probe taken individually, with errors on cosmological parameters being reduced by up to an order of magnitude. In particular, the best improvements are for the parameters defining the dynamical evolution of dark energy, where cluster counts break degeneracies. Moreover, the resulting error on neutrino mass is at the level of σ(M_ν)\sim 0.9 eV, comparable with that derived from present Ly-alpha forest measurements and Cosmic Microwave background (CMB) data in the framework of a non-flat Universe. Further adopting Planck priors and reducing the number of free parameters to a LCDM+nu cosmology allows to place constraints on the total neutrino mass of σ(M_ν) \sim 0.08 eV, close to the lower bound enforced by neutrino oscillation experiments. [abridged]

preprint2012arXiv

MOKA: a new tool for Strong Lensing Studies

Strong gravitational lensing is a powerful tool for probing the matter distribution in the cores of massive dark matter haloes. Recent and ongoing analyses of galaxy cluster surveys (MACS, CFHTLS, SDSS, SGAS, CLASH, LoCuSS) have adressed the question of the nature of the dark matter distribution in clusters. N-body simulations of cold dark-matter haloes consistently find that haloes should be characterized by a concentration-mass relation that decreases monotonically with halo mass, and populated by a large amount of substructures, representing the cores of accreted progenitor halos. It is important for our understanding of dark matter to test these predictions. We present MOKA, a new algorithm for simulating the gravitational lensing signal from cluster-sized haloes. It implements the most recent results from numerical simulations to create realistic cluster-scale lenses with properties independent of numerical resolution. We perform systematic studies of the strong lensing cross section as a function of halo structures. We find that the strong lensing cross sections depend most strongly on the concentration and on the inner slope of the density profile of a halo, followed in order of importance by halo triaxiality and the presence of a bright central galaxy.

preprint2012arXiv

On the modelling of the excesses of galaxy clusters over high-mass thresholds

In this work we present for the first time an application of the Pareto approach to the modelling of the excesses of galaxy clusters over high-mass thresholds. The distribution of those excesses can be described by the generalized Pareto distribution (GPD), which is closely related to the generalized extreme value (GEV) distribution. After introducing the formalism, we study the impact of different thresholds and redshift ranges on the distributions, as well as the influence of the survey area on the mean excess above a given mass threshold. We also show that both the GPD and the GEV approach lead to identical results for rare, thus high-mass and high-redshift, clusters. As an example, we apply the Pareto approach to ACT-CL J0102-4915 and SPT-CL J2106-5844 and derive the respective cumulative distribution functions of the exceedance over different mass thresholds. We also study the possibility to use the GPD as a cosmological probe. Since in the maximum likelihood estimation of the distribution parameters all the information from clusters above the mass threshold is used, the GPD might offer an interesting alternative to GEV-based methods that use only the maxima in patches. When comparing the accuracy with which the parameters can be estimated, it turns out that the patch-based modelling of maxima is superior to the Pareto approach. In an ideal case, the GEV approach is capable to estimate the location parameter with a percent level precision for less than 100 patches. This result makes the GEV based approach potentially also interesting for cluster surveys with a smaller area.

preprint2012arXiv

Predicting the number of giant arcs expected in the next generation wide-field surveys from space

In this paper we estimate the number of gravitational arcs detectable in a wide-field survey such as that which will be operated by the Euclid space mission, assuming a ΛCDM cosmology. We use the publicly available code MOKA to obtain realistic deflection angle maps of mock gravitational lenses. The maps are processed by a ray-tracing code to estimate the strong lensing cross sections of each lens. Our procedure involves 1) the generation of a light-cone which is populated with lenses drawn from a theoretical mass-function; 2) the modeling of each single lens using a triaxial halo with a NFW (Navarro-Frenk-White) density profile and theoretical concentration-mass relation, including substructures, 3) the determination of the lensing cross section as a function of redshift for each lens in the light-cone, 4) the simulation of mock observations to characterize the redshift distribution of sources that will be detectable in the Euclid images. We focus on the so-called giant arcs, i.e. gravitational arcs characterized by large length-to-width ratios (l/w > 5, 7.5 and 10). We quantify the arc detectability at different significances above the level of the background. Performing 128 different realizations of a 15,000 sq. degree survey, we find that the number of arcs detectable at 1σ above the local background will be 8912,2914, and 1275 for l/w>5, 7.5 and 10, respectively. The expected arc numbers decrease to 2409, 790, and 346 for a detection limit at 3σ above the background level. From our analysis, we find that most of the lenses which contribute to the lensing optical depth are located at redshifts 0.4<zl<0.7 and that the 50% of the arcs are images of sources at zs > 3. This is the first step towards the full characterization of the population of strong lenses that will be observed by Euclid. [abridged]

preprint2011arXiv

An application of extreme value statistics to the most massive galaxy clusters at low and high redshifts

In this work we present an application of general extreme value statistics (GEV) to very massive single clusters at high and low redshifts. After introducing the formalism, we apply this statistics to four very massive high redshift clusters. Those clusters comprise ACT-CL J0102-4915 with a mass of M_200m=(2.16+/-0.32)x10^{15} M_sun at a redshift of z=0.87, SPT-CL J2106-5844 with a mass of M_200m=(1.27+/-0.21)x10^{15} M_sun at z=1.132 and two clusters found by the XMM-Newton Distant Cluster Project survey: XMMU J2235.32557 with a mass of M_200c= (7.3+/-1.3)x10^{14} M_sun located at a redshift of z=1.4 and XMMU J0044.0-2033 having a mass in the range of M_200c= (3.5-5.0)x10^{14} M_sun at z=1.579. By relating those systems to their corresponding distribution functions of being the most massive system in a given survey area, we find that none of the systems alone is in extreme tension with LCDM. We confront these results with a GEV analysis of four very massive low redshift clusters: A2163, A370, RXJ1347-1145 and 1E0657-558, finding no tendency of the high-z systems to be more extreme than the low-z ones. In addition, we study the extreme quantiles of single clusters at high-z and present contour plots for fixed quantiles in the mass vs. survey area plane for four redshift intervals, finding that, in order to be significantly in conflict with LCDM, cluster masses would have to be substantially higher than the currently observed ones.

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.

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.

preprint2011arXiv

The clustering of galaxies and galaxy clusters: constraints on primordial non-Gaussianity from future wide-field surveys

We investigate the constraints on primordial non-Gaussianity with varied bispectrum shapes that can be derived from the power spectrum of galaxies and clusters of galaxies detected in future wide field optical/near-infrared surveys. Having in mind the proposed ESA space mission \emph{Euclid} as a specific example, we combine the spatial distribution of spectroscopically selected galaxies with that of weak lensing selected clusters. We use the physically motivated halo model in order to represent the correlation function of arbitrary tracers of the Large Scale Structure in the Universe. As naively expected, we find that galaxies are much more effective in jointly constrain the level of primordial non-Gaussianity $f_\mathrm{NL}$ and the amplitude of the matter power spectrum $σ_8$ than clusters of galaxies, due to the much lower abundance of the latter that is not adequately compensated by the larger effect on the power spectrum. Nevertheless, combination of the galaxy power spectrum with the cluster-galaxy cross spectrum can decrease the error on the determination of $f_\mathrm{NL}$ by up to a factor of $\sim 2$. This decrement is particularly evident for the less studied non-Gaussian bispectrum shapes, the so-called enfolded and the orthogonal ones. Setting constraints on these models can shed new light on various aspects of the physics of the early Universe, and it is hence of extreme importance. By combining the power spectra of clusters and galaxies with the cluster-galaxy cross spectrum we find constraints on primordial non-Gaussianity of the order $Δf_\mathrm{NL} \sim $ a few, competitive and possibly superior to future CMB experiments.

preprint2010arXiv

Non-Gaussian halo assembly bias

The strong dependence of the large-scale dark matter halo bias on the (local) non-Gaussianity parameter, f_NL, offers a promising avenue towards constraining primordial non-Gaussianity with large-scale structure surveys. In this paper, we present the first detection of the dependence of the non-Gaussian halo bias on halo formation history using N-body simulations. We also present an analytic derivation of the expected signal based on the extended Press-Schechter formalism. In excellent agreement with our analytic prediction, we find that the halo formation history-dependent contribution to the non-Gaussian halo bias (which we call non-Gaussian halo assembly bias) can be factorized in a form approximately independent of redshift and halo mass. The correction to the non-Gaussian halo bias due to the halo formation history can be as large as 100%, with a suppression of the signal for recently formed halos and enhancement for old halos. This could in principle be a problem for realistic galaxy surveys if observational selection effects were to pick galaxies occupying only recently formed halos. Current semi-analytic galaxy formation models, for example, imply an enhancement in the expected signal of ~23% and ~48% for galaxies at z=1 selected by stellar mass and star formation rate, respectively.

preprint2010arXiv

Optimal filtering of optical and weak lensing data to search for galaxy clusters: application to the COSMOS field

Galaxy clusters are usually detected in blind optical surveys via suitable filtering methods. We present an optimal matched filter which maximizes their signal-to-noise ratio by taking advantage of the knowledge we have of their intrinsic physical properties and of the data noise properties. In this paper we restrict our application to galaxy magnitudes, positions and photometric redshifts if available, and we also apply the filter separately to weak lensing data. The method is suitable to be naturally extended to a multi-band approach which could include not only additional optical bands but also observables with different nature such as X-rays. For each detection, the filter provides its significance, an estimate for the richness and for the redshift even if photo-z are not given. The provided analytical error estimate is tested against numerical simulations. We finally apply our method to the COSMOS field and compare the results with previous cluster detections obtained with different methods. Our catalogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including both optical and weak-lensing information.

preprint2010arXiv

The distribution of Lyman-limit absorption systems during and after reionization

Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of reionization can drive the large-scale distribution of absorption systems. Specifically, we find that absorbers are highly correlated with the density field on small scales, and then become anti-correlated with the UVB on large scales. After reionization, the large-scale power spectrum of the absorbers traces the UVB power spectrum, which can be predicted with a simple analytic extension of the halo-model. During reionization, absorbers tend to preferentially lie inside overdensities (i.e. filaments) of the recently-ionized intergalactic medium (IGM). Absorbers may also dominate the small-scale (k > 1/Mpc) 21-cm power during and after reionization. Conversely, they smooth the contrast on moderate scales. Once the HII regions grow to surpass the mfp, the absorbers add to the large-scale 21-cm power. Our results should prove useful in interpreting future observations of the reionization epoch.

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

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.

preprint2009arXiv

X-ray and strong lensing mass estimate of MS2137.3-2353

We present new mass estimates of the galaxy cluster MS2137.3-2353, inferred from X-ray and strong lensing analyses. This cluster exhibits an outstanding strong lensing configuration and indicates a well-relaxed dynamical state, being most suitable for a mass reconstruction which combines both techniques. Despite this, several previous studies have claimed a significant discrepancy between the X-ray and the strong lensing mass estimates. The primary aim of this paper is to address and explain this mismatch. For this purpose, we have analysed Chandra observations to recover the profiles of the intra-cluster medium properties and, assuming a functional form for the matter density, the total mass distribution. The notable strong-lensing features of MS2137.3 allow us to reconstruct its projected mass in the central regions with good accuracy, by taking advantage of the lensing inversion code Lenstool. We find a good agreement between our X-ray and strong lensing mass estimates for MS2137, leading to a mean, extrapolated value of M_200 ~ 4.4 +/- 0.3*10^14 M_sun, under the assumption of the Navarro-Frenk-White (NFW) density profile. However, the strong lensing mass estimate is affected by the details of the BCG mass modeling. In particular, we do not find evidence for a high concentration for the NFW density profile, as reported in some earlier works.

preprint2008arXiv

The Effect of Primordial Non--Gaussianity on the Topology of Large-Scale Structure

We study the effect of primordial non-Gaussianity on the development of large-scale cosmic structure using high-resolution N-body simulations. In particular, we focus on the topological properties of the "cosmic web", quantitatively characterized by the Minkowski Functionals, for models with quadratic non-linearities with different values of the usual non-Gaussianity parameter fNL. In the weakly non-linear regime, we find that analytic formulae derived from perturbation theory agree with the numerical results within a few percent of the amplitude of each MF when |fNL|<1000. In the non-linear regime, the detailed behavior of the MFs as functions of threshold density deviates more strongly from the analytical curves, while the overall amplitude of the primordial non-Gaussian effect remains comparable to the perturbative prediction. When smaller-scale information is included, the influence of primordial non-Gaussianity becomes increasingly significant statistically due to decreasing sample variance. We find that the effect of the primordial non-Gaussianity with |fNL|=50 is comparable to the sample variance of mass density fields with a volume of 0.125(Gpc/h)^3 when they are smoothed by Gaussian filter at a scale of 5Mpc/h. The detectability of this effect in actual galaxy surveys will strongly depend upon residual uncertainties in cosmological parameters and galaxy biasing.

preprint2003arXiv

Numerical study of halo concentrations in dark-energy cosmologies

We study the concentration parameters, their mass dependence and redshift evolution, of dark-matter halos in different dark-energy cosmologies with constant and time-variable equation of state, and compare them with "standard'' Lambda-CDM and OCDM models. We find that previously proposed algorithms for predicting halo concentrations can be well adapted to dark-energy models. When centred on the analytically expected values, halo concentrations show a log-normal distribution with a uniform standard deviation of ~0.2. The dependence of averaged halo concentrations on mass and redshift permits a simple fit of the form (1+z) c=c0 (M/M0)^a, with a~-0.1 throughout. We find that the cluster concentration depends on the dark energy equation of state at the cluster formation redshift z_{coll} through the linear growth factor D_+(z_{coll}). As a simple correction accounting for dark-energy cosmologies, we propose scaling c0 from Lambda-CDM with the ratio of linear growth factors, c0 -> c0 D_+(z_{coll})/D_{+,Lambda-CDM}(z_{coll}).

preprint1998arXiv

Large-scale bias in the Universe II: redshift space bispectrum

The determination of the density parameter $Ω_0$ from the large-scale distribution of galaxies is one of the major goals of modern cosmology. However, if galaxies are biased tracers of the underlying mass distribution, linear perturbation theory leads to a degeneracy between $Ω_0$ and the linear bias parameter $b$, and the density parameter cannot be estimated. In Matarrese, Verde & Heavens (1997) we developed a method based on second-order perturbation theory to use the bispectrum to lift this degeneracy by measuring the bias parameter in an $Ω_0$-independent way. The formalism was developed assuming that one has perfect information on the positions of galaxies in three dimensions. In galaxy redshift surveys, the three-dimensional information is imperfect, because of the contaminating effects of peculiar velocities, and the resulting clustering pattern in redshift space is distorted. In this paper, we combine second-order perturbation theory with a model for collapsed, virialised structures, to extend the method to redshift space, and demonstrate that the method should be successful in determining with reasonable accuracy the bias parameter from state-of-the-art surveys such as the Anglo-Australian 2 degree field survey and the Sloan digital sky survey.

preprint1996arXiv

Redshift evolution of clustering

We discuss how the redshift dependence of the observed two-point correlation function of various classes of objects can be related to theoretical predictions. This relation involves first a calculation of the redshift evolution of the underlying matter correlations. The next step is to relate fluctuations in mass to those of any particular class of cosmic objects; in general terms, this means a model for the bias and how it evolves with cosmic epoch. Only after these two effects have been quantified can one perform an appropriate convolution of the non-linearly evolved two-point correlation function of the objects with their redshift distribution to obtain the `observed' correlation function for a given sample. This convolution in itself tends to mask the effect of evolution by mixing amplitudes at different redshifts. We develop a formalism which incorporates these requirements and, in particular, a set of plausible models for the evolution of the bias factor. We apply this formalism to the spatial, angular and projected correlation functions from different samples of high-redshift objects, assuming a simple phenomenological model for the initial power-spectrum and an Einstein-de Sitter cosmological model. We find that our model is roughly consistent with data on the evolution of QSO and galaxy clustering, but only if the effective degree of biasing is small. We discuss the differences between our analysis and other theoretical studies of clustering evolution and argue that the dominant barrier to making definitive predictions is uncertainty about the appropriate form of the bias and its evolution with cosmic epoch.

preprint1996arXiv

The epoch of structure formation in blue mixed dark matter models

Recent data on the high--redshift abundance of damped Ly$α$ systems are compared with theoretical predictions for `blue' (i.e. $n>1$) Mixed Dark Matter models. The results show that decreasing the hot component fraction $Ω_ν$ and/or increasing the primordial spectral index $n$ implies an earlier epoch of cosmic structure formation. However, we also show that varying $Ω_ν$ and $n$ in these directions makes the models barely consistent with the observed abundance of galaxy clusters. Therefore, requiring at the same time observational constraints on damped Ly$α$ systems and cluster abundance to be satisfied represents a challenge for the Mixed Dark Matter class of models.

preprint1995arXiv

COLLISIONAL VERSUS COLLISIONLESS MATTER: A ONE-DIMENSIONAL ANALYSIS OF GRAVITATIONAL CLUSTERING

We present the results of a series of one-dimensional N-body and hydrodynamical simulations which have been used for testing the different clustering properties of baryonic and dark matter in an expanding background. Initial Gaussian random density perturbations with a power-law spectrum $P(k) \propto k^n$ are assumed. We analyse the distribution of density fluctuations and thermodynamical quantities for different spectral indices $n$ and discuss the statistical properties of clustering in the corresponding simulations. At large scales the final distribution of the two components is very similar while at small scales the dark matter presents a lumpiness which is not found in the baryonic matter. The amplitude of density fluctuations in each component depends on the spectral index $n$ and only for $n=-1$ the amplitude of baryonic density fluctuations is larger than that in the dark component. This result is also confirmed by the behaviour of the bias factor, defined as the ratio between the r.m.s of baryonic and dark matter fluctuations at different scales: while for $n=1,\ 3$ it is always less than unity except at very large scales where it tends to one, for $n=-1$ it is above 1.4 at all scales. All simulations show also that there is not an exact correspondence between the positions of largest peaks in dark and baryonic components, as confirmed by a cross-correlation analysis. The final temperatures depend on the initial spectral index: the highest values are obtained for $n=-1$ and are in proximity of high density regions.

preprint1995arXiv

Eulerian Perturbation Theory in Non-Flat Universes: Second-Order Approximation

The problem of solving perturbatively the equations describing the evolution of self-gravitating collisionless matter in an expanding universe considerably simplifies when directly formulated in terms of the gravitational and velocity potentials: the problem can be solved {\it exactly}, rather than approximately, even for cosmological models with arbitrary density parameter $Ω$. The Eulerian approach we present here allows to calculate the higher-order moments of the initially Gaussian density and velocity fields: in particular, we compute the gravitationally induced skewness of the density and velocity-divergence fields for any value of $Ω$, confirming the extremely weak $Ω$-dependence of the skewness previously obtained via Lagrangian perturbation theory. Our results show that the separability assumption of higher-order Eulerian perturbative solutions is restricted to the Einstein-de Sitter case only.

preprint1995arXiv

The Cluster Distribution as a Test of Dark Matter Models. III. the Cluster Velocity Field

We study the large-scale velocity fields traced by galaxy clusters in numerical simulations of a box of side 960 $h^{-1}$ Mpc, and compare them with available data on real clusters. In order to test the reliability of the simulations, which are based on an optimized version of the Zel'dovich approximation, we compare their cluster velocities with those of "exact" N-body simulations, and find a remarkable agreement between the two according to a variety of statistical tests. We analyse Cold Dark Matter models with density parameter in the range $0.2\le Ω_0\le 1$, both with and without the cosmological constant term to provide a flat geometry. We also simulate a Cold+Hot Dark Matter model, with 30% provided by the hot component. Comparison with real data is performed by applying tests based on the cumulative velocity frequency distribution (CVFD) and bulk flow statistics. For the CVFD, we use observational velocity data from different authors. By merging all available data in a combined sample, the CVFD test is able to exclude only open models with $Ω_0\le 0.4$ and flat models with $Ω_0=0.2$. However, the analysis of individual observational samples gives contradictory results; in particular the recent Giovanelli data, which imply much lower velocities (with significantly smaller uncertainties) than the other data, are actually only consistent with the open $Ω_0=0.4$ and the flat $Ω_0=0.2$ models, i.e. the models which are excluded by the other data. Furthermore, CVFD and bulk flow analyses of the Branchini et al. reconstructed velocity data disfavour precisely those models accepted on the grounds of Giovanelli's sample. Finally, we confirm that the Lauer & Postman bulk flow determination is an extremely rare event in the cosmological models we have analysed.

preprint1993arXiv

Testing the Frozen-Flow Approximation

We investigate the accuracy of the frozen--flow approximation (FFA), recently proposed by Matarrese \etal (1992), for following the nonlinear evolution of cosmological density fluctuations under gravitational instability. We compare a number of statistics between results of the FFA and nbody simulations, including those used by Melott, Pellman \& Shandarin (1993) to test the Zel'dovich approximation. The FFA performs reasonably well in a statistical sense, e.g. in reproducing the counts--in--cell distribution, at small scales, but it does poorly in the crosscorrelation with nbody which means it is generally not moving mass to the right place, especially in models with high small--scale power.