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.

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

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.

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

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

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.

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

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

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.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

We use the statistics of the VIPERS survey to investigate the relation between the surface mean stellar mass density Sigma=Mstar/(2*pi*Re^2) of massive passive galaxies (MPGs, Mstar>10^11 Msun) and their environment in the redshift range 0.5<z<0.8. Passive galaxies were selected on the basis of their NUVrK colors (~900 objects), and the environment was defined as the galaxy density contrast, delta, using the fifth nearest-neighbor approach. The analysis of Sigma vs. delta was carried out in two stellar mass bins. In galaxies with Mstar<2*10^11 Msun, no correlation between Sigma and delta is observed. This implies that the accretion of satellite galaxies, which is more frequent in denser environments and efficient in reducing the galaxy Sigma, is not relevant in the formation and evolution of these systems. Conversely, in galaxies with Mstar>2*10^11 Msun, we find an excess of MPGs with low Sigma and a deficit of high-Sigma MPGs in the densest regions wrt other environments. We interpret this result as due to the migration of some high-Sigma MPGs (<1% of the total population of MPGs) into low-Sigma MPGs, probably through mergers or cannibalism of small satellites. In summary, our results imply that the accretion of satellite galaxies has a marginal role in the mass-assembly history of most MPGs. We have previously found that the number density of VIPERS massive star-forming galaxies (MSFGs) declines rapidily from z=0.8 to z=0.5, which mirrors the rapid increase in the number density of MPGs. This indicates that the MSFGs at z>0.8 migrate to the MPG population. Here, we investigate the Sigma-delta relation of MSFGs at z>0.8 and find that it is consistent within 1 sigma with that of low-Sigma MPGs at z<0.8. Thus, the results of this and our previous paper show that MSFGs at z>0.8 are consistent in terms of number and environment with being the progenitors of low-Sigma MPGs at z<0.8.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

Cosmic voids are large underdense regions that, together with galaxy clusters, filaments and walls, build up the large-scale structure of the Universe. The void size function provides a powerful probe to test the cosmological framework. However, to fully exploit this statistics, the void sample has to be properly cleaned from spurious objects. Furthermore, the bias of the mass tracers used to detect these regions has to be taken into account in the size function model. In our work we test a cleaning algorithm and a new void size function model on a set of simulated dark matter halo catalogues, with different mass and redshift selections, to investigate the statistics of voids identified in a biased mass density field. We then investigate how the density field tracers' bias affects the detected size of voids. The main result of this analysis is a new model of the size function, parameterised in terms of the linear effective bias of the tracers used, which is straightforwardly inferred from the large-scale two-point correlation function. This represents a crucial step to exploit the method on real data catalogues. The proposed size function model has been accurately calibrated on mock catalogues, and used to validate the possibility to provide forecasts on the cosmological constraints, namely on the matter density contrast, $\Omega_{\rm M}$, and on the normalisation of the linear matter power spectrum, $\sigma_8$, at different redshifts.

Modified gravity and massive neutrino cosmologies are two of the most interesting scenarios that have been recently explored to account for possible observational deviations from the concordance $\Lambda$-cold dark matter ($\Lambda$CDM) model. In this context, we investigated the large-scale structure of the Universe by exploiting the \dustp simulations that implement, simultaneously, the effects of $f(R)$ gravity and massive neutrinos. To study the possibility of breaking the degeneracy between these two effects, we analysed the redshift-space distortions in the clustering of dark matter haloes at different redshifts. Specifically, we focused on the monopole and quadrupole of the two-point correlation function, both in real and redshift space. The deviations with respect to $\Lambda$CDM model have been quantified in terms of the linear growth rate parameter. We found that redshift-space distortions provide a powerful probe to discriminate between $\Lambda$CDM and modified gravity models, especially at high redshifts ($z \gtrsim 1$), even in the presence of massive neutrinos.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

The lack of power anomaly is an intriguing feature at the largest angular scales of the CMB anisotropy temperature pattern, whose statistical significance is not strong enough to claim any new physics beyond the standard cosmological model. We revisit the former statement by also considering polarisation data. We propose a new one-dimensional estimator which takes jointly into account the information contained in the TT, TE and EE CMB spectra. By employing this estimator on Planck 2015 low-$\ell$ data, we find that a random $\Lambda$CDM realisation is statistically accepted at the level of $3.68 \%$. Even though Planck polarisation contributes a mere $4 \%$ to the total information budget, its use pushes the lower-tail-probability down from the $7.22 \%$ obtained with only temperature data. Forecasts of future CMB polarised measurements, as e.g. the LiteBIRD satellite, can increase the polarisation contribution up to $6$ times with respect to Planck at low-$\ell$. We argue that the large-scale E-mode polarisation may play an important role in analysing CMB temperature anomalies with future mission.

We study the dependency of the concentration on mass and redshift using three large N-body cosmological hydrodynamic simulations carried out by the Magneticum project. We constrain the slope of the mass-concentration relation with an unprecedented mass range for hydrodynamic simulations and find a negative trend on the mass-concentration plane and a slightly negative redshift dependency, in agreement with observations and other numerical works. We also show how the concentration correlates with the fossil parameter, defined as the stellar mass ratio between the central galaxy and the most massive satellite, in agreement with observations. We find that haloes with high fossil parameter have systematically higher concentration and investigate the cause in two different ways. First we study the evolution of haloes that lives unperturbed for a long period of time, where we find that the internal region keeps accreting satellites as the fossil parameter increases and the scale radius decreases (which increases the concentration). We also study the dependency of the concentration on the virial ratio and the energy term from the surface pressure $E_s$. We conclude that fossil objects have higher concentration because they are dynamically relaxed, with no in-fall/out-fall material and had time to accrete their satellites.

Cosmic voids are the major volume component in the matter distribution of the Universe. They posses great potential for constraining dark energy as well as for testing theories of gravity. Nevertheless, in spite of their growing popularity as cosmological probes, a gap of knowledge between cosmic void observations and theory still persists. In particular, the void size function models proposed in literature have been proven unsuccessful in reproducing the results obtained from cosmological simulations in which cosmic voids are detected from biased tracers of the density field, undermining the possibility of using them as cosmological probes. The goal of this work is to cover this gap. In particular, we make use of the findings of a previous work in which we have improved the void selection procedure, presenting an algorithm that redefines the void ridges and, consequently, their radius. By applying this algorithm, we validate the volume conserving model of the void size function on a set of unbiased simulated density field tracers. We highlight the difference in the internal structure between voids selected in this way and those identified by the popular VIDE void finder. We also extend the validation of the model to the case of biased tracers. We find that a relation exists between the tracer used to sample the underlying dark matter density field and its unbiased counterpart. Moreover, we demonstrate that, as long as this relation is accounted for, the size function is a viable approach for studying cosmology with cosmic voids.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

We present the mass calibration for galaxy clusters detected with the AMICO code in KiDS DR3 data. The cluster sample comprises $\sim$ 7000 objects and covers the redshift range 0.1 < $z$ < 0.6. We perform a weak lensing stacked analysis by binning the clusters according to redshift and two different mass proxies provided by AMICO, namely the amplitude $A$ (measure of galaxy abundance through an optimal filter) and the richness $\lambda^*$ (sum of membership probabilities in a consistent radial and magnitude range across redshift). For each bin, we model the data as a truncated NFW profile plus a 2-halo term, taking into account uncertainties related to concentration and miscentring. From the retrieved estimates of the mean halo masses, we construct the $A$-$M_{200}$ and the $\lambda^*$-$M_{200}$ relations. The relations extend over more than one order of magnitude in mass, down to $M_{200} \sim 2 (5) \times 10^{13} M_\odot/h$ at $z$ = 0.2 (0.5), with small evolution in redshift. The logarithmic slope is $\sim 2.0$ for the $A$-mass relation, and $\sim 1.7$ for the $\lambda^*$-mass relation, consistent with previous estimations on mock catalogues and coherent with the different nature of the two observables.

Various galaxy classification schemes have been developed so far to constrain the main physical processes regulating evolution of different galaxy types. In the era of a deluge of astrophysical information and recent progress in machine learning, a new approach to galaxy classification becomes imperative. We employ a Fisher Expectation-Maximization unsupervised algorithm working in a parameter space of 12 rest-frame magnitudes and spectroscopic redshift. The model (DBk) and the number of classes (12) were established based on the joint analysis of standard statistical criteria and confirmed by the analysis of the galaxy distribution with respect to a number of classes and their properties. This new approach allows us to classify galaxies based just on their redshifts and UV-NIR spectral energy distributions. The FEM unsupervised algorithm has automatically distinguished 12 classes: 11 classes of VIPERS galaxies and an additional class of broad-line AGNs. After a first broad division into blue, green and red categories we obtained a further sub-division into three red, three green, and five blue galaxy classes. The FEM classes follow the galaxy sequence from the earliest to the latest types that is reflected in their colours (which are constructed from rest-frame magnitudes used in classification procedure) but also their morphological, physical, and spectroscopic properties (not included in the classification scheme). We demonstrate that the members of each class share similar physical and spectral properties. In particular, we are able to find three different classes of red passive galaxy populations. Thus, we demonstrate the potential of an unsupervised approach to galaxy classification and we retrieve the complexity of galaxy populations at z~0.7, a task that usual simpler colour-based approaches cannot fulfil.

Techniques to classify galaxies solely based on photometry will be necessary for future large cosmology missions, such as Euclid or LSST. However, the precision of classification is always lower in photometric surveys and can be systematically biased with respect to classifications based upon spectroscopic data. We verified how precisely the detailed classification scheme introduced by Siudek et al. (2018, hereafter: S1) for galaxies at z~0.7 could be reproduced if only photometric data are available. We applied the Fisher Expectation-Maximization (FEM) unsupervised clustering algorithm to 54,293 VIPERS galaxies working in a parameter space of reliable photometric redshifts and 12 corresponding rest-frame magnitudes. The FEM algorithm distinguishes four main groups: (1) red, (2) green, (3) blue, and (4) outliers. Each group is further divided into 3, 3, 4, and 2 subclasses, respectively. The accuracy of reproducing galaxy classes using spectroscopic data is high: 92%, 84%, 96% for red, green, and blue classes, respectively, except for dusty star-forming galaxies. The presented verification of the photometric classification demonstrates that large photometric samples can be used to distinguish different galaxy classes at z > 0.5 with an accuracy provided so far only by spectroscopic data except for particular galaxy classes.

We investigate the cosmological information contained in the cross-correlation between the Integrated Sachs-Wolfe (ISW) of the Cosmic Microwave Background (CMB) anisotropy pattern and galaxy clusters from future wide surveys. Future surveys will provide cluster catalogues with a number of objects comparable with galaxy catalogues currently used for the detection of the ISW signal by cross-correlation with the CMB anisotropy pattern. By computing the angular power spectra of clusters and the corresponding cross-correlation with CMB, we perform a signal-to-noise ratio (SNR) analysis for the ISW detection as expected from the eROSITA and the Euclid space missions. We discuss the dependence of the SNR of the ISW-cluster cross-correlation on the specifications of the catalogues and on the reference cosmology. We forecast that the SNRs for ISW-cluster cross-correlation are alightly smaller compared to those which can be obtained from future galaxy surveys but the signal is expected to be detected at high significance, i.e. more than $> 3\,\sigma$. We also forecast the joint constraints on parameters of model extensions of the concordance $\Lambda$CDM cosmology by combining CMB and the ISW-cluster cross-correlation.

We present the first catalogue of galaxy cluster candidates derived from the third data release of the Kilo Degree Survey (KiDS-DR3). The sample of clusters has been produced using the Adaptive Matched Identifier of Clustered Objects (AMICO) algorithm. In this analysis AMICO takes advantage of the luminosity and spatial distribution of galaxies only, not considering colours. In this way, we prevent any selection effect related to the presence or absence of the red-sequence in the clusters. The catalogue contains 7988 candidate galaxy clusters in the redshift range 0.1<z<0.8 down to S/N>3.5 with a purity approaching 95% over the entire redshift range. In addition to the catalogue of galaxy clusters we also provide a catalogue of galaxies with their probabilistic association to galaxy clusters. We quantify the sample purity, completeness and the uncertainties of the detection properties, such as richness, redshift, and position, by means of mock galaxy catalogues derived directly from the data. This preserves their statistical properties including photo-z uncertainties, unknown absorption across the survey, missing data, spatial correlation of galaxies and galaxy clusters. Being based on the real data, such mock catalogues do not have to rely on the assumptions on which numerical simulations and semi-analytic models are based on. This paper is the first of a series of papers in which we discuss the details and physical properties of the sample presented in this work.

Context. We present an estimation of cosmological parameters with clusters of galaxies. Aims. We constrain the $\Omega_m$, $\sigma_8$, and $w$ parameters from a stand-alone sample of X-ray clusters detected in the 50 deg$^2$ XMM-XXL survey with a well-defined selection function. Methods. We analyse the redshift distribution of a sample comprising 178 high S/N clusters out to a redshift of unity. The cluster sample scaling relations are determined in a self-consistent manner. Results. In a lambda cold dark matter ($\Lambda$CDM) model, the cosmology favoured by the XXL clusters compares well with results derived from the Planck S-Z clusters for a totally different sample (mass/redshift range, selection biases, and scaling relations). However, with this preliminary sample and current mass calibration uncertainty, we find no inconsistency with the Planck CMB cosmology. If we relax the $w$ parameter, the Planck CMB uncertainties increase by a factor of $\sim$10 and become comparable with those from XXL clusters. Combining the two probes allows us to put constraints on $\Omega_m$=0.316$\pm$0.060, $\sigma_8$=0.814$\pm$0.054, and $w$=-1.02$\pm$0.20. Conclusions. This first self-consistent cosmological analysis of a sample of serendipitous XMM clusters already provides interesting insights into the constraining power of the XXL survey. Subsequent analysis will use a larger sample extending to lower confidence detections and include additional observable information, potentially improving posterior uncertainties by roughly a factor of 3.

We present a novel suite of cosmological N-body simulations called the DUSTGRAIN-pathfinder, implementing simultaneously the effects of an extension to General Relativity in the form of $f(R)$ gravity and of a non-negligible fraction of massive neutrinos. We describe the generation of simulated weak lensing and cluster counts observables within a past light-cone extracted from these simulations. The simulations have been performed by means of a combination of the MG-GADGET code and a particle-based implementation of massive neutrinos, while the light-cones have been generated using the MapSim pipeline allowing us to compute weak lensing maps through a ray-tracing algorithm for different values of the source plane redshift. The mock observables extracted from our simulations will be employed for a series of papers focussed on understanding and possibly breaking the well-known observational degeneracy between $f(R)$ gravity and massive neutrinos, i.e. the fact that some specific combinations of the characteristic parameters for these two phenomena (the $f_{R0}$ scalar amplitude and the total neutrino mass $\Sigma m_{\nu}$) may result indistinguishable from the standard $\mathrm{\Lambda CDM}$ cosmology through several standard observational probes. In particular, in the present work we show how a tomographic approach to weak lensing statistics could allow - especially for the next generation of wide-field surveys - to disentangle some of the models that appear statistically indistinguishable through standard single-redshift weak lensing probe.

The VIPERS galaxy survey has measured the clustering of $0.5<z<1.2$ galaxies, enabling a number of measurements of galaxy properties and cosmological redshift-space distortions (RSD). Because the measurements were made using one-pass of the VIMOS instrument on the Very Large Telescope (VLT), the galaxies observed only represent approximately 47\% of the parent target sample, with a distribution imprinted with the pattern of the VIMOS slitmask. Correcting for the effect on clustering has previously been achieved using an approximate approach developed using mock catalogues. Pairwise inverse probability (PIP) weighting has recently been proposed by Bianchi & Percival to correct for missing galaxies, and we apply it to mock VIPERS catalogues to show that it accurately corrects the clustering for the VIMOS effects, matching the clustering measured from the observed sample to that of the parent. We then apply PIP-weighting to the VIPERS data, and fit the resulting monopole and quadrupole moments of the galaxy two-point correlation function with respect to the line-of-sight, making measurements of RSD. The results are close to previous measurements, showing that the previous approximate methods used by the VIPERS team are sufficient given the errors obtained on the RSD parameter.

Galaxy clusters trace the highest density peaks in the large-scale structure of the Universe. Their clustering provides a powerful probe that can be exploited in combination with cluster mass measurements to strengthen the cosmological constraints provided by cluster number counts. We investigate the spatial properties of a homogeneous sample of X-ray selected galaxy clusters from the XXL survey, the largest programme carried out by the XMM-Newton satellite. The measurements are compared to $\Lambda$-cold dark matter predictions, and used in combination with self-calibrated mass scaling relations to constrain the effective bias of the sample, $b_{eff}$, and the matter density contrast, $\Omega_{\rm M}$. We measured the angle-averaged two-point correlation function of the XXL cluster sample. The analysed catalogue consists of $182$ X-ray selected clusters from the XXL second data release, with median redshift $\langle z \rangle=0.317$ and median mass $\langle M_{500} \rangle\simeq1.3\cdot10^{14} M_\odot$. A Markov chain Monte Carlo analysis is performed to extract cosmological constraints using a likelihood function constructed to be independent of the cluster selection function. Modelling the redshift-space clustering in the scale range $10<r\,[$\Mpch$]<40$, we obtain $\Omega_{\rm M}=0.27_{-0.04}^{+0.06}$ and $b_{eff}=2.73_{-0.20}^{+0.18}$. This is the first time the two-point correlation function of an X-ray selected cluster catalogue at such relatively high redshifts and low masses has been measured. The XXL cluster clustering appears fully consistent with standard cosmological predictions. The analysis presented in this work demonstrates the feasibility of a cosmological exploitation of the XXL cluster clustering, paving the way for a combined analysis of XXL cluster number counts and clustering.

Galaxy clusters form at the highest density nodes of the cosmic web. The clustering of massive halos is enhanced relative to the general mass distribution and matter beyond the virial region is strongly correlated to the halo mass (halo bias). Clustering can be further enhanced depending on halo properties other than mass (secondary bias). The questions of how much and why the regions surrounding rich clusters are over-dense are still unanswered. Here, we report the analysis of the environment bias in a sample of very massive clusters, selected through the Sunyaev-Zel'dovich effect by the Planck mission. We present the first detection of the correlated dark matter associated to a single cluster, PSZ2 G099.86+58.45. The system is extremely rare in the current paradigm of structure formation. The gravitational lensing signal was traced up to 30 megaparsecs with high signal-to-noise ratio ~3.4. The measured shear is very large and points at environment matter density in notable excess of the cosmological mean. The boosting of the correlated dark matter density around high mass halos can be very effective. Together with ensemble studies of the large scale structure, lensing surveys can picture the surroundings of single haloes.

In this paper, we study the statistical properties of weak lensing peaks in light-cones generated from cosmological simulations. In order to assess the prospects of such observable as a cosmological probe, we consider simulations that include interacting Dark Energy (hereafter DE) models with coupling term between DE and Dark Matter. Cosmological models that produce a larger population of massive clusters have more numerous high signal-to-noise peaks; among models with comparable numbers of clusters those with more concentrated haloes produce more peaks. The most extreme model under investigation shows a difference in peak counts of about $20\%$ with respect to the reference $\mathrm{\Lambda}$CDM model. We find that peak statistics can be used to distinguish a coupling DE model from a reference one with the same power spectrum normalisation. The differences in the expansion history and the growth rate of structure formation are reflected in their halo counts, non-linear scale features and, through them, in the properties of the lensing peaks. For a source redshift distribution consistent with the expectations of future space-based wide field surveys, we find that typically seventy percent of the cluster population contributes to weak-lensing peaks with signal-to-noise ratios larger than two, and that the fraction of clusters in peaks approaches one-hundred percent for haloes with redshift z$\leq$0.5. Our analysis demonstrates that peak statistics are an important tool for disentangling DE models by accurately tracing the structure formation processes as a function of the cosmic time.

The XMM-XXL Survey spans two fields of $\rm 25$ deg$^2$ each observed for more than 6Ms with XMM, which provided a sample of tens of thousands of point sources with a flux limit of $\sim 2.2 \times 10^{-15}$ and $\sim 1.4 \times 10^{-14}$ erg s$^{-1}$ cm$^{2}$, corresponding to 50% of the area curve, in the soft band and hard band, respectively. In this paper we present the spatial clustering properties of $\sim 3100$ and $\sim 1900$ X-ray active galactic nuclei (AGNs) in the 0.5-2 and 2-10 keV bands, respectively, which have been spectroscopically observed with the AAOmega facility. This sample is 90% redshift complete down to an optical magnitude limit of $r\lesssim 21.8$. The sources span the redshift interval $0<z<5.2$, although in the current analysis we limit our samples to $z\le 3$, with corresponding sample median values of $\bar{z}\simeq 0.96$ and 0.79 for the soft band and hard band, respectively. We employ the projected two-point correlation function to infer the spatial clustering and find a correlation length $r_0=7.0 (\pm 0.34)$ and 6.42$(\pm 0.42)$ $h^{-1}$ Mpc, respectively, for the soft- and hard-band detected sources with a slope for both cases of $\gamma=1.44 (\pm 0.1)$. The power-law clustering was detected within comoving separations of 1 and $\sim$25 $h^{-1}$ Mpc. These results, as well as those derived in two separate redshift ranges, provide bias factors of the corresponding AGN host dark matter halos that are consistent with a halo mass of $\log_{10} [M_h/(h^{-1} M_{\odot})]=13.04\pm 0.06$, confirming the results of most recent studies based on smaller X-ray AGN samples.

We study the capability of future measurements of the galaxy clustering power spectrum to probe departures from a power-law spectrum for primordial fluctuations. On considering the information from the galaxy clustering power spectrum up to quasi-linear scales, i.e. $k<0.1$ h Mpc$^{-1}$, we present forecasts for DESI, Euclid and SPHEREx in combination with CMB measurements. As examples of departures in the primordial power spectrum from a simple power-law, we consider four $Planck$ 2015 best-fits motivated by inflationary models with different breaking of the slow-roll approximation. These four representative models provide an improved fit to CMB temperature anisotropies, although not at statistical significant level. As for other extensions in the matter content of the simplest $\Lambda$CDM model, the complementarity of the information in the resulting matter power spectrum expected from these galaxy surveys and in the primordial power spectrum from CMB anisotropies can be effective in constraining cosmological models. We find that the three galaxy surveys can add significant information to CMB to better constrain the extra parameters of the four models considered.

We investigate the possibility of using future photometric and radio surveys to constrain the power spectrum of primordial fluctuations that is predicted by inflationary models with a violation of the slow-roll phase. We forecast constraints with a Fisher analysis on the amplitude of the parametrized features on ultra-large scales, in order to assess whether these could be distinguishable over the cosmic variance. We find that the next generation of photometric and radio surveys has the potential to test these models at a sensitivity better than current CMB experiments and that the synergy between galaxy and CMB observations is able to constrain models with many extra parameters. In particular, an SKA continuum survey with a huge sky coverage and a flux threshold of a few $\mu$Jy could confirm the presence of a new phase in the early Universe at more than 3$\sigma$.

Using an unconventional single line diagnostic that unambiguously identifies AGNs in composite galaxies we report statistical differences in the properties (stellar age, [OII] luminosity, colour) between active and inactive galaxies at 0.62<z<1.2 extracted from the VIMOS Public Extragalactic Redshift Survey (VIPERS). The nuclear activity is probed by the high-ionization [NeV] emission line and along with their parent samples, the galaxies are properly selected according to their stellar mass, redshift, and colour distributions. We report younger underlying stellar ages and higher [OII] luminosities of active galaxies in the green valley and in the blue cloud compared to control samples. We observe higher fractions of green galaxies hosting AGN activity at progressively bluer (r-K) colours. Depending on the location of the host galaxy in the NUVrK colour diagram we find higher AGN fractions in massive blue galaxies and in the least massive red galaxies, in agreement with the picture that black holes vary their properties when hosted in either star-forming or passive galaxies. Exactly where the fast quenching processes are expected to play a role, we identify a novel class of active galaxies in the blue cloud with signatures typical for a suddenly suppression of their star formation activity after a burst happening in the recent past. Their optical spectra resemble those of post-starburst galaxies, that would never be identified in a spectroscopic search using classical selection techniques. Broadly, these active galaxies selected on the [NeV] line are not commonly represented in shallow X-ray, mid-IR, or classical line diagnostics. If we consider that our results are limited by the shallow observational limits and rapid AGN variability, the impact of AGN feedback on galaxy formation and evolution may represent an important channel of fast-transiting galaxies moving to the red sequence.

We use the VIPERS final data release to investigate the performance of colour-selected populations of galaxies as tracers of linear large-scale motions. We empirically select volume-limited samples of blue and red galaxies as to minimise the systematic error on the estimate of the growth rate $f\sigma_8$ from the anisotropy of the two-point correlation function. To this end, rather than rigidly splitting the sample into two colour classes we define the red/blue fractional contribution of each object through a weight based on the $(U-V)$ colour distribution. Using mock surveys that are designed to reproduce the observed properties of VIPERS galaxies, we find the systematic error in recovering the fiducial value of $f\sigma_8$ to be minimized when using a volume-limited sample of luminous blue galaxies. We model non-linear corrections via the Scoccimarro extension of the Kaiser model, finding systematic errors on $f\sigma_8$ of below $1-2\%$, using scales as small as 5 $h^{-1}\mathrm{Mpc}$. We interpret this result as indicating that selection of luminous blue galaxies maximises the fraction that are central objects in their dark matter haloes; this in turn minimises the contribution to the measured $\xi(r_p,\pi)$ from the 1-halo term, which is dominated by non-linear motions. The gain is inferior if one uses the full magnitude-limited sample of blue objects, consistent with the presence of a significant fraction of blue, fainter satellites dominated by non-streaming, orbital velocities. We measure a value of $f\sigma_8=0.45 \pm 0.11$ over the single redshift range $0.6\le z\le 1.0$, corresponding to an effective redshift for the blue galaxies $\left<z\right>=0.85$. Including in the likelihood the potential extra information contained in the blue-red galaxy cross-correlation function does not lead to an appreciable improvement in the error bars, while it increases the systematic error.

We present AMICO (Adaptive Matched Identifier of Clustered Objects), a new algorithm for the detection of galaxy clusters in photometric surveys. AMICO is based on the Optimal Filtering technique, which allows to maximise the signal-to-noise ratio of the clusters. In this work we focus on the new iterative approach to the extraction of cluster candidates from the map produced by the filter. In particular, we provide a definition of membership probability for the galaxies close to any cluster candidate, which allows us to remove its imprint from the map, allowing the detection of smaller structures. As demonstrated in our tests, this method allows the deblending of close-by and aligned structures in more than $50\%$ of the cases for objects at radial distance equal to $0.5 \times R_{200}$ or redshift distance equal to $2 \times \sigma_z$, being $\sigma_z$ the typical uncertainty of photometric redshifts. Running AMICO on mocks derived from N-body simulations and semi-analytical modelling of the galaxy evolution, we obtain a consistent mass-amplitude relation through the redshift range $0.3 < z < 1$, with a logarithmic slope $\sim 0.55$ and a logarithmic scatter $\sim 0.14$. The fraction of false detections is steeply decreasing with S/N, and negligible at S/N > 5.

We carry out a joint analysis of redshift-space distortions and galaxy-galaxy lensing, with the aim of measuring the growth rate of structure; this is a key quantity for understanding the nature of gravity on cosmological scales and late-time cosmic acceleration. We make use of the final VIPERS redshift survey dataset, which maps a portion of the Universe at a redshift of $z \simeq 0.8$, and the lensing data from the CFHTLenS survey over the same area of the sky. We build a consistent theoretical model that combines non-linear galaxy biasing and redshift-space distortion models, and confront it with observations. The two probes are combined in a Bayesian maximum likelihood analysis to determine the growth rate of structure at two redshifts $z=0.6$ and $z=0.86$. We obtain measurements of $f\sigma_8(0.6) = 0.48 \pm 0.12$ and $f\sigma_8(0.86) = 0.48 \pm 0.10$. The additional galaxy-galaxylensing constraint alleviates galaxy bias and $\sigma_8$ degeneracies, providing direct measurements of $[f(0.6),\sigma_8(0.6)] = [0.93 \pm 0.22, 0.52 \pm 0.06]$ and $f(0.86),\sigma_8(0.86)] = [0.99 \pm 0.19, 0.48 \pm 0.04]$. These measurements are statistically consistent with a Universe where the gravitational interactions can be described by General Relativity, although they are not yet accurate enough to rule out some commonly considered alternatives. Finally, as a complementary test we measure the gravitational slip parameter, $E_G$ , for the first time at $z>0.6$. We find values of $\smash{\overline{E}_G}(0.6) = 0.16 \pm 0.09$ and $\smash{\overline{E}_G}(0.86) = 0.09 \pm 0.07$, when $E_G$ is averaged over scales above $3 h^{-1} \rm{Mpc}$. We find that our $E_G$ measurements exhibit slightly lower values than expected for standard relativistic gravity in a {\Lambda}CDM background, although the results are consistent within $1-2\sigma$.

The large-scale structure of the universe can only be observed via luminous tracers of the dark matter. However, the clustering statistics of tracers are biased and depend on various properties, such as their host-halo mass and assembly history. On very large scales this tracer bias results in a constant offset in the clustering amplitude, known as linear bias. Towards smaller nonlinear scales, this is no longer the case and tracer bias becomes a complicated function of scale and time. We focus on tracer bias centred on cosmic voids, depressions of the density field that spatially dominate the universe. We consider three types of tracers: galaxies, galaxy clusters and AGN, extracted from the hydrodynamical simulation Magneticum Pathfinder. In contrast to common clustering statistics that focus on auto-correlations of tracers, we find that void-tracer cross-correlations are successfully described by a linear-bias relation. The tracer-density profile of voids can thus be related to their matter-density profile by a single number. We show that it coincides with the linear tracer bias extracted from the large-scale auto-correlation function and expectations from theory, if sufficiently large voids are considered. For smaller voids we observe a shift towards higher values. This has important consequences on cosmological parameter inference, as the problem of unknown tracer bias is alleviated up to a constant number. The smallest scales in existing datasets become accessible to simpler models, providing numerous modes of the density field that have been disregarded so far, but may help to further reduce statistical errors in constraining cosmology.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We use the final catalogue of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to measure the power spectrum of the galaxy distribution at high redshift, presenting results that extend beyond $z=1$ for the first time. We apply an FFT technique to four independent sub-volumes comprising a total of $51,728$ galaxies at $0.6<z<1.1$ (out of the nearly $90,000$ included in the whole survey). We concentrate here on the shape of the direction-averaged power spectrum in redshift space, explaining the level of modelling of redshift-space anisotropies and the anisotropic survey window function that are needed to deduce this in a robust fashion. We then use covariance matrices derived from a large ensemble of mock datasets in order to fit the spectral data. The results are well matched by a standard $\Lambda$CDM model, with density parameter $\Omega_M h =\smash{0.227^{+0.063}_{-0.050}}$ and baryon fraction $\smash{f_B=\Omega_B/\Omega_M=0.220^{+0.058}_{-0.072}}$. These inferences from the high-$z$ galaxy distribution are consistent with results from local galaxy surveys, and also with the Cosmic Microwave Background. Thus the $\Lambda$CDM model gives a good match to cosmic structure at all redshifts so far accessible to observational study.

Full ray-tracing maps of gravitational lensing, constructed from N-Body simulations, represent a fundamental tool to interpret present and future weak lensing data. However the limitation of computational resources and storage capabilities severely restrict the number of realizations that can be performed in order to accurately sample both the cosmic shear models and covariance matrices. In this paper we present a halo model formalism for weak gravitational lensing that alleviates these issues by producing weak-lensing mocks at a reduced computational cost. Our model takes as input the halo population within a desired light-cone and the linear power spectrum of the underlined cosmological model. We examine the contribution given by the presence of substructures within haloes to the cosmic shear power spectrum and quantify it to the percent level. Our method allows us to reconstruct high-resolution convergence maps, for any desired source redshifts, of light-cones that realistically trace the matter density distribution in the universe, account for masked area and sample selections. We compare our analysis on the same large scale structures constructed using ray-tracing techniques and find very good agreements both in the linear and non-linear regimes up to few percent levels. The accuracy and speed of our method demonstrate the potential of our halo model for weak lensing statistics and the possibility to generate a large sample of convergence maps for different cosmological models as needed for the analysis of large galaxy redshift surveys.

The three-point correlation function (3PCF) is a powerful probe to investigate the clustering of matter in the Universe in a complementary way with respect to lower-order statistics, providing additional information with respect to the two-point correlation function and allowing us to shed light on biasing, nonlinear processes, and deviations from Gaussian statistics. In this paper, we analyse the first data release of the VIMOS Public Extragalactic Redshift Survey (VIPERS), determining the dependence of the three-point correlation function on luminosity and stellar mass at $z=[0.5,1.1]$. We exploit the VIPERS Public Data Release 1, consisting of more than 50,000 galaxies with B-band magnitudes in the range $-21.6\lesssim M_{\rm B}-5\log(h)\lesssim-19.9$ and stellar masses in the range $9.8\lesssim\log(M_\star[h^{-2}\,M_\odot])\lesssim 10.7$. We measure both the connected 3PCF and the reduced 3PCF in redshift space, probing different configurations and scales, in the range $2.5<r\,$[Mpc/h]$<20$. We find a significant dependence of the reduced 3PCF on scales and triangle shapes, with stronger anisotropy at larger scales ($r\sim10$ Mpc/h) and an almost flat trend at smaller scales, $r\sim2.5$ Mpc/h. Massive and luminous galaxies present a larger connected 3PCF, while the reduced 3PCF is remarkably insensitive to magnitude and stellar masses in the range we explored. These trends, already observed at low redshifts, are confirmed for the first time to be still valid up to $z=1.1$, providing support to the hierarchical scenario for which massive and bright systems are expected to be more clustered. The possibility of using the measured 3PCF to provide independent constraints on the linear galaxy bias $b$ has also been explored, showing promising results in agreement with other probes.

Redshift-space clustering anisotropies caused by cosmic peculiar velocities provide a powerful probe to test the gravity theory on large scales. However, to extract unbiased physical constraints, the clustering pattern has to be modelled accurately, taking into account the effects of non-linear dynamics at small scales, and properly describing the link between the selected cosmic tracers and the underlying dark matter field. We use a large hydrodynamic simulation to investigate how the systematic error on the linear growth rate, $f$, caused by model uncertainties, depends on sample selections and comoving scales. Specifically, we measure the redshift-space two-point correlation function of mock samples of galaxies, galaxy clusters and Active Galactic Nuclei, extracted from the Magneticum simulation, in the redshift range 0.2 < z < 2, and adopting different sample selections. We estimate $f\sigma_8$ by modelling both the monopole and the full two-dimensional anisotropic clustering, using the dispersion model. We find that the systematic error on $f\sigma_8$ depends significantly on the range of scales considered for the fit. If the latter is kept fixed, the error depends on both redshift and sample selection, due to the scale-dependent impact of non-linearities, if not properly modelled. On the other hand, we show that it is possible to get unbiased constraints on $f\sigma_8$ provided that the analysis is restricted to a proper range of scales, that depends non trivially on the properties of the sample. This can have a strong impact on multiple tracers analyses, and when combining catalogues selected at different redshifts.

[Abridged] We use the final data of the VIMOS Public Extragalactic Redshift Survey (VIPERS) to investigate the effect of environment on the evolution of galaxies between $z=0.5$ and $z=0.9$. We characterise local environment in terms of the density contrast smoothed over a cylindrical kernel, the scale of which is defined by the distance to the $5^{th}$ nearest neighbour. We find that more massive galaxies tend to reside in higher-density environments over the full redshift range explored. Defining star-forming and passive galaxies through their (NUV$-r$) vs ($r-K$) colours, we then quantify the fraction of star-forming over passive galaxies, $f_{\rm ap}$, as a function of environment at fixed stellar mass. $f_{\rm ap}$ is higher in low-density regions for galaxies with masses ranging from $\log(\mathcal{M}/\mathcal{M}_\odot)=10.38$ (the lowest value explored) to at least $\log(\mathcal{M}/\mathcal{M}_\odot)\sim11.3$, although with decreasing significance going from lower to higher masses. This is the first time that environmental effects on high-mass galaxies are clearly detected at redshifts as high as $z\sim0.9$. We compared these results to VIPERS-like galaxy mock catalogues based on the galaxy formation model of De Lucia & Blaizot. The model correctly reproduces $f_{\rm ap}$ in low-density environments, but underpredicts it at high densities. The discrepancy is particularly strong for the lowest-mass bins. We find that this discrepancy is driven by an excess of low-mass passive satellite galaxies in the model. Looking at the accretion history of these model galaxies, i.e. the times when they become satellites, a better (yet not perfect) agreement with observations can be obtained in high density regions by assuming either that a not-negligible fraction of satellites is destroyed, or that their quenching time-scale is longer than $\sim 2$ Gyr.

We present measurements of the growth rate of cosmological structure from the modelling of the anisotropic galaxy clustering measured in the final data release of the VIPERS survey. The analysis is carried out in configuration space and based on measurements of the first two even multipole moments of the anisotropic galaxy auto-correlation function, in two redshift bins spanning the range $0.5 < z < 1.2$. We provide robust and cosmology-independent corrections for the VIPERS angular selection function, allowing recovery of the underlying clustering amplitude at the percent level down to the Mpc scale. We discuss several improvements on the non-linear modelling of redshift-space distortions (RSD) and perform detailed tests of a variety of approaches against a set of realistic VIPERS-like mock realisations. This includes using novel fitting functions to describe the velocity divergence and density power spectra $P_{\theta\theta}$ and $P_{\delta\theta}$ that appear in RSD models. These tests show that we are able to measure the growth rate with negligible bias down to separations of $5h^{-1}Mpc$. Interestingly, the application to real data shows a weaker sensitivity to the details of non-linear RSD corrections compared to mock results. We obtain consistent values for the growth rate times the matter power spectrum normalisation parameter of $f\sigma_8=0.55\pm 0.12$ and $0.40\pm0.11$ at effective redshifts of $z = 0.6$ and $z=0.86$ respectively. These results are in agreement with standard cosmology predictions assuming Einstein gravity in a $\Lambda \rm{CDM}$ background.

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

Identifying spurious reduction artefacts in galaxy spectra is a challenge for large surveys. We present an algorithm for identifying and repairing residual spurious features in sky-subtracted galaxy spectra with application to the VIPERS survey. The algorithm uses principal component analysis (PCA) applied to the galaxy spectra in the observed frame to identify sky line residuals imprinted at characteristic wavelengths. We further model the galaxy spectra in the rest-frame using PCA to estimate the most probable continuum in the corrupted spectral regions, which are then repaired. We apply the method to 90,000 spectra from the VIPERS survey and compare the results with a subset where careful editing was performed by hand. We find that the automatic technique does an extremely good job in reproducing the time-consuming manual cleaning and does it in a uniform and objective manner across a large data sample. The mask data products produced in this work are released together with the VIPERS second public data release (PDR-2).

We use the full VIPERS redshift survey in combination with SDSS-DR7 to explore the relationships between star-formation history (using d4000), stellar mass and galaxy structure, and how these relationships have evolved since z~1. We trace the extents and evolutions of both the blue cloud and red sequence, by fitting double Gaussians to the d4000 distribution of galaxies in narrow stellar mass bins, for four redshift intervals over 0<z<1. This reveals downsizing in star formation, as the high-mass limit of the blue cloud retreats steadily with time from M*~10^11.2 M_sun at z~0.9 to M*~10^10.7 M_sun by the present day. The number density of massive blue-cloud galaxies (M*>10^11 M_sun, d4000<1.55) drops sharply by a factor five between z~0.8 and z~0.5. These galaxies are becoming quiescent at a rate that largely matches the increase in the numbers of massive passive galaxies seen over this period. We examine the size-mass relation of blue cloud galaxies, finding that its high-mass boundary runs along lines of constant M*/r_e or equivalently inferred velocity dispersion. Larger galaxies can continue to form stars to higher stellar masses than smaller galaxies. As blue cloud galaxies approach this high-mass limit, they start to be quenched, their d4000 values increasing to push them towards the green valley. In parallel, their structures change, showing higher Sersic indices and central stellar mass densities. For these galaxies, bulge growth is necessary for them to reach the high-mass limit of the blue cloud and be quenched by internal mechanisms. The blue cloud galaxies that are being quenched at z~0.8 lie along the same size-mass relation as present day quiescent galaxies, and seem the likely progenitors of today's S0s.

We use the final data from the VIPERS redshift survey to extract an unparalleled sample of more than 2000 massive M > 10^11 M_sun passive galaxies (MPGs) at redshift 0.5 < z < 1.0, based on their NUVrK colours. This enables us to investigate how the population of these objects was built up over cosmic time. We find that the evolution of the number density depends on the galaxy mean surface stellar mass density, Sigma. In particular, dense (Sigma > 2000 M_sun pc^-2) MPGs show a constant comoving number density over this redshift range, whilst this increases by a factor ~ 4 for the least dense objects, defined as having Sigma < 1000 M_sun pc^-2. We estimate stellar ages for the MPG population both fitting the Spectral Energy Distribution (SED) and through the D4000_n index, obtaining results in good agreement. Our findings are consistent with passive ageing of the stellar content of dense MPGs. We show that at any redshift the less dense MPGs are younger than dense ones and that their stellar populations evolve at a slower rate than predicted by passive evolution. This points to a scenario in which the overall population of MPGs was built up over the cosmic time by continuous addition of less dense galaxies: on top of an initial population of dense objects that passively evolves, new, larger, and younger MPGs continuously join the population at later epochs. Finally, we demonstrate that the observed increase in the number density of MPGs is totally accounted for by the observed decrease in the number density of correspondingly massive star forming galaxies (i.e. all the non-passive M > 10^11 M_sun objects). Such systems observed at z ~ 1 in VIPERS, therefore, represent the most plausible progenitors of the subsequent emerging class of larger MPGs.

We identified voids in the completed VIMOS Public Extragalactic Redshift Survey (VIPERS), using an algorithm based on searching for empty spheres. We measured the cross-correlation between the centres of voids and the complete galaxy catalogue. The cross-correlation function exhibits a clear anisotropy in both VIPERS fields (W1 and W4), which is characteristic of linear redshift space distortions. By measuring the projected cross-correlation and then deprojecting it we are able to estimate the undistorted cross-correlation function. We propose that given a sufficiently well measured cross-correlation function one should be able to measure the linear growth rate of structure by applying a simple linear Gaussian streaming model for the redshift space distortions (RSD). Our study of voids in 306 mock galaxy catalogues mimicking the VIPERS fields would suggest that VIPERS is capable of measuring $\beta$ with an error of around $25\%$. Applying our method to the VIPERS data, we find a value for the redshift space distortion parameter, $\beta = 0.423^{+0.104}_{-0.108}$, which given the bias of the galaxy population we use gives a linear growth rate of $f\sigma_8 = 0.296^{+0.075}_{-0.078}$ at $z = 0.727$. These results are consistent with values observed in parallel VIPERS analysis using standard techniques.

We present the full public data release (PDR-2) of the VIMOS Public Extragalactic Redshift Survey (VIPERS), performed at the ESO VLT. We release redshifts, spectra, CFHTLS magnitudes and ancillary information (as masks and weights) for a complete sample of 86,775 galaxies (plus 4,732 other objects, including stars and serendipitous galaxies); we also include their full photometrically-selected parent catalogue. The sample is magnitude limited to i_AB < 22.5, with an additional colour-colour pre-selection devised as to exclude galaxies at z<0.5. This practically doubles the effective sampling of the VIMOS spectrograph over the range 0.5<z<1.2 (reaching 47% on average), yielding a final median local galaxy density close to 5 10^{-3} h^3 Mpc}^{-3}. The total area spanned by the final data set is ~ 23.5 deg^2, corresponding to 288 VIMOS fields with marginal overlaps, split over two regions within the CFHTLS-Wide W1 and W4 equatorial fields (at R.A. ~2 and ~22 hours, respectively). Spectra were observed at a resolution R=220, covering a wavelength range 5500-9500 Angstrom. Data reduction and redshift measurements were performed through a fully automated pipeline; all redshift determinations were then visually validated and assigned a quality flag. Measurements with a quality flag >= 2 are shown to have a confidence level of 96% or larger and make up 88% of all measured galaxy redshifts (76,552 out of 86,775), constituting the VIPERS prime catalogue for statistical investigations. For this sample the rms redshift error, estimated using repeated measurements of about 3,000 galaxies, is found to be sigma_z = 0.00054(1+z). All data are available at http://vipers.inaf.it and on the ESO Archive.

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.

We trace the evolution and the star formation history of passive galaxies, using a subset of the VIMOS Public Extragalactic Redshift Survey (VIPERS). We extracted from the VIPERS survey a sample of passive galaxies in the redshift range 0.4<z<1.0 and stellar mass range 10<$log(M_{star}/M_{\odot})$<12. The sample was selected using an evolving cut in the rest-frame U-V color distribution and additional quality-ensuring cuts. We use the stacked spectra to measure the 4000$\AA$ break (D4000) and the $H\delta$ Lick index ($H\delta_{A}$) with high precision. We compare the results with a grid of synthetic spectra to constrain the star formation epochs of these galaxies. We characterize the formation redshift-stellar mass relation for intermediate-redshift passive galaxies. We find that at $z\sim1$ stellar populations in low-mass passive galaxies are younger than in high-mass passive galaxies, similarly to what is observed at the present epoch. Over the full analyzed redshift and stellar mass range, the $D4000$ index increases with redshift, while $H\delta_{A}$ gets lower. This implies that the stellar populations are getting older with increasing stellar mass. Comparison to the spectra of passive galaxies in the SDSS survey shows that the shape of the relations of $D4000$, and $H\delta_{A}$ with stellar mass has not changed significantly with redshift. Assuming a single burst formation, this implies that high-mass passive galaxies formed their stars at $z_{form}\sim2$, while low-mass galaxies formed their main stellar population more recently, at $z_{form}\sim1$. The consistency of these results, obtained using two independent estimator of the formation redshift ($D4000$ and $H\delta_{A}$), further strengthens a scenario in which star formation proceeds from higher- to lower-mass systems as time passes, i.e. what has become known as the 'downsizing' picture.

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

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.

We explore the accuracy of the clustering-based redshift estimation proposed by M\'enard 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.

We explore the evolution of the statistical distribution of galaxy morphological properties and colours over the redshift range $0.5<z<1$, combining high-quality imaging data from the CFHT Legacy Survey with the large number of redshifts and extended photometry from the VIPERS survey. Galaxy structural parameters are measured by fitting S\'ersic profiles to $i$-band images and then combined with absolute magnitudes, colours and redshifts, to trace the evolution in a multi-parameter space. We analyse, using a new method, the combination of colours and structural parameters of early- and late-type galaxies in luminosity--redshift space. We found that both the rest-frame colour distributions in the (U-B) vs. (B-V) plane and the S\'ersic index distributions are well fitted by a sum of two Gaussians, with a remarkable consistency of red-spheroidal and blue-disky galaxy populations, over the explored redshift ($0.5<z<1$) and luminosity ($-1.5<B-B_*<1.0$) ranges. The combination of the UBV rest-frame colour and S\'ersic index $n$ as a function of redshift and luminosity allows us to present the structure of early- and late-type galaxies and their evolution. We found that early type galaxies display only a slow change of their concentrations since $z\sim1$; it is already established by $z\sim1$ and depends much more strongly on their luminosities. In contrast, late-type galaxies get clearly more concentrated with cosmic time since $z\sim1$, with only little evolution in colour, which remains dependent mainly on their luminosity. This flipped luminosity (mass) and redshift dependence likely reflects different evolutionary tracks of early- and late-type galaxies before and after $z\sim1$. The combination of rest-frame colours and S\'ersic index $n$ as a function of redshift and luminosity leads to a precise statistical description of the structure of galaxies and their evolution.

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.

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

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

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

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

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 $\sigma_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]

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

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

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

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

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\sigma_8 = 0.38$ with 18% uncertainty at redshift 0.7.

We investigate how galaxies in VIPERS (the VIMOS Public Extragalactic Redshift Survey) inhabit the cosmological density field by examining the correlations across the observable parameter space of galaxy properties and clustering strength. The high-dimensional analysis is made manageable by the use of group-finding and regression tools. We find that the major trends in galaxy properties can be explained by a single parameter related to stellar mass. After subtracting this trend, residual correlations remain between galaxy properties and the local environment pointing to complex formation dependencies. As a specific application of this work we build subsamples of galaxies with specific clustering properties for use in cosmological tests.

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 $\sigma_{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{\sigma}-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=\sigma_{8g}/\sigma_{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.

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

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 \Delta f_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter \gamma, which signals deviations from GR, will be constrained to Delta \gamma=0.02, and the neutrino density parameter to Delta Omega_\nu=0.0013 (or Delta \sum m_\nu=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]

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.

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

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

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

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

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.

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

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

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

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.

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

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 \sigma_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.

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

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.

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

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

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

We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $\Lambda$CDM cosmological scenario - namely an $f(R)$ theory of Modified Gravity (MG) and a cosmological background of massive neutrinos - with the aim to investigate their possible observational degeneracies. We focus on three basic statistics of the large-scale matter distribution, more specifically the nonlinear matter power spectrum, the halo mass function, and the halo bias, for which we determine the deviation with respect to the fiducial $\Lambda$CDM cosmology in the context of both separate and combined simulations of $f(R)$ MG and massive neutrinos scenarios. Our results show that while these two extended models separately determine very prominent and potentially detectable features in all the three statistics, when we allow them to be simultaneously at work these features are strongly suppressed, resulting in much weaker deviations from the standard model's predictions. In particular, when an $f(R)$ gravity model with $f_{R0}=-1\times 10^{-4}$ is combined with a total neutrino mass of $\Sigma _{i}m_{\nu _{i}}=0.4$ eV, the resulting matter power spectrum, halo mass function, and bias at z=0 are found to be consistent with the standard model's predictions at the 10%, 20%, and 5% accuracy levels, respectively. Therefore, our results imply an intrinsic theoretical limit to the effective discriminating power of present and future observational data sets with respect to these widely considered extensions of the standard cosmological scenario in the absence of independent measurements of the neutrino masses from laboratory experiments, even though the high-redshift evolution might still allow to partially break the degeneracy [Abridged].

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.

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} $\eta_{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, $\xi(nR)$, and its variance $\sigma^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 $\eta_{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 $\eta_{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 $\Omega_{m,0}=0.270_{-0.025}^{+0.029}$. In addition to the great precision achieved on our estimation of $\Omega_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 $\Omega_{m,0}=0.274\pm0.017$ which highlights the great consistency between our estimation and other cosmological probes such as BAOs, CMB and Supernovae.

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.

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

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]

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.

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.

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

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\sigma_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.

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)

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.

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]

We investigate the impact of statistical and systematic errors on measurements of linear redshift-space distortions (RSD) in future cosmological surveys, analyzing large catalogues of dark-matter halos from the BASICC simulation. These allow us to estimate the dependence of errors on typical survey properties, as volume, galaxy density and mass (i.e. bias factor) of the adopted tracer. We find that measures of the specific growth rate \beta=f/b using the Hamilton/Kaiser harmonic expansion of the redshift-space correlation function \xi(r_p,\pi) on scales larger than 3/h Mpc are typically under-estimated by up to 10% for galaxy sized halos. This is significantly larger than the corresponding statistical errors, which amount to a few percent, indicating the importance of non-linear improvements to the Kaiser model to obtain accurate measurements of the growth rate. We compare the statistical errors to predictions obtained with the Fisher information matrix, based on the usual FKP prescription for the errors on the power spectrum. We show that this produces parameter errors fairly similar to the standard deviations from the halo catalogues, but only if applied to strictly linear scales in Fourier space (k<0.2 h/Mpc). Finally, we present an accurate scaling formula describing the relative error on {\beta} as a function of the survey parameters, which closely matches the simulation results in all explored regimes. This provides a handy and plausibly more realistic alternative to the Fisher matrix approach, to quickly and accurately predict RSD statistical errors expected from future surveys.

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

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

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

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

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

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

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.

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)

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 $\sigma_8$, are at the level of $\Delta f_\mathrm{NL} \sim 100$. We consider three additional bispectrum shapes, for which the cosmic shear constraints range from $\Delta f_\mathrm{NL}\sim 340$ (equilateral shape) up to $\Delta 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 $\sigma_8$ and $f_\mathrm{NL}$ being improved by a factor of $\sim 3$ upon their combination. (abridged)

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

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

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.

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

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

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

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

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

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

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

Cosmological neutrinos strongly affect the evolution of the largest structures in the Universe, i.e. galaxies and galaxy clusters. We use large box-size full hydrodynamic simulations to investigate the non-linear effects that massive neutrinos have on the spatial properties of cold dark matter (CDM) haloes. We quantify the difference with respect to the concordance LambdaCDM model of the halo mass function and of the halo two-point correlation function. We model the redshift-space distortions and compute the errors on the linear distortion parameter beta introduced if cosmological neutrinos are assumed to be massless. We find that, if not taken correctly into account and depending on the total neutrino mass, these effects could lead to a potentially fake signature of modified gravity. Future nearly all-sky spectroscopic galaxy surveys will be able to constrain the neutrino mass if it is larger than 0.6 eV, using beta measurements alone and independently of the value of the matter power spectrum normalisation. In combination with other cosmological probes, this will strengthen neutrino mass constraints and help breaking parameter degeneracies.

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

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.

(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 $\Lambda$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 $\Lambda$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.

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.

(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 $\Delta f_\mathrm{NL}\sim 300$, while at 90% Confidence Level it is of $\Delta 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 $\Delta f_\mathrm{NL} \sim$ a few tens at 68.3% Confidence Level for an all-sky weak lensing survey.

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.

We assess the possibility to detect the warm-hot intergalactic medium (WHIM) in emission and to characterize its physical conditions and spatial distribution through spatially resolved X-ray spectroscopy, in the framework of the recently proposed DIOS, EDGE, Xenia, and ORIGIN missions, all of which are equipped with microcalorimeter-based detectors. For this purpose we analyze a large set of mock emission spectra, extracted from a cosmological hydrodynamical simulation. These mock X-ray spectra are searched for emission features showing both the OVII K alpha triplet and OVIII Ly alpha line, which constitute a typical signature of the warm hot gas. Our analysis shows that 1 Ms long exposures and energy resolution of 2.5 eV will allow us to detect about 400 such features per deg^2 with a significance >5 sigma and reveals that these emission systems are typically associated with density ~100 above the mean. The temperature can be estimated from the line ratio with a precision of ~20%. The combined effect of contamination from other lines, variation in the level of the continuum, and degradation of the energy resolution reduces these estimates. Yet, with an energy resolution of 7 eV and all these effects taken into account, one still expects about 160 detections per deg^2. These line systems are sufficient to trace the spatial distribution of the line-emitting gas, which constitute an additional information, independent from line statistics, to constrain the poorly known cosmic chemical enrichment history and the stellar feedback processes.

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.

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

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

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 $\Lambda$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).

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.

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

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

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

(abridged) We explore the statistical properties of non-linear cosmic structures in a flat $\Lambda$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 $\alpha_{\mathrm{S},0}$ allowed by an updated combination of CMB anisotropy data and the 2dF Galaxy Redshift Survey. Our analysis constrains $\alpha_{\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.

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, \sigma_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 \Delta f_nl~10 at 1 \sigma 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.

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.

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.

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.

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.

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.

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

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

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

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.

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 --> \delta_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.

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

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.

In the study of the process of cosmic structure formation numerical simulations are crucial tools to interface observational data to theoretical models and to investigate issues otherwise unexplored. Enormous advances have been achieved in the last years thanks to the availability of sophisticated codes, now allowing to tackle the problem of cosmic structure formation and subsequent evolution by covering larger and larger dynamical ranges. Moreover, computational cosmology is the ideal interpretative framework for the overwhelming amount of new data from extragalactic surveys and from large sample of individual objects. The Workshop Novicosmo 2008 "The Impact of Simulations in Cosmology and Galaxy Formation' held in SISSA was aimed at providing the state-of-the-art on the latest numerical simulations in Cosmology and in Galaxy Formation. Particular emphasis was given to the implementation of new physical processes in simulation codes, to the comparison between different codes and numerical schemes and how to use best supercomputing facilities of the next generation. Finally, the impact on our knowledge on the Physics of the Universe brought by this new channel of investigation has also been focused. The Workshop was divided in three sections corresponding (roughly) to three main areas of study: Reionization and Intergalactic medium; Dark and Luminous matter in galaxies; Clusters of galaxies and Large scale Structures. This paper will provide i) a short resume' of the scientific results of the Workshop ii) the complete list of the talks and the instructions on how to retrieve the .pdf of the related (powerpoint) presentations iii) a brief presentation of the associated Exhibition "Space Art"

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.

We model the cosmological co-evolution of galaxies and their central supermassive black holes (BHs) within a semi-analytical framework developed on the outputs of the Millennium Simulation (Croton et al., 2006; De Lucia & Blaizot, 2007). In this work, we analyze the model BH scaling relations, fundamental plane and mass function, and compare them with the most recent observational data. Furthermore, we extend the original code developed by Croton et al. (2006) to follow the evolution of the BH mass accretion and its conversion into radiation, and compare the derived AGN bolometric luminosity function with the observed one. We find, for the most part, a very good agreement between predicted and observed BH properties. Moreover, the model is in good agreement with the observed AGN number density in 0<z<5, provided it is assumed that the cold gas fraction accreted by BHs at high redshifts is larger than at low redshifts (Marulli et al., 2008).

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.

We use semi-analytic modelling of the galaxy-cluster population and its strong lensing efficiency to explore how the expected abundance of large gravitational arcs on the sky depends on $\sigma_8$. Our models take all effects into account that have been shown to affect strong cluster lensing substantially, in particular cluster asymmetry, substructure, merging, and variations in the central density concentrations. We show that the optical depth for long and thin arcs increases by approximately one order of magnitude when $\sigma_8$ increases from 0.7 to 0.9, owing to a constructive combination of several effects. Models with high $\sigma_8$ are also several orders of magnitude more efficient in producing arcs at intermediate and high redshifts. Finally, we use realistic source number counts to quantitatively predict the total number of arcs brighter than several magnitude limits in the R and I bands. We confirm that, while $\sigma_8\sim0.9$ may come close to the known abundance of arcs, even $\sigma_8\sim0.8$ falls short by almost an order of magnitude in reproducing known counts. We conclude that, should $\sigma_8\sim0.8$ be confirmed, we would fail to understand the strong-lensing efficiency of the galaxy cluster population, and in particular the abundance of arcs in high-redshift clusters. We argue that early-dark energy or non-Gaussian density fluctuations may indicate one way out of this problem.

Observations of distant supernovae indicate that the Universe is now in a phase of accelerated expansion the physical cause of which is a mystery. Formally, this requires the inclusion of a term acting as a negative pressure in the equations of cosmic expansion, accounting for about 75 per cent of the total energy density in the Universe. The simplest option for this "dark energy" corresponds to a cosmological constant, perhaps related to the quantum vacuum energy. Physically viable alternatives invoke either the presence of a scalar field with an evolving equation of state, or extensions of general relativity involving higher-order curvature terms or extra dimensions. Although they produce similar expansion rates, different models predict measurable differences in the growth rate of large-scale structure with cosmic time. A fingerprint of this growth is provided by coherent galaxy motions, which introduce a radial anisotropy in the clustering pattern reconstructed by galaxy redshift surveys. Here we report a measurement of this effect at a redshift of 0.8. Using a new survey of more than 10,000 faint galaxies, we measure the anisotropy parameter b = 0.70 +/- 0.26, which corresponds to a growth rate of structure at that time of f = 0.91 +/- 0.36. This is consistent with the standard cosmological-constant model with low matter density and flat geometry, although the error bars are still too large to distinguish among alternative origins for the accelerated expansion. This could be achieved with a further factor-of-ten increase in the sampled volume at similar redshift.

The redshifts of all cosmologically distant sources are expected to experience a small, systematic drift as a function of time due to the evolution of the Universe's expansion rate. A measurement of this effect would represent a direct and entirely model-independent determination of the expansion history of the Universe over a redshift range that is inaccessible to other methods. Here we investigate the impact of the next generation of Extremely Large Telescopes on the feasibility of detecting and characterising the cosmological redshift drift. We consider the Lyman alpha forest in the redshift range 2 < z < 5 and other absorption lines in the spectra of high redshift QSOs as the most suitable targets for a redshift drift experiment. Assuming photon-noise limited observations and using extensive Monte Carlo simulations we determine the accuracy to which the redshift drift can be measured from the Ly alpha forest as a function of signal-to-noise and redshift. Based on this relation and using the brightness and redshift distributions of known QSOs we find that a 42-m telescope is capable of unambiguously detecting the redshift drift over a period of ~20 yr using 4000 h of observing time. Such an experiment would provide independent evidence for the existence of dark energy without assuming spatial flatness, using any other cosmological constraints or making any other astrophysical assumption.

In this paper we investigate the effects that a dynamic dark energy component dominant in the universe at late epochs has on reionization. We follow the evolution of HII regions with the analytic approach of Furlanetto and Oh (2005) in two different universes for which we assume the Peebles and Ratra (2003) and Brax and Martin (2000) quintessence models and we compare our results to the LCDM scenario. We show that, for a fixed ionization efficiency, at the same cosmological epoch the topology of bubbles is dominated by high-mass objects and the characteristic size of the ionized regions is slightly smaller than in the LCDM model, especially at the latest stages of reionization, due to the higher recombination efficiency. As a consequence, the bubbles' `epoch of overlap' happens earlier than in LCDM. Finally, we show how the different evolution of the HII regions affects the transmission of the high-z QSO spectra, reducing the Lyman flux absorption at small optical depths.

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

We present a newly developed code that allows simulations of optical observations of galaxy fields with a variety of instruments. The code incorporates gravitational lensing effects and is targetted at simulating lensing by galaxy clusters. Our goal is to create the tools required for comparing theoretical expectations with observations to obtain a better understanding of how observational noise affects lensing applications such as mass estimates, studies on the internal properties of galaxy clusters and arc statistics. Starting from a set of input parameters, characterizing both the instruments and the observational conditions, the simulator provides a virtual observation of a patch of the sky. It includes several sources of noise such as photon-noise, sky background, seeing, and instrumental noise. Ray-tracing through simulated mass distributions accounts for gravitational lensing. Source morphologies are realistically simulated based on shapelet decompositions of galaxy images retrieved from the GOODS-ACS archive. According to their morphological class, spectral-energy-distributions are assigned to the source galaxies in order to reproduce observations of each galaxy in arbitrary photometric bands. We illustrate our techniques showing virtual observations of a galaxy-cluster core as it would be observed with the space telescope DUNE, which was recently proposed to ESA within its "Cosmic vision" programme. (Abridged)

We model the cosmological co-evolution of galaxies and their central supermassive black holes (BHs) within a semi-analytical framework developed on the outputs of the Millennium Simulation. This model, described in detail in Croton et al. (2006) and De Lucia & Blaizot (2007), introduces a `radio mode' feedback from Active Galactic Nuclei (AGN) at the centre of X-ray emitting atmospheres in galaxy groups and clusters. Thanks to this mechanism, the model can simultaneously explain: (i) the low observed mass drop-out rate in cooling flows; (ii) the exponential cut-off in the bright end of the galaxy luminosity function; and (iii) the bulge-dominated morphologies and old stellar ages of the most massive galaxies in clusters. This paper is the first of a series in which we investigate how well this model can also reproduce the physical properties of BHs and AGN. Here we analyze the scaling relations, the fundamental plane and the mass function of BHs, and compare them with the most recent observational data. Moreover, we extend the semi-analytic model to follow the evolution of the BH mass accretion and its conversion into radiation, and compare the derived AGN bolometric luminosity function with the observed one. While we find for the most part a very good agreement between predicted and observed BH properties, the semi-analytic model underestimates the number density of luminous AGN at high redshifts, independently of the adopted Eddington factor and accretion efficiency. However, an agreement with the observations is possible within the framework of our model, provided it is assumed that the cold gas fraction accreted by BHs at high redshifts is larger than at low redshifts.

(Abridged)The X-ray measurements of the ICM metallicity are becoming more frequent due to the availability of powerful X-ray telescope with excellent spatial and spectral resolutions. The information which can be extracted from the measurements of the alpha-elements, like Oxygen, Magnesium and Silicon with respect to the Iron abundance is extremely important to better understand the stellar formation and its evolutionary history. In this paper we investigate possible source of bias connected to the plasma physics when recovering metal abundances from X-ray spectra. To do this we analyze 6 simulated galaxy clusters processed through the new version of our X-ray MAp Simulator, which allows to create mock XMM-Newton EPIC MOS1 and MOS2 observations. By comparing the spectroscopic results to the input values we find that: i) Fe is recovered with high accuracy for both hot (T>3 keV) and cold (T<2 keV) systems; at intermediate temperatures, however, we find a systematic overestimate which depends on the number counts; ii) O is well recovered in cold clusters, while in hot systems its measure may overestimate by a factor up to 2-3; iii) Being a weak line, the measurement of Mg is always difficult; despite of this, for cold systems (T<2 keV) we do not find any systematic behavior, while for very hot systems (T>5 keV) the spectroscopic measurement may be strongly overestimated up to a factor of 4; iv) Si is well recovered for all the clusters in our sample. We investigate in detail the nature of the systematic effects and biases found. We conclude that they are mainly connected with the multi-temperature nature of the projected observed spectra and to the intrinsic limitation of the XMM-Newton EPIC spectral resolution that does not always allow to disentangle among the emission lines produced by different elements.

We study the reliability of dark-matter halo detections with three different linear filters applied to weak-lensing data. We use ray-tracing in the multiple lens-plane approximation through a large cosmological simulation to construct realizations of cosmic lensing by large-scale structures between redshifts zero and two. We apply the filters mentioned above to detect peaks in the weak-lensing signal and compare them with the true population of dark matter halos present in the simulation. We confirm the stability and performance of a filter optimized for suppressing the contamination by large-scale structure. It allows the reliable detection of dark-matter halos with masses above a few times 1e13 M_sun/h with a fraction of spurious detections below ~10%. For sources at redshift two, 50% of the halos more massive than ~7e13 M_sun/h are detected, and completeness is reached at ~2e14 M_sun/h.

We use simulated merger trees of galaxy-cluster halos to study the effect of the halo concentration distribution on strong lensing and X-ray emission. Its log-normal shape typically found in simulations favors outliers with high concentration. Since, at fixed mass, more concentrated halos tend to be more efficient lenses, the scatter in the concentration increases the strong-lensing optical depth by $\lesssim50%$. Within cluster samples, mass and concentration have counteracting effects on strong lensing and X-ray emission because the concentration decreases for increasing mass. Selecting clusters by concentration thus has no effect on the lensing cross section. The most efficiently lensing and hottest clusters are typically the \textit{least} concentrated in samples with a broad mass range. Among cluster samples with a narrow mass range, however, the most strongly lensing and X-ray brightest clusters are typically 10% to 25% more concentrated.

We present a precise estimate of the bulk virial scaling relation of halos formed via hierarchical clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to cosmological parameters, the presence of a trace, dissipationless gas component, and numerical resolution down to a limit of ~1000 particles. The dark matter velocity dispersion scales with total mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026) log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly log-normal, with scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher dispersions containing 10% of the population that are merger transients. We combine this relation with the halo mass function in LCDM models, and show that a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69, favored by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with the current generation of direct simulations of cluster formation. A higher normalization, S_8 = 0.80, alleviates this tension and implies that the hot gas fraction within r_{500} is (0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent Sunyaev-Zel'dovich observations.

We have performed high-resolution cosmological N-body simulations of a concordance LCDM model to study the evolution of virialized, dark matter haloes in the presence of primordial non-Gaussianity. Following a standard procedure, departures from Gaussianity are modeled through a quadratic Gaussian term in the primordial gravitational potential, characterized by a dimensionless non-linearity strength parameter f_NL. We find that the halo mass function and its redshift evolution closely follow the analytic predictions of Matarrese et al.(2000). The existence of precise analytic predictions makes the observation of rare, massive objects at large redshift an even more attractive test to detect primordial non-Gaussian features in the large scale structure of the universe.

[Abridged] We present an analysis of the scaling relations between X-ray properties and Sunyaev-Zel'dovich (SZ) parameters for a sample of 24 X-ray luminous galaxy clusters observed with Chandra and with measured SZ effect. These objects are in the redshift range 0.14--0.82 and have X-ray bolometric luminosity L>10^45 erg/s. We perform a spatially resolved spectral analysis and recover the density, temperature and pressure profiles of the ICM, just relying on the spherical symmetry of the cluster and the hydrostatic equilibrium hypothesis. We observe that the correlations among X-ray quantities only are in agreement with previous results obtained for samples of high-z X-ray luminous galaxy clusters. On the relations involving SZ quantities, we obtain that they correlate with the gas temperature with a logarithmic slope significantly larger than the predicted value from the self-similar model. The measured scatter indicates, however, that the central Compton parameter y_0 is a proxy of the gas temperature at the same level of other X-ray quantities like luminosity. Our results on the X-ray and SZ scaling relations show a tension between the quantities more related to the global energy of the system (e.g. gas temperature, gravitating mass) and the indicators of the structure of the ICM (e.g. gas density profile, central Compton parameter y_0), showing the most significant deviations from the values of the slope predicted from the self-similar model in the L-T, L-M_{tot}, M_{gas}-T, y_0-T relations. When the slope is fixed to the self-similar value, these relations consistently show a negative evolution suggesting a scenario in which the ICM at higher redshift has lower both X-ray luminosity and pressure in the central regions than the expectations from self-similar model.

The polarization sensitivity of the upcoming millimetric observatories will open new possibilities for studying the properties of galaxy clusters and for using them as powerful cosmological probes. For this reason it is necessary to investigate in detail the characteristics of the polarization signals produced by their highly ionized intra-cluster medium (ICM). This work is focussed on the polarization effect induced by the ICM bulk motions, the so-called kpSZ signal, which has an amplitude proportional to the optical depth and to the square of the tangential velocity. In particular we study how this polarization signal is affected by the internal dynamics of galaxy clusters and what is its dependence on the physical modelling adopted to describe the baryonic component. This is done by producing realistic kpSZ maps starting from the outputs of two different sets of high-resolution hydrodynamical N-body simulations. The first set (17 objects) follows only non-radiative hydrodynamics, while for each of 9 objects of the second set we implement four different kinds of physical processes. Our results shows that the kpSZ signal turns out to be a very sensitive probe of the dynamical status of galaxy clusters. We find that major merger events can amplify the signal up to one order of magnitude with respect to relaxed clusters, reaching amplitude up to about 100 nuK. This result implies that the internal ICM dynamics must be taken into account when evaluating this signal because simplicistic models, based on spherical rigid bodies, may provide wrong estimates. Finally we find that the dependence on the physical modelling of the baryonic component is relevant only in the very inner regions of clusters.

Using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model, we study the global properties of the Sunyaev-Zel'dovich (SZ) effects, both considering the thermal (tSZ) and the kinetic (kSZ) component. The simulation follows gravitation and gas dynamics and includes also several physical processes that affect the baryonic component, like a simple reionization scenario, radiative cooling, star formation and supernova feedback. Starting from the outputs of the simulation we create mock maps of the SZ signals due to the large structures of the Universe integrated in the range 0 < z < 6. We predict that the Compton y-parameter has an average value of (1.19 +/- 0.32) 10^-6 and is lognormally distributed in the sky; half of the whole signal comes from z < 1 and about 10 per cent from z > 2. The Doppler b-parameter shows approximately a normal distribution with vanishing mean value and a standard deviation of 1.6 10^-6, with a significant contribution from high-redshift (z > 3) gas. We find that the tSZ effect is expected to dominate the primary CMB anisotropies for l >~ 3000 in the Rayleigh-Jeans limit, while interestingly the kSZ effect dominates at all frequencies at very high multipoles (l >~ 7 10^4). We also analyse the cross-correlation between the two SZ effects and the soft (0.5-2 keV) X-ray emission from the intergalactic medium and we obtain a strong correlation between the three signals, especially between X-ray emission and tSZ effect (r_l ~ 0.8-0.9) at all angular scales.

We use high-resolution hydrodynamic re-simulations to investigate the properties of the thermal Sunyaev-Zel'dovich (SZ) effect from galaxy clusters. We compare results obtained using different physical models for the intracluster medium (ICM), and show how they modify the SZ emission in terms of cluster profiles and scaling relations. We also produce realistic mock observations to verify whether the results from hydrodynamic simulations can be confirmed. We find that SZ profiles depend marginally on the modelled physical processes, while they exhibit a strong dependence on cluster mass. The central and total SZ emission strongly correlate with the cluster X-ray luminosity and temperature. The logarithmic slopes of these scaling relations differ from the self-similar predictions by less than 0.2; the normalization of the relations is lower for simulations including radiative cooling. The observational test suggests that SZ cluster profiles are unlikely to be able to probe the ICM physics. The total SZ decrement appears to be an observable much more robust than the central intensity, and we suggest using the former to investigate scaling relations.

Using a set of hydrodynamical simulations of 9 galaxy clusters with masses in the range 1.5 10^{14} M_sun < M_vir < 3.4 10^{15} M_sun, we have studied the density, temperature and X-ray surface brightness profiles of the intracluster medium in the regions around the virial radius. We have analyzed the profiles in the radial range well above the cluster core, the physics of which are still unclear and matter of tension between simulated and observed properties, and up to the virial radius and beyond, where present observations are unable to provide any constraints. We have modeled the radial profiles between 0.3 R_200 and 3 R_200 with power laws with one index, two indexes and a rolling index. The simulated temperature and [0.5-2] keV surface brightness profiles well reproduce the observed behaviours outside the core. The shape of all these profiles in the radial range considered depends mainly on the activity of the gravitational collapse, with no significant difference among models including extraphysics. The profiles steepen in the outskirts, with the slope of the power-law fit that changes from -2.5 to -3.4 in the gas density, from -0.5 to -1.8 in the gas temperature, and from -3.5 to -5.0 in the X-ray soft surface brightness. We predict that the gas density, temperature and [0.5-2] keV surface brightness values at R_200 are, on average, 0.05, 0.60, 0.008 times the measured values at 0.3 R_200. At 2 R_200, these values decrease by an order of magnitude in the gas density and surface brightness, by a factor of 2 in the temperature, putting stringent limits on the detectable properties of the intracluster-medium (ICM) in the virial regions.

We consider simple semi-analytic models that relate the active galactic nuclei (AGN) evolution to the merging history of their host dark matter haloes and quantify their ability of matching the AGN luminosity function and its spatial clustering at low and intermediate redshifts. In particular, we focus on the recent determinations of the AGN luminosity function in the hard X-ray band at z~0 which constitutes the most stringent observational test for our models. Indeed, while we find an acceptable agreement between the model bolometric luminosity function and the data at 1<z<2 and for luminosities larger than 1e10 solar bolometric luminosities, no semi-analytic model is capable of reproducing the number density of faint X-ray sources in the local universe. Some improvement can be obtained by advocating energy feedback that we model through a time-dependent Eddington ratio. Even in this case, however, the number density of faint AGNs is significantly below observations. This failure indicates that major mergers cannot constitute the only trigger to accretion episodes in the local AGN population.

We present the first hydrodynamic N-body simulations of primordial gas clouds responsible for the reionisation process in dark energy cosmologies. We compare the cosmological constant scenario with a SUGRA quintessence model with marked dynamics in order to highlight effects due to the different acceleration histories imposed by the dark energy. We show that both the number density of gas clouds and their clumpiness keep a record of the expansion rate during evolution, similar to the non-linear dark matter profile at virialisation, as was recently demonstrated by Dolag et al. (2004). Varying the shape of the primordial power spectrum, we show how this effect is mitigated by a running spectral index decreasing the power at small scales. Our results demonstrate that, in order to constrain the dark energy from large scale structures, one must track its effects down to the distribution of luminous matter.

We propose a method for recovering the two-dimensional gravitational potential of galaxy clusters which combines data from weak and strong gravitational lensing. A first estimate of the potential from weak lensing is improved at the approximate locations of critical curves. The method can be fully linearised and does not rely on the existence and identification of multiple images. We use simulations to show that it recovers the surface-mass density profiles and distributions very accurately, even if critical curves are only partially known and if their location is realistically uncertain. We further describe how arcs at different redshifts can be combined, and how deviations from weak lensing can be included.

We investigate how ellipticity, asymmetries and substructures separately affect the ability of galaxy clusters to produce strong lensing events, i.e. gravitational arcs, and how they influence the arc morphologies and fluxes. This is important for those studies aiming, for example, at constraining cosmological parameters from statistical lensing, or at determining the inner structure of galaxy clusters through gravitational arcs. We do so by creating two-dimensional gradually smoothed, differently elliptical and asymmetric versions of some numerical models. On average, we find that the contributions of ellipticity, asymmetries and substructures amount to ~40%, ~10% and ~30% of the total strong lensing cross section, respectively. However, our analysis shows that substructures play a more important role in less elliptical and asymmetric clusters, even if located at large distances from the cluster centers (~1Mpc/h). Conversely, their effect is less important in highly asymmetric lenses. The morphology, position and flux of individual arcs are strongly affected by the presence of substructures in the clusters. Removing substructures on spatial scales <~50kpc/h, roughly corresponding to mass scales <~5 10^{10}M_\odot/h, alters the image multiplicity of ~35% of the sources used in the simulations and causes position shifts larger than 5'' for ~40% of the arcs longer than 5''. We conclude that any model for cluster lens cannot neglect the effects of ellipticity, asymmetries and substructures. On the other hand, the high sensitivity of gravitational arcs to deviations from regular, smooth and symmetric mass distributions suggests that strong gravitational lensing is potentially a powerfull tool to measure the level of substructures and asymmetries in clusters.

We apply the linear filter for the weak-lensing signal of dark-matter halos developed in Maturi et al. (2005) to the cosmic-shear data extracted from the Garching-Bonn-Deep-Survey (GaBoDS). We wish to search for dark-matter halos through weak-lensing signatures which are significantly above the random and systematic noise level caused by intervening large-scale structures. We employ a linear matched filter which maximises the signal-to-noise ratio by minimising the number of spurious detections caused by the superposition of large-scale structures (LSS). This is achieved by suppressing those spatial frequencies dominated by the LSS contamination. We confirm the improved stability and reliability of the detections achieved with our new filter compared to the commonly-used aperture mass (Schneider, 1996; Schneider et al., 1998) and to the aperture mass based on the shear profile expected for NFW haloes (see e.g. Schirmer et al., 2004; Hennawi & Spergel, 2005). Schirmer et al.~(2006) achieved results comparable to our filter, but probably only because of the low average redshift of the background sources in GaBoDS, which keeps the LSS contamination low. For deeper data, the difference will be more important, as shown by Maturi et al. (2005). We detect fourteen halos on about eighteen square degrees selected from the survey. Five are known clusters, two are associated with over-densities of galaxies visible in the GaBoDS image, and seven have no known optical or X-ray counterparts.

Using the results of a high-resolution, cosmological hydrodynamical re-simulation of a supercluster-like region we investigate the physical properties of the gas located along the filaments and bridges which constitute the so-called cosmic web. First we analyze the main characteristics of the density, temperature and velocity fields, which have quite different distributions, reflecting the complex dynamics of the structure formation process. Then we quantify the signals which originate from the matter in the filaments by considering different observables. Inside the cosmic web, we find that the halo density is about 4-6 times larger than in the neighbouring region; the bremsstrahlung X-ray surface brightness reaches at most 10e-16 erg/s/cm^2/armin^2; the Compton-y parameter due to the thermal Sunyaev-Zel'dovich effect is about 10e-6; the reduced shear produced by the weak lensing effect is ~ 0.01-0.02. These results confirm the difficulty of an observational detection of the cosmic web. Finally we find that projection effects of the filamentary network can affect the estimates of the properties of single clusters, increasing their X-ray luminosity by less than 10% and their central Compton-y parameter by up to 30% cent.

We use the GOODS-MUSIC sample, a catalog of ~3000 Ks-selected galaxies based on VLT and HST observation of the GOODS-South field with extended multi-wavelength coverage (from 0.3 to 8 micron) and accurate estimates of the photometric redshifts to select 179 DRGs with J-Ks>1.3 in an area of 135 sq. arcmin. We first show that the J-Ks>1.3 criterion selects a rather heterogeneous sample of galaxies, going from the targeted high-redshift luminous evolved systems, to a significant fraction of lower redshift (1<z<2) and less luminous dusty starbursts. These low-redshift DRGs are significantly less clustered than higher-z DRGs. With the aid of extreme and simplified theoretical models of clustering evolution we show that it is unlikely that the two samples are drawn from the same population observed at two different stages of evolution. High-z DRGs likely represent the progenitors of the more massive and more luminous galaxies in the local Universe and might mark the regions that will later evolve into structures of intermediate mass, like groups or small galaxy clusters. Low-z DRGs, on the other hand, will likely evolve into slightly less massive field galaxies.

We examine the systematics affecting the X-ray mass estimators applied to a set of five simulated galaxy clusters. They have been processed through the X-ray Map Simulator, X-MAS, to provide Chandra-like long exposures that are analyzed to reconstruct the gas temperature, density, and mass profiles used as input. We find that at R_2500 the mass profile obtained via a direct application of the hydrostatic equilibrium equation is consistent within 1 sigma with the actual mass; although we notice this estimator shows high statistical errors due to high level of Chandra background. Instead, the poorness of the beta-model in describing the gas density profile makes the evaluated masses to be underestimated by \sim 40 per cent with respect to the true mass, both with an isothermal and a polytropic temperature profile. We also test ways to recover the mass by adopting an analytic mass model, such as those proposed by Navarro et al. (1997) and Rasia et al. (2004), and fitting the temperature profile expected from the hydrostatic equilibrium equation to the observed one. These methods and the one of the hydrostatic equilibrium equation provide a more robust mass estimation than the ones based on the beta-model. In the present work the main limitation for a precise mass reconstruction is to ascribe to the relatively high level of the background chosen to reproduce the Chandra one. After artificially reducing it by a factor of 100, we find that the estimated mass significantly underestimates the true mass profiles. This is manly due (i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to the bias towards lower values of the X-ray temperature measurements because of the complex thermal structure of the emitting plasma.

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

We study the properties of the diffuse X-ray background by using the results of a cosmological hydrodynamical simulation of the concordance LambdaCDM model. The simulation follows gravitational and gas dynamics and includes a treatment of physical processes like radiative cooling, star formation and supernova feedback. From the simulation outputs, we produce a set of two-dimensional maps of the intergalactic medium X-ray emission integrated over redshift. We find that the signal in the soft (0.5-2 keV) band is lognormally distributed with a mean intensity of about 4 10^-12 erg s^-1 cm^-2 deg^-2; approximately 40 per cent of the emission originates from warm-hot gas (defined as baryons with 10^5<T<10^7 K), and 90 per cent comes from structures at z<0.9. Since the spectrum is soft, being mostly provided by the intergalactic medium at low temperature, the total mean intensity in the hard (2-10 keV) X-ray band is smaller by a factor of about 4. In order to constrain the physical processes included in our simulation, we compare our results with the observed upper limit (1.2 +/- 0.3) 10^-12 erg s^-1 cm^-2 deg^-2 of the soft X-ray emission due to diffuse gas. To this purpose, we remove the contributions of observable extended objects (groups and clusters of galaxies) from the simulated maps by adopting different detectability criteria which are calibrated on the properties of systems at intermediate redshifts observed by Chandra. We show that the simulated diffuse soft X-ray emission is consistent with the present observed upper limit. However, if future measurements will decrease the level of the unresolved X-ray background by a factor of two, a more efficient feedback mechanism should be required to suppress the soft emission of the gas residing within filaments and group-size haloes.

We discuss an extended set of Tree+SPH simulations of galaxy clusters, with the goal of investigating the interplay between numerical resolution effects and star-formation/feedback processes. The simulated clusters span the mass range (0.1-2.3) 10^{15}Msun/h, with mass resolution varying by several decades. At the highest achieved resolution, we resolve the virial region of a Virgo-like cluster with more than 2 million gas particles and with at least as many dark-matter (DM) particles. Our resolution study confirms that, in the absence of an efficient feedback mechanism, runaway cooling leads to about 35 per cent of baryons in clusters to be locked up in long lived stars at our highest resolution, with no evidence of convergence. However, including feedback causes the fraction of cooled baryons to converge at about 15 per cent already at modest resolution. Feedback also stabilizes other gas-related quantities, such as radial profiles of entropy, gas density and temperature, against variations due to changes in resolution. We also investigate the influence of the gravitational force softening length, and that of numerical heating of the gas induced by two-body encounters between DM and lighter gas particles. We show that simulations where more DM than gas particles are used, show a significantly enhanced efficiency of star formation at z>3. Our results are important for establishing and delineating the regime of numerical reliability of the present generation of hydrodynamical simulations of galaxy clusters.

We propose a fast and efficient bispectrum statistic for Cosmic Microwave Background (CMB) temperature anisotropies to constrain the amplitude of the primordial non-Gaussian signal measured in terms of the non-linear coupling parameter f_NL. We show how the method can achieve a remarkable computational advantage by focussing on subsets of the multipole configurations, where the non-Gaussian signal is more concentrated. The detection power of the test, increases roughly linearly with the maximum multipole, as shown in the ideal case of an experiment without noise and gaps. The CPU-time scales as l_{max}^3 instead of l_{max}^5 for the full bispectrum which for Planck resolution l_{max} \sim 3000 means an improvement in speed of a factor 10^7 compared to the full bispectrum analysis with minor loss in precision. We find that the introduction of a galactic cut partially destroys the optimality of the configuration, which will then need to be dealt with in the future. We find for an ideal experiment with l_{max}=2000 that upper limits of f_{NL}<8 can be obtained at 1 sigma. For the case of the WMAP experiment, we would be able to put limits of |f_{NL}|<40 if no galactic cut were present. Using the real data with galactic cut, we obtain an estimate of -80<f_{NL}<80 and -160<f_{NL}<160 at 1 and 2 sigma respectively.

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.

We use high-resolution hydrodynamical simulations of large-scale structure formation to study the imprints of the local superclusters onto the full-sky Sunyaev-Zel'dovich (SZ) signals. Following (Mathis et al. 2002), the initial conditions have been statistically constrained to reproduce the density field within a sphere of 110 Mpc around the Milky Way, as observed in the IRAS 1.2-Jy all-sky redshift survey. As a result, the positions and masses of prominent galaxy clusters and superclusters in our simulations coincide closely with their real counterparts in the local universe. We present the results of two different runs, one with adiabatic gas physics only, and one also including cooling, star formation and feedback. By analysing the full-sky maps for the thermal and kinetic SZ signals extracted from these simulations, we find that for multipoles with l<100 the power spectrum is dominated by the prominent local superclusters, and its amplitude at these scales is a factor of two higher than that obtained from unconstrained simulations; at lower multipoles (l<20) this factor can even reach one order of magnitude. We check the influence of the SZ effect from local superclusters on the cosmic microwave background (CMB) power spectrum at small multipoles and find it negligible and with no signs of quadrupole-octopole alignment. However, performing simulations of the CMB radiation including the experimental noise at the frequencies which will be observed by the Planck satellite, we find results suggesting that an estimate of the SZ power spectrum at large scales can be extracted.

We test the consistency of estimates of the non-linear coupling constant f_{NL} using non-Gaussian CMB maps generated by the method described in (Liguori, Matarrese and Moscardini 2003). This procedure to obtain non-Gaussian maps differs significantly from the method used in previous works on estimation of f_{NL}. Nevertheless, using spherical wavelets, we find results in very good agreement with (Mukherjee and Wang 2004), showing that the two ways of generating primordial non-Gaussian maps give equivalent results. Moreover, we introduce a new method for estimating the non-linear coupling constant from CMB observations by using the local curvature of the temperature fluctuation field. We present both Bayesian credible regions (assuming a flat prior) and proper (frequentist) confidence intervals on f_{NL}, and discuss the relation between the two approaches. The Bayesian approach tends to yield lower error bars than the frequentist approach, suggesting that a careful analysis of the different interpretations is needed. Using this method, we estimate f_{NL}=-10^{+270}_{-260} at the 2\sigma level (Bayesian) and f_{NL}=-10^{+310}_{-270} (frequentist). Moreover, we find that the wavelet and the local curvature approaches, which provide similar error bars, yield approximately uncorrelated estimates of f_{NL} and therefore, as advocated in (Cabella et al. 2004), the estimates may be combined to reduce the error bars. In this way, we obtain f_{NL}=-5\pm 85 and f_{NL}=-5\pm 175 at the 1\sigma and 2\sigma level respectively using the frequentist approach.

In the presence of strong clustering, low-resolution surveys measure the summed contributions of groups of sources within the beam. The counts of bright intensity peaks are therefore shifted to higher flux levels compared to the counts of individual sources detected with high-resolution instruments. If the beam-width corresponds a sizable fraction of the clustering size, as in the case of Planck/HFI, one actually detects the fluxes of clumps of sources. We argue that the distribution of clump luminosities can be modelled in terms of the two- and three-point correlation functions, and apply our formalism to the Planck/HFI 850 micron surveys. The effect on counts is found to be large and sensitive also to the evolution of the three-point correlation function; in the extreme case that the latter function is redshift-independent, the source confusion due to clustering keeps being important above the canonical 5*sigma detection limit. Detailed simulations confirm the reliability of our approach. As the ratio of the beam-width to the clustering angular size decreases, the observed fluxes approach those of the brightest sources in the beam and the clump formalism no longer applies. However, simulations show that also in the case of the Herschel/SPIRE 500 micron survey the enhancement of the bright source counts due to clustering is important.

We construct a linear filter optimised for detecting dark-matter halos in weak-lensing data. The filter assumes a mean radial profile of the halo shear pattern and modifies that shape by the noise power spectrum. Aiming at separating dark-matter halos from spurious peaks caused by large-scale structure lensing, we model the noise as being composed of weak lensing by large-scale structures and Poisson noise from random galaxy positions and intrinsic ellipticities. Optimal filtering against the noise requires the optimal filter scale to be smaller than typical halo sizes. Although a perfect separation of halos from spurious large-scale structure peaks is strictly impossible, we use numerical simulations to demonstrate that our filter produces substantially more sensitive, reliable and stable results than the conventionally used aperture-mass statistic.

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. Here we show that, if the intra-cluster medium is thermally complex, T_ew fails at reproducing T_spec. We propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we also find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and beta-model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X--ray temperature function could bias low sigma_8 by 10-20 per cent. This may alleviate the tension between the value of sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

The fluctuations in the Cosmic Microwave Background (CMB) intensity due to the Sunyaev-Zeldovich (SZ) effect are the sum of a thermal and a kinetic contribution. Separating the two components to measure the peculiar velocity of galaxy clusters requires radio and microwave observations at three or more frequencies, and knowledge of the temperature T_e of the intracluster medium weighted by the electron number density. To quantify the systematics of this procedure, we extract a sample of 117 massive clusters at redshift z=0 from an N-body hydrodynamical simulation, with 2x480^3 particles, of a cosmological volume 192 Mpc/h on a side of a flat Cold Dark Matter model with Omega_0=0.3 and Lambda=0.7. Our simulation includes radiative cooling, star formation and the effect of feedback and galactic winds from supernovae. We find that (1) our simulated clusters reproduce the observed scaling relations between X-ray and SZ properties; (2) bulk flows internal to the intracluster medium affect the velocity estimate by less than 200 km/s in 93 per cent of the cases; (3) using the X-ray emission weighted temperature, as an estimate of T_e, can overestimate the peculiar velocity by 20-50 per cent, if the microwave observations do not spatially resolve the cluster. For spatially resolved clusters, the assumptions on the spatial distribution of the ICM, required to separate the two SZ components, still produce a velocity overestimate of 10-20 per cent, even with an unbiased measure of T_e. Thanks to the large size of our cluster samples, these results set a robust lower limit of 200 km/s to the systematic errors that will affect upcoming measures of cluster peculiar velocities with the SZ effect.

The detection of excess of soft X-ray or Extreme Ultraviolet (EUV) radiation, above the thermal contribution from the hot intracluster medium (ICM), has been a controversial subject ever since the initial discovery of this phenomenon. We use a large--scale hydrodynamical simulation of a concordance $\Lambda$CDM model, to investigate the possible thermal origin for such an excess in a set of 20 simulated clusters having temperatures in the range 1--7 keV. Simulated clusters are analysed by mimicking the observational procedure applied to ROSAT--PSPC data, which for the first time showed evidences for the soft X-ray excess. For cluster--centric distances $0.4< R/R_{\rm vir}< 0.7$ we detect a significant excess in most of the simulated clusters, whose relative amount changes from cluster to cluster and, for the same cluster, by changing the projection direction. In about 30 per cent of the cases, the soft X-ray flux is measured to be at least 50 per cent larger than predicted by the one--temperature plasma model. We find that this excess is generated in most cases within the cluster virialized regions. It is mainly contributed by low--entropy and high--density gas associated with merging sub--halos, rather than to diffuse warm gas. Only in a few cases the excess arises from fore/background groups observed in projection, while no evidence is found for a significant contribution from gas lying within large--scale filaments. We compute the distribution of the relative soft excess, as a function of the cluster--centric distance, and compare it with the observational result by Bonamente et al. (2003) for the Coma cluster. Similar to observations, we find that the relative excess increases with the distance from the cluster center, with no significant excess detected for $R<0.4R_{\rm vir}$. (abridged)

Theoretical studies of the physical processes in clusters of galaxies are mainly based on the results of numerical simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We show that the projected spectroscopic temperature of clusters obtained from X-ray observations is always lower than the emission-weighed temperature. This bias is related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. In real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps. We propose a new formula, the spectroscopic-like temperature function that better approximates the spectroscopic temperature, making simulations more directly comparable to observations

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is thermally complex, T_ew fails at reproducing T_spec, and propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T_sl relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and \beta--model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X-ray temperature function could bias low \sigma_8 by 10-20 per cent. This may alleviate the tension between the value of \sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

We adapt a non-linear filter proposed by Hu 2001 for detecting lensing of the CMB by large-scale structures to recover surface-density profiles of galaxy clusters from their localised, weak gravitational lensing effect on CMB fields. Shifting the band-pass of the filter to smaller scales, and normalising it such as to reproduce the convergence rather than the deflection angle, we find that the mean density profile of a sample of 100 clusters can be recovered to better than 10% from well within the scale radius to almost the virial radius. The kinetic Sunyaev-Zel'dovich effect is shown to be a negligible source of error. We test the filter applying it to data simulated using the characteristics of the Atacama Cosmology Telescope (ACT), showing that it will be possible to recover mean cluster profiles outside a radius of 1' corresponding to ACT's angular resolution.

(ABRIDGED) We analyze the X-ray properties of a sample of local and high redshift galaxy clusters extracted from a large cosmological hydrodynamical simulation. This simulation has been realized using the Tree+SPH code GADGET-2 for a LambdaCDM model. In our analysis, we consider only objects with T_ew >2 keV and adopt an approach that mimics observations, associating with each measurement an error comparable with recent observations and providing best-fit results via robust techniques. Within the clusters, baryons are distributed among (i) a cold neutral phase, with a relative contribution that increases from less than 1 to 3 per cent at higher redshift, (ii) stars which contribute with about 20 per cent and (iii) the X-ray emitting plasma that contributes by 80 (76) per cent at z=0 (1) to the total baryonic budget. A depletion of the cosmic baryon fraction of ~7 (at z=0) and 5 (at z=1) per cent is measured at the virial radius, R_vir, in good agreement with adiabatic hydrodynamical simulations. We confirm that, also at redshift >0.5, power-law relations hold between gas temperature, T, bolometric luminosity, L, central entropy, S, gas mass, M_gas, and total gravitating mass, M_tot and that these relations are steeper than predicted by simple gravitational collapse. A significant, negative evolution in the L-T and L-M_tot relations and positive evolution in the S-T relation are detected at 0.5 < z < 1 in this set of simulated galaxy clusters. This is partially consistent with recent analyses of the observed properties of z>0.5 X-ray galaxy clusters. By fixing the slope to the values predicted by simple gravitational collapse, we measure at high redshift normalizations lower by 10-40 per cent in the L-T, M_tot-T, M_gas-T, f_gas-T and L-M_tot relations than the observed estimates.

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.

We study the properties of the diffuse light in galaxy clusters forming in a large hydrodynamical cosmological simulation of the Lambda-CDM cosmology. The simulation includes a model for radiative cooling, star formation in dense cold gas, and feedback by SN-II explosions. We select clusters having mass M>10^(14) h^(-1) Msun and study the spatial distribution of their star particles. While most stellar light is concentrated in gravitationally bound galaxies orbiting in the cluster potential, we find evidence for a substantial diffuse component, which may account for the extended halos of light observed around central cD galaxies. We find that more massive simulated clusters have a larger fraction of stars in the diffuse light than the less massive ones. The intracluster light is more centrally concentrated than the galaxy light, and the stars in the diffuse component are on average older than the stars in cluster galaxies, supporting the view that the diffuse light is not a random sampling of the stellar population in the cluster galaxies. We thus expect that at least ~10% of the stars in a cluster may be distributed as intracluster light, largely hidden thus far due to its very low surface brightness.

We study the efficiency of numerically simulated galaxy clusters for producing strong gravitational lensing events in different dark-energy cosmologies with constant and time-variable equation of state and we compare it with ``standard'' $\Lambda$CDM and OCDM models. Our main results are that: (1) the expected abundance of gravitational arcs with large length-to-width ratio depends on the equation of state of dark energy at the epoch of formation of the halo; (2) as found in previous studies, the strong lensing cross section of galaxy clusters is very sensitive to dynamical processes like mergers, whose impact however is different for different cosmologies, being stronger for models in which halos are less concentrated. As expected, the largest differences in the lensing optical depth occur at moderate and high redshift.

We present detailed predictions for the confusion noise due to extragalactic sources in the far-IR/(sub)-millimeter channels of ESA/ISO, NASA/Spitzer, ESA/Herschel and ESA/Planck satellites, including the contribution from clustering of unresolved SCUBA galaxies. Clustering is found to increase the confusion noise, compared to the case of purely Poisson fluctuations, by 10-15% for the lowest frequency (i.e. lowest angular resolution) Spitzer and Herschel channels, by 25-35% for the 175 micron ISOPHOT channel, and to dominate in the case of Planck/HFI channels at nu>143GHz. Although our calculations make use of a specific evolutionary model (Granato et al. 2004), the results are strongly constrained by the observed counts and by data on the redshift distribution of SCUBA sources, and therefore are not expected to be heavily model dependent. The main uncertainty arises from the poor observational definition of the source clustering properties. Two models have been used for the latter: a power-law with constant slope and a redshift-independent comoving correlation length,r_0, and the standard theoretical model for clustering evolution in a LambdaCDM universe, with a redshift-dependent bias factor. In both cases, the clustering amplitude has been normalized to yield a unit angular correlation function at theta_0=1-2 arcsec for 850 micron sources fainter than 2 mJy, consistent with the results by Peacock et al. (2000). This normalization yields, for the first model, r_0=8.3$ Mpc/h, and, for the second model, an effective mass of dark matter haloes in which these sources reside of M_halo=1.8*10^{13} M_sun/h. These results are consistent with independent estimates for SCUBA galaxies and for other, likely related, sources.

High-resolution, numerical simulations of 17 cluster-sized dark-matter haloes in eight different cosmologies with and without dynamical dark energy confirm the picture that core halo densities are imprinted early during their formation by the mean cosmological density. Quite independent of cosmology, halo concentrations have a log-normal distribution with a scatter of ~0.2 about the mean. We propose a simple scaling relation for halo concentrations in dark-energy cosmologies.

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.

We analyze the dynamical and thermal evolution of dark matter and ICM in hydrodynamical Tree-SPH simulations of galaxy clusters. Starting from a sample of 17 high-resolution objects, with virial mass ranging from 3 10^{14} to 1.7 10^{15} M_sun, we follow the build-up of the systems in dark matter and hot gas through the repeated merging of satellites along their merging history trees. We measure the self-bound mass fraction of subhaloes as a function of time after the merging, estimate the satellite mean orbital properties as a function of the mass ratio with the main cluster at merging time, and study the evolution of their internal velocity dispersion, gas temperature and entropy as the substructure is disrupted by various dynamical processes, reaching eventually thermo-dynamic equilibrium in the gravitational potential of the main cluster. We model some relevant properties of subhalo orbits, as the time of the first pericentric and apocentric passages, and the typical distances and velocities at the corresponding times. This survival study can be used to interpret the dynamics of observed merging clusters; as an example we apply our results to the system 1E0657-56. We show that, in the light of our results, the most likely interpretation of the data for this cluster points to the merger of a small group with mass M\approx 1 10^{13} M_sun with a massive cluster with M\approx 1.3 10^{15} M_sun.

Using high-resolution non-radiative hydrodynamic simulations of galaxy clusters we obtain simple analytic formulae for DM and gas distribution in the spherical approximation. We derive fits for the DM density, velocity dispersion and velocity anisotropy. We use these models to test the dynamical equilibrium hypothesis through the Jeans equation: we find that this is satisfied to good accuracy. This result show that our fits constitute a self-consistent dynamical model. We then extend our analysis to the gas, studying its density, temperature and velocity structure, with no further hypothesis on the dynamical status or state equation. Gas and DM show self-similar density profiles down to 0.06 R_vir, while at smaller radii the gas produces a central core. Gas temperatures are almost isothermal out to 0.2 R_vir, then steeply decrease reaching at R_vir a value a factor of 2 lower. We find that the gas is not at rest inside R_vir: velocity dispersions are increasing functions of the radius, motions are isotropic to slightly tangential, and contribute to the total pressure support. We test our model using a generalization of the hydrostatic equilibrium equation, where the gas motion is properly taken into account. We find that our fits provide an accurate description of the system: the gas is in equilibrium and is a good tracer the overall cluster potential if all terms are taken into account, while simpler assumptions, as the beta-model, cause systematic mass underestimates. We also find that, if gas velocities are neglected, then a simple isothermal model fares better at large radii than a non-isothermal one. (abridged)

We obtain constraints on quintessence models from magnitude-redshift measurements of 176 type Ia Supernovae. The considered quintessence models are ordinary quintessence, with Ratra-Peebles and SUGRA potentials, and extended quintessence with a Ratra-Peebles potential. We compute confidence regions in the $\Omega_{m0}-\alpha$ plane and find that for SUGRA potentials it is not possible to obtain useful constraints on these parameters; for the Ratra-Peebles case, both for the extended and ordinary quintessence we find $\alpha\mincir 0.8$, at the $1\sigma$ level. We also consider simulated dataset for the SNAP satellite for the same models: again, for a SUGRA potential it will not be possible to obtain constraints on $\alpha$, while with a Ratra-Peebles potential its value will be determined with an error $\mincir 0.6$. We evaluate the inaccuracy made by approximating the time evolution of the equation of state with a linear or constant $w\diz$, instead of using its exact redshift evolution. Finally we discuss the effects of different systematic errors in the determination of quintessence parameters.

We study the impact of merger events on the strong lensing properties of galaxy clusters. Previous lensing simulations were not able to resolve dynamical time scales of cluster lenses, which arise on time scales which are of order a Gyr. In this case study, we first describe qualitatively with an analytic model how some of the lensing properties of clusters are expected to change during merging events. We then analyse a numerically simulated lens model for the variation in its efficiency for producing both tangential and radial arcs while a massive substructure falls onto the main cluster body. We find that: (1) during the merger, the shape of the critical lines and caustics changes substantially; (2) the lensing cross sections for long and thin arcs can grow by one order of magnitude and reach their maxima when the extent of the critical curves is largest; (3) the cross section for radial arcs also grows, but the cluster can efficiently produce this kind of arcs only while the merging substructure crosses the main cluster centre; (4) while the arc cross sections pass through their maxima as the merger proceeds, the cluster's X-ray emission increases by a factor of $\sim5$. Thus, we conclude that accounting for these dynamical processes is very important for arc statistics studies. In particular, they may provide a possible explanation for the arc statistics problem.

Numerical hydro-N-body simulations are very important tools for making theoretical predictions for the formation of galaxy clusters. They show that the atmospheres of clusters of galaxies have quite complex angular and thermal structures. The full understanding of the physical processes behind these features can be only achieved by direct comparison of observations to hydro-N-body simulations. Although simple in principle, these comparisons are not always trivial. In fact, real data are convolved with the instrument response which may substantially influence the apparent properties of the studied features. To overcome this problem we build the software package X-MAS devoted to simulate X-ray observations of galaxy clusters obtained from hydro-N-body simulations. In this paper we present how this software package works and discuss its application to the simulation of Chandra ACIS-S3 observations. We compare some of the main physical properties of the input data to the ones derived from simulated observations after performing a standard imaging and spectral analysis. We show that, if the thermal structure of the cluster along a particular line of sight is quite complex, the projected spectroscopic temperature obtained from the observation is significantly lower than the emission-weighed value inferred directly from hydrodynamical simulation. This implies that much attention must be paid in the theoretical interpretation of observational temperatures.

We present results on the X-ray properties of clusters and groups of galaxies, extracted from a large hydrodynamical simulation. We used the GADGET code to simulate a LambdaCDM model within a box of 192 Mpc/h on a side, with 480^3 dark matter particles and as many gas particles. The simulation includes radiative cooling, star formation and supernova feedback. The simulated M-T relation is consistent with observations once we mimic the procedure for mass estimates applied to real clusters. Also, with the adopted choices of Omega_m=0.3 and sigma_8=0.8 the resulting XTF agrees with observational determinations. The L-T relation also agrees with observations for clusters with T>2 keV, with no change of slope at the scale of groups. The entropy in central cluster regions is higher than predicted by gravitational heating alone, the excess being almost the same for clusters and groups. The simulated clusters appear to have suffered some overcooling, with f*~0.2, thus about twice as large as the value observed. Interestingly, temperature profiles are found to steadily increase toward cluster centers. They decrease in the outer regions, much like observational data do at r>0.2r_vir, while not showing an isothermal regime followed by a smooth temperature decline in the innermost regions.

This is the third paper of a series describing the Asiago-ESO/RASS QSO survey (AERQS), a project aimed at the construction of an all-sky statistically well-defined sample of relatively bright QSOs (B<15) at z<0.3. We present here the clustering analysis of the full spectroscopically identified database (392 AGN). The clustering signal at 0.02<z<0.22 is detected at a 3-4 sigma level and its amplitude is measured to be r_0=8.6\pm 2.0 h^{-1} Mpc (in a LambdaCDM model). The comparison with other classes of objects shows that low-redshift QSOs are clustered in a similar way to Radio Galaxies, EROs and early-type galaxies in general, although with a marginally smaller amplitude. The comparison with recent results from the 2QZ shows that the correlation function of QSOs is constant in redshift or marginally increasing toward low redshift. We discuss this behavior with physically motivated models, deriving interesting constraints on the typical mass of the dark matter halos hosting QSOs, M_DMH= 10^{12.7} h^{-1} M_sun (10^{12.0}-10^{13.5}h^{-1} M_sun at 1 sigma confidence level). Finally, we use the clustering data to infer the physical properties of local AGN, obtaining M_BH=2 10^8 h^{-1} M_sun (10^7-3 10^9 h^{-1} M_sun) for the mass of the active black holes, tau_{AGN}= 8 10^6 yr (2 10^{6}-5 10^{7} yr) for their life-time and eta = 0.14 for their efficiency (always for a LambdaCDM model).

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

We perform ray-tracing simulations evaluating the effect of a cD galaxy on the strong lensing properties of five galaxy cluster halos obtained from N-body simulations. The cD galaxy is modelled using both axially symmetric and elliptical models and assuming several masses for its dark matter halo. The effect of the cD orientation with respect to the mass distribution of the host galaxy cluster is also investigated. The simulations are carried out in an open and a flat model universe with cosmological constant. We find that the enhancement of the cluster lensing cross sections for long and thin arcs due to the presence of a massive cD at the cluster centre is typically less than 100%, depending on the model used for the cD galaxy and its orientation. The impact of the cD on the cluster efficiency for producing radially magnified images is larger only for those clusters whose lensing cross section for radial arcs is very small in absence of the central galaxy. We conclude that the presence of a cD galaxy at the cluster centre can only moderately influence the cluster efficiency for strong lensing and in particular fails to explain the discrepancy between the observed number of giant arcs in galaxy clusters and their abundance predicted from lensing simulations in the currently most favoured $\Lambda$CDM model.

We investigate how the probability of the formation of giant arcs in galaxy clusters is expected to change in cosmological models dominated by dark energy with an equation of state p=w rho c^2 compared to cosmological-constant or open models. To do so, we use a simple analytic model for arc cross sections based on the Navarro-Frenk-White density profile which we demonstrate reproduces essential features of numerically determined arc cross sections. Since analytic lens models are known to be inadequate for accurate absolute quantifications of arc probabilities, we use them only for studying changes relative to cosmological-constant models. Our main results are (1) the order of magnitude difference between the arc probabilities in low density, spatially flat and open CDM models found numerically is reproduced by our analytic model, and (2) dark-energy cosmologies with w>-1 increase the arc optical depth by at most a factor of two and are thus unlikely to reconcile arc statistics with spatially flat cosmological models with low matter density.

We describe a new numerical algorithm to obtain high-resolution simulated maps of the Cosmic Microwave Background (CMB), for a broad class of non-Gaussian models. The kind of non-Gaussianity we account for is based on the simple idea that the primordial gravitational potential is obtained by a non-linear but local mapping from an underlying Gaussian random field, as resulting from a variety of inflationary models. Our technique, which is based on a direct realization of the potential in spherical coordinates and fully accounts for the radiation transfer function, allows to simulate non-Gaussian CMB maps down to the Planck resolution ($\ell_{\rm max} \sim 3,000$), with reasonable memory storage and computational time.

Arc statistics is known to be a powerful cosmological tool. Numerical lensing simulations show that orders of magnitude differences in the number of {\em giant} arcs on the whole sky are expected in different cosmological models. In this paper, we discuss the analytic and numerical methods in arc statistics and show that analytic models fail to reproduce the efficiency for strong lensing of more realistic numerical cluster models. Then, we discuss two recent extensions of the lensing simulations, i.e. the effects of cD galaxies in the lensing clusters and the impact of cluster mergers on arc statistics. We show that cD galaxies can increase the lensing cross sections for long and thin arcs by perhaps up to $\sim50%$, while major mergers can change the cluster efficiency for producing such arcs by up to one order of magnitude.

The statistics of gravitationally lensed arcs was recognised earlier as a potentially powerful cosmological probe. However, while fully numerical models find orders of magnitude difference between the arc probabilities in different cosmological models, analytic models tend to find markedly different results. We introduce in this paper an analytic cluster lens model which improves upon existing analytic models in four ways. (1) We use the more realistic Navarro-Frenk-White profile instead of singular isothermal spheres, (2) we include the effect of cosmology on the compactness of the lenses, (3) we use elliptical instead of axially symmetric lenses, and (4) we take the intrinsic ellipticity of sources into account. While these improvements to the analytic model lead to a pronounced increase of the arc probability, comparisons with numerical models of the same virial mass demonstrate that the analytic models still fall short by a substantial margin of reproducing the results obtained with numerical models. Using multipole expansions of cluster mass distributions, we show that the remaining discrepancy can be attributed to substructure inside clusters and tidal fields contributed by the cluster surroundings, effects that cannot reasonably and reliably be mimicked in analytic models.

We study the clustering properties of the galaxy clusters detectable for the Planck satellite due to their thermal Sunyaev-Zel'dovich effect. We take the past light-cone effect and the redshift evolution of both the underlying dark matter correlation function and the cluster bias factor into account. A theoretical mass-temperature relation allows us to convert the sensitivity limit of a catalogue into a minimum mass for the dark matter haloes hosting the clusters. We confirm that the correlation length is an increasing function of the sensitivity limits defining the survey. Using the expected characteristics of the Planck cluster catalogue, which will be a quite large and unbiased sample, we predict the two-point correlation function and power spectrum for different cosmological models. We show that the wide redshift distribution of the Planck survey, will allow to constrain the cluster clustering properties up to z=1. The dependence of our results on the main cosmological parameters (the matter density parameter, the cosmological constant and the normalisation of the density power-spectrum) is extensively discussed. We find that the future Planck clustering data place only mild constraints on the cosmological parameters, because the results depend on the physical characteristics of the intracluster medium, like the baryon fraction and the mass-temperature relation. Once the cosmological model and the Hubble constant are determined, the clustering data will allow a determination of the baryon fraction with an accuracy of few per cent.

We present an analysis of the evolution of galaxy clustering in the redshift interval 0<z<4.5 in the HDF-S. The HST optical data are combined with infrared ISAAC/VLT observations, and photometric redshifts are used for all the galaxies brighter than I_AB<27.5. The clustering signal is obtained in different redshift bins using two different approaches: a standard one, which uses the best redshift estimate of each object, and a second one, which takes into account the redshift probability function of each object. This second method makes it possible to improve the information in the redshift intervals where contamination from objects with insecure redshifts is important. With both methods, we find that the clustering strength up to z~3.5 in the HDF-S is consistent with the previous results in the HDF-N. While at redshift lower than z~1 the HDF galaxy population is un/anti-biased (b<1) with respect to the underlying dark matter, at high redshift the bias increases up to b~2-3, depending on the cosmological model. These results support previous claims that, at high redshift, galaxies are preferentially located in massive haloes, as predicted by the biased galaxy formation scenario. The impact of cosmic errors on our analyses has been quantified, showing that errors in the clustering measurements in the HDF surveys are indeed dominated by shot-noise in most regimes. Future observations with instruments like the ACS on HST will improve the S/N by at least a factor of two and more detailed analyses of the errors will be required. In fact, pure shot-noise will give a smaller contribution with respect to other sources of errors, such as finite volume effects or non-Poissonian discreteness effects.

Upcoming wide-area surveys in the submillimetre regime will allow the construction of complete galaxy cluster samples through their thermal Sunyaev-Zel'dovich effect. We propose an analytic method to predict the number of gravitationally lensed giant arcs produced by the cluster sample expected to be detectable for the Planck satellite. The statistics of lensed arcs can be used to constrain cosmological parameters. This computation implies the choice of an analytic cluster model which has to reproduce the lensing properties of realistic galaxy clusters. Comparing the results of analytic computations and numerical ray-tracing simulations, we show that spherical models with the Navarro-Frenk-White density profile fail by far to produce as many giant arcs as realistic galaxy clusters. This confirms the great importance of asymmetries and substructures in the lensing matter distribution. By elliptically distorting the lensing potential of spherical NFW density profiles, we show that the discrepancy between analytic and numerical computations can almost completely be removed.

The clustering properties of SCUBA-selected galaxies are investigated within the framework of a unifying scheme relating the formation of QSOs and spheroids. The theoretical angular correlation function is derived for different bias functions, corresponding to different values of the ratio $M_{halo}/M_{sph}$ between the mass of the dark halos hosting such galaxies and the mass in stars produced at the end of the major star-formation burst. SCUBA sources are predicted to be strongly clustered, with a clustering strength increasing with mass. Comparisons with the best available measurements show better fits for $M_{halo}/M_{sph}\simeq 100$. The model can also account for the clustering of Lyman-break galaxies, seen as the optical counterpart of low- to intermediate-mass primeval spheroidal galaxies. Best agreement is once again obtained for high values of the $M_{halo}/M_{sph}$ ratio. We also discuss implications for small scale fluctuations observed at different wavelengths by forthcoming experiments such as the Planck mission aimed at mapping the Cosmic Microwave Background (CMB).

(abridged) We use a theoretical model to predict the clustering properties of galaxy clusters. Our technique accounts for past light-cone effects on the observed clustering and follows the non-linear evolution of the dark matter correlation function and cluster bias factor. A linear treatment of redshift-space distortions is also included. We perform a maximum-likelihood analysis by comparing the theoretical predictions to a set of observational data, both in the optical (APM and EDCC catalogues) and X-ray band (RASS1 Bright Sample, BCS, XBACs, REFLEX). In the framework of CDM models, we compute the constraints on cosmological parameters. Our results show that X-ray data are more powerful than optical ones, allowing smaller regions in the parameter space. We find that the clustering properties of clusters are almost independent of the matter density parameter and of the presence of a cosmological constant, while they appear to be strongly dependent on the shape parameter. Using the X-ray data only, we obtain Gamma~0.1 and 0.4<sigma_8<1.1 for the Einstein-de Sitter model, while 0.14<Gamma<0.22 and 0.6<sigma_8<1.3 for open and flat models with Omega_0m=0.3. Finally, we use our model to make predictions on the correlation length of galaxy clusters expected in future surveys. In particular, we show the results for an optical catalogue with characteristics similar to the EIS project and for a very deep X-ray catalogue with the characteristics of the XMM/LSS survey. We find that clusters at high redshifts are expected to have larger a correlation length than local ones.

We carry out ray tracing through five high resolution simulations of a galaxy cluster to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long (``collisionless'') and short (``fluid'') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing ``extreme'' arcs of various sizes. Even a few collisions per particle modify the core structure enough to destroy the cluster's ability to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the ``fluid'' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm^2 g^-1, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.

This paper investigates the clustering properties of SCUBA-selected galaxies in the framework of a unifying scheme relating the formation of QSOs and spheroids (Granato et al. 2000). The theoretical angular correlation function is derived for different bias functions, corresponding to different values of the ratio $M_{\rm halo}/M_{\rm sph}$ between the mass of the dark halo and the final mass in stars. SCUBA sources are predicted to be strongly clustered, with a clustering strength increasing with mass. We show that the model accounts for the clustering of Lyman-break Galaxies, seen as the optical counterpart of low- to intermediate-mass primeval spheroidal galaxies and is also consistent with the observed angular correlation function of Extremely Red Objects. Best agreement is obtained for $M_{\rm halo}/M_{\rm sph}=100$. We also consider the implications for small scale fluctuations observed at sub-mm wavelengths by current or forthcoming experiments aimed at mapping the Cosmic Microwave Background (CMB). The predicted amplitude of the clustering signal in the 350 GHz channel of the Planck mission strongly depends on the halo-to-bulge mass ratio and may be of comparable amplitude to primary CMB anisotropies for multipole numbers $l\simgt 50$.

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.

We present the results of a set of numerical simulations aiming at evaluating the effects of cluster galaxies on the arc statistics. At this goal we use nine different galaxy clusters obtained from N-body simulations. We mimic the presence of a population of galaxies inside each cluster, trying to reproduce the observed luminosity function and the spatial distribution. We compare the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of the properties of long arcs.

(abridged) We present a model to predict the clustering properties of X-ray clusters in flux-limited surveys. Our technique correctly accounts for past light-cone effects on the observed clustering and follows the non-linear evolution in redshift of the underlying DM correlation function and cluster bias factor. The conversion of the limiting flux of a survey into the corresponding minimum mass of the hosting DM haloes is obtained by using theoretical and empirical relations between mass, temperature and X-ray luminosity of clusters. Finally, our model is calibrated to reproduce the observed cluster counts adopting a temperature-luminosity relation moderately evolving with redshift. We apply our technique to three existing catalogues: BCS, XBACs and REFLEX samples. Moreover, we consider an example of possible future space missions with fainter limiting flux. In general, we find that the amplitude of the spatial correlation function is a decreasing function of the limiting flux and that the EdS models always give smaller correlation amplitudes than open or flat models with low matter density parameter. In the case of XBACs, the comparison with previous estimates of the observational spatial correlation shows that only the predictions of models with Omega_0m=0.3 are in good agreement with the data, while the EdS models have too low a correlation strength. Finally, we use our technique to discuss the best strategy for future surveys. Our results show that the choice of a wide area catalogue, even with a brighter limiting flux, is preferable to a deeper, but with smaller area, survey.

We analyse the spatial clustering properties of the RASS1 Bright Sample, an X-ray flux-limited catalogue of galaxy clusters selected from the southern part of the $ROSAT$ All-Sky Survey. The two-point correlation function $\xi(r)$ of the whole sample is well fitted (in an Einstein-de Sitter model) by the power-law $\xi=(r/r_0)^{-\gamma}$, with $r_0= 21.5^{+3.4}_{-4.4} h^{-1}$ Mpc and $\gamma=2.11^{+0.53}_{-0.56}$ (95.4 per cent confidence level with one fitting parameter). We use the RASS1 Bright Sample as a first application of a theoretical model which aims at predicting the clustering properties of X-ray clusters in flux-limited surveys for different cosmological scenarios. The model uses the theoretical and empirical relations between mass, temperature and X-ray cluster luminosity, and fully accounts for the redshift evolution of the underlying dark matter clustering and cluster bias factor. The comparison between observational results and theoretical predictions shows that the Einstein-de Sitter models display too low a correlation length, while models with a matter density parameter $\Omega_{\rm 0m}=0.3$ (with or without a cosmological constant) are successful in reproducing the observed clustering. The dependence of the correlation length $r_0$ on the X-ray limiting flux and luminosity of the sample is generally consistent with the predictions of all our models. Quantitative agreement is however only reached for $\Omega_{\rm 0m} = 0.3$ models. The model presented here can be reliably applied to future deeper X-ray cluster surveys: the study of their clustering properties will provide a useful complementary tool to the traditional cluster abundance analyses to constrain the cosmological parameters.

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics. We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N-body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White (1997) density profile. We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12'', while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4''.

We present a theoretical model which aims at predicting the clustering properties of X-ray clusters in flux-limited surveys for different cosmological scenarios. The model uses the theoretical and empirical relations between mass, temperature and X-ray luminosity, and fully accounts for the redshift evolution of the underlying dark matter clustering and cluster bias factor. We apply the model to the RASS1 Bright Sample and to the XBACs catalogue. The results show that the Einstein-de Sitter models display too low a correlation length, while models with a matter density parameter Omega_0m=0.3 (with or without a cosmological constant) are successful in reproducing the observed clustering.

(abridged) We present the results of a redshift survey of intercluster galaxies in the central region of the Shapley Concentration supercluster, aimed at determining the distribution of galaxies in between obvious overdensities. Our sample is formed by 442 new redshifts, mainly in the b_J magnitude range 17-18.8. Adding the data from our redshift surveys on the A3558 and A3528 complexes, which are close to the geometrical centre of this supercluster, we obtain a total sample of ~2000 radial velocities. Using the 1440 galaxies of our total sample in the magnitude range 17 - 18.8, we reconstructed the density profile in the central part of the Shapley Concentration; moreover we detected another significant overdensity at ~30000 km/s (dubbed S300). We estimate the total overdensity in galaxies, the mass and the dynamical state of these structures, and discuss the effect of considering a bias between the galaxy distribution and the underlying matter. We find an indication that the value of the bias between clusters and galaxies in the Shapley Concentration is higher that that reported in literature, confirming the impression that this supercluster is very rich in clusters. Finally, from the comparison with some theoretical scenarios, we find that the existence of the Shapley Concentration is more consistent with the predictions of the models with a matter density parameter <1, such as open CDM and Lambda CDM.

(abridged) The evolution of galaxy clustering from z=0 to z=4.5 is analyzed using the angular correlation function and the photometric redshift distribution of galaxies brighter than I_{AB}\le 28.5 in the HDF North. The reliability of the photometric redshift estimates is discussed on the basis of the available spectroscopic redshifts, comparing different codes and investigating the effects of photometric errors. The redshift bins in which the clustering properties are measured are then optimized to take into account the uncertainties of the photometric redshifts. The results show that the comoving correlation length has a small decrease in the range 0<z<1 followed by an increase at higher z. We compare these results with the theoretical predictions of a variety of cosmological models belonging to the general class of CDM. The comparison with the expected mass clustering evolution indicates that the observed high-redshift galaxies are biased tracers of the dark matter with an effective bias b strongly increasing with redshift. Assuming an Einstein-de Sitter universe, we obtain b\simeq 2 at z=2 and b\simeq 5 at z=4. A comparison of the clustering amplitudes that we measured at z=3 with those reported for LBG suggests that the clustering depends on the abundance of the objects: more abundant objects are less clustered, as expected in the paradigm of hierarchical galaxy formation. The strong clustering and high bias measured at z=3 are consistent with the expected density of massive haloes predicted for the various cosmologies here considered. At z=4, the strong clustering observed in the HDF requires a significant fraction of massive haloes to be already formed by that epoch. This feature could be a discriminant test for the cosmological parameters if confirmed by future observations.

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

We calculate the statistical clustering of Lyman-break galaxies predicted in a selection of currently-fashionable structure formation scenarios. These models are all based on the cold dark matter model, but vary in the amount of dark matter, the initial perturbation spectrum, the background cosmology and in the presence or absence of a cosmological constant term. If Lyman-break galaxies form as a result of hierarchical merging, the amplitude of clustering depends quite sensitively on the minimum halo mass that can host such a galaxy. Interpretation of the recent observations by Giavalisco et al. (1998) would therefore be considerably clarified by a direct determination of the relevant halo properties. For a typical halo mass around $10^{11} h^{-1} M_\odot$ the observations do not discriminate strongly between cosmological models, but if the appropriate mass is larger, say $10^{12} ~h^{-1} M_\odot$ (which seems likely on theoretical grounds), then the data strongly favour models with a low matter-density.

We discuss the theoretical interpretation of observational data concerning the clustering of galaxies at high redshifts. Building on the theoretical machinery developed by Matarrese et al. (1997), we make detailed quantitative predictions of galaxy clustering statistics for a variety of cosmological models, taking into account differences in spatial geometry and initial fluctuation spectra and exploring the role of bias as a complicating factor in these calculations. We demonstrate that the usual description of evolution (in terms of the parameters $\epsilon$ and $r_0$) is not useful for realistic galaxy clustering models. We compare the detailed predictions of the variation of correlation functions with redshift against current observational data to constrain available models of structure formation. Theories that fit the present-day abundance of rich clusters are generally compatible with the observed redshift evolution of galaxy clustering if galaxies are no more than slightly biased at $z\sim 1$. We also discuss the interpretation of a concentration of Lyman-break galaxies found by Steidel et al. (1998), coming to the conclusion that such concentrations are not unexpected in `standard' models of structure formation.

The observed clustering of Lyman-$\alpha$ lines is reviewed and compared with the clustering of CIV systems. We argue that a continuity of properties exists between Lyman-$\alpha$ and metal systems and show that the small-scale clustering of the absorbers is consistent with a scenario of gravitationally induced correlations. At large scales statistically significant over and under-densities (including voids) are found on scales of tens of Mpc.

We study the properties of X-ray galaxy clusters in four cold dark matter models with different baryon fraction $\Omega_{BM}$ ranging from 5 to 20 per cent. By using an original three-dimensional hydrodynamic code based on the piecewise parabolic method, we run simulations on a box with size $64 h^{-1}$ Mpc and we identify the clusters by selecting the peaks in the X-ray luminosity field. We analyse these mock catalogues by computing the mass function, the luminosity function, the temperature distribution and the luminosity-temperature relation. By comparing the predictions of the different models to a series of recent observational results, we find that only the models with low baryonic content agree with the data, while models with larger baryon fraction are well outside the 1-$\sigma$ errorbars. In particular, the analysis of the luminosity functions, both bolometric and in the energy band [0.5--2] keV, requires $\Omega_{BM}\mincir 0.05$ when we fix the values $h=0.5$ and $n=0.8$ for the Hubble parameter and the primordial spectral index, respectively. Moreover we find that, independently of the cosmological scenario, all the considered quantities have a very little redshift evolution, particularly between $z=0.5$ and $z=0$.

We study the geometry and topology of the large-scale structure traced by galaxy clusters in numerical simulations of a box of side 320 $h^{-1}$ Mpc, and compare them with available data on real clusters. The simulations we use are generated by the Zel'dovich approximation, using the same methods as we have used in the first three papers in this series. We consider the following models to see if there are measurable differences in the topology and geometry of the superclustering they produce: (i) the standard CDM model (SCDM); (ii) a CDM model with $\Omega_0=0.2$ (OCDM); (iii) a CDM model with a `tilted' power spectrum having $n=0.7$ (TCDM); (iv) a CDM model with a very low Hubble constant, $h=0.3$ (LOWH); (v) a model with mixed CDM and HDM (CHDM); (vi) a flat low-density CDM model with $\Omega_0=0.2$ and a non-zero cosmological $\Lambda$ term ($\Lambda$CDM). We analyse these models using a variety of statistical tests based on the analysis of: (i) the Euler-Poincar\'{e} characteristic; (ii) percolation properties; (iii) the Minimal Spanning Tree construction. Taking all these tests together we find that the best fitting model is $\Lambda$CDM and, indeed, the others do not appear to be consistent with the data. Our results demonstrate that despite their biased and extremely sparse sampling of the cosmological density field, it is possible to use clusters to probe subtle statistical diagnostics of models which go far beyond the low-order correlation functions usually applied to study superclustering.

We present a hydrodynamical code for cosmological simulations which uses the Piecewise Parabolic Method (PPM) to follow the dynamics of gas component and an N-body Particle-Mesh algorithm for the evolution of collisionless component. The gravitational interaction between the two components is regulated by the Poisson equation which is solved by a standard FFT procedure. In order to simulate cosmological flows we have introduced several modifications to the original PPM scheme which we describe in detail. Various tests of the code are presented including adiabatic expansion, single and multiple pancake formation and three-dimensional cosmological simulations with initial conditions based on the cold dark matter scenario.

(Shortened Abstract). We analyze a redshift sample of Abell/ACO clusters and compare them with numerical simulations based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, on which we estimate cosmic variance effects. We analyse correlation statistics, the probability density function, and supercluster properties from percolation analysis. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 < \Gamma < 0.25 and an effective spectral index at 20Mpc/h in the range [-1.1,-0.9] are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press--Schechter approach, whose reliability is explicitly tested against N--body simulations. We conclude that, of the cosmological models considered here, the only viable models are either Cold+Hot DM ones with \Omega_\nu = [0.2-0.3], better if shared between two massive neutrinos, and flat low-density CDM models with \Omega_0 = [0.3-0.5].

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

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

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

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.

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 $\Omega_0=1$, while the other two have $\Omega_0=0.2$, one with and one without a cosmological constant term. For the Abell/ACO sample, we evaluate the parameter $\beta=\Omega_0^{0.6}/b_{ cl}$, which reaches its asymptotic value at $R=160 Mpc/h$. We find $\beta=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 $\beta$, 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.

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 \Omega_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.

Recently, there has been growing interest on primordial ``blue" ($n>1$) perturbation spectra, motivated both by a composite set of observational data on large scales and from the point of view of theoretical model building. After reviewing the theoretical (inflationary) motivations for these blue spectra, we consider various observational constraints, within both Cold Dark Matter and Mixed Dark Matter scenarios. In particular, using linear theory, we discuss large--scale bulk flows and the X--ray cluster abundance. We also perform various N--body simulations of these models to study the clustering properties of the matter distribution and the peculiar velocity field.

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

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

We estimate the variance and the skewness of the cluster distribution in several dark matter (DM) models. The cluster simulations are based on the Zel'dovich approximation, the low computational cost of which allows us to run 50 random realizations of each model. We compare our results with those coming from a similar analysis of a redshift sample of Abell/ACO clusters. Within the list of the considered models, we find that only the model based on Cold+Hot DM (with $\Omega_{\rm hot}=0.3$) provides a good fit to the data. The standard CDM model and the low-density ($\Omega_{\circ}=0.2$) CDM models, both with and without a cosmological constant term ($\Omega_\Lambda =0.8$) are ruled out. The tilted CDM model with primordial spectral index $n=0.7$ and a low Hubble constant ($h=0.3$) CDM model are only marginally consistent with the data.

We present extended simulations of the distribution of galaxy clusters in different dark matter models, using an optimized version of the Zel'dovich approximation. Six different dark matter models are studied: Standard CDM, Open CDM, Tilted CDM, Cold + Hot DM, low-density CDM with cosmological constant, and low Hubble constant ($h=0.3$) CDM. We compare simulations with an Abell/ACO redshift sample. We find that the models that best reproduce the clustering of the real data sets are the Cold + Hot DM model and the CDM model with a cosmological constant. The probability density function of all models and data is always well approximated by a lognormal distribution. We also find that the linear biasing parameter for the simulated clusters is nearly constant over a large range of scales, but its value depends on the model. We also note that the abundances of clusters predicted by the Press-Schechter theory provide strong constraints on these models: only the Cold + Hot DM model and the low-H CDM model appear to produce the correct number-density of clusters. Taking this constraint together with the cluster clustering statistics leads one to conclude that the best of our models is the Cold + Hot DM scenario.

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.

We analyze the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We perform a statistical analysis, by calculating the bulk flow, Cosmic Mach Number and velocity correlation function for this sample and for mock catalogs extracted from a set of N--body simulations. We run four cold dark matter (CDM) simulations: two tilted models (with spectral index $n=0.6$ and $n=0.8$), the standard model ($n=1$) and a ``blue" one ($n=1.2$), with different values of the linear bias parameter $b$. By means of a Maximum Likelihood analysis we estimate the ability of our models to fit the observations, as measured by the above statistics, and to reproduce the Local Group properties. On the basis of this analysis we conclude that the best model is the unbiased standard model $(n,~b)=(1,~1)$, even though the overall flatness of the joint likelihood function implies that one cannot strongly discriminate models in the range $0.8 \le n \le 1$, and $1 \le b \le 1.5$. Models with $b \geq 2.5$ are rejected at the $95\%$ confidence level. For $n=0.8$ the values of $b$ preferred by the present analysis, together with the {\em COBE} data, require a negligible contribution to $\Delta T/T$ by gravitational waves. Finally, the blue model, normalized to {\em COBE}, does not provide a good fit to the velocity data.

We discuss the non--linear growth of the kurtosis of the smoothed peculiar velocity field (along an arbitrary direction), in an Einstein--de Sitter universe, induced by Gaussian primordial density fluctuations. Applying the perturbative theory, we show that, for different cosmological models, a departure from the original Gaussian distribution is gravitationally induced only on small scales ($\simlt 20-30$ Mpc). Models with scale--free power spectrum $P(k)\propto k^{\,n}$, with $-1< n \leq 1$, are also considered. When the fluid particles move according to the Zel'dovich approximation the probability distribution of the peculiar velocity field remains unaltered during the evolution.

{}From three quasar samples with a total of 1038 objects in the redshift range $1.0 \div 2.2$ we measure the variance $\sigma^2$ of counts in cells of volume $V_u$. By a maximum likelihood analysis applied separately on these samples we obtain estimates of $\sigma^2(\ell)$, with $\ell \equiv V_u^{1/3}$. The analysis from a single catalog for $\ell = ~40~h^{-1}$ Mpc and from a suitable average over the three catalogs for $\ell = ~60,~80$ and $100~h^{-1}$ Mpc, gives $\sigma^2(\ell) = 0.46^{+0.27}_{-0.27}$, $0.18^{+0.14}_{-0.15}$, $0.05^{+0.14}_{-0.05}$ and $0.12^{+0.13}_{-0.12}$, respectively, where the $70\%$ confidence ranges account for both sampling errors and statistical fluctuations in the counts. This allows a comparison of QSO clustering on large scales with analogous data recently obtained both for optical and IRAS galaxies: QSOs seem to be more clustered than these galaxies by a biasing factor $b_{QSO}/b_{gal} \sim 1.4 - 2.3$.

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.

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.

Some recently proposed approximations to follow the non--linear evolution of collisionless matter perturbations in the universe are reviewed. The first one, called frozen--flow approximation, is an Eulerian method within Newtonian theory, and is based on neglecting the role of particle inertia compared to the damping implied by the Hubble drag. The second approach is General Relativistic and Lagrangian; it is based on following the evolution of fluid and geometric observables in the rest frame of each fluid element, under the only approximation of neglecting the back--reaction on the system of the gravitational radiation emitted during non--linear collapse.

We study the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We carry out a statistical analysis, by calculating bulk flows and velocity correlation functions for this sample and for mock catalogs which we extract from N--body simulations. For the simulations we consider tilted (i.e. with spectral index $n\leq 1$) CDM models with different values of the linear bias parameter $b$. By mean of a maximum likelihood analysis we estimate the ability of the models in fitting the observations as measured by the above statistics and in reproducing the Local Group properties.

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

We discuss the non--linear growth of the excess kurtosis parameter of the smoothed density fluctuation field $\delta$, $S_4\equiv[\lan\delta^{\,4}\ran-3\lan\delta^{\,2}\ran^2]/ \lan\delta^{\,2}\ran^3$ in an Einstein--de Sitter universe. We assume Gaussian primordial density fluctuations with scale--free power spectrum $P(k)\propto k^{\,n}$ and analyze the dependence of $S_4$ on primordial spectral index $n$, after smoothing with a Gaussian filter. As already known for the skewness ratio $S_3$, the kurtosis parameter is a {\it decreasing function} of $n$, both in exact perturbative theory and in the Zel'dovich approximation. The parameter $S_4$ provides a powerful statistics to test different cosmological scenarios.

We calculate reduced moments $\overline \xi_q$ of the matter density fluctuations, up to order $q=5$, from counts in cells produced by Particle--Mesh numerical simulations with scale--free Gaussian initial conditions. We use power--law spectra $P(k) \propto k^n$ with indices $n=-3,~-2,~-1,~0,~1$. Due to the supposed absence of characteristic times or scales in our models, all quantities are expected to depend on a single scaling variable. For each model, the moments at all times can be expressed in terms of the variance $\overline \xi_2$, alone. We look for agreement with the hierarchical scaling ansatz, according to which $\overline \xi_q \propto \overline \xi_2^{~q-1}$. For $n\leq -2$ models we find strong deviations from the hierarchy, which are mostly due to the presence of boundary problems in the simulations. A small, residual signal of deviation from the hierarchical scaling is however also found in $n \geq -1$ models. For the first time, due to our large dynamic range and careful checks of scaling and shot--noise effects, we are able to detect evolution away from the perturbation theory result.

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.

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.

We generate artificial Lick maps using $N$--body simulations and compare the angular correlation function, $w(\vartheta)$, measured from the simulations with the APM correlation. For the Gaussian CDM model, neither the standard biassed model nor a more evolved model (as suggested by the COBE data), reproduce the correlations on large angular scales. We come to a similar conclusion about CDM models with positively skewed initial fluctuation distributions. In contrast, models with initially negatively skewed fluctuations produce a $w(\vartheta)$ that declines much more gently on large scales. Such models are therefore in principle capable of reconciling the lack of large--scale power of the CDM spectrum with the observed clustering of APM galaxies.