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.

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

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

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

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.

We test the history of structure formation from redshift 1 to today by matching galaxies from VIPERS and SDSS with dark matter haloes in the MultiDark SMDPL N-body simulation. We first show that the standard subhalo abundance matching (SHAM) recipe implemented with MultiDark characterizes the clustering of galaxies well both at redshift 0 for SDSS and at redshift 1 for VIPERS. This is an important validation of the SHAM model at high redshift. We then remap the simulation timesteps to test alternative growth histories and infer the growth index $\gamma=0.6\pm0.3$. This analysis demonstrates the power of using N-body simulations to forward model galaxy surveys for cosmological inference.

We investigate the clustering properties of radio sources in the Alternative Data Release 1 of the TIFR GMRT Sky Survey (TGSS), focusing on large angular scales, where previous analyses have detected a large clustering signal. After appropriate data selection, the TGSS sample we use contains ~110,000 sources selected at 150 MHz over ~70% of the sky. The survey footprint is largely superimposed on that of the NRAO VLA Sky Survey (NVSS) with the majority of TGSS sources having a counterpart in the NVSS sample. These characteristics make TGSS suitable for large-scale clustering analyses and facilitate the comparison with the results of previous studies. In this analysis we focus on the angular power spectrum, although the angular correlation function is also computed to quantify the contribution of multiple-component radio sources. We find that on large angular scales, corresponding to multipoles $2 \leq \ell \leq 30$, the amplitude of the TGSS angular power spectrum is significantly larger than that of the NVSS. We do not identify any observational systematic effects that may explain this mismatch. We have produced a number of physically motivated models for the TGSS angular power spectrum and found that all of them fail to match observations, even when taking into account observational and theoretical uncertainties. The same models provide a good fit to the angular spectrum of the NVSS sources. These results confirm the anomalous nature of the TGSS large-scale power, which has no obvious physical origin and seems to indicate that unknown systematic errors are present in the TGSS dataset.

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.

A new fully non-linear reconstruction algorithm for the accurate recovery of the Baryonic Acoustic Oscillations (BAO) scale in two-point correlation functions is proposed, based on the least-action principle and extending the Fast Action Minimisation method by Nusser & Branchini (2000). Especially designed for massive spectroscopic surveys, it is tested on dark-matter halo catalogues extracted from the DEUS-FUR $\Lambda$CDM simulation to trace the trajectories of up to $\sim207,000$ haloes backward-in-time, well beyond the first-order Lagrangian approximation. The new algorithm successfully recovers the BAO feature in real and redshift-space in both the monopole and the anisotropic two-point correlation function, also for anomalous samples showing misplaced or absent signature of BAO. In redshift-space the non-linear displacement parameter $\Sigma_\mathrm{NL}$ is reduced from $11.8\pm0.3h^{-1}$Mpc at redshift $z=0$ (pre-reconstruction) to $4.0\pm0.5h^{-1}$Mpc at $z\simeq37$ after reconstruction. A comparison with the first-order Lagrangian reconstruction is presented, showing that this techniques outperforms the linear approximation in recovering an unbiased measurement of the acoustic scale.

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.

It has been known for decades that the observed number of baryons in the local universe falls about 30-40% short of the total number of baryons predicted by Big-Bang Nucleosynthesis, as inferred from density fluctuations of the Cosmic Microwave Background and seen during the first 2-3 billion years of the universe in the so called Lyman-alpha Forest. A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas - hydrogen - is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios. Indeed, despite the large observational efforts, only a few marginal claims of detection have been made so far. Here we report observations of two absorbers of highly ionized oxygen (OVII) in the high signal-to-noise-ratio X-ray spectrum of a quasar at redshift >0.4. These absorbers show no variability over a 2-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The OVII systems lie in regions characterized by large (x4 compared to average) galaxy over-densities and their number (down to the sensitivity threshold of our data), agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium (WHIM). We conclude that the missing baryons have been found.

Using the almost all-sky 2MASS Photometric Redshift catalogue (2MPZ) we perform for the first time a tomographic analysis of galaxy angular clustering in the local Universe ($z<0.24$). We estimate the angular auto- and cross-power spectra of 2MPZ galaxies in three photometric redshift bins, and use dedicated mock catalogues to assess their errors. We measure a subset of cosmological parameters, having fixed the others at their Planck values, namely the baryon fraction $f_{b}=0.14^{+0.09}_{-0.06}$, the total matter density parameter $\Omega_{\rm m}=0.30\pm 0.06$, and the effective linear bias of 2MPZ galaxies $b_{\rm eff}$, which grows from $1.1^{+0.3}_{-0.4}$ at $\langle z \rangle=0.05$ up to $2.1^{+0.3}_{-0.5}$ at $\langle z\rangle=0.2$, largely because of the flux-limited nature of the dataset. The results obtained here for the local Universe agree with those derived with the same methodology at higher redshifts, and confirm the importance of the tomographic technique for next-generation photometric surveys such as Euclid or LSST.

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

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

We investigate the nature of the extragalactic unresolved gamma-ray background (UGRB) by cross-correlating several galaxy catalogs with sky-maps of the UGRB built from 78 months of Pass 8 Fermi-Large Area Telescope data. This study updates and improves similar previous analyses in several aspects. Firstly, the use of a larger gamma-ray dataset allows us to investigate the energy dependence of the cross-correlation in more detail, using up to 8 energy bins over a wide energy range of [0.25-500] GeV. Secondly, we consider larger and deeper catalogs (2MASS Photometric-Redshift catalog, 2MPZ; WISE x SuperCOSMOS, WISC; and SDSS-DR12 photometric-redshift dataset) in addition to the ones employed in the previous studies (NVSS and SDSS-QSOs). Thirdly, we exploit the redshift information available for the above catalogs to divide them into redshift bins and perform the cross-correlation separately in each of them. Our results confirm, with higher statistical significance, the detection of cross-correlation signal between the UGRB maps and all the catalogs considered, on angular scales smaller than 1 degree. Significances range from 16.3 sigma for NVSS, 7 sigma for SDSS-DR12 and WISC, 5 sigma for 2MPZ and 4 sigma for SDSS-QSOs. Furthermore, including redshift tomography, the significance of the SDSS-DR12 signal strikingly rises up to 12 sigma and the one of WISC to 10.6 sigma. We offer a simple interpretation of the signal in the framework of the halo model. The precise redshift and energy information allows us to clearly detect a change over redshift in the spectral and clustering behavior of the gamma-ray sources contributing to the UGRB.

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

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.

Observed galaxy luminosities (derived from redshifts) hold information on the large-scale peculiar velocity field in the form of spatially correlated scatter, which allows for bounds on bulk flows and the growth rate of matter density perturbations using large galaxy redshift surveys. We apply this luminosity approach to galaxies from the recent SDSS Data Release 13. Our goal is twofold. First, we take advantage of the recalibrated photometry to identify possible systematic errors relevant to our previous analysis of earlier data. Second, we seek improved constraints on the bulk flow and the normalized growth rate f$\sigma_{8}$ at z ~ 0.1. Our results confirm the robustness of our method. Bulk flow amplitudes, estimated in two redshift bins with 0.02 < z$_{1}$ < 0.07 < z$_{2}$ < 0.22, are generally smaller than in previous measurements, consistent with both the updated photometry and expectations for the $\Lambda$CDM model. The obtained growth rate, f$\sigma_{8}$ = 0.48 +/- 0.16, is larger than, but still compatible with, its previous estimate, and closer to the reference value of Planck. Rather than precision, the importance of these results is due to the fact that they follow from an independent method that relies on accurate photometry, which is a top requirement for next-generation photometric catalogs.

[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 report the detection of a cross-correlation signal between {\it Fermi} Large Area Telescope diffuse gamma-ray maps and catalogs of clusters. In our analysis, we considered three different catalogs: WHL12, redMaPPer and PlanckSZ. They all show a positive correlation with different amplitudes, related to the average mass of the objects in each catalog, which also sets the catalog bias. The signal detection is confirmed by the results of a stacking analysis. The cross-correlation signal extends to rather large angular scales, around 1 degree, that correspond, at the typical redshift of the clusters in these catalogs, to a few to tens of Mpc, i.e. the typical scale-length of the large scale structures in the Universe. Most likely this signal is contributed by the cumulative emission from AGNs associated to the filamentary structures that converge toward the high peaks of the matter density field in which galaxy clusters reside. In addition, our analysis reveals the presence of a second component, more compact in size and compatible with a point-like emission from within individual clusters. At present, we cannot distinguish between the two most likely interpretations for such a signal, i.e. whether it is produced by AGNs inside clusters or if it is a diffuse gamma-ray emission from the intra-cluster medium. We argue that this latter, intriguing, hypothesis might be tested by applying this technique to a low redshift large mass cluster sample.

We forecast the future constraints on scale-dependent parametrizations of galaxy bias and their impact on the estimate of cosmological parameters from the power spectrum of galaxies measured in a spectroscopic redshift survey. For the latter we assume a wide survey at relatively large redshifts, similar to the planned Euclid survey, as baseline for future experiments. To assess the impact of the bias we perform a Fisher matrix analysis and we adopt two different parametrizations of scale-dependent bias. The fiducial models for galaxy bias are calibrated using a mock catalogs of H$\alpha$ emitting galaxies mimicking the expected properties of the objects that will be targeted by the Euclid survey. In our analysis we have obtained two main results. First of all, allowing for a scale-dependent bias does not significantly increase the errors on the other cosmological parameters apart from the rms amplitude of density fluctuations, $\sigma_{8}$, and the growth index $\gamma$, whose uncertainties increase by a factor up to two, depending on the bias model adopted. Second, we find that the accuracy in the linear bias parameter $b_{0}$ can be estimated to within 1-2\% at various redshifts regardless of the fiducial model. The non-linear bias parameters have significantly large errors that depend on the model adopted. Despite of this, in the more realistic scenarios departures from the simple linear bias prescription can be detected with a $\sim2\,\sigma$ significance at each redshift explored. Finally, we use the Fisher Matrix formalism to assess the impact of assuming an incorrect bias model and found that the systematic errors induced on the cosmological parameters are similar or even larger than the statistical ones.

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 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 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 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 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 compare the measured angular cross-correlation between the Fermi-LAT gamma-ray sky and catalogues of extra-galactic objects with the expected signal induced by weakly interacting massive particle (WIMP) dark matter (DM). We include a detailed description of the contribution of astrophysical gamma-ray emitters such as blazars, misaligned AGN and star forming galaxies, and perform a global fit to the measured cross-correlation. Five catalogues are considered: SDSS-DR6 quasars, 2MASS galaxies, NVSS radio galaxies, SDSS-DR8 Luminous Red Galaxies and SDSS-DR8 main galaxy sample. To model the cross-correlation signal we use the halo occupation distribution formalism to estimate the number of galaxies of a given catalogue in DM halos and their spatial correlation properties. We discuss uncertainties in the predicted cross-correlation signal arising from the DM clustering and WIMP microscopic properties, which set the DM gamma-ray emission. The use of different catalogues probing objects at different redshifts reduces significantly, though not completely, the degeneracy among the different gamma-ray components. We find that the presence of a significant WIMP DM signal is allowed by the data but not significantly preferred by the fit, although this is mainly due to a degeneracy with the misaligned AGN component. With modest substructure boost, the sensitivity of this method excludes thermal annihilation cross sections at 95% C.L. for WIMP masses up to few tens of GeV. Constraining the low-redshift properties of astrophysical populations with future data will further improve the sensitivity to DM.

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

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.

The observed distribution of galaxies has local transverse isotropy around the line-of- sight (LOS) with respect to the observer. The difference in the statistical clustering signal along and across the line-of-sight encodes important information about the ge- ometry of the Universe, its expansion rate and the rate of growth of structure within it. Because the LOS varies across a survey, the standard Fast Fourier Transform (FFT) based methods of measuring the Anisotropic Power-Spectrum (APS) cannot be used for surveys with wide observational footprint, other than to measure the monopole moment. We derive a simple analytic formula to quantify the bias for higher-order Legendre moments and we demonstrate that it is scale independent for a simple sur- vey model, and depends only on the observed area. We derive a similar numerical correction formula for recently proposed alternative estimators of the APS that are based on summing over galaxies rather than using an FFT, and can therefore in- corporate a varying LOS. We demonstrate that their bias depends on scale but not on the observed area. For a quadrupole the bias is always less than 1 per cent for k > 0.01h/Mpc at z > 0.32. For a hexadecapole the bias is below 5 per cent for k>0.05h/Mpc at z>0.32.

We propose a new approach to the missing baryons problem. Building on the common assumption that the missing baryons are in the form of the Warm Hot Intergalactic Medium (WHIM), we further assumed here that the galaxy luminosity density can be used as a tracer of the WHIM. The latter assumption is supported by our finding of a significant correlation between the WHIM density and the galaxy luminosity density in the hydrodynamical simulations of Cui et al. (2012). We further found that the fraction of the gas mass in the WHIM phase is substantially (by a factor of $\sim$1.6) higher within the large scale galactic filaments, i.e. $\sim$70\%, compared to the average in the full simulation volume of $\sim$0.1\,Gpc$^3$. The relation between the WHIM overdensity and the galaxy luminosity overdensity within the galactic filaments is consistent with linear: $\delta_{\rm whim}\,=\,0.7\,\pm\,0.1\,\times\,\delta_\mathrm{LD}^{0.9 \pm 0.2}$. We applied our procedure to the line of sight to the blazar H2356-309 and found evidence for the WHIM in correspondence of the Sculptor Wall (z $\sim$0.03 and $\log{N_H}$ = $19.9^{+0.1}_{-0.3}$) and Pisces-Cetus superclusters (z $\sim$0.06 and $\log{N_H}$ = $19.7^{+0.2}_{-0.3}$), in agreement with the redshifts and column densities of the X-ray absorbers identified and studied by Fang et al. (2010) and Zappacosta et al. (2010). This agreement indicates that the galaxy luminosity density and galactic filaments are reliable signposts for the WHIM and that our method is robust in estimating the WHIM density. The signal that we detected cannot originate from the halos of the nearby galaxies since they cannot account for the large WHIM column densities that our method and X-ray analysis consistently find in the Sculptor Wall and Pisces-Cetus superclusters.

Spatial variations in the distribution of galaxy luminosities, estimated from redshifts as distance proxies, are correlated with the peculiar velocity field. Comparing these variations with the peculiar velocities inferred from galaxy redshift surveys is a powerful test of gravity and dark energy theories on cosmological scales. Using ~ 2 $\times$ 10$^{5}$ galaxies from the SDSS Data Release 7, we perform this test in the framework of gravitational instability to estimate the normalized growth rate of density perturbations f$\sigma_{8}$ = 0.37 +/- 0.13 at z ~ 0.1, which is in agreement with the $\Lambda$CDM scenario. This unique measurement is complementary to those obtained with more traditional methods, including clustering analysis. The estimated accuracy at z ~ 0.1 is competitive with other methods when applied to similar datasets.

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.

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.

If dark matter (DM) is composed by particles which are non-gravitationally coupled to ordinary matter, their annihilations or decays in cosmic structures can result in detectable radiation. We show that the most powerful technique to detect a particle DM signal outside the Local Group is to study the angular cross-correlation of non-gravitational signals with low-redshift gravitational probes. This method allows to enhance signal-to-noise from the regions of the Universe where the DM-induced emission is preferentially generated. We demonstrate the power of this approach by focusing on GeV-TeV DM and on the recent cross-correlation analysis between the 2MASS galaxy catalogue and the Fermi-LAT gamma-ray maps. We show that this technique is more sensitive than other extragalactic gamma-ray probes, such as the energy spectrum and angular autocorrelation of the extragalactic background, and emission from clusters of galaxies. Intriguingly, we find that the measured cross-correlation can be well fitted by a DM component, with thermal annihilation cross section and mass between 10 and 100 GeV, depending on the small-scale DM properties and gamma-ray production mechanism. This solicits further data collection and dedicated analyses.

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.

Building on our previous cross-correlation analysis (Xia et al. 2011) between the isotropic gamma-ray background (IGRB) and different tracers of the large-scale structure of the universe, we update our results using 60-months of data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. We perform a cross-correlation analysis between the IGRB and objects that may trace the astrophysical sources of the IGRB: SDSS-DR6 QSOs, the SDSS-DR8 Main Galaxy Sample, Luminous Red Galaxies (LRGs) in the SDSS catalog, 2MASS galaxies, and radio NVSS galaxies. The benefit of correlating the Fermi-LAT signal with catalogs of objects at various redshifts is to provide tomographic information on the IGRB which is crucial to separate the various contributions and to clarify its origin. We observe a significant (>3.5 sigma) cross-correlation signal on angular scales smaller than 1 deg in the NVSS, 2MASS and QSO cases and, at lower statistical significance (~3.0 sigma), with SDSS galaxies. These results are robust against the choice of the statistical estimator, estimate of errors, map cleaning procedure and instrumental effects. Finally, we test the hypothesis that the IGRB observed by Fermi-LAT originates from the summed contributions of three types of unresolved extragalactic sources: BL Lacs, FSRQs and Star-Forming Galaxies (SFGs). We find that a model in which the IGRB is mainly produced by SFGs ($72^{+23}_{-37}$% with 2 sigma errors), with BL Lacs and FSRQs giving a minor contribution, provides a good fit to the data. We also consider a possible contribution from Misaligned Active Galactic Nuclei, and we find that, depending on the details of the model and its uncertainty, they can also provide a substantial contribution, partly degenerate with the SFG one. (abridged)

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

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

We have developed a new method to approach the missing baryons problem. We assume that the missing baryons reside in a form of Warm Hot Intergalactic Medium, i.e. the WHIM. Our method consists of (a) detecting the coherent large scale structure in the spatial distribution of galaxies that traces the Cosmic Web and that in hydrodynamical simulations is associated to the WHIM, (b) map its luminosity into a galaxy luminosity density field, (c) use numerical simulations to relate the luminosity density to the density of the WHIM, (d) apply this relation to real data to trace the WHIM using the observed galaxy luminosities in the Sloan Digital Sky Survey and 2dF redshift surveys. In our application we find evidence for the WHIM along the line of sight to the Sculptor Wall, at redshifts consistent with the recently reported X-ray absorption line detections. Our indirect WHIM detection technique complements the standard method based on the detection of characteristic X-ray absorption lines, showing that the galaxy luminosity density is a reliable signpost for the WHIM. For this reason, our method could be applied to current galaxy surveys to optimise the observational strategies for detecting and studying the WHIM and its properties. Our estimates of the WHIM hydrogen column density in Sculptor agree with those obtained via the X-ray analysis. Due to the additional column density estimate, our method has potential for improving the constrains of the physical parameters of the WHIM as derived with X-ray absorption, and thus for improving the understanding of the missing baryons problem.

Peculiar motion introduces systematic variations in the observed luminosity distribution of galaxies. This allows one to constrain the cosmic peculiar velocity field from large galaxy redshift surveys. Using around half a million galaxies from the SDSS Data Release 7 at z ~ 0.1, we demonstrate the applicability of this approach to large datasets and obtain bounds on peculiar velocity moments and $\sigma_{8}$, the amplitude of the linear matter power spectrum. Our results are in good agreement with the $\Lambda$CDM model and consistent with the previously reported ~ 1% zero-point tilt in the SDSS photometry. Finally, we discuss the prospects of constraining the growth rate of density perturbations by reconstructing the full linear velocity field from the observed galaxy clustering in redshift space.

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.

Spatial modulations in the distribution of observed luminosities (computed using redshifts) of ~ 5 $\times$ 10$^5$ galaxies from the SDSS Data Release 7, probe the cosmic peculiar velocity field out to z ~ 0.1. Allowing for luminosity evolution, the r-band luminosity function, determined via a spline-based estimator, is well represented by a Schechter form with M$^{\star}$(z) - 5log$_{10}$h = -20.52 - 1.6(z - 0.1) $\pm$ 0.05 and $\alpha^{\star}$ = -1.1 $\pm$ 0.03. Bulk flows and higher velocity moments in two redshift bins, 0.02 < z < 0.07 and 0.07 < z < 0.22, agree with the predictions of the $\Lambda$CDM model, as obtained from mock galaxy catalogs designed to match the observations. Assuming a $\Lambda$CDM model, we estimate $\sigma_{8}$ $\approx$ 1.1 $\pm$ 0.4 for the amplitude of the linear matter power spectrum, where the low accuracy is due to the limited number of galaxies. While the low-z bin is robust against coherent photometric uncertainties, the bias of results from the second bin is consistent with the ~ 1% magnitude tilt reported by the SDSS collaboration. The systematics are expected to have a significantly lower impact in future datasets with larger sky coverage and better photometric calibration.

We report the discovery of a luminosity distance estimator using Active Galactic Nuclei (AGN). We combine the correlation between the X-ray variability amplitude and the Black Hole (BH) mass with the single epoch spectra BH mass estimates which depend on the AGN luminosity and the line width emitted by the broad line region. We demonstrate that significant correlations do exist which allows one to predict the AGN (optical or X-ray) luminosity as a function of the AGN X-ray variability and either the HBeta or the PaBeta line widths. In the best case, when the PaBeta is used, the relationship has an intrinsic dispersion of ~0.6 dex. Although intrinsically more disperse than Supernovae Ia, this relation constitutes an alternative distance indicator potentially able to probe, in an independent way, the expansion history of the Universe. With this respect, we show that the new mission concept Athena should be able to measure the X-ray variability of hundreds of AGN and then constrain the distance modulus with uncertainties of 0.1 mag up to z~0.6. We also discuss how, using a new dedicated wide field X-ray telescope able to measure the variability of thousands of AGNs, our estimator has the prospect to become a cosmological probe even more sensitive than current Supernovae Ia samples.

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

There is a $\sim 150 km s^{-1}$ discrepancy between the measured motion of the Local Group of galaxies (LG) with respect to the CMB and the linear theory prediction based on the gravitational force field of the large scale structure in full-sky redshift surveys. We perform a variety of tests which show that the LG motion cannot be recovered to better than $150-200 km s^{-1}$ in amplitude and within a $\approx10^\circ$ in direction. The tests rely on catalogs of mock galaxies identified in the Millennium simulation using semi-analytic galaxy formation models. We compare these results to the $K_s=11.75$ Two-Mass Galaxy Redshift Survey, which provides the deepest, widest and most complete spatial distribution of galaxies available so far. In our analysis we use a new, concise relation for deriving the LG motion and bulk flow from the true distribution of galaxies in redshift space. Our results show that the main source of uncertainty is the small effective depth of surveys like the 2MRS that prevents a proper sampling of the large scale structure beyond $\sim100 h^{-1} Mpc$. Deeper redshift surveys are needed to reach the "convergence scale" of $\approx 250 h^{-1}Mpc $ in a $\Lambda$CDM universe. Deeper survey would also mitigate the impact of the "Kaiser rocket" which, in a survey like 2MRS, remains a significant source of uncertainty. Thanks to the quiet and moderate density environment of the LG, purely dynamical uncertainties of the linear predictions are subdominant at the level of $\sim 90 km s^{-1}$. Finally, we show that deviations from linear galaxy biasing and shot noise errors provide a minor contribution to the total error budget.

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 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 backbone of the large-scale structure of the Universe is determined by processes on a cosmological scale and by the gravitational interaction of the dominant dark matter. However, the mobile baryon population shapes the appearance of these structures. Theory predicts that most of the baryons reside in vast unvirialized filamentary structures that connect galaxy groups and clusters, but the observational evidence is currently lacking. Because the majority of the baryons are supposed to exist in a large-scale, hot and dilute gaseous phase, X-rays provide the ideal tool to progress our understanding. Observations with the Athena+ X-ray Integral Field Unit will reveal the location, chemical composition, physical state and dynamics of the active population of baryons.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

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 present the first results from our pilot 500 ks Chandra-LETG Large Program observation of the soft X-ray brightest source in the z>=0.4 sky, the blazar 1ES 1553+113, aimed to secure the first uncontroversial detections of the missing baryons in the X-rays. We identify a total of 11 possible absorption lines, with single-line statistical significances between 2.2-4.1\sigma. Six of these lines are detected at high single-line statistical significance (3.6 <= sigma <= 4.1), while the remaining five are regarded as marginal detections in association with either other X-ray lines detected at higher significance and/or Far-Ultraviolet (FUV) signposts. In particular, five of these possible intervening absorption lines, are identified as CV and CVI K\alpha absorbers belonging to three WHIM systems at z_X = 0.312, z_X = 0.237 and <z_X> = 0.133, which also produce broad HI (and OVI for the z_X = 0.312 system) absorption in the FUV. For two of these systems (z_X = 0.312 and 0.237), the Chandra X-ray data led the a-posteriori discovery of physically consistent broad HI associations in the FUV, so confirming the power of the X-ray-FUV synergy for WHIM studies. The true statistical significances of these three X-ray absorption systems, after properly accounting for the number of redshift trials, are 5.8 sigma (z_X = 0.312; 6.3 sigma if the low-significance OV and CV K-beta associations are considered), 3.9 sigma (z_X = 0.237), and 3.8 sigma (\langle z_X \rangle = 0.133), respectively.

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

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

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

The cosmological peculiar velocity field (deviations from the pure Hubble flow) of matter carries significant information on dark energy, dark matter and the underlying theory of gravity on large scales. Peculiar motions of galaxies introduce systematic deviations between the observed galaxy redshifts z and the corresponding cosmological redshifts z_cos. A novel method for estimating the angular power spectrum of the peculiar velocity field based on observations of galaxy redshifts and apparent magnitudes m (or equivalently fluxes) is presented. This method exploits the fact that a mean relation between z_cos and m of galaxies can be derived from all galaxies in a redshift-magnitude survey. Given a galaxy magnitude, it is shown that the z_cos(m) relation yields its cosmological redshift with a 1-sigma error of sigma_z~0.3 for a survey like Euclid (~10^9 galaxies at z<~2), and can be used to constrain the angular power spectrum of z-z_cos(m) with a high signal-to-noise ratio. At large angular separations corresponding to l<~15, we obtain significant constraints on the power spectrum of the peculiar velocity field. At 15<~l<~60, magnitude shifts in the z_cos(m) relation caused by gravitational lensing magnification dominate, allowing us to probe the line-of-sight integral of the gravitational potential. Effects related to the environmental dependence in the luminosity function can easily be computed and their contamination removed from the estimated power spectra. The amplitude of the combined velocity and lensing power spectra at z~1 can be measured with <~5% accuracy.

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.

In this Letter we point out that redshift surveys can break the degeneracy between the galaxy bias, the power spectrum normalization, \sigma_{8,0} and the growth factor, without the need for external information by using a simple and rather general parametrization for the growth rate, the well known \gamma-parametrization and measuring the power spectrum at least at two different redshifts. We find that in next-generation surveys like Euclid, \sigma_{8,0} and \gamma can be measured to within 1% and 5%, respectively, while the bias b(z) can be measured to within 1-2% in each of 14 equal-width redshift bins spanning 0.7<=z<=2.

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 constructed a Bayesian hyper-parameter statistical method to quantify the difference between predicted velocities derived from the observed galaxy distribution in the \textit{IRAS}-PSC$z$ redshift survey and peculiar velocities measured using different distance indicators. In our analysis we find that the model--data comparison becomes unreliable beyond $70 \hmpc$ because of the inadequate sampling by \textit{IRAS} survey of prominent, distant superclusters, like the Shapley Concentration. On the other hand, the analysis of the velocity residuals show that the PSC$z$ gravity field provides an adequate model to the local, $\le 70 \hmpc$, peculiar velocity field. The hyper-parameter combination of ENEAR, SN, A1SN and SFI++ catalogues in the Bayesian framework constrains the amplitude of the linear flow to be $\beta=0.53 \pm 0.014$. For an rms density fluctuations in the PSC$z$ galaxy number density $\sigma_8^{\rm gal}=0.42\pm0.03$, we obtain an estimate of the growth rate of density fluctuations $f\sigma_{8}(z\sim0) = 0.42 \pm 0.033$, which is in excellent agreement with independent estimates based on different techniques.

Using redshifts as a proxy for galaxy distances, estimates of the 2D transverse peculiar velocities of distant galaxies could be obtained from future measurements of proper motions. We provide the mathematical framework for analyzing 2D transverse motions and show that they offer several advantages over traditional probes of large scale motions. They are completely independent of any intrinsic relations between galaxy properties, hence they are essentially free of selection biases. They are free from homogeneous and inhomogeneous Malmquist biases that typically plague distance indicator catalogs. They provide additional information to traditional probes which yield line-of-sight peculiar velocities only. Further, because of their 2D nature, fundamental questions regarding vorticity of large scale flows can be addressed. Gaia for example is expected to provide proper motions of at least bright galaxies with high central surface brightness, making proper motions a likely contender traditional probes based on current and future distance indicator measurements.

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

Using the nearly full sky Ks=11.75 2MASS Redshift Survey [2MRS]of ~45,000 galaxies we reconstruct the underlying peculiar velocity field and constrain the cosmological bulk flow within ~100. These results are obtained by maximizing the probability to estimate the absolute magnitude of a galaxy given its observed apparent magnitude and redshift. At a depth of ~60 Mpc/h we find a bulk flow Vb=(90\pm65,-230\pm65,50\pm65) km/s in agreement with the theoretical predictions of the LCDM model. The reconstructed peculiar velocity field that maximizes the likelihood is characterized by the parameter beta=0.323 +/- 0.08. Both results are in agreement with those obtained previously using the ~23,000 galaxies of the shallower Ks=11.25 2MRS survey. In our analysis we find that the luminosity function of 2MRS galaxies is poorly fitted by the Schechter form and that luminosity evolves such that objects become fainter with increasing redshift according to L(z)=L(z=0)(1+z)^(+2.7 +/-0.15).

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.

In this work, starting from 21 months of data from the Fermi-Large Area Telescope, we derive maps of the residual isotropic gamma-ray emission, a relevant fraction of which is expected to be contributed by the extragalactic diffuse gamma-ray background. We compute the angular two-point auto-correlation function of the residual Fermi-LAT maps at energies E>1GeV, E>3GeV and E>30GeV well above the Galactic plane and find no significant correlation signal. This is, indeed, what is expected if the EGB were contributed by BL Lacertae, Flat Spectrum Radio Quasars or star-forming galaxies, since, in this case, the predicted signal is very weak. Then, we search for the Integrated Sachs-Wolfe signature by cross-correlating the Fermi-LAT maps with the WMAP7-Cosmic Microwave Background map. We find a cross-correlation consistent with zero, even though the expected signal is larger than that of the EGB auto-correlation. Finally, in an attempt to constrain the nature of the gamma-ray background we cross-correlate the Fermi-LAT maps with the angular distributions of objects that may contribute to the EGB: QSOs in the SDSS-DR6 catalog, NVSS galaxies, 2MASS galaxies and LRG in the SDSS catalog. The cross-correlation is always consistent with zero, in agreement with theoretical expectations, but we find (with low statistical significance) some interesting features that may indicate that some specific classes of objects contribute to the EGB. A chi2 analysis confirms that the correlation properties of the 21-month data do not provide strong constraints of the EGB origin. However, the results suggest that the situation will significantly improve with the 5- and 10-year Fermi-LAT data. The future EGB analysis will then allow placing significant constraints on the nature of the EGB and might provide in addition a detection of the ISW signal. (Abridged)

Given a redshift survey of galaxies with measurements of apparent magnitudes, we present a novel method for measuring the growth rate $f(\Omega)$ of cosmological linear perturbations. We use the galaxy distribution within the survey to solve for the peculiar velocity field which depends in linear perturbation theory on $\beta=f(\Omega)/b$, where $b$ is the bias factor of the galaxy distribution. The recovered line-of-sight peculiar velocities are subtracted from the redshifts to derive the distances, which thus allows an estimate of the absolute magnitude of each galaxy. A constraint on $\beta$ is then found by minimizing the spread of the estimated magnitudes from their distribution function. We apply the method to the all sky $K = 11.25$ Two-MASS Redhsift Survey (2MRS) and derive $\beta=0.35\pm 0.1$ at $z\sim 0$, remarkably consistent with our previous estimate from the velocity-velocity comparison. The method could easily be applied to subvolumes extracted from the SDSS survey to derive the growth rate at $z \sim 0.1$. Further, it should also be applicable to ongoing and future spectroscopic redshift surveys to trace the evolution of $f(\Omega)$ to $z\sim1$. Constraints obtained from this method are entirely independent from those obtained from the two-dimensional distortion of $\xi(s)$ and provide an important check on $f(\Omega)$, as alternative gravity models predict observable differences.

We present a simple method for measuring cosmological bulk flows from large redshift surveys, based on the apparent dimming or brightening of galaxies due to their peculiar motion. It is aimed at estimating bulk flows of cosmological volumes containing large numbers of galaxies. Constraints on the bulk flow are obtained by minimizing systematic variations in galaxy luminosities with respect to a reference luminosity function measured from the whole survey. This method offers two advantages over more popular bulk flow estimators: it is independent of error-prone distance indicators and of the poorly-known galaxy bias. We apply the method to the 2MASS redshift survey (2MRS) to measure the local bulk flows of spherical shells centered on the Milky Way (MW). The result is consistent with that obtained by Nusser and Davis (2011) using the SFI++ catalogue of Tully-Fisher distance indicators. We also make an assessment of the ability of the method to constrain bulk flows at larger redshifts ($z=0.1-0.5$) from next generation datasets. As a case study we consider the planned EUCLID survey. Using this method we will be able to measure a bulk motion of $ \sim 200 \kms$ of $10^6$ galaxies with photometric redshifts, at the $3\sigma$ level for both $z\sim 0.15$ and $z\sim 0.5$. Thus the method will allow us to put strong constraints on dark energy models as well as alternative theories for structure formation.

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

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

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.

Many experiments in the near future will test dark energy through its effects on the linear growth of matter perturbations. In this paper we discuss the constraints that future large-scale redshift surveys can put on three different parameterizations of the linear growth factor and how these constraints will help ruling out different classes of dark energy and modified gravity models. We show that a scale-independent bias can be estimated to a few percent per redshift slice by combining redshift distortions with power spectrum amplitude, without the need of an external estimation. We find that the growth rate can be constrained to within 2-4% for each $\Delta z=0.2$ redshift slice, while the equation of state $w$ and the index $\gamma$ can be simultaneously estimated both to within 0.02. We also find that a constant dimensionless coupling between dark energy and dark matter can be constrained to be smaller than 0.14.

We study the prospects for detecting the annihilation products of Dark Matter [DM] in the framework of the two highest-resolution numerical simulations currently available, i.e. {\it Via Lactea II} and {\it Aquarius}. We propose a strategy to determine the shape and size of the region around the Galactic center that maximizes the probability of observing a DM signal, and we show that although the predicted flux can differ by a factor of 10 for a given DM candidate in the two simulation setups, the search strategy remains actually unchanged, since it relies on the angular profile of the annihilation flux, not on its normalization. We present mock gamma-ray maps that keep into account the diffuse emission produced by unresolved halos in the Galaxy, and we estimate that in an optimistic DM scenario a few individual clumps can be resolved above the background with the Fermi-LAT. Finally we calculate the energy-dependent boost factors for positrons and antiprotons, and show that they are always of $\cal O$(1), and therefore they cannot lead to the large enhancements of the antimatter fluxes required to explain the recent PAMELA, ATIC, Fermi and HESS data.

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

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

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

We study the flux and the angular power spectrum of gamma-rays produced by Dark Matter (DM) annihilations in the Milky Way (MW) and in extra-galactic halos. The annihilation signal receives contributions from: a) the smooth MW halo, b) resolved and unresolved substructures in the MW, c) external DM halos at all redshifts, including d) their substructures. Adopting a self-consistent description of local and extra-galactic substructures, we show that the annihilation flux from substructures in the MW dominates over all the other components for angles larger than O(1) degrees from the Galactic Center, unless an extreme prescription is adopted for the substructures concentration. We also compute the angular power spectrum of gamma-ray anisotropies and find that, for an optimistic choice of the particle physics parameters, an interesting signature of DM annihilations could soon be discovered by the Fermi LAT satellite at low multipoles, l<100, where the dominant contribution comes from MW substructures with mass M>10^4 solar masses. For the substructures models we have adopted, we find that the contribution of extra-galactic annihilations is instead negligible at all scales.

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

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

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

We use high-resolution hydrodynamical simulations to investigate the spatial correlation between weak (N(HI) <10^{15} cm^{-2}) Ly-alpha absorbers and gas-rich galaxies in the local universe. We confirm that Ly-alpha absorbers are preferentially expected near gas-rich galaxies and that the degree of correlation increases with the column density of the absorber. The real-space galaxy auto-correlation is stronger than the cross-correlation (correlation lengths r_{0,gg}=3.1 \pm 0.1 Mpc/h and r_{0,ag}=1.4 \pm 0.1 Mpc/h, respectively), in contrast with the recent results of Ryan-Weber (2006, RW06), and the auto-correlation of absorbers is very weak. These results are robust to the presence of strong galactic winds in the hydrodynamical simulations. In redshift-space a further mismatch arises since at small separations the distortion pattern of the simulated galaxy-absorber cross-correlation function is different from the one measured by RW06. However, when sampling the intergalactic medium along a limited number of lines-of-sight, as in the real data, uncertainties in the cross correlation estimates are large enough to account for these discrepancies. Our analysis suggests that the statistical significance of difference between the cross-correlation and auto-correlation signal in current datasets is ~ 1-sigma only.

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

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

Upcoming $\gamma$-ray satellites will search for Dark Matter annihilations in Milky Way substructures (or 'clumps'). The prospects for detecting these objects strongly depend on the assumptions made on the distribution of Dark Matter in substructures, and on the distribution of substructures in the Milky Way halo. By adopting simplified, yet rather extreme, prescriptions for these quantities, we compute the number of sources that can be detected with upcoming experiments such as GLAST, and show that, for the most optimistic particle physics setup ($m_\chi=40$ GeV and annihilation cross section $\sigma v = 3 \times 10^{-26}$ cm$^3$ s$^{-1}$), the result ranges from zero to $\sim$ hundred sources, all with mass above $10^{5}M\odot$. However, for a fiducial DM candidate with mass $m_\chi=100$ GeV and $\sigma v = 10^{-26}$ cm$^3$ s$^{-1}$, at most a handful of large mass substructures can be detected at $5 \sigma$, with a 1-year exposure time, by a GLAST-like experiment. Scenarios where micro-clumps (i.e. clumps with mass as small as $10^{-6}M\odot$) can be detected are severely constrained by the diffuse $\gamma$-ray background detected by EGRET.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

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

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.

Analytical calculations and recent numerical experiments have shown that a sizable of the mass in our Galaxy is in a form of clumpy, virialized substructures that, according to \cite{dms:05}, can be as light as $10^{-6} \msun$. In this work we estimate the gamma-rays flux expected from dark matter annihilation occurring within these minihalos, under the hypothesis that the bulk of dark matter is composed by neutralinos. We generate mock sky maps showing the angular distribution of the expected gamma-ray signal. We compare them with the sensitivities of satellite-borne experiments such as GLAST and find that a possible detection of minihalos is indeed very challenging.

We assess the extent to which the PSCz survey encapsulates the complete or major share of matter inhomogeneities responsible for the external tidal forces affecting the peculiar velocity flow within the Local Supercluster [LS]. We investigate this issue on the basis of mock galaxy catalogs drawn from two large Nbody simulations of cosmic structure formation in two different cosmological scenarios. From these a variety of datasets is selected imitating the observational conditions of either the volume-limited LS mimicking NBG catalog or the deeper magnitude-limited PSCz catalog. The mildly nonlinear dynamics in the mock NBG and PSCz velocities are analyzed by means of an optimized implementation of the Least Action Principle, the Fast Action Method [FAM]. By comparing the modeled and true velocities of the corresponding mock catalogs we assess the extent and nature of the external tidal influence on the LS volume. We find that the dynamics in the inner 30Mpc/h volume is strongly affected by the external forces. Most of the external forces can be traced back to a depth of no more than 100Mpc/h. This is concluded from the fact that the FAM reconstructions of the 100Mpc/h PSCz volume appear to have included most gravitational influences. In addition, we demonstrate that for all considered cosmological models the bulk flow and shear components of the tidal velocity field generated by the external distribution of PSCz galaxies provides sufficient information for representing the full external tidal force field.

In the frame of indirect dark matter searches we investigate the flux of high-energy $\gamma$-ray photons produced by annihilation of dark matter in caustics within our Galaxy under the hypothesis that the bulk of dark matter is composed of the lightest supersymmetric particles. Unfortunately, the detection of the caustics annihilation signal with currently available instruments is rather challenging. Indeed, with realistic assumptions concerning particle physics and cosmology, the $\gamma $-ray signal from caustics is below the detection threshold of both $\check {\rm C}$erenkov telescopes and satellite-borne experiments. Nevertheless, we find that this signal is more prominent than that expected if annihilation only occurs in the smoothed Galactic halo, with the possible exception of a $\sim 15^{\circ}$ circle around the Galactic center if the mass density profile of our Galaxy exhibits a sharp cusp there. We show that the angular distribution of this $\gamma$-ray flux changes significantly if DM annihilation preferentially occurs within virialized sub-halos populating our Galaxy rather than in caustics.

We present a simple method for tracing the spatial distribution and predicting the physical properties of the Warm-Hot Intergalactic Medium (WHIM), from the map of galaxy light in the local universe. Under the assumption that biasing is local and monotonic we map the ~ 2 Mpc/h smoothed density field of galaxy light into the mass density field from which we infer the spatial distribution of the WHIM in the local supercluster. Taking into account the scatter in the WHIM density-temperature and density-metallicity relation, extracted from the z=0 outputs of high-resolution and large box size hydro-dynamical cosmological simulations, we are able to quantify the probability of detecting WHIM signatures in the form of absorption features in the X-ray spectra, along arbitrary directions in the sky. To illustrate the usefulness of this semi-analytical method we focus on the WHIM properties in the Virgo Cluster region.

We investigate the clustering properties of galaxies in the recently completed ELAIS-S1 redshift survey through their spatial two point autocorrelation function. We used a sub-sample of the ELAIS-S1 catalog covering approximately 4 deg^2 and consisting of 148 objects selected at 15 micron with a flux >0.5 mJy and redshift z<0.5. We detected a positive signal in the correlation function that, in the range of separations 1-10 h mpc is well approximated by a power law with a slope gamma = 1.4 and a correlation length s =5.4 h mpc. This result is in good agreement with the redshift-space correlation function measured in more local samples of mid infrared selected galaxies like the IRAS PSC$z$ redshift survey. This suggests a lack of significant clustering evolution of infrared selected objects out to $z=0.5$ that is further confirmed by the consistency found between the correlation functions measured in a local (z<0.2) and a distant (0.2<z<0.5) subsample of ELAIS-S1 galaxies. We also confirm that optically selected galaxies in the local redshift surveys, especially those of the SDSS sample, are significantly more clustered than infrared objects.

We use the PSCz galaxy redshift catalogue combined with constrained simulations based on the IRAS 1.2 Jy galaxy density field to estimate the contribution of hot gas in the local universe to the SZ-effect on large scales. We produce a full sky Healpix map predicting the SZ-effect from clusters as well as diffuse hot gas within 80Mpc/h. Performing cross-correlation tests between this map and the WMAP data in pixel, harmonic and wavelet space we can put an upper limit on the effect. We conclude that the SZ effect from diffuse gas in the local universe cannot be detected in current CMB data and is not a contaminating factor on large scales (l<60) in studies of the CMB angular anisotropies. However, for future high sensitivity experiments observing at a wider range of frequencies, the predicted large scale SZ effect could be of importance.

We introduce a simple linear equation relating the line-of-sight peculiar velocity and density contrast correlation functions. The relation, which we call the "Gaussian cell two-point energy-like equation", is valid in the distant-observer-limit and requires Gaussian smoothed fields. In the variance case, i.e., at zero lag, the equation is similar in its mathematical form to the Layzer-Irvine cosmic energy equation. \beta estimation with this equation from the PSC$z$ redshift galaxy survey and the SEcat catalogue of peculiar velocities is carried out, returning a value of \beta = 0.44+-0.08. The applicability of the method for the 6dF galaxy redshift and peculiar motions survey is demonstrated with mock data where it is shown that beta could be determined with ~5% accuracy. The prospects for constraining the dark energy equation of state with this method from the kinematic and thermal Sunyaev-Zel'dovich cluster surveys are discussed. The equation is also used to construct a nonparametric mass density power spectrum estimator from peculiar velocity data.

We analyze in detail the peculiar velocity field traced by 56 clusters within 120 h^-1 Mpc in the Streaming Motions of Abell Clusters (SMAC) sample. The bulk flow of the SMAC sample is 687 +- 203 km/s, toward l = 260 +- 13, b = 0 +- 11. We discuss possible systematic errors and show that no systematic effect is larger than half of the random error. The flow does not drop off significantly with depth, which suggests that it is generated by structures on large scales. In particular, a Great Attractor as originally proposed by Lynden-Bell et al. cannot be responsible for the SMAC bulk flow. The SMAC data suggest infall into an attractor at the location of the Shapley Concentration, but the detection is marginal (at the 90% confidence level). We find that distant attractors in addition to the Shapley Concentration are required to explain the SMAC bulk flow. A comparison with peculiar velocities predicted from the IRAS PSCz redshift survey shows good agreement with a best fit value of Beta_I = Omega^0.6/b_I = 0.39 +- 0.17. However, the PSCz density field is not sufficient to acount for all of the SMAC bulk motion. We also detect, at the 98% confidence level, a residual bulk flow of 372 +- 127 km/s toward l = 273, b = 6 which must be generated by sources not included in the PSCz catalogue, that is, either beyond 200 h^-1 Mpc, in the Zone of Avoidance or in superclusters undersampled by IRAS. Finally, we compare the SMAC bulk flow with other recent measurements. We suggest that a bulk flow of 225 km/s toward l = 300, b = 10 at depths greater than 60 h^-1 Mpc is consistent with all peculiar velocity surveys, when allowance is made for sparse sampling effects. (abridged)

Under the hypothesis of a Dark Matter composed by supersymmetric particles like neutralinos, we investigate the possibility that their annihilation in the haloes of nearby galaxies could produce detectable fluxes of $\gamma$-photons. Expected fluxes depend on several, poorly known quantities such as the density profiles of Dark Matter haloes, the existence and prominence of central density cusps and the presence of a population of sub-haloes. We find that, for all reasonable choices of Dark Matter halo models, the intensity of the $\gamma$-ray flux from some of the nearest extragalactic objects, like M31, is comparable or higher than the diffuse Galactic foreground. We show that next generation ground-based experiments could have the sensitivity to reveal such fluxes which could help us unveiling the nature of Dark Matter particles.

We have used the two IRAS redshift surveys, 1.2-Jy (Fisher et al 1995) and PSCz (Saunders et al 2000), to model the linear velocity fields within a redshift of 8,000 km/s and have compared them in redshift space. The two velocity fields only differ significantly in their monopole components. The monopole discrepancy cannot be solely ascribed to shot-noise errors and incomplete sky coverage. The mismatch seems to arise from incompleteness of the PSCz catalog at fluxes $\le$ 1.2 Jy. The 1.2-Jy and PSCz higher order velocity multipoles, particularly the dipole and quadrupole components, appear to be consistent, suggesting that the dipole residuals found by Davis, Nusser and Willick (1996) when comparing 1.2-Jy and MarkIII velocity fields probably originates from the MarkIII velocities calibration procedure rather than from uncertainties in the model velocity field. Our results illustrate the efficiency of the spherical harmonics decomposition techniques in detecting possible differences between real and model velocity fields. Similar analyses shall prove to be very useful in the future to test the reliability of next generation model velocity fields derived from new redshift catalogs like 2dFGRS (Colless et al 2001), SDSS (York et al 2000), 6dF and 2MRS.

Recent studies have suggested the possibility that the lightest supersymmetric particle is a suitable dark matter candidate. In this theoretical framework, annihilations in high density environments like the center of dark matter haloes may produce an intense flux of gamma-rays. In this paper we discuss the possibility of detecting the signatures of neutralino annihilation in nearby galaxies with next generation ground-based detectors.

We propose a possible way to detect baryons at low redshifts from the analysis of X-ray absorption spectra of bright AGN pairs. A simple semi-analytical model to simulate the spectra is presented. We model the diffuse warm-hot intergalactic medium (WHIM) component, responsible for the X-ray absorption, using inputs from high-resolution hydro-dynamical simulations and analytical prescriptions. We show that the number of OVII absorbers per unit redshift with column density larger than $10^{13.5}$ cm$^{-2}$ - corresponding to an equivalent width of $\sim$ 1 km/s - which will be possibly detectable by {\it XEUS}, is $\magcir 30$ per unit redshift. {\it Constellation-X} will detect $\sim 6$ OVII absorptions per unit redshift with an equivalent width of 10 km/s. Our results show that, in a $\Lambda$CDM Universe, the characteristic size of these absorbers at $z\sim 0.1$ is $\sim 1$ $h^{-1}$ Mpc. The filamentary structure of WHIM can be probed by finding coincident absorption lines in the spectra of background AGN pairs. We estimate that at least 20 AGN pairs at separation $\mincir 20$ arcmin are needed to detect this filamentary structure at a 3$\sigma$ level. Assuming observations of distant sources using {\it XEUS} for exposure times of 500 ksec, we find that the minimum source flux to probe the filamentary structure is $\sim 2\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$, in the 0.1-2.4 keV energy band. Thus, most pairs of these extragalactic X-ray bright sources have already been identified in the {\it ROSAT} All-Sky Survey. Re-observation of these objects by future missions could be a powerful way to search for baryons in the low redshift Universe.

We apply a new algorithm, called the Unbiased Minimal Variance (hereafter UMV) estimator, to reconstruct the cosmic density and peculiar velocity fields in our local universe from the SEcat catalog of peculiar velocities comprising both early (ENEAR) and late type (SFI) galaxies. The reconstructed fields are compared with those predicted from the IRAS PSCz galaxy redshift survey to constrain the value of beta = Omega_m^{0.6}/b, where Omega_m and b are the mass density and the bias parameters. The comparison of the density and velocity fields is carried out within the same methodological framework, and leads, for the first time, to consistent values of beta, yielding beta = 0.57_{-0.13}^{+0.11} and beta = 0.51 +/- {0.06}, respectively. We find that the distribution of the density and velocity residuals, relative to their respective errors, is consistent with a Gaussian distribution with sigma approximately 1, indicating that the density field predicted from the PSCz is an acceptable fit to that deduced from the peculiar velocities of the SEcat galaxies.

We present extensive tests of the Fast Action Method (FAM) for recovering the past orbits of mass tracers in an expanding universe from their redshift-space coordinates at the present epoch. The tests focus on the reconstruction of present-day peculiar velocities using mock catalogs extracted from high resolution $N$-body simulations. The method allows for a self-consistent treatment of redshift-space distortions by direct minimization of a modified action for a cosmological gravitating system. When applied to ideal, volume limited catalogs, FAM recovers unbiased peculiar velocities with a 1-D, 1\sigma error of ~220 km/s, if velocities are smoothed on a scale of 5 Mpc/h. Alternatively, when no smoothing is applied, FAM predicts nearly unbiased velocities for objects residing outside the highest density regions. In this second case the 1\sigma$error decreases to a level of ~150 km/s. The correlation properties of the peculiar velocity fields are also correctly recovered on scales larger than 5 Mpc/h. Similar results are obtained when FAM is applied to flux limited catalogs mimicking the IRAS PSCz survey. In this case FAM reconstructs peculiar velocities with similar intrinsic random errors, while velocity-velocity correlation properties are well reproduced beyond scales of ~8 Mpc/h. We also show that FAM provides better velocity predictions than other, competing methods based on linear theory or Zel'dovich approximation. These results indicate that FAM can be successfully applied to presently available galaxy redshift surveys such as IRAS PSCz.

We present the results of three independent analyses in which we show that IRAS galaxies faithfully trace the underlying distribution the mass in the local universe. In the first analysis we compare the mass and the galaxy density fields and show that they are consistent within 60 $h^{-1}$Mpc. In the second one we show that the tidal velocity field of the Mark III catalog is consistent with the tidal velocity field predicted from the distribution of IRAS galaxies, hence indicating that IRAS galaxies trace the mass distribution well beyond 60 $h^{-1}$Mpc. Finally, a third analysis aimed at determining the mean biasing relation of IRAS galaxies showed no appreciable deviations from the linear biasing prescription.

We assess the extent to which the flux-limited PSCz redshift sample is capable of accounting for the major share of mass concentrations inducing the external tidal forces affecting the peculiar velocities within the Local Supercluster (LS). The investigation is based upon a comparison of the ``true'' velocities in 2 large N-body simulations and their reconstruction from ``observation-mimicking'' mock catalogues. The mildly nonlinear ``mock'' LS and PSCz velocities are analyzed by means of the Least Action Principle technique in its highly optimized implementation of Nusser & Branchini's Fast Action Method (FAM). For both model N-body Universes, we conclude that the dipolar and quadrupolar force field implied by the PSCz galaxy distribution would indeed be sufficiently representing the full external tidal force field.

We compare the peculiar velocities derived from the I-band Tully-Fisher (TF) relation for 989 field spiral galaxies in the SFI catalogue with the predicted velocity field derived from the IRAS PSCz galaxy redshift survey. We assume linear gravitational instability theory and apply the maximum likelihood technique, VELMOD, to SFI galaxies within cz<6000 km/s. The resulting calibration of the TF relation is consistent with a previous, independent calibration for a similar sample of spirals residing in clusters. Our analysis provides an accurate estimate of the quantity beta. Using the forward TF relation and smoothing the predicted velocity field with a Gaussian filter of radius 300 km/s, we obtain beta = 0.42 +/- 0.04 (1 sigma uncertainty). This value, as well as other parameters in the fit, are robust to varying the smoothing radius to 500 km/s and splitting the sample into spherical shells in redshift space. The one exception is the small-scale velocity dispersion, sigma_v, which varies from 200 km/s within cz= 4000 km/s 500 km/s at larger distance. For beta = 0.42, the residuals between the TF data and the PSCz gravity field are uncorrelated, indicating that the model provides a good fit to the data. More generally, a chi^2 statistic indicates that the PSCz model velocity field provides an acceptable (3 sigma) fit to the data for 0.3 < beta < 0.5.

We measure the autocorrelation function of galaxies in the IRAS Point Source Catalogue galaxy redshift (PSCz) survey and investigate its dependence on the far-infrared colour and absolute luminosity of the galaxies. We find that the PSCz survey correlation function can be modelled out to a scale of 10 h^{-1}Mpc as a power law of slope 1.30 +/- 0.04 and correlation length 4.77 +/- 0.20. At a scale of 75 h^{-1}Mpc we find the value of J_3 to be 1500 +/- 400. We also find that galaxies with cooler dust temperatures, are more strongly clustered than warmer galaxies. Splitting the survey into three colour subsamples, we find that, between 1 and 10 h^{-1}Mpc, the ratio of the correlation function is factor of 1.5 higher for the cooler galaxies compared to the hotter galaxies. This is consistent with the suggestion that hotter galaxies have higher star-formation rates, and correspond to later-type galaxies which are less clustered than earlier types. Using volume limited sub-samples, we find a weak variation of the correlation function as a function of absolute luminosity, in the sense that more luminous galaxies are less clustered than fainter galaxies. The trend is consistent with the colour dependence of the correlation function and the observed colour-luminosity correlation, but the large uncertainties mean that it has a low statistical significance.

We use the Least Action Principle to predict the peculiar velocities of PSC-z galaxies inside $cz=2000$\km. Linear theory is used to account for tidal effects to $cz=15000$\km, and we iterate galaxy positions to account for redshift distortions. As the Least Action Principle is valid beyond Linear theory, we can predict reliable peculiar velocities even for very nearby galaxies (ie $cz \leq 500$\km). These predicted peculiar velocities are then compared with the observed velocities of 12 galaxies with Cepheid distances. The combination of the PSC-z galaxy survey (with its large sky coverage and uniform selection), with the accurate Cepheid distances, makes this comparison relatively free from systematic effects. We find that galaxies are good tracers of the mass, even at small ($\leq 10$\mpc) scales; and under the assumption of no biasing, $0.25 \leq \beta \leq 0.75$ (at 90% confidence). We use the reliable predicted peculiar velocities to estimate the Hubble constant \ho from the local volume without ``stepping up'' the distance ladder, finding a confidence range of 65-75\kmM (at 90% confidence).

We investigate the spatial clustering of galaxies in the PSCz galaxy redshift survey, as revealed by the two-point correlation function, the luminosity mark correlations, and the moments of counts-in-cells. We construct volume-limited subsamples at different depths, and search for a luminosity dependence of the clustering pattern. We find no statistically significant effect in either the two-point correlation function or the mark correlations and so we take each subsample (of different characteristic luminosity) as representing the same statistical process. We then carry out a counts-in-cells analysis of the volume-limited subsamples, including a rigorous error calculation based on the recent theory of Szapudi, Colombi and Bernardeau. In this way, we derive the best estimates to date of the skewness and kurtosis of IRAS galaxies in redshift space. Our results agree well with previous measurements in both the parent angular catalogue, and in the derived redshift surveys. This is in contrast with smaller, optically selected surveys, were there is a discrepancy between the redshift space and projected measurements. Predictions from cold dark matter theory, obtained using the recent semi-analytical model of galaxy formation of Benson {\it et al}, provide an excellent description of our clustering data.

The algorithm ZTRACE of Monaco & Efstathiou (1999) is applied to the IRAS PSCz catalogue to reconstruct the initial conditions of our local Universe with a resolution down to ~5 Mpc/h. The 1-point PDF of the reconstructed initial conditions is consistent with the assumptions that (i) IRAS galaxies trace mass on scales of ~5 Mpc/h, and (ii) the statistics of primordial density fluctuations is Gaussian. We use simulated PSCz catalogues, constructed from N-body simulations with Gaussian initial conditions, to show that local non-linear bias can cause the recovered initial PDF (assuming no bias) to be non-Gaussian. However, for plausible bias models, the distortions of the recovered PDF would be difficult to detect using the volume finely sampled by the PSCz catalogue. So, for Gaussian initial conditions, a range of bias models remain compatible with our PSCz reconstruction results.

We present a comparison between the peculiar velocity field measured from the ENEAR all-sky $D_n-\sigma$ catalog and that derived from the galaxy distribution of the IRAS PSCz redshift survey. The analysis is based on a modal expansion of these data in redshift space by means of spherical harmonics and Bessel functions. The effective smoothing scale of the expansion is almost linear with redshift reaching 1500km/s at 3000km/s. The general flow patterns in the filtered ENEAR and PSCz velocity fields agree well within 6000km/s, assuming a linear biasing relation between the mass and the PSCz galaxies. The comparison allows us to determine the parameter $\beta=\Omega^{0.6}/b$, where $\Omega$ is the cosmological density parameter and $b$ is the linear biasing factor. A likelihood analysis of the ENEAR and PSCz modes yields $\beta=0.5 +- 0.1$, in good agreement with values obtained from Tully-Fisher surveys.

We present a simple method for evaluating the nonlinear biasing function of galaxies from a redshift survey. The nonlinear biasing is characterized by the conditional mean of the galaxy density fluctuation given the underlying mass density fluctuation, or by the associated parameters of mean biasing and nonlinearity (following Dekel & Lahav 1999). Using the distribution of galaxies in cosmological simulations, at smoothing of a few Mpc, we find that the mean biasing can be recovered to a good accuracy from the cumulative distribution functions (CDFs) of galaxies and mass, despite the biasing scatter. Then, using a suite of simulations of different cosmological models, we demonstrate that the matter CDF is robust compared to the difference between it and the galaxy CDF, and can be approximated for our purpose by a cumulative log-normal distribution of 1+\delta with a single parameter \sigma. Finally, we show how the nonlinear biasing function can be obtained with adequate accuracy directly from the observed galaxy CDF in redshift space. Thus, the biasing function can be obtained from counts in cells once the rms mass fluctuation at the appropriate scale is assumed a priori. The relative biasing function between different galaxy types is measurable in a similar way. The main source of error is sparse sampling, which requires that the mean galaxy separation be smaller than the smoothing scale. Once applied to redshift surveys such as PSCz, 2dF, SDSS, or DEEP, the biasing function can provide valuable constraints on galaxy formation and structure evolution.

The density and velocity fields as extracted from the Abell/ACO clusters are compared to the corresponding fields recovered by the POTENT method from the Mark~III peculiar velocities of galaxies. In order to minimize non-linear effects and to deal with ill-sampled regions we smooth both fields using a Gaussian window with radii ranging between $12 - 20\hmpc$. The density and velocity fields within $70\hmpc$ exhibit similarities, qualitatively consistent with gravitational instability theory and a linear biasing relation between clusters and mass. The random and systematic errors are evaluated with the help of mock catalogs. Quantitative comparisons within a volume containing $\sim\!12$ independent samples yield $\betac\equiv\Omega^{0.6}/b_c=0.22\pm0.08$, where $b_c$ is the cluster biasing parameter at $15\hmpc$. If $b_c \sim 4.5$, as indicated by the cluster correlation function, our result is consistent with $\Omega \sim 1$.

We present the results of reconstruction analysis of the galaxy distribution in a spherical region of radius 50 Mpc/h centered on the Local Group, as mapped by the IRAS Point Source Catalog Redshift Survey (PSCz). We reconstruct it using 15 different models for structure formation in the universe, each consisting of a set of assumptions regarding the value of the density parameter Omega_m, and the amplitude and nature of the biasing between IRAS galaxies and mass. For every model, we also reconstruct ten mock PSCz catalogs derived from the outputs of numerical simulations that have the appropriate assumptions about Omega_m and bias. We use a variety of statistics to compare the accuracy of each reconstruction of the PSCz catalog to the accuracy expected based on the corresponding mock catalog reconstructions. We find that gravitational instability of Gaussian mass density fluctuations can account for the galaxy distribution in the PSCz catalog, at least for some plausible assumptions about the value of Omega_m and the bias. Unbiased models in which IRAS galaxies trace mass fail to reconstruct the PSCz catalog accurately, both for Omega_m=0.4 and for Omega_m=1. Low Omega_m models in which IRAS galaxies are antibiased with respect to the mass are the most successful in reconstructing the PSCz catalog. In particular, a model with Omega_m=0.4 and IRAS galaxies related to the mass according to the predictions of a semi-analytic galaxy formation model is very successful in reconstructing the PSCz galaxy distribution.

We present the results for the predicted density and peculiar velocity fields and the dipole from the PSCz survey of 15,000 IRAS galaxies over 84% of the sky. We find a significant component to the dipole arising between 6000 and 15,000 km/s, but no significant component from greater distances. The misalignment with the CMB is 20 degrees. The most remarkable feature of the PSCz model velocity field is a coherent large-scale flow along the baseline connecting Perseus-Pisces, the Local Supercluster, Great Attractor and the Shapley Concentration. We have measured the parameter beta using the amplitude of the dipole, bulk flow and point by point comparisons between the individual velocities of galaxies in the MarkIII and SFI datasets, and the large-scale clustering distortion in redshift space.All our results are consistent with beta = 0.6 +- 0.1.

We perform a detailed comparison of the IRAS PSCz and 1.2-Jy spherical harmonic coefficients of the density and velocity fields in redshift space. The monopole terms predicted from the two surveys show some differences. The mismatch between the velocity monopoles arises from faint sources and disappears when extracting a PSCz subsample of galaxies with fluxes larger than 1.2 Jy. The analysis of PSCz dipole moments confirms the same inconsistencies found by Davis, Nusser and Willick (1996) between the IRAS 1.2-Jy gravity field and MARK III peculiar velocities. We conclude that shot-noise, which is greatly reduced in our PSCz gravity field, cannot be responsible for the observed mismatch.

Given the present distribution of mass tracing objects in an expanding universe, we develop and test a fast method for recovering their past orbits using the least action principle. In this method, termed FAM for Fast Action Minimization, the orbits are expanded in a set of orthogonal time-base functions satisfying the appropriate boundary conditions at the initial and final times. The conjugate gradient method is applied to locate the extremum of the action in the space of the expansion coefficients of the orbits. The TREECODE gravity solver routine is used for computing the gravitational fields appearing in the action and its gradients. The time integration of the Lagrangian is done using Gaussian quadratures. FAM allows us to increase the number of galaxies used in previous numerical action principle implementations by more than one order of magnitude. For example, orbits for the $\sim 15,000$ IRAS PSC$z$ galaxies can be recovered in $\sim 12,000$ CPU seconds on a 400MHz DEC-Alpha machine. FAM can recover the present peculiar velocities of particles and the initial fluctuations field. It successfully recovers the flow field down to clusters scales, where deviations of the flow from the Zel'dovich solution are significant. We also show how to recover orbits from the present distribution of objects as in redshift space by direct minimization of a modified action, without iterating the solution between real and redshift spaces.

We present a version of the Fourier Bessel method first introduced by Fisher et al (1994) and Zaroubi et al (1994) with two extensions: (a) we amend the formalism to allow a generic galaxy weight which can be constant rather than the more conventional overweighting of galaxies at high distances, and (b) we correct for the masked zones by extrapolation of Fourier Bessel modes rather than by cloning from the galaxy distribution in neighbouring regions. We test the procedure extensively on N-body simulations and find that it gives generally unbiased results but that the reconstructed velocities tend to be overpredicted in high-density regions. Applying the formalism to the PSZz redshift catalog, we find that beta = 0.7 +/- 0.5 from a comparison of the reconstructed Local Group velocity to the CMB dipole. From an anisotropy test of the velocity field, we find that beta = 1 CDM models models normalized to the current cluster abundance can be excluded with 90% confidence. The density and velocity fields reconstructed agree with the fields found by Branchini et al (1998) in most points. We find a back-infall into the Great Attractor region (Hydra-Centaurus region) but tests suggest that this may be an artifact. We identify all the major clusters in our density field and confirm the existence of some previously identified possible ones.

We present a self consistent nonparametric model of the local cosmic velocity field based on the density distribution in the PSCz redshift survey of IRAS galaxies. The error analysis, carried out on mock PSCz catalogues constructed from N-body simulations, reveals uncertainties of ~70 km/sec. The denser sampling provided by the PSCz survey compared to previous IRAS galaxy surveys allows us to reconstruct the velocity field out to larger distances. The most striking feature of the model velocity field is a coherent large-scale streaming motion along a baseline connecting Perseus-Pisces, the Local Supercluster, the Great Attractor, and the Shapley Concentration. We find no evidence for back-infall onto the Great Attractor. Instead, material behind and around the Great Attractor in inferred to be streaming towards the Shapley Concentration The PSCz model velocities compare well with those predicted from the 1.2Jy redshift survey of IRAS galaxies and with those predicted from the distribution of Abell/ACO clusters, out to 140 Mpc/h. Comparison of the real-space density fields (or, alternatively, the peculiar velocity fields) inferred from the PSC$z$ and cluster catalogues gives a relative (linear) bias parameter between clusters and IRAS galaxies of 4.4 +/- 0.6. Finally, we compare the cumulative bulk flows predicted from the PSCz gravity field with those measured from the Mark III and SFI catalogues of peculiar velocities. A conservative estimate of beta=Omega**0.6/b, where b is the bias parameter for IRAS galaxies, gives beta =0.6 + 0.22 -0.15, in agreement with other recent determinations.

We compute the acceleration on the Local Group using 11206 IRAS galaxies from the recently completed all-sky PSCz redshift survey. Measuring the acceleration vector in redshift space generates systematic uncertainties due to the redshift space distortions in the density field. We therefore assign galaxies to their real space positions by adopting a non-parametric model for the velocity field that solely relies on the linear gravitational instability and linear biasing hypotheses. Remaining systematic contributions to the measured acceleration vector are corrected for by using PSCz mock catalogues from N-body experiments. The resulting acceleration vector points approx. 15 degrees away from the CMB dipole apex, with a remarkable alignment between small and large scale contributions. A considerable fraction of the measured acceleration is generated within 40 h-1 Mpc with a non-negligible contribution from scales between 90 and 140 h-1 Mpc after which the acceleration amplitude seems to have converged. The local group acceleration from PSCz appears to be consistent with the one determined from the IRAS 1.2 Jy galaxy catalogue once the different contributions from shot noise have been taken into account. The results are consistent with the gravitational instability hypothesis and do not indicate any strong deviations from the linear biasing relation on large scales. A maximum-likelihood analysis of the comulative PSCz dipole is performed within a radius of 150 h-1 Mpc in which we account for nonlinear effects, shot noise and finite sample size. We obtain \beta = 0.70 (+ 0.35)(-0.2) at 1 \sigma confidence level.

We present a self-consistent nonparametric model of the cosmic velocity field based on the spatial distribution of IRAS galaxies in the recently completed all-sky PSC$z$ redshift survey. The most streaking feature of the PSCz model velocity field is a coherent large-scale streaming motion along the Perseus Pisces, Local Supercluster, Great Attractor and Shapley Concentration baseline, with no evidence for a backinfall into the Great Attractor region. A likelihood analysis has been performed to measure $\beta$. We have obtained $\beta =0.6^{+0.22}_{-0.15}$ (1-$\sigma$), in agreement with other recent determinations.

We compare the predicted Abell/ACO cluster density and velocity fields, using linear theory, with that based on the MarkIII catalog and the POTENT machinery within 7000 km/sec. Both density & velocity fields are in good qualitative agreement indicating that apart from a linear biasing factor the optical Abell/ACO clusters do trace the underlying mass distribution quite well.

Using a dynamical 3-D reconstruction procedure we estimate the peculiar velocities of $R\ge0$ Abell/ACO galaxy clusters from their measured redshift within 25000 km/sec. The reconstruction algorithm relies on the linear gravitational instability hypothesis, assumes linear biasing and requires an input value of the cluster $\beta$-parameter ($\beta_c \equiv \Omega_{\circ}^{0.6}/b_c$), which we estimated in Branchini \& Plionis (1995) to be $\beta_c\simeq 0.21$. The resulting cluster velocity field is dominated by a large scale streaming motion along the Perseus Pisces--Great Attractor base-line directed towards the Shapley concentration, in qualitative agreement with the galaxy velocity field on smaller scales. Fitting the predicted cluster peculiar velocities to a dipole term, in the local group frame and within a distance of $\sim 18000$ km/sec, we recover extremely well both the local group velocity and direction, in disagreement with the Lauer \& Postman (1994) observation. However, we find a $\sim 6\%$ probability that their observed velocity field could be a realization of our corresponding one, if the latter is convolved with their large distance dependent errors. Our predicted cluster bulk velocity amplitude agrees well with that deduced by the POTENT and the da Costa et al. (1995) analyses of observed galaxy motions at $\sim 5000 - 6000$ km/sec; it decreases thereafter while at the Lauer \& Postman limiting depth ($\sim 15000$ km/sec) its amplitude is $\sim 150$ km/sec, in comfortable agreement with most cosmological models.

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

Starting from the observed redshift distribution of a volume limited Abell/ACO cluster sample we use a 2-step procedure to recover their distances and peculiar velocities (assuming linear theory and linear biasing). The resulting 3D dipole is found to be 23\% less than the corresponding z-space one which leads to a value of $\beta \approx 0.21$ consistent with an $\Omega_{\circ}=1$ Universe. The cluster velocity field is dominated by a large-scale streaming motion towards the Shapley concentration and the resulting bulk velocity agrees extremely well with that deduced by POTENT, on smaller scales. Our results are in strong disagreement with those of Lauer \& Postman.

Starting from the observed distribution of galaxy clusters in redshift space we use a two--step procedure to recover their distances and peculiar velocities. After statistically correcting for the unobserved cluster distribution in the zone of avoidance ($|b|\le 20^{\circ}$) and also for a smooth absorption at higher $|b|$'s, we use a dynamical iterative algorithm to recover the real--space cluster positions by minimizing the redshift space distortions. The whole procedure assumes that clusters trace the mass, that peculiar velocities are caused by gravity and that linear perturbation theory applies. The amplitude of the cluster dipole measured in the 3D space turns out to be $\sim 23\%$ less than that measured in redshift space. In both cases the dipole direction is aligned with the Cosmic Microwave Background dipole within $\sim 10^{\circ}$, taking into account the Virgocentric infall component of the Local Group motion. Observational errors, limitations in the reconstruction procedure and the intrinsic cosmological variance, which is the dominant source of uncertainty, render a stringent determination of the $\beta$ parameter whose central value turns out to be $\beta \approx 0.2$ while its total uncertainty is $\pm 0.1$. This implies that for a cluster-mass bias parameter of $\sim 5$, a flat Universe is not excluded, contrary to previous cluster-dipole $z$-space analysis.

The least action principle (LAP) is a dynamically rigorous method for deriving the history of galaxy orbits. In particular it is an Omega_0 test, predicting current epoch galaxy velocities as a function of position and of the cosmological background. It is most usefully applied to in--falling structures, such as the local group, where its application indicates that the preferred cosmological model is Omega_0 = 0.1 and h=0.75 (h is the Hubble parameter in units of 100 Km s^-1 Mpc^-1). The method assumes that all the mass acts as if it were distributed as the visible galaxies. We test the reliability of the LAP to Local Group-like systems extracted from Omega_0=1 N--body simulations. While the orbits of the galaxies are qualitatively well reconstructed, the LAP systematically underestimates the mass of the system. This failure is attributed to the presence of extended halos weakly clustered around visible galaxies which prevent a large fraction of the group mass from being detected by the LAP technique. We conclude that the LAP method cannot rule out an Omega_0=1 value on the Local Group scale. Better constraints on Omega_0 may be obtained by applying this technique to in--falling systems, such as clusters, containing objects with separations large compared to galaxy sizes.

We show how, based on considerations on the observed form of the galaxy 2-point spatial correlation function xi(r), a very simplified -- yet surprisingly effective -- model for the linear density fluctuations power spectrum can be constructed. We first relate the observed large-scale shape of xi(r) to a power-law form for the power spectrum, P(k)\propto k^{-2.2}. For a plausible value of the bias parameter b = 1/sigma_8 ~ 1.8, one has (delta_rho / rho)_{rms} ~ 1 r ~ 3.5/h Mpc, suggesting that the change of slope observed in xi(r) around this scale marks the transition between the linear and nonlinear gravitational regimes. Under this working hypothesis, we use a simple analytical form to fit the large-scale correlations constraints together with the COBE CMB anisotropy measurement, thus constructing a simple phenomenological model for the linear power spectrum. Despite its simplicity, the model fits remarkably well directly estimated power spectra from different optical galaxy samples, and when evolved through an N-body simulation it provides a good match to the observed galaxy correlations. One of the most interesting features of the model is the small-scale one-dimensional velocity dispersion produced: sigma_{1d} = 450 Km s^{-1} at 0.5/h Mpc and sigma_{1d} = 350 Km s^{-1} for separations larger than ~ 2/h Mpc.