Source author record

L. Moscardini

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

94works
4topics
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

94 published item(s)

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid's unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10's stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10's orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1).

preprint2023arXiv

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

preprint2022arXiv

AMICO galaxy clusters in KiDS-DR3: Constraints on cosmological parameters and on the normalisation of the mass-richness relation from clustering

We analysed the clustering of a photometric sample of galaxy clusters selected from the Third Data Release of the Kilo-Degree Survey, focusing on the redshift-space two-point correlation function (2PCF). We compared our measurements to theoretical predictions of the standard $Λ$ cold dark matter ($Λ$CDM) cosmological model. We measured the 2PCF of the sample in the cluster centric radial range $r\in[5,80]$ $h^{-1}$Mpc, considering 4934 galaxy clusters with richness $λ^*\geq15$ in the redshift range $z\in[0.1,0.6]$. A Markov chain Monte Carlo analysis has been performed to constrain the cosmological parameters $Ω_{\rm m}$, $σ_8$, and $S_8 \equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$, assuming Gaussian priors on the mass-richness relation given by the posteriors obtained from a joint analysis of cluster counts and weak lensing. In addition, we constrained the normalisation of the mass-richness relation, $α$, with fixed cosmological parameters. We obtained $Ω_{\rm m}=0.28^{+0.05}_{-0.04}$, $σ_8=0.82^{+0.14}_{-0.12}$, and $S_8=0.80^{+0.08}_{-0.08}$. The constraint on $S_8$ is consistent within 1$σ$ with the results from WMAP and Planck. Furthermore, by fixing the cosmological parameters to those provided by Planck, we obtained $α=0.12^{+0.06}_{-0.06}$, which is fully consistent with the result obtained from the joint analysis of cluster counts and weak lensing performed for this sample.

preprint2022arXiv

Euclid preparation: XIII. Forecasts for galaxy morphology with the Euclid Survey using Deep Generative Models

We present a machine learning framework to simulate realistic galaxies for the Euclid Survey. The proposed method combines a control on galaxy shape parameters offered by analytic models with realistic surface brightness distributions learned from real Hubble Space Telescope observations by deep generative models. We simulate a galaxy field of $0.4\,\rm{deg}^2$ as it will be seen by the Euclid visible imager VIS and show that galaxy structural parameters are recovered with similar accuracy as for pure analytic Sérsic profiles. Based on these simulations, we estimate that the Euclid Wide Survey will be able to resolve the internal morphological structure of galaxies down to a surface brightness of $22.5\,\rm{mag}\,\rm{arcsec}^{-2}$, and $24.9\,\rm{mag}\,\rm{arcsec}^{-2}$ for the Euclid Deep Survey. This corresponds to approximately $250$ million galaxies at the end of the mission and a $50\,\%$ complete sample for stellar masses above $10^{10.6}\,\rm{M}_\odot$ (resp. $10^{9.6}\,\rm{M}_\odot$) at a redshift $z\sim0.5$ for the wide (resp. deep) survey. The approach presented in this work can contribute to improving the preparation of future high-precision cosmological imaging surveys by allowing simulations to incorporate more realistic galaxies.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Euclid preparation. XVIII. The NISP photometric system

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μ$m range, to a 5$σ$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime.

preprint2022arXiv

Euclid: Constraining ensemble photometric redshift distributions with stacked spectroscopy

The ESA Euclid mission will produce photometric galaxy samples over 15000 square degrees of the sky that will be rich for clustering and weak lensing statistics. The accuracy of the cosmological constraints derived from these measurements will depend on the knowledge of the underlying redshift distributions based on photometric redshift calibrations. A new approach is proposed to use the stacked spectra from Euclid slitless spectroscopy to augment broad-band photometric information to constrain the redshift distribution with spectral energy distribution fitting. The high spectral resolution available in the stacked spectra complements the photometry and helps to break the colour-redshift degeneracy and constrain the redshift distribution of galaxy samples. We modelled the stacked spectra as a linear mixture of spectral templates. The mixture may be inverted to infer the underlying redshift distribution using constrained regression algorithms. We demonstrate the method on simulated Vera C. Rubin Observatory and Euclid mock survey data sets based on the Euclid Flagship mock galaxy catalogue. We assess the accuracy of the reconstruction by considering the inference of the baryon acoustic scale from angular two-point correlation function measurements. We selected mock photometric galaxy samples at redshift z>1 using the self-organising map algorithm. Considering the idealised case without dust attenuation, we find that the redshift distributions of these samples can be recovered with 0.5% accuracy on the baryon acoustic scale. The estimates are not significantly degraded by the spectroscopic measurement noise due to the large sample size. However, the error degrades to 2% when the dust attenuation model is left free. We find that the colour degeneracies introduced by attenuation limit the accuracy considering the wavelength coverage of Euclid near-infrared spectroscopy.

preprint2022arXiv

Euclid: Covariance of weak lensing pseudo-$C_\ell$ estimates. Calculation, comparison to simulations, and dependence on survey geometry

An accurate covariance matrix is essential for obtaining reliable cosmological results when using a Gaussian likelihood. In this paper we study the covariance of pseudo-$C_\ell$ estimates of tomographic cosmic shear power spectra. Using two existing publicly available codes in combination, we calculate the full covariance matrix, including mode-coupling contributions arising from both partial sky coverage and non-linear structure growth. For three different sky masks, we compare the theoretical covariance matrix to that estimated from publicly available N-body weak lensing simulations, finding good agreement. We find that as a more extreme sky cut is applied, a corresponding increase in both Gaussian off-diagonal covariance and non-Gaussian super-sample covariance is observed in both theory and simulations, in accordance with expectations. Studying the different contributions to the covariance in detail, we find that the Gaussian covariance dominates along the main diagonal and the closest off-diagonals, but further away from the main diagonal the super-sample covariance is dominant. Forming mock constraints in parameters describing matter clustering and dark energy, we find that neglecting non-Gaussian contributions to the covariance can lead to underestimating the true size of confidence regions by up to 70 per cent. The dominant non-Gaussian covariance component is the super-sample covariance, but neglecting the smaller connected non-Gaussian covariance can still lead to the underestimation of uncertainties by 10--20 per cent. A real cosmological analysis will require marginalisation over many nuisance parameters, which will decrease the relative importance of all cosmological contributions to the covariance, so these values should be taken as upper limits on the importance of each component.

preprint2022arXiv

Euclid: Forecast constraints on consistency tests of the $Λ$CDM model

The standard cosmological model is based on the fundamental assumptions of a spatially homogeneous and isotropic universe on large scales. An observational detection of a violation of these assumptions at any redshift would immediately indicate the presence of new physics. We quantify the ability of the Euclid mission, together with contemporary surveys, to improve the current sensitivity of null tests of the canonical cosmological constant $Λ$ and the cold dark matter (LCDM) model in the redshift range $0<z<1.8$. We considered both currently available data and simulated Euclid and external data products based on a LCDM fiducial model, an evolving dark energy model assuming the Chevallier-Polarski-Linder (CPL) parameterization or an inhomogeneous Lema\^ıtre-Tolman-Bondi model with a cosmological constant $Λ$ (LLTB), and carried out two separate but complementary analyses: a machine learning reconstruction of the null tests based on genetic algorithms, and a theory-agnostic parametric approach based on Taylor expansion and binning of the data, in order to avoid assumptions about any particular model. We find that in combination with external probes, Euclid can improve current constraints on null tests of the LCDM by approximately a factor of three when using the machine learning approach and by a further factor of two in the case of the parametric approach. However, we also find that in certain cases, the parametric approach may be biased against or missing some features of models far from LCDM. Our analysis highlights the importance of synergies between Euclid and other surveys. These synergies are crucial for providing tighter constraints over an extended redshift range for a plethora of different consistency tests of some of the main assumptions of the current cosmological paradigm.

preprint2022arXiv

KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using eight bands in the multipole range $76<\ell<1500$ for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-$C_{\ell}$ estimates by forward modelling the effects of the mask. Assuming a spatially flat $Λ$CDM cosmology, we constrained the structure growth parameter $S_8 = σ_8(Ω_{\rm m}/0.3)^{1/2} = 0.754_{-0.029}^{+0.027}$. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-$α$ forest and its cross-correlation with quasars, we tightened these constraints to $S_8 = 0.771^{+0.006}_{-0.032}$. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a $\sim 3σ$ tension with early-Universe constraints from cosmic microwave background experiments.

preprint2022arXiv

Testing General Relativity: new measurements of gravitational redshift in galaxy clusters

The peculiar velocity distribution of cluster member galaxies provides a powerful tool to directly investigate the gravitational potentials within galaxy clusters and to test the gravity theory on megaparsec scales. We exploit spectroscopic galaxy and galaxy cluster samples extracted from the latest releases of the Sloan Digital Sky Survey (SDSS) to derive new constraints on the gravity theory. We consider a spectroscopic sample of $3058$ galaxy clusters, with a maximum redshift of $0.5$. We analyse the velocity distribution of the cluster member galaxies to make new measurements of the gravitational redshift effect inside galaxy clusters. We accurately estimate the cluster centres, computing them as the average of angular positions and redshifts of the closest galaxies to the brightest cluster galaxies. We find that this centre definition provides a better estimation of the centre of the cluster gravitational potential wells, relative to simply assuming the brightest cluster galaxies as the cluster centres, as done in the past literature. We compare our measurements with the theoretical predictions of three different gravity theories: general relativity (GR), the $f(R)$ model, and the Dvali-Gabadadze-Porrati (DGP) model. A new statistical procedure is used to fit the measured gravitational redshift signal and thus to discriminate among the considered gravity theories. Finally, we investigate the systematic uncertainties possibly affecting the analysis. We clearly detect the gravitational redshift effect in the exploited cluster member catalogue. We recover an integrated gravitational redshift signal of $-11.4 \pm 3.3$ km s$^{-1}$, which is in agreement, within the errors, with past literature works. Overall, our results are consistent with both GR and DGP predictions, while they are in marginal disagreement with the predictions of the considered $f(R)$ strong field model.

preprint2021arXiv

$Euclid$ preparation: XV. Forecasting cosmological constraints for the $Euclid$ and CMB joint analysis

The combination and cross-correlation of the upcoming $Euclid$ data with cosmic microwave background (CMB) measurements is a source of great expectation since it will provide the largest lever arm of epochs, ranging from recombination to structure formation across the entire past light cone. In this work, we present forecasts for the joint analysis of $Euclid$ and CMB data on the cosmological parameters of the standard cosmological model and some of its extensions. This work expands and complements the recently published forecasts based on $Euclid$-specific probes, namely galaxy clustering, weak lensing, and their cross-correlation. With some assumptions on the specifications of current and future CMB experiments, the predicted constraints are obtained from both a standard Fisher formalism and a posterior-fitting approach based on actual CMB data. Compared to a $Euclid$-only analysis, the addition of CMB data leads to a substantial impact on constraints for all cosmological parameters of the standard $Λ$-cold-dark-matter model, with improvements reaching up to a factor of ten. For the parameters of extended models, which include a redshift-dependent dark energy equation of state, non-zero curvature, and a phenomenological modification of gravity, improvements can be of the order of two to three, reaching higher than ten in some cases. The results highlight the crucial importance for cosmological constraints of the combination and cross-correlation of $Euclid$ probes with CMB data.

preprint2021arXiv

AMICO galaxy clusters in KiDS-DR3: cosmological constraints from counts and stacked weak-lensing

We present a cosmological analysis of abundances and stacked weak-lensing profiles of galaxy clusters, exploiting the AMICO KiDS-DR3 catalogue. The sample consists of 3652 galaxy clusters with intrinsic richness $λ^*\geq20$, over an effective area of 377 deg$^2$, in the redshift range $z\in[0.1,\,0.6]$. We quantified the purity and completeness of the sample through simulations. The statistical analysis has been performed by simultaneously modelling the comoving number density of galaxy clusters and the scaling relation between the intrinsic richnesses and the cluster masses, assessed through a stacked weak-lensing profile modelling. The fluctuations of the matter background density, caused by super-survey modes, have been taken into account in the likelihood. Assuming a flat $Λ$CDM model, we constrained $Ω_{\rm m}$, $σ_8$, $S_8 \equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$, and the parameters of the mass-richness scaling relation. We obtained $Ω_{\rm m}=0.24^{+0.03}_{-0.04}$, $σ_8=0.86^{+0.07}_{-0.07}$, $S_8=0.78^{+0.04}_{-0.04}$. The constraint on $S_8$ is consistent within 1$σ$ with the results from WMAP and Planck. Furthermore, we got constraints on the cluster mass scaling relation in agreement with those obtained from a previous weak-lensing only analysis.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

Euclid preparation: XVI. Exploring the ultra low-surface brightness Universe with Euclid/VIS

While Euclid is an ESA mission specifically designed to investigate the nature of Dark Energy and Dark Matter, the planned unprecedented combination of survey area ($\sim15\,000$ deg$^2$), spatial resolution, low sky-background, and depth also make Euclid an excellent space observatory for the study of the low surface brightness Universe. Scientific exploitation of the extended low surface brightness structures requires dedicated calibration procedures yet to be tested. We investigate the capabilities of Euclid to detect extended low surface brightness structure by identifying and quantifying sky background sources and stray-light contamination. We test the feasibility of generating sky flat-fields to reduce large-scale residual gradients in order to reveal the extended emission of galaxies observed in the Euclid Survey. We simulate a realistic set of Euclid/VIS observations, taking into account both instrumental and astronomical sources of contamination, including cosmic rays, stray-light, zodiacal light, ISM, and the CIB, while simulating the effects of the presence of background sources in the FOV. We demonstrate that a combination of calibration lamps, sky flats and self-calibration would enable recovery of emission at a limiting surface brightness magnitude of $μ=29.5^{+0.08}_{-0.27} $ mag arcsec$^{-2}$ ($3σ$, $10\times10$ arcsec$^2$) in the Wide Survey, reaching regions 2 magnitudes deeper in the Deep Surveys. Euclid/VIS has the potential to be an excellent low surface brightness observatory. Covering the gap between pixel-to-pixel calibration lamp flats and self-calibration observations for large scales, the application of sky flat-fielding will enhance the sensitivity of the VIS detector at scales of larger than 1 degree, up to the size of the FOV, enabling Euclid to detect extended surface brightness structures below $μ=31$ mag arcsec$^{-2}$ and beyond.

preprint2021arXiv

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2021arXiv

Euclid: Forecast constraints on the cosmic distance duality relation with complementary external probes

In metric theories of gravity with photon number conservation, the luminosity and angular diameter distances are related via the Etherington relation, also known as the distance-duality relation (DDR). A violation of this relation would rule out the standard cosmological paradigm and point at the presence of new physics. We quantify the ability of Euclid, in combination with contemporary surveys, to improve the current constraints on deviations from the DDR in the redshift range $0<z<1.6$. We start by an analysis of the latest available data, improving previously reported constraints by a factor of 2.5. We then present a detailed analysis of simulated Euclid and external data products, using both standard parametric methods (relying on phenomenological descriptions of possible DDR violations) and a machine learning reconstruction using Genetic Algorithms. We find that for parametric methods Euclid can (in combination with external probes) improve current constraints by approximately a factor of six, while for non-parametric methods Euclid can improve current constraints by a factor of three. Our results highlight the importance of surveys like Euclid in accurately testing the pillars of the current cosmological paradigm and constraining physics beyond the standard cosmological model.

preprint2021arXiv

Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids

Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmology. In particular, we study the imprints of dynamic and geometric distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11000 voids and estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on $f/b$ and $D_M H$, where $f$ is the linear growth rate of density fluctuations, $b$ the galaxy bias, $D_M$ the comoving angular diameter distance, and $H$ the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of $f/b$ and 0.5% on $D_M H$ in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter $w$ for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts.

preprint2020arXiv

Euclid preparation: VIII. The Complete Calibration of the Colour-Redshift Relation survey: VLT/KMOS observations and data release

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

preprint2020arXiv

Euclid: The reduced shear approximation and magnification bias for Stage IV cosmic shear experiments

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

preprint2020arXiv

Nature versus nurture: relic nature and environment of the most massive passive galaxies at $z < 0.5$

Relic galaxies are thought to be the progenitors of high-redshift red nuggets that for some reason missed the channels of size growth and evolved passively and undisturbed since the first star formation burst (at $z>2$). These local ultracompact old galaxies are unique laboratories for studying the star formation processes at high redshift and thus the early stage of galaxy formation scenarios. Counterintuitively, theoretical and observational studies indicate that relics are more common in denser environments, where merging events predominate. To verify this scenario, we compared the number counts of a sample of ultracompact massive galaxies (UCMGs) selected within the third data release of the Kilo Degree Survey, that is, systems with sizes $R_{\rm e} < 1.5 \, \rm kpc$ and stellar masses $M_{\rm \star} > 8 \times 10^{10}\, \rm M_{\odot}$, with the number counts of galaxies with the same masses but normal sizes in field and cluster environments. Based on their optical and near-infrared colors, these UCMGs are likely to be mainly old, and hence representative of the relic population. We find that both UCMGs and normal-size galaxies are more abundant in clusters and their relative fraction depends only mildly on the global environment, with denser environments penalizing the survival of relics. Hence, UCMGs (and likely relics overall) are not special because of the environment effect on their nurture, but rather they are just a product of the stochasticity of the merging processes regardless of the global environment in which they live.

preprint2019arXiv

Euclid preparation: VI. Verifying the Performance of Cosmic Shear Experiments

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

preprint2016arXiv

Clustering-based redshift estimation: application to VIPERS/CFHTLS

We explore the accuracy of the clustering-based redshift estimation proposed by Ménard et al. (2013) when applied to VIPERS and CFHTLS real data. This method enables us to reconstruct redshift distributions from measurement of the angular clus- tering of objects using a set of secure spectroscopic redshifts. We use state of the art spectroscopic measurements with iAB < 22.5 from the VIMOS Public Extragalactic Redshift Survey (VIPERS) as reference population to infer the redshift distribution of galaxies from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) T0007 release. VIPERS provides a nearly representative sample to the flux limit iAB < 22.5 at redshift > 0.5 which allows us to test the accuracy of the clustering-based red- shift distributions. We show that this method enables us to reproduce the true mean color-redshift relation when both populations have the same magnitude limit. We also show that this technique allows the inference of redshift distributions for a population fainter than the one of reference and we give an estimate of the color-redshift mapping in this case. This last point is of great interest for future large redshift surveys which suffer from the need of a complete faint spectroscopic sample.

preprint2016arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS): Galaxy segregation inside filaments at $z \simeq 0.7$

We present the first quantitative detection of large-scale filamentary structure at $z \simeq 0.7$ in the large cosmological volume probed by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We use simulations to show the capability of VIPERS to recover robust topological features in the galaxy distribution, in particular the filamentary network. We then investigate how galaxies with different stellar masses and stellar activities are distributed around the filaments and find a significant segregation, with the most massive or quiescent galaxies being closer to the filament axis than less massive or active galaxies. The signal persists even after down-weighting the contribution of peak regions. Our results suggest that massive and quiescent galaxies assemble their stellar mass through successive mergers during their migration along filaments towards the nodes of the cosmic web. On the other hand, low-mass star-forming galaxies prefer the outer edge of filaments, a vorticity rich region dominated by smooth accretion, as predicted by the recent spin alignment theory. This emphasizes the role of large scale cosmic flows in shaping galaxy properties.

preprint2016arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Measuring nonlinear galaxy bias at z~0.8

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

preprint2016arXiv

Zoomed high-resolution simulations of Multi-coupled Dark Energy: cored galaxy density profiles at high redshift

We perform for the first time high-resolution zoom-in re-simulations of individual halos in the context of the Multi-coupled Dark Energy (McDE) scenario, which is characterised by the existence of two distinct dark matter particle species with opposite couplings to a Dark Energy scalar field. We compare the structural properties of the simulated halos to the standard Lambda-CDM results. The zoomed-in initial conditions are set up using a specifically designed code called ZInCo that we publicly release along with the present paper. Our numerical results allow to investigate in detail and with unprecedented resolution the halo segregation process that characterises McDE cosmologies from its very early stages. In particular, we find that in contrast to what could be inferred from previous numerical analysis at lower resolution, the segregation process is already in place at redshifts as high as z ~ 7. Most remarkably, we find that the subsequent evolution of the segregation leads to the formation of cored total matter density profiles with a core size that progressively increases in time. The shape of the cored profiles can be accurately predicted as the superposition of two NFW profiles with an increasing offset, thereby confirming the interpretation of the simulations results in terms of the segregation of the two dark matter components of the halo as a consequence of their different coupling to the Dark Energy field.

preprint2015arXiv

A weak lensing analysis of the PLCK G100.2-30.4 cluster

We present a mass estimate of the Planck-discovered cluster PLCK G100.2-30.4, derived from a weak lensing analysis of deep SUBARU griz images. We perform a careful selection of the background galaxies using the multi-band imaging data, and undertake the weak lensing analysis on the deep (1hr) r-band image. The shape measurement is based on the KSB algorithm; we adopt the PSFex software to model the Point Spread Function (PSF) across the field and correct for this in the shape measurement. The weak lensing analysis is validated through extensive image simulations. We compare the resulting weak lensing mass profile and total mass estimate to those obtained from our re-analysis of XMM-Newton observations, derived under the hypothesis of hydrostatic equilibrium. The total integrated mass profiles are in remarkably good agreement, agreeing within 1$σ$ across their common radial range. A mass $M_{500} \sim 7 \times 10^{14} M_\odot$ is derived for the cluster from our weak lensing analysis. Comparing this value to that obtained from our reanalysis of XMM-Newton data, we obtain a bias factor of (1-b) = 0.8 $\pm$ 0.1. This is compatible within 1$σ$ with the value of (1-b) obtained by Planck Collaboration XXIV from their calibration of the bias factor using newly-available weak lensing reconstructed masses.

preprint2015arXiv

Next Generation Cosmology: Constraints from the Euclid Galaxy Cluster Survey

We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to Δf_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter γ, which signals deviations from GR, will be constrained to Delta γ=0.02, and the neutrino density parameter to Delta Omega_ν=0.0013 (or Delta \sum m_ν=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]

preprint2015arXiv

The galaxy-halo connection from a joint lensing, clustering and abundance analysis in the CFHTLenS/VIPERS field

We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by $\sim60\,000$ secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy lensing and the stellar mass function. We interpret our measurements within the halo occupation distribution (HOD) framework, separating the contributions from central and satellite galaxies. We find that the SHMR for the central galaxies peaks at $M_{\rm h, peak} = 1.9^{+0.2}_{-0.1}\times10^{12} M_{\odot}$ with an amplitude of $0.025$, which decreases to $\sim0.001$ for massive halos ($M_{\rm h} > 10^{14} M_{\odot}$). Compared to central galaxies only, the total SHMR (including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size halos), a result consistent with cluster analyses from the literature based on fully independent methods. After properly accounting for differences in modelling, we have compared our results with a large number of results from the literature up to $z=1$: we find good general agreement, independently of the method used, within the typical stellar-mass systematic errors at low to intermediate mass (${M}_{\star} < 10^{11} M_{\odot}$) and the statistical errors above. We have also compared our SHMR results to semi-analytic simulations and found that the SHMR is tilted compared to our measurements in such a way that they over- (under-) predict star formation efficiency in central (satellite) galaxies.

preprint2015arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS): On the correct recovery of the count-in-cell probability distribution function

We compare three methods to measure the count-in-cell probability density function of galaxies in a spectroscopic redshift survey. From this comparison we found that when the sampling is low (the average number of object per cell is around unity) it is necessary to use a parametric method to model the galaxy distribution. We used a set of mock catalogues of VIPERS, in order to verify if we were able to reconstruct the cell-count probability distribution once the observational strategy is applied. We find that in the simulated catalogues, the probability distribution of galaxies is better represented by a Gamma expansion than a Skewed Log-Normal. Finally, we correct the cell-count probability distribution function from the angular selection effect of the VIMOS instrument and study the redshift and absolute magnitude dependency of the underlying galaxy density function in VIPERS from redshift $0.5$ to $1.1$. We found very weak evolution of the probability density distribution function and that it is well approximated, independently from the chosen tracers, by a Gamma distribution.

preprint2015arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field

Aims. Using the VIMOS Public Extragalactic Redshift Survey (VIPERS) we aim to jointly estimate the key parameters that describe the galaxy density field and its spatial correlations in redshift space. Methods. We use the Bayesian formalism to jointly reconstruct the redshift-space galaxy density field, power spectrum, galaxy bias and galaxy luminosity function given the observations and survey selection function. The high-dimensional posterior distribution is explored using the Wiener filter within a Gibbs sampler. We validate the analysis using simulated catalogues and apply it to VIPERS data taking into consideration the inhomogeneous selection function. Results. We present joint constraints on the anisotropic power spectrum as well as the bias and number density of red and blue galaxy classes in luminosity and redshift bins as well as the measurement covariances of these quantities. We find that the inferred galaxy bias and number density parameters are strongly correlated although these are only weakly correlated with the galaxy power spectrum. The power spectrum and redshift-space distortion parameters are in agreement with previous VIPERS results with the value of the growth rate $fσ_8 = 0.38$ with 18% uncertainty at redshift 0.7.

preprint2015arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Environmental effects shaping the galaxy stellar mass function

We exploit the first public data release of VIPERS to investigate environmental effects in galaxy evolution between $z\sim0.5$ and $0.9$. The large number of spectroscopic redshifts over an area of about $10\,\mathrm{deg}^2$ provides a galaxy sample with high statistical power. The accurate redshift measurements, with $σ_z = 0.00047(1+z_\mathrm{spec})$, allow us to robustly isolate galaxies living in the lowest- and highest-density environments, as defined in terms of spatial 3D density contrast. We estimate the stellar mass function (SMF) of galaxies residing in these two environments, and constrain its high-mass end with unprecedented precision. We find that the galaxy SMF in the densest regions has a different shape than that measured at low densities, with an enhancement of massive galaxies and a hint of a flatter (less negative) slope at $z<0.8$. We normalise each SMF to the comoving volume occupied by the corresponding environment, and relate estimates from different redshift bins. We observe an evolution of the SMF of VIPERS galaxies in high densities, while the low-density one is nearly constant. We compare these results to semi-analytical models and find consistent environmental signatures. We discuss how the halo mass function and fraction of central/satellite galaxies depend on the environments considered, making intrinsic and environmental properties of galaxies physically coupled, and therefore difficult to disentangle. The evolution of our low-density regions is well described by the formalism introduced by Peng et al.~(2010), and is consistent with the idea that galaxies become progressively passive because of internal physical processes. The same formalism could also describe the evolution of the SMF in the high density regions, but only if a significant contribution from dry mergers is considered. [Abridged]

preprint2015arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Hierarchical scaling and biasing

We investigate the higher-order correlation properties of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to test the hierarchical scaling hypothesis at z~1 and the dependence on galaxy luminosity, stellar mass, and redshift. We also aim to assess deviations from the linearity of galaxy bias independently from a previously performed analysis of our survey (Di Porto et al. 2014). We have measured the count probability distribution function in cells of radii 3 < R < 10 Mpc/h, deriving $σ_{8g}$, the volume-averaged two-,three-,and four-point correlation functions and the normalized skewness $S_{3g}$ and kurtosis $S_{4g}$ for volume-limited subsamples covering the ranges $-19.5 \le M_B(z=1.1)-5log(h) \le -21.0$, $9.0 < log(M*/M_{\odot} h^{-2}) \le 11.0$, $0.5 \le z < 1.1$. We have thus performed the first measurement of high-order correlations at z~1 in a spectroscopic redshift survey. Our main results are the following. 1) The hierarchical scaling holds throughout the whole range of scale and z. 2) We do not find a significant dependence of $S_{3g}$ on luminosity (below z=0.9 $S_{3g}$ decreases with luminosity but only at 1σ-level). 3) We do not detect a significant dependence of $S_{3g}$ and $S_{4g}$ on scale, except beyond z~0.9, where the dependence can be explained as a consequence of sample variance. 4) We do not detect an evolution of $S_{3g}$ and $S_{4g}$ with z. 5) The linear bias factor $b=σ_{8g}/σ_{8m}$ increases with z, in agreement with previous results. 6) We quantify deviations from the linear bias by means of the Taylor expansion parameter $b_2$. Our results are compatible with a null non-linear bias term, but taking into account other available data we argue that there is evidence for a small non-linear bias term.

preprint2014arXiv

An improved measurement of baryon acoustic oscillations from the correlation function of galaxy clusters at $z \sim 0.3$

We detect the peak of baryon acoustic oscillations (BAO) in the two-point correlation function of a spectroscopic sample of $25226$ clusters selected from the Sloan Digital Sky Survey. Galaxy clusters, as tracers of massive dark matter haloes, are highly biased structures. The linear bias $b$ of the sample considered in this work, that we estimate from the projected correlation function, is $b σ_8 = 1.72 \pm 0.03$. Thanks to the high signal in the cluster correlation function and to the accurate spectroscopic redshift measurements, we can clearly detect the BAO peak and determine its position, $s_p$, with high accuracy, despite the relative paucity of the sample. Our measurement, $s_p = 104 \pm 7 \, \mathrm{Mpc} \, h^{-1}$, is in good agreement with previous estimates from large galaxy surveys, and has a similar uncertainty. The BAO measurement presented in this work thus provides a new strong confirmation of the concordance cosmological model and demonstrates the power and promise of galaxy clusters as key probes for cosmological applications based on large scale structures.

preprint2014arXiv

The MUSIC of CLASH: predictions on the concentration-mass relation

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

preprint2014arXiv

The VIMOS Public Extragalactic Redshift Survey - Searching for Cosmic Voids

The characterisation of cosmic voids gives unique information about the large-scale distribution of galaxies, their evolution and the cosmological model. We identify and characterise cosmic voids in the VIMOS Public Extragalactic Redshift Survey (VIPERS) at redshift 0.55 < z < 0.9. A new void search method is developed based upon the identification of empty spheres that fit between galaxies. The method can be used to characterise the cosmic voids despite the presence of complex survey boundaries and internal gaps. We investigate the impact of systematic observational effects and validate the method against mock catalogues. We measure the void size distribution and the void-galaxy correlation function. We construct a catalogue of voids in VIPERS. The distribution of voids is found to agree well with the distribution of voids found in mock catalogues. The void-galaxy correlation function shows indications of outflow velocity from the voids.

preprint2014arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS): A quiescent formation of massive red-sequence galaxies over the past 9 Gyr

We explore the evolution of the Colour-Magnitude Relation (CMR) and Luminosity Function (LF) at 0.4<z<1.3 from the VIMOS Public Extragalactic Redshift Survey (VIPERS) using ~45,000 galaxies with precise spectroscopic redshifts down to i'_AB<22.5 over ~10.32 deg^2 in two fields. From z=0.5 to z=1.3 the LF and CMR are well defined for different galaxy populations and M^*_B evolves by ~1.04(1.09)+/-0.06(0.10) mag for the total (red) galaxy sample. We compare different criteria for selecting early-type galaxies (ETGs): (1) fixed cut in rest-frame (U-V) colours, (2) evolving cut in (U-V) colours, (3) rest-frame (NUV-r')-(r'-K) colour selection, and (4) SED classification. Regardless of the method we measure a consistent evolution of the red-sequence (RS). Between 0.4<z<1.3 we find a moderate evolution of the RS intercept of Delta(U-V)=0.28+/-0.14 mag, favouring exponentially declining star formation (SF) histories with SF truncation at 1.7<=z<=2.3. Together with the rise in the ETG number density by 0.64 dex since z=1, this suggests a rapid build-up of massive galaxies (M>10^11 M_sun) and expeditious RS formation over a short period of ~1.5 Gyr starting before z=1. This is supported by the detection of ongoing SF in ETGs at 0.9<z<1.0, in contrast with the quiescent red stellar populations of ETGs at 0.5<z<0.6. There is an increase in the observed CMR scatter with redshift, two times larger than in galaxy clusters and at variance with theoretical models. We discuss possible physical mechanisms that support the observed evolution of the red galaxy population. Our findings point out that massive galaxies have experienced a sharp SF quenching at z~1 with only limited additional merging. In contrast, less-massive galaxies experience a mix of SF truncation and minor mergers which build-up the low- and intermediate-mass end of the CMR.

preprint2014arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Never mind the gaps: comparing techniques to restore homogeneous sky coverage

[Abridged] Non-uniform sampling and gaps in sky coverage are common in galaxy redshift surveys, but these effects can degrade galaxy counts-in-cells and density estimates. We carry out a comparison of methods that aim to fill the gaps to correct for the systematic effects. Our study is motivated by the analysis of the VIMOS Extragalactic Redshift Survey (VIPERS), a flux-limited survey (i<22.5) based on one-pass observations with VIMOS, with gaps covering 25% of the surveyed area and a mean sampling rate of 35%. Our findings are applicable to other surveys with similar observing strategies. We compare 1) two algorithms based on photometric redshift, that assign redshifts to galaxies based on the spectroscopic redshifts of the nearest neighbours, 2) two Bayesian methods, the Wiener filter and the Poisson-Lognormal filter. Using galaxy mock catalogues we quantify the accuracy of the counts-in-cells measurements on scales of R=5 and 8 Mpc/h after applying each of these methods. We also study how they perform to account for spectroscopic redshift error and inhomogeneous and sparse sampling rate. We find that in VIPERS the errors in counts-in-cells measurements on R<10 Mpc/h scales are dominated by the sparseness of the sample. All methods underpredict by 20-35% the counts at high densities. This systematic bias is of the same order as random errors. No method outperforms the others. Random and systematic errors decrease for larger cells. We show that it is possible to separate the lowest and highest densities on scales of 5 Mpc/h at redshifts 0.5<z<1.1, over a large volume such as in VIPERS survey. This is vital for the characterisation of cosmic variance and rare populations (e.g, brightest galaxies) in environmental studies at these redshifts.

preprint2013arXiv

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

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

preprint2013arXiv

Large-scale inhomogeneities of the intracluster medium: improving mass estimates using the observed azimuthal scatter

Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study the ICM of inhomogeneities, focusing on the ones on the large scale that, unlike clumps, are the most difficult to identify. To this purpose we introduce the concept of residual clumpiness, C_R, that quantifies the large-scale inhomogeneity of the ICM. After showing that this quantity can be robustly defined for relaxed systems, we characterize how it varies with radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces an overestimate in the determination of the density profile from the X-ray emission, which translates into a systematic overestimate of 6 (12)% in the measurement of M_gas at R_200 for our relaxed (perturbed) cluster sample. At the same time, the increase of C_R with radius introduces also a ~2% systematic underestimate in the measurement of the hydrostatic-equilibrium mass (M_he), which adds to the previous one generating a systematic ~8.5% overestimate in f_gas in our relaxed sample. Since the residual clumpiness of the ICM is not directly observable, we study its correlation with the azimuthal scatter in the X-ray surface brightness of the halo and in the y-parameter profiles. We find that their correlation is highly significant (r_S = 0.6-0.7), allowing to define the azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as robust proxies of C_R. After providing a function that connects the two quantities, we obtain that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias in the measurement of M_gas for relaxed objects, which becomes (0+/-2)% up to 2R_200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of M_he and f_gas, although a significant halo to halo scatter remains. (abridged)

preprint2013arXiv

Structure formation in cosmologies with oscillating dark energy

{abridged} We study the imprints on the formation and evolution of cosmic structures of dynamical dark energy models, characterized by an oscillating equation of state. The redshift evolution of the equation of state parameter w(z) for dark energy is characterized by two parameters, describing the amplitude and the frequency of the oscillations. We consider six different oscillating dark energy models, each characterized by a different set of parameter values. Under the common assumption that dark energy is not clustering on the scales of interest, we study different aspects of cosmic structure formation. In particular, we self-consistently solve the spherical collapse problem. We then estimate the behavior of several cosmological observables, such as the linear growth factor, the Integrated Sachs-Wolfe (ISW) effect, the number counts of massive structures, and the matter and cosmic shear power spectra. We show that, independently of the amplitude and the frequency of the dark energy oscillations, none of the aforementioned observables show an oscillating behavior as a function of redshift. This is a consequence of the said observables' being integrals over some functions of the expansion rate over cosmic history. We also notice that deviations with respect to the expectations for a fiducial LambdaCDM cosmology are generically small, and in the majority of the cases distinguishing an oscillating dark energy model would be difficult. Exceptions to this conclusion are provided by the cosmic shear power spectrum, which for some of the models shows a difference at the level of \sim 10% over a wide range of angular scales, and the abundance of galaxy clusters, which is modified at the $\sim 10-20%$ level at $z \gtrsim 0.6$ for future wide weak lensing surveys.

preprint2013arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS): $Ω_{\rm m_0}$ from the galaxy clustering ratio measured at $z \sim 1$

We use a sample of about 22,000 galaxies at $0.65<z<1.2$ from the VIPERS PDR-1 catalogue, to constrain the cosmological model through a measurement of the galaxy {\it clustering ratio} $η_{g,R}$. This statistic has favourable properties, being defined as the ratio of two quantities characterizing the smoothed density field in spheres of given radius $R$: the value of its correlation function on a multiple of this scale, $ξ(nR)$, and its variance $σ^2(R)$. For sufficiently large values of $R$, this is a universal number, capturing 2-point clustering information independently of the linear bias and linear redshift-space distortions of the specific galaxy tracers. In this paper we discuss in detail how to extend the application of $η_{g,R}$ to quasi-linear scales and how to control and remove observational selection effects which are typical of redshift surveys as VIPERS. We verify the accuracy and efficiency of these procedures using mock catalogues that match the survey selection process. These results evidence the robustness of $η_{g,R}$ to non-linearities and observational effects, which is related to its very definition as a ratio of quantities that are similarly affected. We measure $Ω_{m,0}=0.270_{-0.025}^{+0.029}$. In addition to the great precision achieved on our estimation of $Ω_m$ using VIPERS PDR-1, this result is remarkable because it appears to be in good agreement with a recent estimate $z\simeq 0.3$, obtained applying the same technique to the SDSS-LRG catalogue. It, therefore, suports the robustness of the present analysis. Moreover, the combination of these two measurements at $z\sim 0.3$ and $z\sim 0.9$ provides us with a very precise estimate $Ω_{m,0}=0.274\pm0.017$ which highlights the great consistency between our estimation and other cosmological probes such as BAOs, CMB and Supernovae.

preprint2013arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). A precise measurement of the galaxy stellar mass function and the abundance of massive galaxies at redshifts 0.5<z<1.3

We measure the evolution of the galaxy stellar mass function from z=1.3 to z=0.5 using the first 53,608 redshifts of the ongoing VIMOS Public Extragalactic Survey (VIPERS). We estimate the galaxy stellar mass function at several epochs discussing in detail the amount of cosmic variance affecting our estimate. We find that Poisson noise and cosmic variance of the galaxy mass function in the VIPERS survey are comparable with the statistical uncertainties of large surveys in the local universe. VIPERS data allow us to determine with unprecedented accuracy the high-mass tail of the galaxy stellar mass function, which includes a significant number of galaxies that are usually too rare to detect with any of the past spectroscopic surveys. At the epochs sampled by VIPERS, massive galaxies had already assembled most of their stellar mass. We apply a photometric classification in the (U-V) rest-frame colour to compute the mass function of blue and red galaxies, finding evidence for the evolution of their contribution to the total number density budget: the transition mass above which red galaxies dominate is found to be about 10^10.4 M_sun at z=0.55 and evolves proportionally to (1+z)^3. We are able to trace separately the evolution of the number density of blue and red galaxies with masses above 10^11.4 M_sun, in a mass range barely studied in previous work. We find that for such large masses, red galaxies show a milder evolution with redshift, when compared to objects at lower masses. At the same time, we detect a population of similarly massive blue galaxies, which are no longer detectable below z=0.7. These results show the improved statistical power of VIPERS data, and give initial promising indications of mass-dependent quenching of galaxies at z~1. [Abridged]

preprint2013arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). A Support Vector Machine classification of galaxies, stars and AGNs

The aim of this work is to develop a comprehensive method for classifying sources in large sky surveys and we apply the techniques to the VIMOS Public Extragalactic Redshift Survey (VIPERS). Using the optical (u*, g', r', i') and NIR data (z', Ks), we develop a classifier, based on broad-band photometry, for identifying stars, AGNs and galaxies improving the purity of the VIPERS sample. Support Vector Machine (SVM) supervised learning algorithms allow the automatic classification of objects into two or more classes based on a multidimensional parameter space. In this work, we tailored the SVM for classifying stars, AGNs and galaxies, and applied this classification to the VIPERS data. We train the SVM using spectroscopically confirmed sources from the VIPERS and VVDS surveys. We tested two SVM classifiers and concluded that including NIR data can significantly improve the efficiency of the classifier. The self-check of the best optical + NIR classifier has shown a 97% accuracy in the classification of galaxies, 97 for stars, and 95 for AGNs in the 5-dimensional colour space. In the test on VIPERS sources with 99% redshift confidence, the classifier gives an accuracy equal to 94% for galaxies, 93% for stars, and 82% for AGNs. The method was applied to sources with low quality spectra to verify their classification, and thus increasing the security of measurements for almost 4 900 objects. We conclude that the SVM algorithm trained on a carefully selected sample of galaxies, AGNs, and stars outperforms simple colour-colour selection methods, and can be regarded as a very efficient classification method particularly suitable for modern large surveys.

preprint2013arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). An unprecedented view of galaxies and large-scale structure at 0.5<z<1.2

We describe the construction and general features of VIPERS, the VIMOS Public Extragalactic Redshift Survey. This `Large Programme' has been using the ESO VLT with the aim of building a spectroscopic sample of ~100,000 galaxies with i_{AB}<22.5 and 0.5<z<1.5. The survey covers a total area of ~24 deg^2 within the CFHTLS-Wide W1 and W4 fields. VIPERS is designed to address a broad range of problems in large-scale structure and galaxy evolution, thanks to a unique combination of volume (~ 5 x 10^7 h^{-3} Mpc^3) and sampling rate (~ 40%), comparable to state-of-the-art surveys of the local Universe, together with extensive multi-band optical and near-infrared photometry. Here we present the survey design, the selection of the source catalogue and the development of the spectroscopic observations. We discuss in detail the overall selection function that results from the combination of the different constituents of the project. This includes the masks arising from the parent photometric sample and the spectroscopic instrumental footprint, together with the weights needed to account for the sampling and the success rates of the observations. Using the catalogue of 53,608 galaxy redshifts composing the forthcoming VIPERS Public Data Release 1 (PDR-1), we provide a first assessment of the quality of the spectroscopic data. Benefiting from the combination of size and detailed sampling of this dataset, we conclude by presenting a map showing in unprecedented detail the large-scale distribution of galaxies between 5 and 8 billion years ago. [abridged]

preprint2013arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Galaxy clustering and redshift-space distortions at z=0.8 in the first data release

We present in this paper the general real- and redshift-space clustering properties of galaxies as measured in the first data release of the VIPERS survey. VIPERS is a large redshift survey designed to probe the distant Universe and its large-scale structure at 0.5 < z < 1.2. We describe in this analysis the global properties of the sample and discuss the survey completeness and associated corrections. This sample allows us to measure the galaxy clustering with an unprecedented accuracy at these redshifts. From the redshift-space distortions observed in the galaxy clustering pattern we provide a first measurement of the growth rate of structure at z = 0.8: fσ_8 = 0.47 +/- 0.08. This is completely consistent with the predictions of standard cosmological models based on Einstein gravity, although this measurement alone does not discriminate between different gravity models.

preprint2013arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS). Luminosity and stellar mass dependence of galaxy clustering at 0.5<z<1.1

We investigate the dependence of galaxy clustering on luminosity and stellar mass in the redshift range 0.5<z<1.1, using the first ~55000 redshifts from the VIMOS Public Extragalactic Redshift Survey (VIPERS). We measured the redshift-space two-point correlation functions (2PCF), and the projected correlation function, in samples covering different ranges of B-band absolute magnitudes and stellar masses. We considered both threshold and binned galaxy samples, with median B-band absolute magnitudes -21.6<MB-5log(h)<-19.5 and median stellar masses 9.8<log(M*[Msun/h^2])<10.7. We assessed the real-space clustering in the data from the projected correlation function, which we model as a power law in the range 0.2<r_p[Mpc/h]<20. Finally, we estimated the galaxy bias as a function of luminosity, stellar mass, and redshift, assuming a flat LCDM model to derive the dark matter 2PCF. We provide the best-fit parameters of the power-law model assumed for the real-space 2PCF -- the correlation length and the slope -- as well as the linear bias parameter, as a function of the B-band absolute magnitude, stellar mass, and redshift. We confirm and provide the tightest constraints on the dependence of clustering on luminosity at 0.5<z<1.1. We prove the complexity of comparing the clustering dependence on stellar mass from samples that are originally flux-limited and discuss the possible origin of the observed discrepancies. Overall, our measurements provide stronger constraints on galaxy formation models, which are now required to match, in addition to local observations, the clustering evolution measured by VIPERS galaxies between z=0.5 and z=1.1 for a broad range of luminosities and stellar masses.

preprint2013arXiv

The VIMOS Public Extragalactic Survey (VIPERS): First Data Release of 57 204 spectroscopic measurements

We present the first Public Data Release (PDR-1) of the VIMOS Public Extragalactic Survey (VIPERS). It comprises 57 204 spectroscopic measurements together with all additional information necessary for optimal scientific exploitation of the data, in particular the associated photometric measurements and quantification of the photometric and survey completeness. VIPERS is an ESO Large Programme designed to build a spectroscopic sample of ' 100 000 galaxies with iAB < 22.5 and 0.5 < z < 1.5 with high sampling rate (~45%). The survey spectroscopic targets are selected from the CFHTLS-Wide five-band catalogues in the W1 and W4 fields. The final survey will cover a total area of nearly 24 deg2, for a total comoving volume between z = 0.5 and 1.2 of ~4x10^7 h^(-3)Mpc^3 and a median galaxy redshift of z~0.8. The release presented in this paper includes data from virtually the entire W4 field and nearly half of the W1 area, thus representing 64% of the final dataset. We provide a detailed description of sample selection, observations and data reduction procedures; we summarise the global properties of the spectroscopic catalogue and explain the associated data products and their use, and provide all the details for accessing the data through the survey database (http://vipers.inaf.it) where all information can be queried interactively.

preprint2013arXiv

VIPERS: An Unprecedented View of Galaxies and Large-Scale Structure Halfway Back in the Life of the Universe

The VIMOS Public Extragalactic Redshift Survey (VIPERS) is an ongoing ESO Large Programme to map in detail the large-scale distribution of galaxies at 0.5 < z <1.2. With a combination of volume and sampling density that is unique for these redshifts, it focuses on measuring galaxy clustering and related cosmological quantities as part of the grand challenge of understanding the origin of cosmic acceleration. VIPERS has also been designed to guarantee a broader legacy, allowing detailed investigations of the properties and evolutionary trends of z~1 galaxies. The survey strategy exploits the specific advantages of the VIMOS spectrograph at the VLT, aiming at a final sample of nearly 100,000 galaxy redshifts to iAB = 22.5 mag, which represents the largest redshift survey ever performed with ESO telescopes. In this introductory article we describe the survey construction, together with early results based on a first sample of ~55,000 galaxies.

preprint2012arXiv

Constraining Primordial Magnetic Fields with Future Cosmic Shear Surveys

The origin of astrophysical magnetic fields observed in galaxies and clusters of galaxies is still unclear. One possibility is that primordial magnetic fields generated in the early Universe provide seeds that grow through compression and turbulence during structure formation. A cosmological magnetic field present prior to recombination would produce substantial matter clustering at intermediate/small scales, on top of the standard inflationary power spectrum. In this work we study the effect of this alteration on one particular cosmological observable, cosmic shear. We adopt the semi-analytic halo model in order to describe the non-linear clustering of matter, and feed it with the altered mass variance induced by primordial magnetic fields. We find that the convergence power spectrum is, as expected, substantially enhanced at intermediate/small angular scales, with the exact amplitude of the enhancement depending on the magnitude and power-law index of the magnetic field power spectrum. We use the predicted statistical errors for a future wide-field cosmic shear survey, on the model of the ESA Cosmic Vision mission \emph{Euclid}, in order to forecast constraints on the amplitude of primordial magnetic fields as a function of the spectral index. We find that the amplitude will be constrained at the level of $\sim 0.1$ nG for $n_B\sim -3$, and at the level of $\sim 10^{-7}$ nG for $n_B\sim 3$. The latter is at the same level of lower bounds coming from the secondary emission of gamma-ray sources, implying that for high spectral indices \emph{Euclid} will certainly be able to detect primordial magnetic fields, if they exist. The present study shows how large-scale structure surveys can be used for both understanding the origins of astrophysical magnetic fields and shedding new light on the physics of the pre-recombination Universe. (abridged)

preprint2012arXiv

Constraining primordial non-Gaussianity with cosmological weak lensing: shear and flexion

We examine the cosmological constraining power of future large-scale weak lensing surveys on the model of \emph{Euclid}, with particular reference to primordial non-Gaussianity. Our analysis considers several different estimators of the projected matter power spectrum, based on both shear and flexion, for which we review the covariances and Fisher matrices. The bounds provided by cosmic shear alone for the local bispectrum shape, marginalized over $σ_8$, are at the level of $Δf_\mathrm{NL} \sim 100$. We consider three additional bispectrum shapes, for which the cosmic shear constraints range from $Δf_\mathrm{NL}\sim 340$ (equilateral shape) up to $Δf_\mathrm{NL}\sim 500$ (orthogonal shape). The competitiveness of cosmic flexion constraints against cosmic shear ones depends on the galaxy intrinsic flexion noise, that is still virtually unconstrained. Adopting the very high value that has been occasionally used in the literature results in the flexion contribution being basically negligible with respect to the shear one, and for realistic configurations the former does not improve significantly the constraining power of the latter. Since the flexion noise decreases with decreasing scale, by extending the analysis up to $\ell_\mathrm{max} = 20,000$ cosmic flexion, while being still subdominant, improves the shear constraints by $\sim 10%$ when added. However on such small scales the highly non-linear clustering of matter and the impact of baryonic physics make any error estimation uncertain. By considering lower, and possibly more realistic, values of the flexion intrinsic shape noise results in flexion constraining power being a factor of $\sim 2$ better than that of shear, and the bounds on $σ_8$ and $f_\mathrm{NL}$ being improved by a factor of $\sim 3$ upon their combination. (abridged)

preprint2012arXiv

The effect of feedback on the emission properties of the Warm-Hot Intergalactic Medium

At present, 30-40 per cent of the baryons in the local Universe is still undetected. According to theoretical predictions, this gas should reside in filaments filling the large-scale structure (LSS) in the form of a Warm-Hot Intergalactic Medium (WHIM), at temperatures of 10^5 - 10^7 K, thus emitting in the soft X-ray energies via free-free interaction and line emission from heavy elements. In this work we characterize the properties of the X-ray emission of the WHIM, and the LSS in general, focusing on the influence of different physical mechanisms, namely galactic winds (GWs), black-hole feedback and star-formation, and providing estimates of possible observational constraints. To this purpose we use a set of cosmological hydrodynamical simulations that include a self-consistent treatment of star-formation and chemical enrichment of the intergalactic medium, that allows us to follow the evolution of different metal species. We construct a set of simulated light-cones to make predictions of the emission in the 0.3-10 keV energy range. We obtain that GWs increase by a factor of 2 the emission of both galaxy clusters and WHIM. The amount of oxygen at average temperature and, consequently, the amount of expected bright Ovii and Oviii lines is increased by a factor of 3 due to GWs and by 20 per cent when assuming a top-heavy IMF. We compare our results with current observational constraints and find that the emission from faint groups and WHIM should account from half to all of the unresolved X-ray background in the 1-2 keV band.

preprint2012arXiv

The VIMOS Public Extragalactic Redshift Survey (VIPERS): spectral classification through Principal Component Analysis

We develop a Principal Component Analysis aimed at classifying a sub-set of 27,350 spectra of galaxies in the range 0.4 < z < 1.0 collected by the VIMOS Public Extragalactic Redshift Survey (VIPERS). We apply an iterative algorithm to simultaneously repair parts of spectra affected by noise and/or sky residuals, and reconstruct gaps due to rest-frame transformation, and obtain a set of orthogonal spectral templates that span the diversity of galaxy types. By taking the three most significant components, we find that we can describe the whole sample without contamination from noise. We produce a catalogue of eigen-coefficients and template spectra that will be part of future VIPERS data releases. Our templates effectively condense the spectral information into two coefficients that can be related to the age and star formation rate of the galaxies. We examine the spectrophotometric types in this space and identify early, intermediate, late and starburst galaxies.

preprint2011arXiv

Abell 611. II. X-ray and strong lensing analyses

We present the results of our analyses of the X-ray emission and of the strong lensing systems in the relaxed galaxy cluster Abell 611 (z=0.288). We infer the X-ray mass estimate deriving the density and temperature profiles of the intra-cluster medium within the radius r ~ 700 kpc through a non-parametric approach; assuming that the cluster is in hydrostatic equilibrium and adopting a matter density profile, we can recover the total mass distribution of Abell 611 via the X-ray data. We derive the total projected mass in the central region of Abell 611 performing a parametric analysis of its strong lensing features through the publicly available analysis software Lenstool. As a final step we compare the results obtained with both methods. We derive a good agreement between the X-ray and strong lensing total mass estimates in the central region (i.e. within the radius r ~100 kpc), while a marginal disagreement is found between the two mass estimates when extrapolating the strong lensing results in the outer spatial range. We suggest that in this case the X-ray/strong lensing mass disagreement can be explained by an incorrect estimate of the relative contributions of the baryonic component and of the dark matter, caused by the intrinsic degeneracy between the different mass components in the strong lensing analysis. We discuss the effect of some possible systematic errors that influence both mass estimates. We find a slight dependence of the measurements of the X-ray temperatures (and therefore of the X-ray total masses) on the background adopted in the spectral analysis, with total deviations on the value of M_200 of the order of the 1-sigma statistical error. The strong lensing mass results are instead sensitive to the parameterisation of the galactic halo mass in the central regions, in particular to the modelling of the Brightest Cluster Galaxy (BCG) baryonic component.

preprint2011arXiv

Comparison of an X-ray selected sample of massive lensing clusters with the MareNostrum Universe LCDM simulation

A long-standing problem of strong lensing by galaxy clusters regards the observed high rate of giant gravitational arcs as compared to the predictions in the framework of the "standard" cosmological model. Recently, few other inconsistencies between theoretical expectations and observations have been claimed which regard the large size of the Einstein rings and the high concentrations of few clusters with strong lensing features. All of these problems consistently indicate that observed galaxy clusters may be gravitational lenses stronger than expected. We use clusters extracted from the MareNostrum Universe to build up mock catalogs of galaxy clusters selected through their X-ray flux. We use these objects to estimate the probability distributions of lensing cross sections, Einstein rings, and concentrations for the sample of 12 MACS clusters at $z>0.5$ presented in Ebeling et al. (2007) and discussed in Zitrin et al. (2010). We find that simulated clusters produce $\sim 50%$ less arcs than observed clusters do. The medians of the distributions of the Einstein ring sizes differ by $\sim 25%$ between simulations and observations. We estimate that, due to cluster triaxiality and orientation biases affecting the lenses with the largest cross sections, the concentrations of the individual MACS clusters inferred from the lensing analysis should be up to a factor of $\sim 2$ larger than expected from the $Λ$CDM model. The arc statistics, the Einstein ring, and the concentration problems in strong lensing clusters are mitigated but not solved on the basis of our analysis. Nevertheless, due to the lack of redshifts for most of the multiple image systems used for modeling the MACS clusters, the results of this work will need to be verified with additional data. The upcoming CLASH program will provide an ideal sample for extending our comparison (abridged).

preprint2011arXiv

Euclid Definition Study Report

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

preprint2011arXiv

Imprints of primordial non-Gaussianity on the number counts of cosmic shear peaks

We studied the effect of primordial non-Gaussianity with varied bispectrum shapes on the number counts of signal-to-noise peaks in wide field cosmic shear maps. The two cosmological contributions to this particular weak lensing statistic, namely the chance projection of Large Scale Structure and the occurrence of real, cluster-sized dark matter halos, have been modeled semi-analytically, thus allowing to easily introduce the effect of non-Gaussian initial conditions. We performed a Fisher matrix analysis by taking into account the full covariance of the peak counts in order to forecast the joint constraints on the level of primordial non-Gaussianity and the amplitude of the matter power spectrum that are expected by future wide field imaging surveys. We find that positive-skewed non-Gaussianity increases the number counts of cosmic shear peaks, more so at high signal-to-noise values, where the signal is mostly dominated by massive clusters as expected. The increment is at the level of ~1 for f_NL=10 and ~10 for f_NL=100 for a local shape of the primordial bispectrum, while different bispectrum shapes give generically a smaller effect. For a future survey on the model of the proposed ESA space mission Euclid and by avoiding the strong assumption of being capable to distinguish the weak lensing signal of galaxy clusters from chance projection of Large Scale Structures we forecasted a 1-sigma error on the level of non-Gaussianity of ~30-40 for the local and equilateral models, and of ~100-200 for the less explored enfolded and orthogonal bispectrum shapes.

preprint2011arXiv

Matter power spectra in dynamical-Dark Energy cosmologies

(abridged) We used a suite of numerical cosmological simulations in order to investigate the effect of gas cooling and star formation on the large scale matter distribution. The simulations follow the formation of cosmic structures in five different Dark Energy models: the fiducial $Λ$CDM cosmology and four models where the Dark Energy density is allowed to have a non-trivial redshift evolution. For each cosmology we have a control run with dark matter only, in order to allow a direct assessment of the impact of baryonic processes. We found that the power spectra of gas and stars, as well as the total matter power spectrum, are in qualitative agreement with the results of previous works in the framework of the fiducial model, although several quantitative differences exist. We used the halo model in order to investigate the backreaction of gas and stars on the dark matter distribution, finding that it is very well reproduced by increasing the average dark matter halo concentration by 17%, irrespective of the mass. Moving to model universes dominated by dynamical Dark Energy, it turns out that they introduce a specific signature on the power spectra of the various matter components, that is qualitatively independent of the exact cosmology considered. This generic shape is well captured by the halo model, however the finer details of the dark matter power spectrum can be precisely captured only at the cost of a few slight modifications to the ingredients entering the model. The backreaction of baryons onto the dark matter distribution works pretty much in the same way as in the reference $Λ$CDM model. Nonetheless, the increment in average concentration is less pronounced than in the fiducial model (only $\sim 10%$), in agreement with a series of other clues pointing toward the fact that star formation is less efficient when Dark Energy displays a dynamical evolution.

preprint2011arXiv

On a novel approach using massive clusters at high redshifts as cosmological probe

In this work we propose a novel method for testing the validity of the fiducial LCDM cosmology by measuring the cumulative distribution function of the most massive haloes in a sample of subvolumes of identical size tiled on the sky at a fixed redshift. The fact that the most massive clusters probe the high-mass tail of the mass function, where the difference between LCDM and alternative cosmological models is strongest, makes our method particularly interesting as a cosmological probe. We utilise general extreme value statistics (GEV) to obtain a cumulative distribution function of the most massive objects in a given volume. We sample this distribution function according to the number of patches covered by the survey area for a range of different "test cosmologies" and for differently accurate mass estimations of the haloes. By fitting this sample with the GEV distribution function, we can study which parameters are the most sensitive with respect to the test cosmologies. We find that the peak of the probability distribution function of the most massive halo is well suited to test the validity of the fiducial LCDM model, once we are able to establish a sufficiently complete large-area survey with M_lim=10^14.5 M_sun/h (M_lim=10^14 M_sun/h) at redshifts above z=1 (z=1.5). Being of cumulative nature the proposed measure is robust and an accuracy of 20-30% in the cluster masses would be sufficient to test for alternative models. Since one only needs the most massive system in each angular patch, this method would be ideally suited as a first fast consistency check before going into a more complex statistical analysis of the observed halo sample.

preprint2011arXiv

ORIGIN: Metal Creation and Evolution from the Cosmic Dawn

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

preprint2011arXiv

The effect of primordial non-Gaussianity on the skeleton of cosmic shear maps

(abridged) We explore the imprints of deviations from Gaussian primordial density fluctuations on the skeleton of the large-scale matter distribution as mapped through cosmological weak lensing. We computed the skeleton length of simulated effective convergence maps covering $\sim 35$ sq. deg each, extracted from a suite of cosmological $n-$body runs with different levels of local primordial non-Gaussianity. The latter is expected to alter the structure formation process with respect to the fiducial Gaussian scenario, and thus to leave a signature on the cosmic web. We found that alterations of the initial conditions consistently modify both the cumulative and the differential skeleton length, although the effect is generically smaller than the cosmic variance and depends on the smoothing of the map prior to the skeleton computation. Nevertheless, the qualitative shape of these deviations is rather similar to their primordial counterparts, implying that skeleton statistics retain good memory of the initial conditions. We performed a statistical analysis in order to find out at what Confidence Level primordial non-Gaussianity could be constrained by the skeleton test on cosmic shear maps of the size we adopted. At 68.3% Confidence Level we found an error on the measured level of primordial non-Gaussianity of $Δf_\mathrm{NL}\sim 300$, while at 90% Confidence Level it is of $Δf_\mathrm{NL}\sim 500$. While these values by themselves are not competitive with the current constraints, weak lensing maps larger than those used here would have a smaller field-to-field variance, and thus would likely lead to tighter constraints. A rough estimate indicates $Δf_\mathrm{NL} \sim$ a few tens at 68.3% Confidence Level for an all-sky weak lensing survey.

preprint2010arXiv

Abell 611: I. weak lensing analysis with LBC

Aims. The Large Binocular Cameras (LBC) are two twin wide field cameras (FOV ~ 23'x 25') mounted at the prime foci of the 8.4m Large Binocular Telescope (LBT). We performed a weak lensing analysis of the z=0.288 cluster Abell 611 on g-band data obtained by the blue-optimized Large Binocular Camera in order to estimate the cluster mass. Methods. Due to the complexity of the PSF of LBC, we decided to use two different approaches, KSB and Shapelets, to measure the shape of background galaxies and to derive the shear signal produced by the cluster. Then we estimated the cluster mass with both aperture densitometry and parametric model fits. Results. The combination of the large aperture of the telescope and the wide field of view allowed us to map a region well beyond the expected virial radius of the cluster and to get a high surface density of background galaxies (23 galaxies/arcmin^2). This made possible to estimate an accurate mass for Abell 611. We find that the mass within 1.5 Mpc is: $(8\pm3)\times 10^{14} M_\odot$ from the aperture mass technique and $(5\pm1)\times 10^{14} M_\odot$ using the model fitting by a NFW mass density profile, for both Shapelets and KSB methods. This analysis demonstrates that LBC is a powerful instrument for weak gravitational lensing studies.

preprint2010arXiv

Cosmic shear statistics in cosmologies with non-Gaussian initial conditions

We computed the power spectrum of weak cosmic shear in models with non-Gaussian primordial density fluctuations. Cosmological initial conditions deviating from Gaussianity have recently attracted much attention in the literature, especially with respect to their effect on the formation of non-linear structures and because of the bounds that they can put on the inflationary epoch. The fully non-linear matter power spectrum was evaluated with the use of the physically motivated, semi-analytic halo model, where the mass function and linear halo bias were suitably corrected for non-Gaussian cosmologies. In agreement with previous work, we found that a level of non-Gaussianity compatible with CMB bounds and with positive skewness produces an increase in power of the order of a few percent at intermediate scales. We then used the matter power spectrum, together with observationally motivated background source redshift distributions in order to compute the cosmological weak lensing power spectrum. We found that the degree of deviation from the power spectrum of the reference Gaussian model is small compared to the statistical error expected from even future weak lensing surveys. However, summing the signal over a large range of multipoles can beat down the noise, bringing to a significant detection of non-Gaussianity at the level of $|f_\mathrm{NL}| \simeq $ few tens, when all other cosmological parameters are held fixed. Finally, we have shown that the constraints on the level of non-Gaussianity can be improved by $\sim 20%$ with the use of weak lensing tomography.

preprint2010arXiv

Primordial density perturbations with running spectral index: impact on non-linear cosmic structures

(abridged) We explore the statistical properties of non-linear cosmic structures in a flat $Λ$CDM cosmology in which the index of the primordial power spectrum for scalar perturbations is allowed to depend on the scale. Within the inflationary paradigm, the running of the scalar spectral index can be related to the properties of the inflaton potential, and it is hence of critical importance to test it with all kinds of observations, which cover the linear and non-linear regime of gravitational instability. We focus on the amount of running $α_{\mathrm{S},0}$ allowed by an updated combination of CMB anisotropy data and the 2dF Galaxy Redshift Survey. Our analysis constrains $α_{\mathrm{S},0} = -0.051^{+0.047}_{-0.053}$ $(-0.034^{+0.039}_{-0.040})$ at 95% Confidence Level when (not) taking into account primordial gravitational waves in a ratio as predicted by canonical single field inflation, in agreement with other works. For the cosmological models best fitting the data both with and without running we studied the abundance of galaxy clusters and of rare objects, the halo bias, the concentration of dark matter halos, the Baryon Acoustic Oscillation, the power spectrum of cosmic shear, and the Integrated Sachs-Wolfe effect. We find that counting galaxy clusters in future X-ray and Sunyaev-Zel'dovich surveys could discriminate between the two models, more so if broad redshift information about the cluster samples will be available. Likewise, measurements of the power spectrum of cosmological weak lensing as performed by planned all-sky optical surveys such as EUCLID could detect a running of the primordial spectral index, provided the uncertainties about the source redshift distribution and the underlying matter power spectrum are well under control.

preprint2010arXiv

Statistical properties of SZ and X-ray cluster detections

Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal SZ signal, and compare them to X-ray cluster detections. Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the simulation volume to identify true and spurious detections. Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations with the characteristics of ACT detect clusters with masses M>6-9e13 M_o/h, quite independent of redshift, reach 50% completeness at ~1e14 M_o/h and 100% completeness at ~2e14 M_o/h. Samples are contaminated by only a few per cent spurious detections. This is broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and X-ray maps confirms that the filters introduce at most a negligible bias.

preprint2010arXiv

The potential of X-ray cluster surveys to constrain primordial non-Gaussianity

We present forecasts for constraints on deviations from Gaussian distribution of primordial density perturbations from future X-ray surveys of galaxy clusters. Our analysis is based on computing the Fisher-Matrix for number counts and large-scale power spectrum of clusters. We consider a survey with high-sensitivity and wide-area to detect about 2.5 x 10^5 extended sources. Based on the self-calibration approach, and including Planck priors in our analysis, we constrain 9 cosmological parameters and 4 nuisance parameters, which define the relation between cluster mass and X-ray flux. Because of the scale dependence of large-scale bias induced by local-shape non-Gaussianity, we find that the power spectrum provides strong constraints on the non-Gaussianity f_nl parameter, which complement the stringent constraints on the power spectrum normalization, σ_8, from the number counts. To quantify the joint constraints on theese two parameters, that specify the timing of structure formation for a fixed background expansion, we define the figure-of-merit FoM_SFT = (det[Cov(sigma_8,f_nl)])^{-1/2}. We find that our surveys constrain deviations from Gaussianity with a precision of Δf_nl~10 at 1 σconfidence level, with FoM_SFT~39. We point out that constraints on f_nl are weakly sensitive to the uncertainties in the knowledge of nuisance parameters. As an application of non-Gaussian constraints from available data, we analyse the impact of positive skewness on the occurrence of XMMU-J2235, a massive distant cluster discovered at z~1.4. We confirm that in a WMAP-7 Gaussian Lambda CDM cosmology, within the survey volume, ~5x10^{-3} objects like this are expected to be found: to increase this probability by a factor of at least 10, one needs to evade either the available constraints on f_nl or sigma_8.

preprint2009arXiv

Imprints of primordial non-Gaussianities in X-ray and SZ signals from galaxy clusters

Several inflationary models predict the possibility that the primordial perturbations of the density field may contain a degree of non-Gaussianity which would influence the subsequent evolution of cosmic structures at large scales. In order to study their impact, we use a set of three cosmological DM-only simulations starting from initial conditions with different levels of non-Gaussianity: f_NL=0,+/-100. More specifically, we focus on the distribution of galaxy clusters at different redshifts and, using suitable scaling relations, we determine their X-ray and SZ signals. Our analysis allows us to estimate the differences in the logN-logS and logN-logY due to the different initial conditions and to predict the cluster counts at different redshifts expected for future surveys (eROSITA and SPT). We also use a second set of simulations assuming a different cosmological scenario to estimate how the dependence on f_NL is degenerate with respect to other parameters. Our results indicate that the effects introduced by a realistic amount of primordial non-Gaussianity are small when compared to the ones connected with current uncertainties in cosmological parameters, particularly with sigma_8. However, if future surveys will be associated with optical follow-up campaigns to determine the cluster redshift, an analysis of the samples at z>1 can provide significant constraints on f_NL. In particular we predict that the SPT cluster survey will be able to detect ~1000 clusters at z>1 for the Gaussian case, with a difference of 15-20 per cent associated to f_NL=+/-100.

preprint2009arXiv

Large-scale non-Gaussian mass function and halo bias: tests on N-body simulations

The description of the abundance and clustering of halos for non-Gaussian initial conditions has recently received renewed interest, motivated by the forthcoming large galaxy and cluster surveys, which can potentially yield constraints of order unity on the non-Gaussianity parameter f_{NL}. We present tests on N-body simulations of analytical formulae describing the halo abundance and clustering for non-Gaussian initial conditions. We calibrate the analytic non-Gaussian mass function of Matarrese et al.(2000) and LoVerde et al.(2008) and the analytic description of clustering of halos for non-Gaussian initial conditions on N-body simulations. We find excellent agreement between the simulations and the analytic predictions if we make the corrections delta_c --> delta_c X sqrt{q} and delta_c --> δ_c X q where q ~ 0.75, in the density threshold for gravitational collapse and in the non-Gaussian fractional correction to the halo bias, respectively. We discuss the implications of this correction on present and forecasted primordial non-Gaussianity constraints. We confirm that the non-Gaussian halo bias offers a robust and highly competitive test of primordial non-Gaussianity.

preprint2009arXiv

Observing the clustering properties of galaxy clusters in dynamical dark-energy cosmologies

We study the clustering properties of galaxy clusters expected to be observed by various forthcoming surveys both in the X-ray and sub-mm regimes by the thermal Sunyaev-Zel'dovich effect. Several different background cosmological models are assumed, including the concordance $Λ$CDM and various cosmologies with dynamical evolution of the dark energy. Particular attention is paid to models with a significant contribution of dark energy at early times which affects the process of structure formation. Past light cone and selection effects in cluster catalogs are carefully modeled by realistic scaling relations between cluster mass and observables and by properly taking into account the selection functions of the different instruments. The results show that early dark-energy models are expected to produce significantly lower values of effective bias and both spatial and angular correlation amplitudes with respect to the standard $Λ$CDM model. Among the cluster catalogues studied in this work, it turns out that those based on \emph{eRosita}, \emph{Planck}, and South Pole Telescope observations are the most promising for distinguishing between various dark-energy models.

preprint2009arXiv

Studying the WHIM with Gamma Ray Bursts

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

preprint2009arXiv

The clustering of galaxy clusters in cosmological models with non-Gaussian initial conditions: Predictions for future surveys

We predict the biasing and clustering properties of galaxy clusters that are expected to be observed in the catalogues produced by two forthcoming X-ray and Sunyaev-Zel'dovich effect surveys. We study a set of flat cosmological models where the primordial density probability distribution shows deviations from Gaussianity in agreement with current observational bounds form the background radiation. We consider both local and equilateral shapes for the primordial bispectrum in non-Gaussian models. The two catalogues investigated are those produced by the \emph{e}ROSITA wide survey and from a survey based on South Pole Telescope observations. It turns out that both the bias and observed power spectrum of galaxy clusters are severely affected in non-Gaussian models with local shape of the primordial bispectrum, especially at large scales. On the other hand, models with equilateral shape of the primordial bispectrum show only a mild effect at all scales, that is difficult to be detected with clustering observations. Between the two catalogues, the one performing better is the \emph{e}ROSITA one, since it contains only the largest masses, that are more sensitive to primordial non-Gaussianity.

preprint2008arXiv

Primordial non-Gaussianities in the Intergalactic Medium

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

preprint2008arXiv

The effects of primordial non-Gaussianity on the cosmological reionization

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

preprint2008arXiv

The mass density field in simulated non-Gaussian scenarios

In this work we study the properties of the mass density field in the non-Gaussian world models simulated by Grossi et al. 2007. In particular we focus on the one-point density probability distribution function of the mass density field in non-Gausian models with quadratic non-linearities quantified by the usual parameter f_NL. We find that the imprints of primordial non-Gaussianity are well preserved in the negative tail of the probability function during the evolution of the density perturbation. The effect is already noticeable at redshifts as large as 4 and can be detected out to the present epoch. At z=0 we find that the fraction of the volume occupied by regions with underdensity delta < -0.9, typical of voids, is about 1.3 per cent in the Gaussian case and increases to ~2.2 per cent if f_NL=-1000 while decreases to ~0.5 per cent if f_NL=+1000. This result suggests that void-based statistics may provide a powerful method to detect non-Gaussianity even at low redshifts which is complementary to the measurements of the higher-order moments of the probability distribution function like the skewness or the kurtosis for which deviations from the Gaussian case are detected at the 25-50 per cent level.

preprint2005arXiv

A fast method for computing strong-lensing cross sections: Application to merging clusters

Strong gravitational lensing by irregular mass distributions, such as galaxy clusters, is generally not well quantified by cross sections of analytic mass models. Computationally expensive ray-tracing methods have so far been necessary for accurate cross-section calculations. We describe a fast, semi-analytic method here which is based on surface integrals over high-magnification regions in the lens plane and demonstrate that it yields reliable cross sections even for complex, asymmetric mass distributions. The method is faster than ray-tracing simulations by factors of $\sim30$ and thus suitable for large cosmological simulations, saving large amounts of computing time. We apply this method to a sample of galaxy cluster-sized dark matter haloes with simulated merger trees and show that cluster mergers approximately double the strong-lensing optical depth for lens redshifts $z_\mathrm{l}\gtrsim0.5$ and sources near $z_\mathrm{s} = 2$. We believe that this result hints at one possibility for understanding the recently detected high arcs abundance in clusters at moderate and high redshifts, and is thus worth further studies.

preprint2004arXiv

Beppo-SAX temperature maps of galaxy clusters in the Corona Borealis supercluster: A2061, A2067 and A2124

In this paper we present the analysis of Beppo-SAX observations of the cluster pairs A2061-A2067 and A2122-A2124, located in the Corona Borealis supercluster, which have been selected as candidate merging clusters. The aim of this work is to study the physics of the intracluster medium and to look for the possible presence of merging signatures. We derived the global temperatures and abundances and the temperature profiles and maps for these clusters. We do not find any significant evidence of interaction between the clusters forming the pairs, but we detect a candidate shock inside A2061. On the basis of the X-ray and optical properties of this cluster we propose a scenario in which a group is falling inside A2061. This interaction is in the phase in which the cores have not encountered yet and in which the formation of a shock is expected.

preprint2004arXiv

Comparing the temperatures of galaxy clusters from hydro-N-body simulations to Chandra and XMM-Newton observations

Theoretical studies of the physical processes guiding the formation and evolution of galaxies and galaxy clusters in the X-ray are mainly based on the results of numerical hydrodynamical N-body simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We demonstrate that the projected spectroscopic temperature of thermally complex clusters obtained from X-ray observations is always lower than the emission-weighed temperature, which is widely used in the analysis of numerical simulations. We show that this temperature bias is mainly related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has important implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. Because of this bias, in real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps, and may even not be detected. This may explain why, although numerical simulations predict that shock fronts are a quite common feature in clusters of galaxies, to date there are very few observations of objects in which they are clearly seen. To fix this problem we propose a new formula, the spectroscopic-like temperature function, and show that, for temperature larger than 3 keV, it approximates the spectroscopic temperature better than few per cent, making simulations more directly comparable to observations.

preprint2000arXiv

Large-scale distribution and spectral properties of galaxies in the Shapley Concentration

We present the results of a redshift survey of both cluster and intercluster galaxies in the central part of the Shapley Concentration, the richest supercluster of clusters in the nearby Universe, consisting of ~2000 radial velocities. We estimate the total overdensity in galaxies of the supercluster and its mass and we discuss the cosmological implications of these results. Moreover, using a Principal Components Analysis technique, we study the influence of the cluster and supercluster dynamics on the galaxy spectral morphology.

preprint1995arXiv

Sampling Effects on Cosmological Dipoles

We use numerical simulations to investigate the behaviour of the dipole moment of the spatial distribution of different kinds of mass tracers. We select density peaks of the simulated matter distribution with mean separations of 38 and 30 $h^{-1}$ Mpc to represent two samples of rich clusters of galaxies. We extract, from the same simulations, samples selected to mimic the full 3D galaxy distribution of IR galaxies, and the flux--limited IRAS and QDOT galaxy samples. We compare the dipole moments of these ``galaxy'' and ``cluster'' samples in order to assess the effects of sampling uncertainties and shot--noise on the relationship between the ``true'' underlying galaxy dipole and the dipoles obtained for clusters and for the flux--limited galaxy samples. The results of this analysis demonstrate that the dipoles of both the IRAS and QDOT--like catalogues should trace the full 3D dipole shape fairly accurately, with the loss however of about 15--20$\%$ of the total 3D dipole amplitude. Furthermore, using a simple argument based on linear perturbation theory, on the linear biasing assumption and on the amplitude of the cluster dipole relative to that of galaxies, we can estimate their relative biasing factors quite accurately and in agreement with results obtained by other methods.

preprint1995arXiv

The Cluster Distribution as a Test for Dark Matter Models. I: Clustering Properties

We present simulations of the cluster distribution in several dark matter models, using an optimized version of the truncated Zel'dovich approximation (TZA). We compare them with N-body cluster simulations and find that the TZA provides a very accurate description of the cluster distribution as long as fluctuations on the cluster mass scale are in the mildly non-linear regime. The simulated dark matter models are: Standard CDM (SCDM), Tilted CDM (TCDM) with n=0.7, Cold+Hot DM (CHDM) with 30% of hot component, low Hubble constant (h=0.3) CDM (LOWH) and a spatially flat low-density CDM model with Ω_0=0.2. We compare the simulations with a redshift sample of Abell/ACO clusters, using the integral of the 2-point correlation function and the probability density function. We find that the best models at reproducing the data are CHDM and LCDM. All the other models are ruled out. The reduced skewness S_3 is fairly constant with S_3=1.9, independent of the DM model and consistent with observational data. The abundances of clusters predicted using the Press--Schechter theory provide strong constraints: only the CHDM, LOWH and LCDM models appear to produce the correct number-density of clusters.

preprint1995arXiv

The Cluster Distribution as a Test for Dark Matter Models. II: The Dipole Structure

We study the dipole for a redshift sample of Abell/ACO clusters. To elucidate the constraints it places on dark matter models, we use numerical simulations based on the Zel'dovich approximation. We run 20 realizations of each of six different dark matter models: four of these have a density parameter $Ω_0=1$, while the other two have $Ω_0=0.2$, one with and one without a cosmological constant term. For the Abell/ACO sample, we evaluate the parameter $β=Ω_0^{0.6}/b_{ cl}$, which reaches its asymptotic value at $R=160 Mpc/h$. We find $β=0.15\pm 0.04$ and $0.25\pm 0.06$ in LG and CMB frames respectively. After identifying in the simulations LG-like observers, we construct mock cluster samples. We find that an alignment between the cluster dipole and observer velocity, such as that observed (about 20 degrees) should not be expected necessarily in the simulations. This, together with the large observer-to-observer variance estimates of $β$, makes it difficult to place any firm constraints on cosmological models. Our results demonstrate that the large amplitude and convergence depth of the cluster dipole cannot be taken as strong evidence either for or against a low--density Universe.

preprint1994arXiv

Velocity Fields in Non--Gaussian Cold Dark Matter Models

We analyse the large--scale velocity field obtained by N--body simulations of cold dark matter (CDM) models with non--Gaussian primordial density fluctuations, considering models with both positive and negative primordial skewness in the density fluctuation distribution. We study the velocity probability distribution and calculate the dependence of the bulk flow, one--point velocity dispersion and Cosmic Mach Number on the filtering size. We find that the sign of the primordial skewness of the density field provides poor discriminatory power on the evolved velocity field. All non--Gaussian models here considered tend to have lower velocity dispersion and bulk flow than the standard Gaussian CDM model, while the Cosmic Mach Number turns out to be a poor statistic in characterizing the models. Next, we compare the large--scale velocity field of a composite sample of optically selected galaxies as described by the Local Group properties, bulk flow, velocity correlation function and Cosmic Mach Number, with the velocity field of mock catalogues extracted from the N--body simulations. The comparison does not clearly permit to single out a best model: the standard Gaussian model is however marginally preferred by the maximum likelihood analysis.

preprint1993arXiv

Angular Distribution of Clustersin Skewed CDM Models

We perform a detailed investigation of the statistical properties of the projected distribution of galaxy clusters obtained in Cold Dark Matter (CDM) models with both Gaussian and skewed primordial density fluctuations. We use N-body simulations to construct a set artificial Lick maps. An objective cluster--finding algorithm is used to identify clusters of different richness. For Gaussian models, the overall number of clusters is too small in the standard CDM case, but a model with higher normalisation fares much better; non--Gaussian models with negative skewness also fit faily well. We apply several statistical tests to compare real and simulated cluster samples, such as the 2-point correlation function, the minimal spanning tree construction, the multifractal analysis and the skewness of cell counts. The emerging picture is that Gaussian models, even with a higher normalization, are in trouble. Skew-positive models are also ruled out, while skew-negative models can reproduce the observed clustering of galaxy clusters in the CDM framework.

preprint1993arXiv

Cluster Correlations in the Zel'dovich Approximation

We show how to simulate the clustering of rich clusters of galaxies using a technique based on the Zel'dovich approximation. This method well reproduces the spatial distribution of clusters obtainable from full N-body simulations at a fraction of the computational cost. We use an ensemble of large--scale simulations to assess the level and statistical significance of cluster clustering in open, tilted and flat versions of the Cold Dark Matter (CDM) model, as well as a model comprising a mixture of Cold and Hot Dark Matter (CHDM). We find the open and flat CDM models are excluded by the data. The tilted CDM model with a slight tilt is in marginal agreement, while larger tilt produces the right amount of clustering; CHDM is the best of all our models at reproducing the observations of cluster clustering. We find that {\em all} our models display a systematically weaker relationship between clustering length and mean cluster separation than seems to be implied by observations. We also note that the cluster bias factor, is not constant in any of the models, showing that one needs to be very careful when relating cluster clustering statistics to primordial density fluctuations.

preprint1993arXiv

Skewness as a Test of Non-Gaussian Primordial Density Fluctuations

We investigate the evolution of the skewness of the distribution of density fluctuations in CDM models with both Gaussian and non--Gaussian initial fluctuations. We show that the method proposed by Coles \& Frenk (1991), which uses the skewness of galaxy counts to test the hypothesis of Gaussian primordial density fluctuations, is a potentially powerful probe of initial conditions. As expected, the mass distribution in models with initially non--Gaussian fluctuations shows systematic departures from the Gaussian behaviour on intermediate to large scales. We investigate the effect of peculiar velocity distortions and normalisation upon the relationship between skewness and variance. These effects are generally small for the models we consider. Comparing our results to the QDOT measurements of the skewness, we find that our initially positive--skew models are clearly excluded by this analysis, but the available data do not rule out the negative--skew models.

preprint1993arXiv

Weighted Bias and Galaxy Clustering

We consider a weighted biasing scheme for galaxy clustering. This differs from previous treatments in the fact that the biased density field coincides with the background mass--density whenever the latter exceeds a given threshold value. There is some physical motivation for this scheme and it is in better accord with intuitive ideas than models based on the Kaiser (1984) analysis of the clustering of rich clusters. We explain how different classes of object could be biased in different ways with respect to the underlying density distribution but still have $b=1$. We also show that if one applies our scheme consistently a weak dependence of $b$ upon density can be implied. This could also be the reason why the correlation function of galaxies in groups does not differ substantially from the correlation function of all galaxies.