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.

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

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

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.

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.

The two-point correlation function of the galaxy distribution is a key cosmological observable that allows us to constrain the dynamical and geometrical state of our Universe. To measure the correlation function we need to know both the galaxy positions and the expected galaxy density field. The expected field is commonly specified using a Monte-Carlo sampling of the volume covered by the survey and, to minimize additional sampling errors, this random catalog has to be much larger than the data catalog. Correlation function estimators compare data-data pair counts to data-random and random-random pair counts, where random-random pairs usually dominate the computational cost. Future redshift surveys will deliver spectroscopic catalogs of tens of millions of galaxies. Given the large number of random objects required to guarantee sub-percent accuracy, it is of paramount importance to improve the efficiency of the algorithm without degrading its precision. We show both analytically and numerically that splitting the random catalog into a number of subcatalogs of the same size as the data catalog when calculating random-random pairs, and excluding pairs across different subcatalogs provides the optimal error at fixed computational cost. For a random catalog fifty times larger than the data catalog, this reduces the computation time by a factor of more than ten without affecting estimator variance or bias.

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.

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.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

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.

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.

In the currently debated context of using clusters of galaxies as cosmological probes, the need for well-defined cluster samples is critical. The XXL Survey has been specifically designed to provide a well characterised sample of some 500 X-ray detected clusters suitable for cosmological studies. The main goal of present article is to make public and describe the properties of the cluster catalogue in its present state, as well as of associated catalogues as super-clusters and fossil groups. We release a sample containing 365 clusters in total. We give the details of the follow-up observations and explain the procedure adopted to validate the cluster spectroscopic redshifts. Considering the whole XXL cluster sample, we have provided two types of selection, both complete in a particular sense: one based on flux-morphology criteria, and an alternative based on the [0.5-2] keV flux within one arcmin of the cluster centre. We have also provided X-ray temperature measurements for 80$\%$ of the clusters having a flux larger than 9$\times$10$^{-15}$$\rm \thinspace erg \, s^{-1} \, cm^{-2}$. Our cluster sample extends from z$\sim$0 to z$\sim$1.2, with one cluster at z$\sim$2. Clusters were identified through a mean number of six spectroscopically confirmed cluster members. Our updated luminosity function and luminosity-temperature relation are compatible with our previous determinations based on the 100 brightest clusters, but show smaller uncertainties. We also present an enlarged list of super-clusters and a sample of 18 possible fossil groups. This intermediate publication is the last before the final release of the complete XXL cluster catalogue when the ongoing C2 cluster spectroscopic follow-up is complete. It provides a unique inventory of medium-mass clusters over a 50~\dd\ area out to z$\sim$1.

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.

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 present new numerical tools to analyse cosmic void catalogues, implemented inside the CosmoBolognaLib, a large set of Open Source C++/Python numerical libraries. The CosmoBolognaLib provides a common numerical environment for cosmological calculations. This work extends these libraries by adding new algorithms for cosmological analyses of cosmic voids, covering the existing gap between theory and observations. We implemented new methods to model the size function of cosmic voids, in both observed and simulated samples of dark matter and biased tracers. Moreover, we provide new numerical tools to construct unambiguous void catalogues. The latter are designed to be independent of the void finder, in order to allow a high versatility in comparing independent results. The implemented Open Source software is available at the GitHub repository https://github.com/federicomarulli/CosmoBolognaLib. We provide also a full doxygen documentation and some example codes that explain how to use these libraries.

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

Combining measurements on the expansion history of the Universe and on the growth rate of cosmic structures is key to discriminate between alternative cosmological frameworks and to test gravity. Recently, Linder (2017) proposed a new diagram to investigate the joint evolutionary track of these two quantities. In this letter, we collect the most recent cosmic growth and expansion rate datasets to provide the state-of-the-art observational constraints on this diagram. By performing a joint statistical analysis of both probes, we test the standard $\Lambda$CDM model, confirming a mild tension between cosmic microwave background predictions from Planck mission and cosmic growth measurements at low redshift ($z<2$). Then we test alternative models allowing the variation of one single cosmological parameter at a time. In particular, we find a larger growth index than the one predicted by general relativity $\gamma=0.65^{+0.05}_{-0.04}$). However, also a standard model with total neutrino mass of $0.26\pm0.10$ eV provides a similarly accurate description of the current data. By simulating an additional dataset consistent with next-generation dark-energy mission forecasts, we show that growth rate constraints at $z>1$ will be crucial to discriminate between alternative models.

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.

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.

We investigate to which accuracy it is possible to recover the real-space two-point correlation function of galaxy clusters from cluster catalogues based on photometric redshifts, and test our ability to measure the redshift and mass evolution of the correlation length and the bias parameter as a function of the redshift uncertainty. We calculate the correlation function for cluster sub-samples covering various mass and redshift bins selected from a light-cone catalogue. To simulate the distribution of clusters in photometric redshift space, we assign to each cluster a redshift randomly extracted from a Gaussian distribution. The dispersion is varied in the range $\sigma_{(z=0)} = 0.001$ to $0.050$. The correlation function in real-space is computed through estimation and deprojection of $w_{p}(r_{p})$. Four mass ranges (from $M_{halo}> 2 \times 10^{13}$ to $M_{halo}> 2 \times 10^{14}$) and six redshift slices covering the redshift range [0,2] are investigated, using cosmological redshifts and photo-z configurations. We find a clear increase of the correlation amplitude as a function of redshift and mass for the $z_{c}$ samples. The evolution of the derived bias parameter is in agreement with theoretical expectations. From our pilot sample limited to $M_{halo}> 5 \times 10^{13} (0.4 < z < 0.7)$, we find that the real-space correlation function can be recovered by deprojection of $w_{p}(r_{p})$ within an accuracy of 5% for $\sigma_{z} = 0.001 \times (1 + z_{c})$ and within 10% for $\sigma_{z} = 0.03 \times (1 + z_{c})$. The evolution of the correlation in redshift and mass is clearly detected for all $\sigma_{z}$ tested. The best-fit parameters $(r_{0}$ and ${\gamma})$ as well as the bias obtained from the deprojection method for all $\sigma_{z}$ are within the $1 \sigma$ uncertainty of the $z_{c}$ sample.

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

We exploit the synergy between low-resolution spectroscopy and photometric redshifts to study environmental effects on galaxy evolution in slitless spectroscopic surveys from space. As a test case, we consider the future Euclid Deep survey (~40deg$^2$), which combines a slitless spectroscopic survey limited at H$\alpha$ flux $\geq5\times 10^{-17}$ erg cm$^{-2}$ s$^{-1}$ and a photometric survey limited in H-band ($H\leq26$). We use Euclid-like galaxy mock catalogues, in which we anchor the photometric redshifts to the 3D galaxy distribution of the available spectroscopic redshifts. We then estimate the local density contrast by counting objects in cylindrical cells with radius from 1 to 10 h$^{-1}$Mpc over the redshift range 0.9<z<1.8. We compare this density field with the one computed in a mock catalogue with the same depth as the Euclid Deep survey (H=26) but without redshift measurement errors. We find that our method successfully separates high from low density environments (the last from the first quintile of the density distribution), with higher efficiency at low redshift and large cell: the fraction of low density regions mistaken by high density peaks is <1% for all scales and redshifts explored, but for scales of 1 h$^{-1}$Mpc for which is a few percent. These results show that we can efficiently study environment in photometric samples if spectroscopic information is available for a smaller sample of objects that sparsely samples the same volume. We demonstrate that these studies are possible in the Euclid Deep survey, i.e. in a redshift range in which environmental effects are different from those observed in the local universe, hence providing new constraints for galaxy evolution models.

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.

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

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

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.

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.

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

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

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

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.

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.

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

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.

By using a suite of large box-size N-body simulations that incorporate massive neutrinos as an extra set of particles, we investigate the impact of neutrino masses on the spatial distribution of dark matter haloes and galaxies. We compute the bias between the spatial distribution of dark matter haloes and the overall matter and cold dark matter distributions using statistical tools such as the power spectrum and the two-point correlation function. Overall we find a scale-dependent bias on large scales for the cosmologies with massive neutrinos. However, our results indicate that the scale-dependence in the bias is reduced if the latter is computed with respect to the cold dark matter distribution only. We find that the value of the bias on large scales is reasonably well reproduced by the Tinker fitting formula once the linear cold dark matter power spectrum is used, instead of the total matter power spectrum. We investigate whether scale-dependent bias really comes from purely neutrino's effect or from nonlinear gravitational collapse of haloes. For this purpose, we address the $\Omega_\nu$-$\sigma_8$ degeneracy and find that such degeneracy is not perfect, implying that neutrinos imprint a slight scale dependence on the large-scale bias. Finally, by using a simple halo occupation distribution (HOD) model, we investigate the impact of massive neutrinos on the distribution of galaxies within dark matter haloes. We use the main galaxy sample in the Sloan Digital Sky Survey II Data Release 7 to investigate if the small-scale galaxy clustering alone can be used to discriminate among different cosmological models with different neutrino masses. Our results suggest that different choices of the HOD parameters can reproduce the observational measurements relatively well, and we quantify the difference between the values of the HOD parameters between massless and massive neutrino cosmologies.

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.

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.

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.

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

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

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 measure the real-space galaxy power spectrum on large scales at redshifts 0.5 to 1.2 using optical colour-selected samples from the CFHT Legacy Survey. With the redshift distributions measured with a preliminary ~14000 spectroscopic redshifts from the VIMOS Public Extragalactic Redshift Survey (VIPERS), we deproject the angular distribution and directly estimate the three-dimensional power spectrum. We use a maximum likelihood estimator that is optimal for a Gaussian random field giving well-defined window functions and error estimates. This measurement presents an initial look at the large-scale structure field probed by the VIPERS survey. We measure the galaxy bias of the VIPERS-like sample to be b_g=1.38 +- 0.05 (sigma_8=0.8) on scales k<0.2h/mpc averaged over 0.5<z<1.2. We further investigate three photometric redshift slices, and marginalising over the bias factors while keeping other LCDM parameters fixed, we find the matter density Omega_m=0.30+-0.06.

There is mounting evidence that a significant fraction of Black Holes (BHs) today live in late-type galaxies, including bulge-less galaxies and those hosting pseudobulges, and are significantly undermassive with respect to the scaling relations followed by their counterpart BHs in classical bulges of similar stellar (or even bulge) mass. Here we discuss the predictions of two state-of-the-art hierarchical galaxy formation models in which BHs grow via mergers and, in one, also via disk instability. Our aim is to understand if the wealth of new data on local BH demography is consistent with standard models. We follow the merger trees of representative subsamples of BHs and compute the fractional contributions of different processes to the final BH mass. We show that the model in which BHs always closely follow the growth of their host bulges, also during late disk instabilities (i.e., bars), produces too narrow a distribution of BHs at fixed stellar mass to account for the numerous low-mass BHs now detected in later-type galaxies. Models with a looser connection between BH growth and bar instability instead predict the existence of a larger number of undermassive BHs, in better agreement with the observations. The scatter in the updated local BH-bulge mass relation (with no restriction on galaxy type) appears to be quite large when including later-type systems, but it can still be managed to be reproduced within current hierarchical models. However, the fuelling of BHs during the late bar-instability mode needs to be better quantified/improved to properly fit the data. We conclude discussing how the possibly large number of BHs in later type galaxies demands for an in-depth revision of the local BH mass function and its modelling.

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.

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.

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.

We discuss some of the main open issues in the evolution of Active Galactic Nuclei which can be solved by the sensitive, wide area surveys to be performed by the proposed Wide Field X-ray Telescope mission.

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

The relation between the mass of supermassive black holes located in the center of the host galaxies and the kinetic energy of random motions of the corresponding bulges can be reinterpreted as an age-temperature diagram for galaxies. This relation fits the experimental data better than the M_bh-M_G, M_bh-L_G, and M_bh-sigma laws. The validity of this statement has been confirmed by using three samples extracted from different catalogues of galaxies. In the framework of the LambdaCDM cosmology our relation has been compared with the predictions of two galaxy formation models based on the Millennium Simulation.

We provide fits to the distribution of galaxy luminosity, size, velocity dispersion and stellar mass as a function of concentration index C_r and morphological type in the SDSS. We also quantify how estimates of the fraction of `early' or `late' type galaxies depend on whether the samples were cut in color, concentration or light profile shape, and compare with similar estimates based on morphology. Our fits show that Es account for about 20% of the r-band luminosity density, rho_Lr, and 25% of the stellar mass density, rho_*; including S0s and Sas increases these numbers to 33% and 40%, and 50% and 60%, respectively. Summed over all galaxy types, we find rho_* ~ 3 * 10^8 M_Sun Mpc^{-3} at z ~ 0. This is in good agreement with expectations based on integrating the star formation history. However, compared to most previous work, we find an excess of objects at large masses, up to a factor of ~ 10 at M_* ~ 5*10^{11} M_Sun. The stellar mass density further increases at large masses if we assume different IMFs for Es and spiral galaxies, as suggested by some recent chemical evolution models, and results in a better agreement with the dynamical mass function. We also show that the trend for ellipticity to decrease with luminosity is primarily because the E/S0 ratio increases at large L. However, the most massive galaxies, M_* > 5 * 10^{11} M_Sun, are less concentrated and not as round as expected if one extrapolates from lower L, and they are not well-fit by pure deVaucouleur laws. This suggests formation histories with recent radial mergers. Finally, we show that the age-size relation is flat for Es of fixed dynamical mass, but, at fixed M_dyn, S0s and Sas with large sizes tend to be younger. Explaining this difference between E and S0 formation is a new challenge for models of early-type galaxy formation.

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

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

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.

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