Modified Gravity (MG) scenarios have been advocated to account for the dark energy phenomenon in the universe. These models predict departures from General Relativity on large cosmic scales that can be tested through a variety of probes such as observations of galaxy clusters among others. Here, we investigate the imprint of MG models on the internal mass distribution of cluster-like halos as probed by the dark matter halo sparsity. To this purpose we perform a comparative analysis of the properties of the halo sparsity using N-body simulation halo catalogs of a standard flat $\Lambda$CDM model and MG scenarios from the DUSTGRAIN-pathfinder simulation suite. We find that the onset of the screening mechanism leaves a distinct signature in the redshift evolution of the ensemble average halos sparsity. Measurements of the sparsity of galaxy clusters from currently available mass estimates are unable to test MG models due to the large uncertainties on the cluster masses. We show that this should be possible in the future provided large cluster samples with cluster masses determined to better than $30\%$ accuracy level.

The late integrated Sachs-Wolfe (ISW) effect correlates the Cosmic Microwave Background (CMB) temperature anisotropies with foreground cosmic large-scale structures. As the correlation depends crucially on the growth history in the era of dark energy, it is a key observational probe for constraining the cosmological model. Here we present a detailed study based on full-sky and deep light cones generated from very large volume numerical N-body simulations, which allow us to avoid the use of standard replica techniques, while capturing the entirety of the late ISW effect on the large scales. We post-process the light cones using an accurate ray-tracing method and construct full-sky maps of the ISW temperature anisotropy for three different dark energy models. We quantify in detail the extent to which the ISW effect can be used to discriminate between different dark energy scenarios when cross-correlated with the matter distribution or the CMB lensing potential. We also investigate the onset of non-linearities, the so-called Rees-Sciama effect which provides a complementary probe of the dark sector. We find the signal of the lensing-lensing and ISW-lensing correlation of the three dark energy models to be consistent with measurements from the Planck satellite. Future surveys of the large-scale structures may provide cross-correlation measurements that are sufficiently precise to distinguish the signal of these models. Our methodology is very general and can be applied to any dark energy or modified gravity scenario as long as the metric seen by photons can still be characterized by a Weyl potential.

The Baryon Acoustic Oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the Purely Geometric BAO methods. To date two different Purely Geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable. We define the linear point and sound horizon and we characterize and compare the two rulers' cosmological-parameter dependence. We find that they are both geometrical within the required accuracy, and they have the same parameter dependence for a wide range of parameter values. We estimate the rulers' best-fit values and errors given the cosmological constraints obtained by the Planck Satellite team from the CMB measurements. We do this for three different cosmological models encompassed by the Purely Geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by $0.5\sigma$ when we do not fix the spatial geometry to be flat in LCDM. This points toward a sensitivity of the rulers to different cosmological models when they are estimated from the CMB.

The dark matter halo sparsity provides a direct observational proxy of the halo mass profile, characterizing halos in terms of the ratio of masses within radii which enclose two different overdensities. Previous numerical simulation analyses have shown that at a given redshift the halo sparsity carries cosmological information encoded in the halo mass profile. Moreover, its ensemble averaged value can be inferred from prior knowledge of the halo mass function at the overdensities of interest. Here, we present a detailed study of the ensemble average properties of the halo sparsity. In particular, using halo catalogs from high-resolution N-body simulations, we show that its ensemble average value can be estimated from the ratio of the averages of the inverse halo masses as well as the ratio of the averages of the halo masses at the overdensity of interests. This can be relevant for galaxy clusters data analyses. As an example, we have estimated the average sparsity properties of galaxy clusters from the LoCuSS and HIFLUGCS datasets respectively. The results suggest that the expected consistency of the different average sparsity estimates can provide a test of the robustness of mass measurements in galaxy cluster samples.

We perform an analysis of the cosmic reionization in the standard cold dark matter (CDM) paradigm and in alternative dark matter scenarios. Building upon the work of Corasaniti et al. (2017), we predict the reionization history for CDM, for warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ADM) models which reproduce state-of-art measurements of the galaxy luminosity function at very high-redshifts $6\le z\le 10$. To this purpose we adopt a reionization model parametrized in terms of the limiting UV-magnitude of galaxies contributing to the reionization $M_{\rm lim}$ and the average effective escape fraction of UV photons reaching the intergalactic medium $\tilde{f}$. For each DM model we compute the redshift evolution of the Thomson scattering optical depth $\tau_e(z)$ and the comoving ionization fraction $Q_{\rm HII}(z)$. We find the different DM models to have similar reionization histories. Differences with respect to the CDM case increase at fainter limiting UV-magnitudes and are degenerate with the effect of varying the reionization model parameters. Using Planck's determination of the integrated optical depth in combination with measurements of the neutral hydrogen fraction at different redshifts, we infer constraints on $\tilde{f}$ and $M_{\rm lim}$. The results are largely independent of the assumed DM scenario, in particular for $M_{\rm lim}\gtrsim -13$ we obtain that the effective escape fraction lies in the range $0.07\lesssim \tilde{f}\lesssim 0.15$ at $2\sigma$.

Leveraging the Baryon Acoustic Oscillations (BAO) feature present in clustering 2-point statistics, we aim to measure cosmological distances independently of the underlying background cosmological model. However this inference is complicated by late-time non-linearities that introduce model and tracer dependencies in the clustering correlation function and power spectrum, which must be properly accounted for. With this in mind, we introduce the "Purely-Geometric-BAO," which provides a rigorous tool to measure cosmological distances without assuming a specific background cosmology. We focus on the 2-point clustering correlation function monopole, and show how to implement such an inference scheme employing two different methodologies: the Linear Point standard ruler (LP) and correlation-function model-fitting (CF-MF). For the first time we demonstrate how, by means of the CF-MF, we can measure very precisely the sound-horizon/isotropic-volume-distance ratio, $r_{d}/D_{V}(\bar{z})$, while correctly propagating all the uncertainties. Using synthetic data, we compare the outcomes of the two methodologies, and find that the LP provides up to $50\%$ more precise measurements than the CF-MF. Finally, we test a procedure widely employed in BAO analyses: fitting the 2-point function while fixing the cosmological and the non-linear-damping parameters at fiducial values. We find that this underestimates the distance errors by nearly a factor of $2$. We thus recommend that this practice be reconsidered, whether for parameter determination or model selection.

We test the distance--duality relation $\eta \equiv d_L / [ (1 + z)^2 d_A ] = 1$ between cosmological luminosity distance ($d_L$) from the JLA SNe Ia compilation (arXiv:1401.4064) and angular-diameter distance ($d_A$) based on Baryon Oscillation Spectroscopic Survey (BOSS; arXiv:1607.03155) and WiggleZ baryon acoustic oscillation measurements (arXiv:1105.2862, arXiv:1204.3674). The $d_L$ measurements are matched to $d_A$ redshift by a statistically consistent compression procedure. With Monte Carlo methods, nontrivial and correlated distributions of $\eta$ can be explored in a straightforward manner without resorting to a particular evolution template $\eta(z)$. Assuming independent constraints on cosmological parameters that are necessary to obtain $d_L$ and $d_A$ values, we find 9% constraints consistent with $\eta = 1$ from the analysis of SNIa + BOSS and an 18% bound results from SNIa + WiggleZ. These results are contrary to previous claims that $\eta < 1$ has been found close to or above the $1 \sigma$ level. We discuss the effect of different cosmological parameter inputs and the use of the apparent deviation from distance--duality as a proxy of systematic effects on cosmic distance measurements. The results suggest possible systematic overestimation of SNIa luminosity distances compared with $d_A$ data when a Planck {\Lambda}CDM cosmological parameter inference (arXiv:1502.01589) is used to enhance the precision. If interpreted as an extinction correction due to a gray dust component, the effect is broadly consistent with independent observational constraints.

We present a new cosmological probe for galaxy clusters, the halo sparsity. This characterises halos in terms of the ratio of halo masses measured at two different radii and carries cosmological information encoded in the halo mass profile. Building upon the work of Balmes et al. (2014) we test the properties of the sparsity using halo catalogs from a numerical N-body simulation of ($2.6$ Gpc/h)$^3$ volume with $4096^3$ particles. We show that at a given redshift the average sparsity can be predicted from prior knowledge of the halo mass function. This provides a quantitative framework to infer cosmological parameter constraints using measurements of the sparsity of galaxy clusters. We show this point by performing a likelihood analysis of synthetic datasets with no systematics, from which we recover the input fiducial cosmology. We also perform a preliminary analysis of potential systematic errors and provide an estimate of the impact of baryonic effects on sparsity measurements. We evaluate the sparsity for a sample of 104 clusters with hydrostatic masses from X-ray observations and derive constraints on the cosmic matter density $\Omega_m$ and the normalisation amplitude of density fluctuations at the $8$ Mpc h$^{-1}$ scale, $\sigma_8$. Assuming no systematics, we find $\Omega_m=0.42\pm 0.17$ and $\sigma_8=0.80\pm 0.31$ at $1\sigma$, corresponding to $S_8\equiv \sigma_8\sqrt{\Omega_m}=0.48\pm 0.11$. Future cluster surveys may provide opportunities for precise measurements of the sparsity. A sample of a few hundreds clusters with mass estimate errors at a few percent level can provide competitive cosmological parameter constraints complementary to those inferred from other cosmic probes.

Aims.We study the dependence of the halo abundance function (AF) on different environments in a whole-sky $\Lambda$CDM light-cone halo catalogue extending to z < 0.65, using a simple and well-defined halo isolation criterion. Methods. The isolation status of each individual dark matter halo is determined by the distance to its nearest neighbour, which defines the maximum spherical region devoid of halos above a threshold mass around it (although the true size of such region may be much larger since it is not necessarily spherical). A versatile double power-law Schechter function is used to fit the dark matter halo AF, and its derived parameters are studied as a function of halo isolation status. Results. (a) Our function fits the halo abundances for all halo isolation statuses extremely well, while the well-established theoretical mass functions, integrated over the volume of the light-cone, provide an adequate but poorer fit than our phenomenological model. (b) As expected, and in agreement with other studies based on snap-shot simulations, we find significant differences of the halo abundance function as a function of halo isolation, indicating different rates of halo formation. The slope of the power law and the characteristic mass of the Schechter-like fitting function decrease with isolation, a result consistent with the formation of less massive haloes in lower density regions. (c) We find an unexpected upturn of the characteristic mass of the most isolated haloes of our sample. This upturn originates and characterises only the higher redshift regime (z > 0.45), which probably implies a significant and recent evolution of the isolation status of the most isolated and most massive haloes.

We show how a characteristic length scale imprinted in the galaxy two-point correlation function, dubbed the "linear point", can serve as a comoving cosmological standard ruler. In contrast to the Baryon Acoustic Oscillation peak location, this scale is constant in redshift and is unaffected by non-linear effects to within $0.5$ percent precision. We measure the location of the linear point in the galaxy correlation function of the LOWZ and CMASS samples from the Twelfth Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS) collaboration. We combine our linear-point measurement with cosmic-microwave-background constraints from the Planck satellite to estimate the isotropic-volume distance $D_{V}(z)$, without relying on a model-template or reconstruction method. We find $D_V(0.32)=1264\pm 28$ Mpc and $D_V(0.57)=2056\pm 22$ Mpc respectively, consistent with the quoted values from the BOSS collaboration. This remarkable result suggests that all the distance information contained in the baryon acoustic oscillations can be conveniently compressed into the single length associated with the linear point.

Due to late time non-linearities, the location of the acoustic peak in the two-point galaxy correlation function is a redshift-dependent quantity, thus it cannot be simply employed as a cosmological standard ruler. This has motivated the recent proposal of a novel ruler, also located in the Baryon Acoustic Oscillation range of scales of the correlation function, dubbed the "linear point". Unlike the peak, it is insensitive at the $0.5\%$ level to many of the non-linear effects that distort the clustering correlation function and shift the peak. However, this is not enough to make the linear point a useful standard ruler. In addition, we require a model-independent method to estimate its value from real data, avoiding the need to deploy a poorly known non-linear model of the correlation function. In this manuscript, we precisely validate a procedure for model-independent estimation of the linear point. We also identify the optimal set-up to estimate the linear point from the correlation function using galaxy catalogs. The methodology developed here is of general validity, and can be applied to any galaxy correlation-function data. As a working example, we apply this procedure to the LOWZ and CMASS galaxy samples of the Twelfth Data Release (DR12) of the Baryon Oscillation Spectroscopic Survey (BOSS), for which the estimates of cosmic distances using the linear point have been presented in Anselmi et al. (2017) [1].

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

We study the imprint of non-standard dark energy (DE) and dark matter (DM) models on the 21cm intensity map power spectra from high-redshift neutral hydrogen (HI) gas. To this purpose we use halo catalogs from N-body simulations of dynamical DE models and DM scenarios which are statistically indistinguishable from the standard Cold Dark Matter model with Cosmological Constant (LCDM) using currently available cosmological observations. We limit our analysis to halo catalogs at redshift z = 1 and 2.3 which are common to all simulations. For each catalog we model the HI distribution by using a simple prescription to associate the HI gas mass to N-body halos. We find that the DE models leave a distinct signature on the HI spectra across a wide range of scales, which correlates with differences in the halo mass function and the onset of the non-linear regime of clustering. In the case of the non-standard DM model significant differences of the HI spectra with respect to the LCDM model only arise from the suppressed abundance of low mass halos. These cosmological model dependent features also appear in the 21cm spectra. In particular, we find that future SKA measurements can distinguish the imprints of DE and DM models at high statistical significance.

We use state-of-art measurements of the galaxy luminosity function (LF) at z=6, 7 and 8 to derive constraints on warm dark matter (WDM), late-forming dark matter (LFDM) and ultra-light axion dark matter (ULADM) models alternative to the cold dark matter (CDM) paradigm. To this purpose we have run a suite of high-resolution N-body simulations to accurately characterise the low mass-end of the halo mass function and derive DM model predictions of the high-z luminosity function. In order to convert halo masses into UV-magnitudes we introduce an empirical approach based on halo abundance matching which allows us to model the LF in terms of the amplitude and scatter of the ensemble average star formation rate halo mass relation of each DM model, $\langle {\rm SFR}({\rm M_{ h}},z)\rangle$. We find that independent of the DM scenario the average SFR at fixed halo mass increases from z=6 to 8, while the scatter remains constant. At halo mass ${\rm M_{h}}\gtrsim 10^{12}\,{\rm M}_\odot$ h$^{-1}$ the average SFR as function of halo mass follows a double power law trend that is common to all models, while differences occur at smaller masses. In particular, we find that models with a suppressed low-mass halo abundance exhibit higher SFR compared to the CDM results. Using deviance statistics we obtain a lower limit on the WDM thermal relic particle mass, $m_{\rm WDM}\gtrsim 1.5$ keV at $2\sigma$. In the case of LFDM models, the phase transition redshift parameter is bounded to $z_t\gtrsim 8\cdot 10^5$ at $2\sigma$. We find ULADM best-fit models with axion mass $m_a\gtrsim 1.6\cdot 10^{-22}$ eV to be well within $2\sigma$ of the deviance statistics. We remark that measurements at $z=6$ slightly favour a flattening of the LF at faint UV-magnitudes. This tends to prefer some of the non-CDM models in our simulation suite, although not at a statistically significant level to distinguish them from CDM.

Bayesian graphical models are an efficient tool for modelling complex data and derive self-consistent expressions of the posterior distribution of model parameters. We apply Bayesian graphs to perform statistical analyses of Type Ia supernova (SN Ia) luminosity distance measurements from the joint light-curve analysis (JLA) data set. In contrast to the $\chi^2$ approach used in previous studies, the Bayesian inference allows us to fully account for the standard-candle parameter dependence of the data covariance matrix. Comparing with $\chi^2$ analysis results, we find a systematic offset of the marginal model parameter bounds. We demonstrate that the bias is statistically significant in the case of the SN Ia standardization parameters with a maximal 6 $\sigma$ shift of the SN light-curve colour correction. In addition, we find that the evidence for a host galaxy correction is now only 2.4 $\sigma$. Systematic offsets on the cosmological parameters remain small, but may increase by combining constraints from complementary cosmological probes. The bias of the $\chi^2$ analysis is due to neglecting the parameter-dependent log-determinant of the data covariance, which gives more statistical weight to larger values of the standardization parameters. We find a similar effect on compressed distance modulus data. To this end, we implement a fully consistent compression method of the JLA data set that uses a Gaussian approximation of the posterior distribution for fast generation of compressed data. Overall, the results of our analysis emphasize the need for a fully consistent Bayesian statistical approach in the analysis of future large SN Ia data sets.

We discuss an extension of the Coyote emulator to predict non-linear matter power spectra of dark energy (DE) models with a scale factor dependent equation of state of the form w = w_0 + ( 1 - a )w_a . The extension is based on the mapping rule between non-linear spectra of DE models with constant equation of state and those with time varying one originally introduced in ref. [40]. Using a series of N-body simulations we show that the spectral equivalence is accurate to sub-percent level across the same range of modes and redshift covered by the Coyote suite. Thus, the extended emulator provides a very efficient and accurate tool to predict non-linear power spectra for DE models with w_0 - w_a parametrization. According to the same criteria we have developed a numerical code, and we have implemented in a dedicated module for the CAMB code, that can be used in combination with the Coyote Emulator in likelihood analyses of non-linear matter power spectrum measurements. All codes can be found at https://github.com/luciano-casarini/PKequal

The covariance matrix of the matter power spectrum is a key element of the statistical analysis of galaxy clustering data. Independent realisations of observational measurements can be used to sample the covariance, nevertheless statistical sampling errors will propagate into the cosmological parameter inference potentially limiting the capabilities of the upcoming generation of galaxy surveys. The impact of these errors as function of the number of independent realisations has been previously evaluated for Gaussian distributed data. However, non-linearities in the late time clustering of matter cause departures from Gaussian statistics. Here, we address the impact of non-Gaussian errors on the sample covariance and precision matrix errors using a large ensemble of numerical N-body simulations. In the range of modes where finite volume effects are negligible ($0.1\lesssim k\,[h\,{\rm Mpc^{-1}}]\lesssim 1.2$) we find deviations of the estimated variance of the sample covariance with respect to Gaussian predictions above $\sim 10\%$ level. These reduce to about $\sim 5\%$ in the case of the precision matrix. Finally, we perform a Fisher analysis to estimate the effect of covariance errors on the cosmological parameter constraints. In particular, assuming Euclid-like survey characteristics we find that a number of independent realisation larger than $\gtrsim 5000$ is necessary to reduce the contribution of sample covariance errors to the cosmological parameter uncertainties at sub-percent level. We also show that restricting the analysis to large scales $k\lesssim0.2\,h\,{\rm Mpc^{-1}}$ results in a considerable loss in constraining power, while using the linear covariance to include smaller scales leads to an underestimation of the errors on the cosmological parameters.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

We perform a study of the nonlinear clustering of matter in the late-forming dark matter (LFDM) scenario in which dark matter results from the transition of a nonminimally coupled scalar field from radiation to collisionless matter. A distinct feature of this model is the presence of a damped oscillatory cutoff in the linear matter power spectrum at small scales. We use a suite of high-resolution N-body simulations to study the imprints of LFDM on the nonlinear matter power spectrum, the halo mass and velocity functions and the halo density profiles. The model largely satisfies high-redshift matter power spectrum constraints from Lyman-$\alpha$ forest measurements, while it predicts suppressed abundance of low-mass halos ($\sim 10^{9}-10^{10}$ h$^{-1}$ M$_\odot$) at all redshifts compared to a vanilla $\Lambda$CDM model. The analysis of the LFDM halo velocity function shows a better agreement than the $\Lambda$CDM prediction with the observed abundance of low-velocity galaxies in the local volume. Halos with mass $M\gtrsim 10^{11}$ h$^{-1}$ M$_\odot$ show minor departures of the density profiles from $\Lambda$CDM expectations, while smaller-mass halos are less dense, consistent with the fact that they form later than their $\Lambda$CDM counterparts.

Artificial fragmentation of the matter density field causes the formation of spurious groups of particles in N-body simulations of non-standard Dark Matter (DM) models which are characterized by a small scale cut-off in the linear matter power spectrum. These spurious halos alter the prediction of the mass function in a range of masses where differences among DM models are most relevant to observational tests. Using a suite of high resolution simulations we show that the contamination of artificial groups of particles significantly affect the statistics of halo spin, shape and virial state parameters. We find that spurious halos have systematically larger spin values, are highly elliptical or prolate and significantly deviate from virial equilibrium. These characteristics allow us to detect the presence of spurious halos even in non-standard DM models for which the low-mass end of the mass function remains well behaved. We show that selecting halos near the virial equilibrium provides a simple and effective method to remove the bulk of spurious halos from numerical halo catalogs and consistently recover the halo mass function at low masses.

The excursion set approach provides a framework for predicting how the abundance of dark matter halos depends on the initial conditions. A key ingredient of this formalism comes from the physics of halo formation: the specification of a critical overdensity threshold (barrier) which protohalos must exceed if they are to form bound virialized halos at a later time. Another ingredient is statistical, as it requires the specification of the appropriate statistical ensemble over which to average when making predictions. The excursion set approach explicitly averages over all initial positions, thus implicitly assuming that the appropriate ensemble is that associated with randomly chosen positions in space, rather than special positions such as peaks of the initial density field. Since halos are known to collapse around special positions, it is not clear that the physical and statistical assumptions which underlie the excursion set approach are self-consistent. We argue that they are at least for low mass halos, and illustrate by comparing our excursion set predictions with numerical data from the DEUS simulations.

The upcoming generation of galaxy surveys will probe the distribution of matter in the universe with unprecedented accuracy. Measurements of the matter power spectrum at different scales and redshifts will provide stringent constraints on the cosmological parameters. However, on non-linear scales this will require an accurate evaluation of the covariance matrix. Here, we compute the covariance matrix of the 3D matter density power spectrum for the concordance $\Lambda$CDM cosmology from an ensemble of N-body simulations of the Dark Energy Universe Simulation - Parallel Universe Runs (DEUS-PUR). This consists of 12288 realisations of a $(656\,h^{-1}\,\textrm{Mpc})^3$ simulation box with $256^3$ particles. We combine this set with an auxiliary sample of 96 simulations of the same volume with $1024^3$ particles. We find N-body mass resolution effect to be an important source of systematic errors on the covariance at high redshift and small intermediate scales. We correct for this effect by introducing an empirical statistical method which provide an accurate determination of the covariance matrix over a wide range of scales including the Baryon Oscillations interval. Contrary to previous studies that used smaller N-body ensembles, we find the power spectrum distribution to significantly deviate from expectations of a Gaussian random density field at $k\gtrsim 0.25\,h\,\textrm{Mpc}^{-1}$ and $z<0.5$. This suggests that in the case of finite volume surveys an unbiased estimate of the ensemble averaged band power at these scales and redshifts may require a careful assessment of non-Gaussian errors more than previously considered.

We show how cross-correlating a high redshift external tracer field, such as the 21cm neutral hydrogen distribution and product maps involving Cosmic Microwave Background (CMB) temperature and polarisation fields, that probe mixed bispectrum involving these fields, can help to determine the reionization history of the Universe, beyond what can be achieved from cross-spectrum analysis. Taking clues from recent studies for the detection of primordial non-Gaussianity (Munshi & Heavens 2010), we develop a set of estimators that can study reionization using a power spectrum associated with the bispectrum (or skew-spectrum). We use the matched filtering inherent in this method to investigate different reionization histories. We check to what extent they can be used to rule out various models of reionization and study cross contamination from different sources such as the lensing of the CMB. The estimators can be fine-tuned to optimize study of a specific reionization history. We consider three different types of tracers in our study, namely: proto-galaxies; 21cm maps of neutral hydrogen; quasars. We also consider four alternative models of reionization. We find that the cumulative signal-to-noise (S/N) for detection at $\ell_{max}=2000$ can reach $\mathcal{O}(70)$ for cosmic variance limited all-sky experiments. Combining 100GHz, 143GHz and 217GHz channels of the Planck experiment, we find that the (S/N) lies in the range $\mathcal{O}(5)-\mathcal{O}(35)$. The (S/N) depends on the specific choice of a tracer field, and multiple tracers can be effectively used to map out the entire reionization history with reasonable S/N. Contamination from weak lensing is investigated and found to be negligible, and the effects of Thomson scattering from patchy reionizationare also considered.

Observations of colliding galaxy clusters with high relative velocity probe the tail of the halo pairwise velocity distribution with the potential of providing a powerful test of cosmology. As an example it has been argued that the discovery of the Bullet Cluster challenges standard $\Lambda$CDM model predictions. Halo catalogs from N-body simulations have been used to estimate the probability of Bullet-like clusters. However, due to simulation volume effects previous studies had to rely on a Gaussian extrapolation of the pairwise velocity distribution to high velocities. Here, we perform a detail analysis using the halo catalogs from the Dark Energy Universe Simulation Full Universe Runs (DEUS-FUR), which enables us to resolve the high-velocity tail of the distribution and study its dependence on the halo mass definition, redshift and cosmology. Building upon these results we estimate the probability of Bullet-like systems in the framework of Extreme Value Statistics. We show that the tail of extreme pairwise velocities significantly deviates from that of a Gaussian, moreover it carries an imprint of the underlying cosmology. We find the Bullet Cluster probability to be two orders of magnitude larger than previous estimates, thus easing the tension with the $\Lambda$CDM model. Finally, the comparison of the inferred probabilities for the different DEUS-FUR cosmologies suggests that observations of extreme interacting clusters can provide constraints on dark energy models complementary to standard cosmological tests.

We investigate the non-linear evolution of Baryon Acoustic Oscillations (BAO) in the low-redshift matter power spectrum from the DEUS-FUR $\Lambda$CDM model simulation. This is the first cosmological N-body simulation encompassing the full observable cosmic volume, thus allowing cosmic variance limited predictions at BAO scales. We control the effect of numerical systematic errors using a series of large volume high-resolution simulations. The combined analysis allows us to measure the matter power spectrum between $z=0$ and $1$ to 1% over the entire BAO range, $0.03<k [\textrm{h Mpc}^{-1}]<0.3$, in bins of size $\Delta k/k\lesssim$ 1%. We define the BAO with respect to a non-linearly evolved wiggle-free spectrum and determine the characteristics of the BAO without recurring to extrapolation from global fitting functions. We quantify the effects of non-linearities on the position and amplitude of the BAO extrema, and the coupling to the broadband slope of the power spectrum. We use these estimates to test non-linear predictions from semi-analytical models. Quite remarkably from the analysis of the redshift evolution of BAO we find that the second dip and third peak remains unaltered by non-linear effects. Furthermore, we find that the square of the damping factor and the shift of the position of BAO extrema scale to good approximation as the square of the growth factor, in agreement with expectations from perturbation theory. This confirms the idea that, besides cosmic distances, an accurate measurement of BAO at different redshifts can directly probe the growth of cosmic structures.

We study the imprint of Dark Energy on the density profile of Dark Matter halos using a set of high-resolution large volume cosmological N-body simulations from the Dark Energy Universe Simulation Series (DEUSS). We first focus on the analysis of the goodness-of-fit of the Navarro-Frenk-White (NFW) profile which we find to vary with halo mass and redshift. We also find that the fraction of halos ill-fitted by NFW varies with cosmology, thus indicating that the mass assembly of halos with perturbed density profiles carries a characteristic signature of Dark Energy. To access this information independently of any parametric profile, we introduce a new observable quantity: the halo sparsity $s_\Delta$. This is defined as the mass ratio $M_{200}/M_\Delta$, i.e. the ratio of mass inside a sphere of radius $r_{200}$ to that contained within a radius $r_\Delta$, enclosing 200 and $\Delta$ times the mean matter density respectively. We find the average sparsity to be nearly independent of the total halo mass, while its value can be inferred to better than a few percent from the ratio of the integrated halo mass functions at overdensities $\Delta$ and 200 respectively. This provides a consistency relation that can validate observational measurements of the halo sparsity. Most importantly, the sparsity significantly varies with the underlying Dark Energy model, thus providing an alternative cosmological probe.

Bayesian model selection methods provide a self-consistent probabilistic framework to test the validity of competing scenarios given a set of data. We present a case study application to strong gravitational lens parametric models. Our goal is to select a homogeneous lens subsample suitable for cosmological parameter inference. To this end we apply a Bayes factor analysis to a synthetic catalog of 500 lenses with power-law potential and external shear. For simplicity we focus on double-image lenses (the largest fraction of lens in the simulated sample) and select a subsample for which astrometry and time-delays provide strong evidence for a simple power-law model description. Through a likelihood analysis we recover the input value of the Hubble constant to within 3\sigma statistical uncertainty. We apply this methodology to a sample of double image lensed quasars. In the case of B1600+434, SBS 1520+530 and SDSS J1650+4251 the Bayes' factor analysis favors a simple power-law model description with high statistical significance. Assuming a flat \LambdaCDM cosmology, the combined likelihood data analysis of such systems gives the Hubble constant H_0=76+15-5 km/s/Mpc having marginalized over the lens model parameters, the cosmic matter density and consistently propagated the observational errors on the angular position of the images. The next generation of cosmic structure surveys will provide larger lens datasets and the method described here can be particularly useful to select homogeneous lens subsamples adapted to perform unbiased cosmological parameter inference

The high-mass end of the halo mass function is a sensitive probe of primordial non-Gaussianity (NG). In a recent study [9] we have computed the NG halo mass function in the context of the Excursion Set theory and shown that the primordial NG imprint is coupled to that induced by the non-linear collapse of dark matter halos. We also found an excellent agreement with N-body simulation results. Here, we perform a more accurate computation which accounts for the interval validity of the bispectrum expansion to next-to-leading order and extend the calculation to the case of a non-vanishing primordial trispectrum.

We have performed the first-ever numerical N- body simulation of the full observable universe (DEUS "Dark Energy Universe Simulation" FUR "Full Universe Run"). This has evolved 550 billion particles on an Adaptive Mesh Refinement grid with more than two trillion computing points along the entire evolutionary history of the universe and across 6 order of magnitudes length scales, from the size of the Milky Way to that of the whole observable universe. To date, this is the largest and most advanced cosmological simulation ever run. It provides unique information on the formation and evolution of the largest structure in the universe and an exceptional support to future observational programs dedicated to mapping the distribution of matter and galaxies in the universe. The simulation has run on 4752 (of 5040) thin nodes of BULL supercomputer CURIE, using more than 300 TB of memory for 10 million hours of computing time. About 50 PBytes of data were generated throughout the run. Using an advanced and innovative reduction workflow the amount of useful stored data has been reduced to 500 TBytes.

The mass distribution of dark matter halos is a sensitive probe of primordial non-Gaussianity (NG). We derive an analytical formula of the halo mass function by perturbatively computing excursion set path-integrals for a non-Gaussian density field with non-vanishing skewness, fnl. We assume a stochastic barrier model which captures the main features of the ellipsoidal collapse of halos. Contrary to previous results based on extensions of the Press-Schechter formalism to NG initial conditions, we find that the non-spherical collapse of halos directly alter the signature of primordial NG. This points toward a potential degeneracy between the effect of primordial non-Gaussianity and that of non-linear halo collapse. The inferred mass function is found to be in remarkable agreement with N-body simulations of NG local type. Deviations are well within numerical uncertainties for all values of -80 < fNL^loc < 300 in the range of validity of the perturbative calculation. Moreover, the comparison with simulation results suggests that for fNL > 150 or fNL < -50 the non-linear collapse of halos, as described by our barrier model, strongly deviates from that of Gaussian initial conditions. This is not surprising since the effect of non-linear gravitational processes may be altered by initially large NG. Hence, in the lack of prior theoretical knowledge, halo collapse model parameters should be included in statistical halo mass function data analysis which aim to constrain the signature of primordial NG

Over the past decade advancements in the understanding of several astrophysical phenomena have allowed us to infer a concordance cosmological model that successfully accounts for most of the observations of our universe. This has opened up the way to studies that aim to better determine the constants of the model and confront its predictions with those of competing scenarios. Here, we use strong gravitational lenses as cosmological probes. Strong lensing, as opposed to weak lensing, produces multiple images of a single source. Extracting cosmologically relevant information requires accurate modeling of the lens mass distribution, the latter being a galaxy or a cluster. To this purpose a variety of models are available, but it is hard to distinguish between them, as the choice is mostly guided by the quality of the fit to the data without accounting for the number of additional parameters introduced. However, this is a model selection problem rather than one of parameter fitting that we address in the Bayesian framework. Using simple test cases, we show that the assumption of more complicate lens models may not be justified given the level of accuracy of the data.

We explore the cosmological constraints expected from wide area XMM-type cluster surveys covering 50-200 deg2, under realistic observing conditions. We perform a Fisher matrix analysis based on cluster number counts in combination with estimates of the 2-point cluster correlation function. The effect of the survey design is implemented through an observationally well tested cluster selection function. Special attention is given to the modelling of the shot noise and sample variance, which we estimate by applying our selection function to numerically simulated surveys. We then infer the constraints on the equation of state of the dark energy considering various survey configurations. We quantitatively investigate the respective impact of the cluster mass measurements, of the correlation function and of the 1<z<2 cluster population. We show that, with some 20 Ms XMM observing time, it is possible to constrain the dark energy parameters at a level which is comparable to that expected from the next generation of cosmic probes. Such a survey has also the power to provide unique insights into the physics of high redshift clusters and AGN properties.

We compute the dark matter halo mass function using the excursion set formalism for a diffusive barrier with linearly drifting average which captures the main features of the ellipsoidal collapse model. We evaluate the non-Markovian corrections due to the sharp filtering of the linear density field in real space with a path-integral method. We find an unprecedented agreement with N-body simulation data with deviations within ~5% level over the range of masses probed by the simulations. This indicates that the Excursion Set in combination with a realistic modelling of the collapse threshold can provide a robust estimation of the halo mass function.

We use the Excursion Set formalism to compute the properties of the halo mass distribution for a stochastic barrier model which encapsulates the main features of the ellipsoidal collapse of dark matter halos. Non-markovian corrections due to the sharp filtering of the linear density field in real space are computed with the path-integral technique introduced by Maggiore & Riotto (2010). Here, we provide a detailed derivation of the results presented in Corasaniti & Achitouv (2011) and extend the mass function analysis to higher redshift. We also derive an analytical expression for the linear halo bias. We find the analytically derived mass function to be in remarkable agreement with N-body simulation data from Tinker et al. (2008) with differences smaller than ~5% over the range of mass probed by the simulations. The excursion set solution from Monte Carlo generated random walks shows the same level of agreement, thus confirming the validity of the path-integral approach for the barrier model considered here. Similarly the analysis of the linear halo bias shows deviations no greater than 20%. Overall these results indicate that the Excursion Set formalism in combination with a realistic modeling of the conditions of halo collapse can provide an accurate description of the halo mass distribution.

The universality of the halo mass function is investigated in the context of dark energy cosmologies. This widely used approximation assumes that the mass function can be expressed as a function of the matter density omega_m and the rms linear density fluctuation sigma only, with no explicit dependence on the properties of dark energy or redshift. In order to test this hypothesis we run a series of 15 high-resolution N-body simulations for different cosmological models. These consists of three LCDM cosmologies best fitting WMAP-1, 3 and 5 years data, and three toy-models characterized by a Ratra-Peebles quintessence potential with different slopes and amounts of dark energy density. These toy models have very different evolutionary histories at the background and linear level, but share the same sigma8 value. For each of these models we measure the mass function from catalogues of halos identified in the simulations using the Friend-of-Friend (FoF) algorithm. We find redshift dependent deviations from a universal behaviour, well above numerical uncertainties and of non-stochastic origin, which are correlated with the linear growth factor of the investigated cosmologies. Using the spherical collapse as guidance, we show that such deviations are caused by the cosmology dependence of the non-linear collapse and virialization process. For practical applications, we provide a fitting formula of the mass function accurate to 5 percents over the all range of investigated cosmologies. We also derive an empirical relation between the FoF linking parameter and the virial overdensity which can account for most of the deviations from an exact universal behavior. Overall these results suggest that the halo mass function contains unique cosmological information since it carries a fossil record of the past cosmic evolution.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

In this "Invisible Universe" proceedings, we introduce the Dark Energy Universe Simulation Series (DEUSS) which aim at investigating the imprints of realistic dark energy models on cosmic structure formation. It represents the largest dynamical dark energy simulation suite to date in term of spatial dynamics. We first present the 3 realistic dark energy models (calibrated on latest SNIa and CMB data): LambdaCDM, quintessence with Ratra-Peebles potential, and quintessence with Sugra potential. We then isolate various contributions for non-linear matter power spectra from a series of pre-DEUSS high-resolution simulations (130 million particles). Finally, we introduce DEUSS which consist in 9 Grand Challenge runs with 1 billion particles each thus probing scales from 4 Gpc down to 3 kpc at z=0. Our goal is to make these simulations available to the community through the "Dark Energy Universe Virtual Observatory" (DEUVO), and the "Dark Energy Universe Simulations" (DEUS) consortium.

We present N-body cosmological numerical simulations including a primordial non-Gaussianity in the density fluctuation field quantified by the non-linear parameter $f_{NL}$. We have used MPGRAFIC code to produce initial conditions and the Adaptive Mesh Refinement (AMR) code RAMSES to evolve the large scale structure formation. We estimated the higher order momenta of the initial distribution of density fluctuations, investigated the redshift evolution of the non-linear power spectrum and estimated the discrepancy introduced by the primordial non-Gaussianity in the non-linear power spectrum.

Non-minimally coupled scalar field models suffer of unstable growing modes at the linear perturbation level. The nature of these instabilities depends on the dynamical state of the scalar field. In particular in systems which admit adiabatic solutions, large scale instabilities are suppressed by the slow-roll dynamics of the field. Here we review these results and present a preliminary likelihood data analysis suggesting that along adiabatic solutions coupled models with coupling of order of gravitational strength can provide viable cosmological scenarios satisfying constraints from SN Ia, CMB and large scale structure data.

Dark energy as a quintessence component causes a typical modification of the background cosmic expansion, which in addition to its clustering properties, can leave a potentially distinctive signature on large scale structures. Many previous studies have investigated this topic, particularly in relation to the non-linear regime of structure formation. However, no careful pre-selection of viable quintessence models with high precision cosmological data was performed. Here we show that this has led to a misinterpretation (and underestimation) of the imprint of quintessence on the distribution of large scale structures. To this purpose we perform a likelihood analysis of the combined Supernova Ia UNION dataset and WMAP5-years data to identify realistic quintessence models. Differences from the vanilla LambdaCDM are especially manifest in the predicted amplitude and shape of the linear matter power spectrum, though these remain within the uncertainties of the SDSS data. We use these models as benchmark for studying the clustering properties of dark matter halos by performing a series of high resolution N-body simulations. We find that realistic quintessence models allow for relevant differences of the dark matter distribution with the respect to the LambdaCDM scenario well into the non-linear regime, with deviations up to 40% in the non-linear power spectrum. Such differences are shown to depend on the nature of DE, as well as the scale and epoch considered. At small scales (k~1-5 h Mpc^{-1}, depending on the redshift) the structure formation process is about 20% more efficient than in LambdaCDM. We show that these imprints are a specific record of the cosmic structure formation history in DE cosmologies and therefore cannot be accounted in standard fitting functions of the non-linear matter power spectrum.

Modifications of the initial-state of the inflaton field can induce a departure from Gaussianity and leave a testable imprint on the higher order correlations of the CMB and large scale structures in the Universe. We focus on the bispectrum statistics of the primordial curvature perturbation and its projection on the CMB. For a canonical single-field action the three-point correlator enhancement is localized, maximizing in the collinear limit, corresponding to enfolded or squashed triangles in comoving momentum space. We show that the available local and equilateral template are very insensitive to this localized enhancement and do not generate noteworthy constraints on initial-state modifications. On the other hand, when considering the addition of a dimension 8 higher order derivative term, we find a dominant rapidly oscillating contribution, which had previously been overlooked and whose significantly enhanced amplitude is independent of the triangle under consideration. Nevertheless, the oscillatory nature of (the sign of) the correlation function implies the signal is nearly orthogonal to currently available observational templates, strongly reducing the sensitivity to the enhancement. Constraints on departures from the standard Bunch-Davies vacuum state can be derived, but also depend on the next-to-leading terms. We emphasize that the construction and application of especially adapted templates could lead to CMB bispectrum constraints on modified initial states already competing with those derived from the power spectrum.

In order to quantitatively test the ability of averaged inhomogeneous cosmologies to correctly describe observations of the large scale properties of the Universe, we introduce a smoothed template metric corresponding to a constant spatial curvature model at any time, but with an evolving curvature parameter. This metric is used to compute quantities along an approximate effective lightcone of the averaged model of the Universe. As opposed to the standard Friedmann model, we parameterize this template metric by exact scaling properties of an averaged inhomogeneous cosmology, and we also motivate this form of the metric by results on a geometrical smoothing of inhomogeneous cosmological hypersurfaces. We test our hypothesis for the template metric against supernova data and the position of the CMB peaks, and infer the goodness--of--fit and parameter uncertainties. We find that averaged inhomogeneous models can reproduce the observations without requiring an additional Dark Energy component (though a volume acceleration is still needed), and that current data do not disfavour our main assumption on the effective lightcone structure. We also show that the experimental uncertainties on the angular diameter distance and the Hubble parameter from Baryon Acoustic Oscillations measurements -- forseen in future surveys like the proposed EUCLID satellite project -- are sufficiently small to distinguish between a FLRW template geometry and the template geometry with consistently evolving curvature. (Abridged)

We study the evolution of linear density perturbations in the context of interacting scalar field-dark matter cosmologies, where the presence of the coupling acts as a stabilization mechanism for the runaway behavior of the scalar self-interaction potential as in the case of the Chameleon model. We show that in the "adiabatic" background regime of the system the rise of unstable growing modes of the perturbations is suppressed by the slow-roll dynamics of the field. Furthermore the coupled system behaves as an inhomogeneous adiabatic fluid. In contrast instabilities may develop for large values of the coupling constant, or along non-adiabatic solutions, characterized by a period of high-frequency dumped oscillations of the scalar field. In the latter case the dynamical instabilities of the field fluctuations, which are typical of oscillatory scalar field regimes, are amplified and transmitted by the coupling to dark matter perturbations.

WMAP observations have accurately determined the position of the first two peaks and dips in the CMB temperature power spectrum. These encode information on the ratio of the distance to the last scattering surface to the sound horizon at decoupling. However pre-recombination processes can contaminate this distance information. In order to assess the amplitude of these effects we use the WMAP data and evaluate the relative differences of the CMB peaks and dips multipoles. We find that the position of the first peak is largely displaced with the respect to the expected position of the sound horizon scale at decoupling. In contrast the relative spacings of the higher extrema are statistically consistent with those expected from perfect harmonic oscillations. This provides evidence for a scale dependent phase shift of the CMB oscillations which is caused by gravitational driving forces affecting the propagation of sound waves before recombination. By accounting for these effects we have performed a MCMC likelihood analysis to constrain in combination with recent BAO data a constant dark energy equation w. For a flat universe we find at 95% upper limit w<-1.10, and including the HST prior w<-1.14, which are only marginally consistent with limits derived from the supernova SNLS sample. Larger limits are obtained for non-flat cosmologies. From the full CMB likelihood analysis we also estimate the values of the shift parameter R and the multipole l_a of the acoustic horizon at decoupling for several cosmologies to test their dependence on model assumptions. Although the analysis of the full CMB spectra should be always preferred, using the position of the CMB peaks and dips provide a simple and consistent method for combining CMB constraints with other datasets.

In astro-ph/0702542v2, Linder and Miquel seek to criticize the use of Bayesian model selection for data analysis and for survey forecasting and design. Their discussion is based on three serious misunderstandings of the conceptual underpinnings and application of model-level Bayesian inference, which invalidate all their main conclusions. Their paper includes numerous further inaccuracies, including an erroneous calculation of the Bayesian Information Criterion. Here we seek to set the record straight.

Standard two-parameter compressions of the infinite dimensional dark energy model space show crippling limitations even with current SN-Ia data. Firstly they cannot cope with rapid evolution - the best-fit to the latest SN-Ia data shows late and very rapid evolution to w_0 = -2.85. However all of the standard parametrisations (incorrectly) claim that this best-fit is ruled out at more than 2-sigma, primarily because they track it well only at very low redshifts, z < 0.2. Further they incorrectly rule out the observationally acceptable region w << -1 for z > 1. Secondly the parametrisations give wildly different estimates for the redshift of acceleration, which vary from z_{acc}=0.14 to z_{acc}=0.59. Although these failings are largely cured by including higher-order terms (3 or 4 parameters) this results in new degeneracies which open up large regions of previously ruled-out parameter space. Finally we test the parametrisations against a suite of theoretical quintessence models. The widely used linear expansion in z is generally the worst, with errors of up to 10% at z=1 and 20% at z > 2. All of this casts serious doubt on the usefulness of the standard two-parameter compressions in the coming era of high-precision dark energy cosmology and emphasises the need for decorrelated compressions with at least three parameters.

We study the prospects for constraining dark energy at very high redshift with the Sandage-Loeb (SL) test - a measurement of the evolution of cosmic redshift obtained by taking quasar spectra at sufficiently separated epochs. This test is unique in its coverage of the ``redshift desert'', corresponding roughly to redshifts between 2 and 5, where other dark energy probes are unable to provide useful information about the cosmic expansion history. Extremely large telescopes planned for construction in the near future, with ultra high resolution spectrographs (such as the proposed CODEX), will indeed be able to measure cosmic redshift variations of quasar Lyman-alpha absorption lines over a period as short as ten years. We find that these measurements can constrain non-standard and dynamical dark energy models with high significance and in a redshift range not accessible with future dark energy surveys. As the cosmic signal increases linearly with time, measurements made over several decades by a generation of patient cosmologists may provide definitive constraints on the expansion history in the era that follows the dark ages but precedes the time when standard candles and rulers come into existence.

Cosmic dust extinction alters the flux of type Ia supernovae. Inhomogeneities in the dust distribution induce correlated fluctuations of the SN fluxes. We find that such correlation can be up to 60% of the signal caused by gravitational lensing magnification, with an opposite sign. Therefore if not corrected, cosmic dust extinction is the dominant source of systematic uncertainty for future SNe Ia lensing measurement limiting the overall S/N to be $\la 10$. On the other hand, SN flux correlation measurements can be used in combination with other lensing data to infer the level of dust extinction. This will provide a viable method to eliminate gray dust contamination from the SN Ia Hubble diagram.

Extinction by intergalactic gray dust introduces a magnitude redshift dependent offset in the standard-candle relation of SN Ia. This leads to overestimated luminosity distances compared to a dust-free universe. Quantifying the amplitude of this systematic effect is crucial for an accurate determination of the dark energy parameters. In this paper we model the gray dust extinction in terms of the star-formation history of the universe and the physical properties of the dust grains. We focus on a class of cosmic dust models which satisfy current observational constraints. These can produce an extinction as large as 0.08 mag at z=1.7 and potentially disrupt the dark energy parameter inference from future SN surveys. In particular depending on the dust model we find that an unaccounted extinction can bias the estimation of a constant dark energy equation of state w by shifting its best fit value up to 20% from its true value. Near-IR broadband photometry will hardly detect this effect, while the induced decrement of the Balmer lines requires high signal-to-noise spectra. Indeed IR-spectroscopy will be needed for high redshift SNe. Cosmic dust extinction may also cause a detectable violation of the distance-duality relation. A more comprehensive knowledge of the physics of the IGM is necessary for an accurate modeling of intergalactic dust. Due to the large magnitude dispersion current luminosity distance measurements are insensitive to such possible extinction effects. In contrast these must be taken into account if we hope to disclose the true nature of dark energy with the upcoming generation of SN Ia surveys.

A key science goal of upcoming dark energy surveys is to seek time evolution of the dark energy. This problem is one of {\em model selection}, where the aim is to differentiate between cosmological models with different numbers of parameters. However, the power of these surveys is traditionally assessed by estimating their ability to constrain parameters, which is a different statistical problem. In this paper we use Bayesian model selection techniques, specifically forecasting of the Bayes factors, to compare the abilities of different proposed surveys in discovering dark energy evolution. We consider six experiments -- supernova luminosity measurements by the Supernova Legacy Survey, SNAP, JEDI, and ALPACA, and baryon acoustic oscillation measurements by WFMOS and JEDI -- and use Bayes factor plots to compare their statistical constraining power. The concept of Bayes factor forecasting has much broader applicability than dark energy surveys.

The Advanced Liquid-Mirror Probe of Asteroids, Cosmology and Astrophysics (ALPACA) is a proposed 8-meter liquid mirror telescope surveying 1000 square degree of the southern-hemisphere sky. It will be a remarkably simple and inexpensive telescope, that nonetheless will deliver a powerful sample of optical data for studying dark energy. The bulk of the cosmological data consists of nightly, high signal-to-noise, multiband light curves of SN Ia. At the end of the three-years run ALPACA is expected to collect 100,000 SNe Ia up to z~1. This will allow us to reduce present systematic uncertainties affecting the standard-candle relation. The survey will also provide several other datasets such as the detection of baryon acoustic oscillations in the matter power spectrum and shear weak lensing measurements. In this preliminary analysis we forecast constraints on dark energy parameters from SN Ia and baryon acoustic oscillations. The combination of these two datasets will provide competitive constraints on the dark energy parameters under minimal prior assumptions. Further studies are needed to address the accuracy of weak lensing measurements.

We show that an interaction between dark matter and dark energy generically results in an effective dark energy equation of state of w<-1. This arises because the interaction alters the redshift-dependence of the matter density. An observer who fits the data treating the dark matter as non-interacting will infer an effective dark energy fluid with w<-1. We argue that the model is consistent with all current observations, the tightest constraint coming from estimates of the matter density at different redshifts. Comparing the luminosity and angular-diameter distance relations with LambdaCDM and phantom models, we find that the three models are degenerate within current uncertainties but likely distinguishable by the next generation of dark energy experiments.

We present an analysis of the constraining power of future measurements of the Integrated Sachs-Wolfe (ISW) effect on models of the equation of state of dark energy as a function of redshift, w(z). To achieve this, we employ a new parameterization of w, which utilizes the mean value of w(z) (<w>) as an explicit parameter. This helps to separate the information contained in the estimation of the distance to the last scattering surface (from the CMB) from the information contained in the ISW effect. We then use Fisher analysis to forecast the expected uncertainties in the measured parameters from future ISW observations for two models of dark energy with very different time evolution properties. For example, we demonstrate that the cross-correlation of Planck CMB data and LSST galaxy catalogs will provide competitive constraints on w(z), compared to a SNAP-like SNe project, for models of dark energy with a rapidly changing equation of state (e.g. 'Kink' models). Our work confirms that, while SNe measurements are more suitable for constraining variations in w(z) at low redshift, the ISW effect can provide important independent constraints on w(z) at high z.

While large scale cosmic microwave background (CMB) anisotropies involve a combination of the scalar and tensor fluctuations, the scalar amplitude can be independently determined through the CMB-galaxy cross-correlation. Using recently measured cross-correlation amplitudes, arising from the cross-correlation between galaxies and the Integrated Sachs Wolfe effect in CMB anisotropies, we obtain a constraint r < 0.5 at 68% confidence level on the tensor-to-scalar fluctuation amplitude ratio. The data also allow us to exclude gravity waves at a level of a few percent, relative to the density field, in a low - Lambda dominated universe(Omega_Lambda~0.5). In future, joining cross-correlation ISW measurements, which captures cosmological parameter information, with independent determinations of the matter density and CMB anisotropy power spectrum, may constrain the tensor-to-scalar ratio to a level above 0.05. This value is the ultimate limit on tensor-to-scalar ratio from temperature anisotropy maps when all other cosmological parameters except for the tensor amplitude are known and the combination with CMB-galaxy correlation allows this limit to be reached easily by accounting for degeneracies in certain cosmological parameters.

One of the fundamental and yet untested predictions of inflationary models is the generation of a very weak cosmic background of gravitational radiation. We investigate the sensitivity required for a space-based gravitational wave laser interferometer with peak sensitivity at $\sim 1$ Hz to observe such signal as a function of the model parameters and compare it with indirect limits that can be set with data from present and future cosmic microwave background missions. We concentrate on signals predicted by slow-roll single field inflationary models and instrumental configurations such as those proposed for the LISA follow-on mission: Big Bang Observer.

We investigate the possibility of constraining dark energy with the Integrated Sachs Wolfe effect recently detected by cross-correlating the WMAP maps with several Large Scale Structure surveys. In agreement with previous works, we found that, under the assumption of a flat universe, the ISW signal is a promising tool for constraining dark energy. Current available data put weak limits on a constant dark energy equation of state w. We also find no constraints on the dark energy sound speed c_e^2. For quintessence-like dark energy (c_e^2=1) we find w<-0.53, while tighter bounds are possible only if the dark energy is ``clustered'' (c_e^2=0), in such a case -1.94<w<-0.63 at 2-sigma. Better measurements of the CMB-LSS correlation will be possible with the next generation of deep redshift surveys. This will provide independent constraints on the dark energy which are alternative to those usually inferred from CMB and SN-Ia data.

Detecting dark energy dynamics is the main quest of current dark energy research. Addressing the issue demands a fully consistent analysis of CMB, large scale structure and SN-Ia data with multi-parameter freedom valid for all redshifts. Here we undertake a ten parameter analysis of general dark energy confronted with the first year WMAP, 2dF galaxy survey and latest SN-Ia data. Despite the huge freedom in dark energy dynamics there are no new degeneracies with standard cosmic parameters apart from a mild degeneracy between reionisation and the redshift of acceleration, both of which effectively suppress small scale power. Breaking this degeneracy will help significantly in detecting dynamics, if it exists. Our best-fit model to the data has significant late-time evolution at z<1.5. Phantom models are also considered and we find that the best-fit crosses w=-1 which, if confirmed, would be a clear signal for radically new physics. Treatment of such rapidly varying models requires careful integration of the dark energy density usually not implemented in standard codes, leading to crucial errors of up to 5%. Nevertheless cosmic variance means that standard LCDM models are still a very good fit to the data and evidence for dynamics is currently very weak. Independent tests of reionisation or the epoch of acceleration (e.g. ISW-LSS correlations) or reduction of cosmic variance at large scales (e.g. cluster polarisation at high redshift) may prove key in the hunt for dynamics.

By combining the recent WMAP measurements of the cosmic microwave background anisotropies and the results of the recent luminosity distance measurements to type-Ia supernovae, we find that the normalization of the matter power spectrum on cluster scales, sigma_8, can be used to discriminate between dynamical models of dark energy (quintessence models) and a conventional cosmological constant model (LCDM).

Dark energy affects the CMB through its perturbations and affects both CMB and Sn-Ia through its background evolution. Using recent CMB and Sn-Ia data sets, together with the most general parameterization of the dark energy equation of state available, we find that today w < -0.8 (2 sigma). We also find that the value of the normalization of the power spectrum on cluster scales, sigma8, can be used to discriminate between dynamical models of dark energy (Quintessence models) and a cosmological constant model (Lambda CDM).

We study the properties of cosmological density perturbations in a multi-component system consisting of a scalar field and a perfect fluid. We discuss the number of degrees of freedom completely describing the system, introduce a full set of dynamical gauge-invariant equations in terms of the curvature and entropy perturbations, and display an efficient formulation of these equations as a first-order system linked by a fairly sparse matrix. Our formalism includes spatial gradients, extending previous formulations restricted to the large-scale limit, and fully accounts for the evolution of an isocurvature mode intrinsic to the scalar field. We then address the issue of the adiabatic condition, in particular demonstrating its preservation on large scales. Finally, we apply our formalism to the quintessence scenario and clearly underline the importance of initial conditions when considering late-time perturbations. In particular, we show that entropy perturbations can still be present when the quintessence field energy density becomes non-negligible.

In this thesis we study some observational consequences of dark energy dominated cosmologies. After reviewing the main characteristics of quintessential models of dark energy, by using the SN Ia data and the position of the CMB peaks we constrain a class of parameterized scalar field potentials. Going beyond constraining specific classes models, we introduce a parameterization of the dark energy equation of state. Such an approach allows us to study in a model independent way the imprint dark energy leaves in the CMB power spectrum. Finally we develop a formalism to model localized CMB anisotropies and compute analytical formulae for the power spectrum and the bispectrum.

The consensus of opinion in cosmology is that the Universe is currently undergoing a period of accelerated expansion. With current and proposed high precision experiments it offers the hope of being able to discriminate between the two competing models that are being suggested to explain the observations, namely a cosmological constant or a time dependent `Quintessence' model. The latter suffers from a plethora of scalar field potentials all leading to similar late time behaviour of the universe, hence to a lack of predictability. In this paper, we develop a model independent approach which simply involves parameterizing the dark energy equation of state in terms of known observables. This allows to analyse the impact dark energy has had on cosmology without the need to refer to particular scalar field models and opens up the possibility that future experiments will be able to constrain the dark energy equation of state in a model independent manner.

We study the integrated Sachs-Wolfe (ISW) effect using a model-independent parameterization of the dark energy equation of state. Cosmic variance severely restricts the class of models distinguishable from Lambda-CDM. In particular if the present value of the equation of state satisfies w_Q^o <-0.8, the only distinguishable models are those with a rapidly varying equation of state, characterized by a late transition occurring at redshifts z<3. Due to the degeneracy with other cosmological parameters, models with a slowly varying behavior cannot be differentiated from each other or from a cosmological constant using even cosmic variance-limited CMB measurements. This may place a fundamental limit on our understanding of the origin of the presently observed acceleration.

Understanding the nature of the dark energy is one of the most important task in cosmology. In principle several cosmological observations can be used to discriminate amid a static cosmological constant contribution and a dynamical quintessence component. In view of the upcoming high resolution CMB experiments, using a model independent approach we study the dark energy imprint in the CMB anisotropy power spectrum and the degeneracy with other cosmological parameters.

Quintessence has been introduced as an alternative to the cosmological constant scenario to account for the current acceleration of the universe. This new dark energy component allows values of the equation of state parameter $w_{Q}^0\geq-1$, and in principle measurements of cosmological distances to Type Ia supernovae can be used to distinguish between these two types of models. Assuming a flat universe, we use the supernovae data and measurements of the position of the acoustic peaks in the Cosmic Microwave Background (CMB) spectra to constrain a rather general class of Quintessence potentials, including inverse power law models and recently proposed Supergravity inspired potentials. In particular we use a likelihood analysis, marginalizing over the dark energy density $\Omega_{Q}$, the physical baryon density $\Omega_{b}h^2$ and the scalar spectral index $n$, to constrain the slopes of our Quintessence potential. Considering only the first Doppler peak the best fit in our range of models gives $w_{Q}^0\sim-0.8$. However, including the SnIa data and the three peaks, we find an upper limit on the present value of the equation of state parameter, $-1\leq w_{Q}^0\leq-0.93$ at $2\sigma$, a result that appears to rule out a class of recently proposed potentials.

The analysis of the CMB anisotropies is a rich source of cosmological informations. In our study, we simulated the signals produced by the relics of a first order phase transition occured during an inflationary epoch in the early Universe. These relics are bubbles of true vacuum that leave a characteristic non-Gaussian imprint on the CMB. We use different statistical estimators in order to evaluate this non-Gaussianity. We obtain some limits on the allowed values of the bubble parameters comparing our results with the experimental data. We also predict the possibility to detect this signal with the next high resolution experiments.

The existence of large scale voids in several galaxy surveys suggests the occurence of an inflationary first order phase transition. This process generates primordial bubbles that, before evolving into the present voids, leave at decoupling a non-Gaussian imprint on the CMB. I this paper we evaluate an analytical expression of the collapsed three point correlation function from the bubble temperature fluctuations. Comparing the results with COBE-DMR measures, we obtain upper limits on the allowed non-Gaussianity and hence on the bubble parameters.