We propose a new method, called MonteCarlo Posterior Fit, to boost the MonteCarlo sampling of likelihood (posterior) functions. The idea is to approximate the posterior function by an analytical multidimensional non-Gaussian fit. The many free parameters of this fit can be obtained by a smaller sampling than is needed to derive the full numerical posterior. In the examples that we consider, based on supernovae and cosmic microwave background data, we find that one needs an order of magnitude smaller sampling than in the standard algorithms to achieve comparable precision. This method can be applied to a variety of situations and is expected to significantly improve the performance of the MonteCarlo routines in all the cases in which sampling is very time-consuming. Finally, it can also be applied to Fisher matrix forecasts, and can help solve various limitations of the standard approach.

We introduce a novel theory of gravity based on the inverse of the Ricci tensor, that we call the anticurvature tensor. We derive the general equations of motion for any Lagrangian function of the curvature and anticurvature scalars. We then demonstrate a no-go theorem: no Lagrangian that contains terms linear in any positive or negative power of the anticurvature scalar can drive an evolution from deceleration to acceleration, as required by observations. This effectively rules out many realizations of this theory, as we illustrate in detail in a particular case. Finally, we speculate on how to circumvent the no-go theorem.

In this work we provide updated constraints on coupled dark energy (CDE) cosmology with Peebles-Ratra (PR) potential and constant coupling strength $\beta$. This modified gravity scenario introduces a fifth force between dark matter particles, mediated by a scalar field that plays the role of dark energy. The mass of the dark matter particles does not remain constant, but changes with time as a function of the scalar field. Here we focus on the phenomenological behavior of the model, and assess its ability to describe updated cosmological data sets that include the Planck 2018 cosmic microwave background (CMB) temperature, polarization and lensing, baryon acoustic oscillations, the Pantheon compilation of supernovae of Type Ia, data on $H(z)$ from cosmic chronometers, and redshift-space distortions. We also study which is the impact of the local measurement of $H_0$ from SH0ES and the strong-lensing time delay data from the H0LICOW collaboration on the parameter that controls the strength of the interaction in the dark sector. We find a peak corresponding to a coupling $\beta > 0$ and to a potential parameter $\alpha > 0$, more or less evident depending on the data set combination. We show separately the impact of each data set and remark that it is especially CMB lensing the one data set that shifts the peak the most towards $\Lambda$CDM. When a model selection criterion based on the full Bayesian evidence is applied, however, $\Lambda$CDM is still preferred in all cases, due to the additional parameters introduced in the CDE model.

We argue that the $\Lambda$CDM tensions of the Hubble-Lemaitre expansion rate $H_0$ and the clustering normalization $\sigma_8$ can be eased, at least in principle, by considering an interaction between dark energy and dark matter in such a way to induce a small and positive early effective equation of state and a weaker gravity. For a dark energy scalar field $\phi$ interacting with dark matter through an exchange of both energy and momentum, we derive a general form of the Lagrangian allowing for the presence of scaling solutions. In a subclass of such interacting theories, we show the existence of a scaling $\phi$-matter-dominated-era ($\phi$MDE) which can potentially alleviate the $H_0$ tension by generating an effective high-redshift equation of state. We also study the evolution of perturbations for a model with $\phi$MDE followed by cosmic acceleration and find that the effective gravitational coupling relevant to the linear growth of large-scale structures can be smaller than the Newton gravitational constant $G$ at low redshifts. The momentum exchange between dark energy and dark matter plays a crucial role for realizing weak gravity, while the energy transfer is also required for the existence of $\phi$MDE.

The 21cm line probes the evolution of matter perturbations over a wide range of redshifts, from the dark ages down to the completion of reionization. Observing the 21cm cosmological signal will extend our understanding of the evolution of the Universe and it is thus important to investigate the predictions of different cosmological models. In this paper we focus on the prospect of constraining coupled quintessence models during the Epoch of Reionization both for global signal experiments and for intensity mapping surveys. To derive the all-sky 21cm signal and fluctuations in coupled quintessence, we simulate cosmological volumes of the 21cm signal including the coupling between dark matter and the quintessence field, where the strength of the coupling is labeled by the parameter $Q$. We show that the coupling between dark matter and quintessence modifies structure formation and expedites the process of reionization. For the upcoming 21cm line surveys like SKA and a fiducial global 21cm signal experiment, we perform a Fisher matrix analysis to constrain the coupling $Q$ and the dark matter density parameter $\Omega_\mathrm{dm}$. The results indicate that SKA will be able to place a 68% upper limit $0.04$ to $|Q|$. At the same time, such a global 21cm detector provides a constraint on the dark matter density parameter $\Omega_\mathrm{dm}$ around $\Delta\Omega_{\rm dm}\approx 0.005$, whereas SKA sets a quite weaker constraint around $\Delta\Omega_{\rm dm}\approx0.1$. These constraints are comparable to those already obtained from the cosmic microwave background, but explore an entirely different redshift range.

The galaxy power spectrum is one of the central quantities in cosmology. It contains information about the primordial inflationary process, the matter clustering, the baryon-photon interaction, the effects of gravity, the galaxy-matter bias, the cosmic expansion, the peculiar velocity field, etc.. Most of this information is however difficult to extract without assuming a specific cosmological model, for instance $\Lambda$CDM and standard gravity. In this paper we explore instead how much information can be obtained that is independent of the cosmological model, both at background and linear perturbation level. We determine the full set of model-independent statistics that can be constructed by combining two redshift bins and two distinct tracers. We focus in particular on the statistics $r(k,z_1,z_2)$, defined as the ratio of $f\sigma_8(z)$ at two redshift shells, and we show how to estimate it with a Fisher matrix approach. Finally, we forecast the constraints on $r$ that can be achieved by future galaxy surveys, and compare it with the standard single-tracer result. We find that $r$ can be measured with a precision from 3 to 11%, depending on the survey. Using two tracers, we find improvements in the constraints up to a factor of two.

Sub-horizon perturbations in scalar-tensor theories have been shown to grow generically faster than in uncoupled models due to a positive, additive Yukawa force. In such cases, the amount of clustering becomes larger than in the standard cosmological model, exacerbating the observed tension in the $\sigma_8$ parameter. Here we show instead that in some simple cases of conformally coupled dark energy one can obtain a transient regime of negative Yukawa force, without introducing ghosts or other instabilities.

We investigate the cosmological perturbations in the context of a Scalar-Tensor-Vector theory of Gravity known as MOG in the literature. Recent investigations show that MOG reproduces a viable background cosmological evolution comparable to $\Lambda$CDM. However, the matter-dominated era is slightly different. In this paper, we study the linear matter perturbations and estimate the relevant modified gravity parameters. We show that MOG reduces the growth rate of the perturbations and comparing with the RSD data reveals that MOG suggests a higher value for $\sigma_8$, compare to $\Lambda$CDM. This point, constitutes a powerful challenge to the cosmological viability of MOG.

We present the new method DALI (Derivative Approximation for LIkelihoods) for reconstructing and forecasting posteriors. DALI extends the Fisher Matrix formalism but allows for a much wider range of posterior shapes. While the Fisher Matrix formalism is limited to yield ellipsoidal confidence contours, our method can reproduce the often observed flexed, deformed or curved shapes of known posteriors. This gain in shape fidelity is obtained by expanding the posterior to higher order in derivatives with respect to parameters, such that non-Gaussianity in the parameter space is taken into account. The resulting expansion is positive definite and normalizable at every order. Here, we present the new technique, highlight its advantages and limitations, and show a representative application to a posterior of dark energy parameters from supernovae measurements.

Supernova Ia magnitude surveys measure the dimensionless luminosity distance $H_{0}D_{L}$. However, from the distances alone one cannot obtain quantities like $H(z)$ or the dark energy equation of state, unless further cosmological assumptions are imposed. Here we show that by measuring the power spectrum of density contrast and of peculiar velocities of supernovae one can estimate also $H(z)/H_{0}$ regardless of background or linearly perturbed cosmology and of galaxy-matter bias. This method, dubbed Clustering of Standard Candles (CSC) also yields the redshift distortion parameter $\beta(k,z)$ and the galaxy power spectrum. We forecast that LSST may be able to constrain $H(z)/H_{0}$ to 3--12\% in redshift bins of $\Delta z=0.1$ up $z=0.8$.

We argue that primordial dark matter halos could be generated during radiation domination by long range attractive forces stronger than gravity. In this paper we derive the conditions under which these structures could dominate the dark matter content of the Universe while passing microlensing constraints and cosmic microwave background energy injection bounds. The dark matter particles would be clumped in objects in the solar mass range with typical sizes of the order of the solar system. Consequences for direct dark matter searches are important.

In this work, we investigate the possibility that the galaxy rotation curves can be explained in the framework of modified gravity models that introduce a Yukawa term in the gravitational potential. We include dark matter and assume that the fifth-force couples differently to dark matter and to baryons. We aim at constraining the modified gravity parameters $\beta$ and $\lambda$, that is, the strength and the range of the Yukawa fifth force, respectively, using a set of 40 galaxy rotation curves data from the SPARC catalogue. We include baryonic gas, disk and bulge components, along with a NFW halo of dark matter. Each galaxy rotation curve is modeled with three free parameters, beside the two global Yukawa parameter. We find that the inclusion of the Yukawa term improves the $\chi^2$ from $680.75$ to $536.23$ for $655$ degrees of freedom. As global best-fit we obtain $\beta = 0.34\pm0.04$ and $\lambda = 5.61\pm0.91$kpc and a dark matter content on average 20\% smaller than without the Yukawa term. The Bayesian evidence in favor of a NFW profile plus Yukawa term is higher than 8$\sigma$ with respect to the standard gravity parametrization.

The notion of self-acceleration has been introduced as a convenient way to theoretically distinguish cosmological models in which acceleration is due to modified gravity from those in which it is due to the properties of matter or fields. In this paper we review the concept of self-acceleration as given, for example, by [1], and highlight two problems. First, that it applies only to universal couplings, and second, that it is too narrow, i.e. it excludes models in which the acceleration can be shown to be induced by a genuine modification of gravity, for instance coupled dark energy with a universal coupling, the Hu-Sawicki f(R) model or, in the context of inflation, the Starobinski model. We then propose two new, more general, concepts in its place: force-acceleration and field-acceleration, which are also applicable in presence of non universal cosmologies. We illustrate their concrete application with two examples, among the modified gravity classes which are still in agreement with current data, i.e. f(R) models and coupled dark energy.

We show how to obtain constraints on $\beta=f/b$, the ratio of the matter growth rate and the bias that quantifies the linear redshift-space distortions, that are independent of the cosmological model, using multiple tracers of large-scale structure. For a single tracer the uncertainties on $\beta$ are constrained by the uncertainties in the amplitude and shape of the power spectrum, which is limited by cosmic variance. However, for two or more tracers this limit does not apply, since taking the ratio of power spectra cosmic variance cancels out, and in the linear (Kaiser) approximation one measures directly the quantity $(1+ \beta_1 \mu^2)^2/(1+ \beta_2 \mu^2)^2$, where $\mu$ is the angle of a given mode with the line of sight. We provide analytic formulae for the Fisher matrix for one and two tracers, and quantify the signal-to-noise ratio needed to make effective use of the multiple-tracer technique. We also forecast the errors on $\beta$ for a survey like Euclid.

This paper is a pedagogical introduction to models of gravity and how to constrain them through cosmological observations. We focus on the Horndeski scalar-tensor theory and on the quantities that can be measured with a minimum of assumptions. Alternatives or extensions of General Relativity have been proposed ever since its early years. Because of Lovelock theorem, modifying gravity in four dimensions typically means adding new degrees of freedom. The simplest way is to include a scalar field coupled to the curvature tensor terms. The most general way of doing so without incurring in the Ostrogradski instability is the Horndeski Lagrangian and its extensions. Testing gravity means therefore, in its simplest term, testing the Horndeski Lagrangian. Since local gravity experiments can always be evaded by assuming some screening mechanism or that baryons are decoupled, or even that the effects of modified gravity are visible only at early times, we need to test gravity with cosmological observations in the late universe (large-scale structure) and in the early universe (cosmic microwave background). In this work we review the basic tools to test gravity at cosmological scales, focusing on model-independent measurements.

We derive new constraints on the Hubble parameter $H_0$ using the available data on $H(z)$ from cosmic chronometers (CCH), and the Hubble rate data points from the supernovae of Type Ia (SnIa) of the Pantheon compilation and the Hubble Space Telescope (HST) CANDELS and CLASH Multy-Cycle Treasury (MCT) programs. We employ two alternative techniques, Gaussian Processes (GPs) and the Weighted Polynomial Regression (WPR) method, to reconstruct the Hubble function, determine the derived values of $H_0$, and compare them with the local HST measurement provided by Riess {\it et al.} (2018), $H_0^{{\rm HST}}=(73.48\pm1.66)$ km/s/Mpc, and with the Planck+$\Lambda$CDM value, $H_0^{{\rm P18}}=(66.88\pm 0.92)$ km/s/Mpc. With GPs we obtain $H_0=(67.99\pm 1.94)$ km/s/Mpc, and with the WPR method $H_0=(68.90\pm 1.96)$ km/s/Mpc. Both are fully compatible at $<1\sigma$ c.l., and also with $H_0^{{\rm P18}}$. In contrast, they are in $\sim 2\sigma$ tension with $H_0^{{\rm HST}}$.

We study the cosmology of a specific class of nonlocal model of modified gravity, the so-called Deser-Woodard (DW) model, modifying the Einstein-Hilbert action by a term $\sim R f(\Box^{-1}R)$, where $f$ is a free function. Choosing $f$ so as to reproduce the $\Lambda{\rm CDM}$ cosmological background expansion history within the nonlocal model, we implement the model in a cosmological linear Einstein--Boltzmann solver and study the deviations to GR the model induces in the scalar and tensor perturbations. We observe that the DW nonlocal model describes a modified propagation for the gravitational waves, as well as a lower linear growth rate and a stronger lensing power as compared to $\Lambda{\rm CDM}$, up to several percents. Such prominent growth and lensing features lead to the inference of a significantly smaller value of $\sigma_8$ with respect to the one in $\Lambda{\rm CDM}$, given \textit{Planck} CMB+lensing data. The prediction for the linear growth rate $f \sigma_8$ within the DW model is therefore significantly smaller than the one in $\Lambda{\rm CDM}$ and the addition of growth rate data $f \sigma_8$ from Redshift-space distortion measurements to \textit{Planck} CMB+lensing, opens a (dominant) tension between Redshift-space distortion data and the reconstructed \textit{Planck} CMB lensing potential. However, model selection issues only result in "weak" evidences for $\Lambda{\rm CDM}$ against the DW model given the data. Such a fact shows that the datasets we consider are not constraining enough for distinguishing between the models. As we discuss, the addition of galaxy WL data or cosmological constraints from future galaxy clustering, weak lensing surveys, but also third generation gravitational wave interferometers, prove to be useful for discriminating modified gravity models such as the DW one from $\Lambda{\rm CDM}$, within the close future.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The gravitational slip parameter is an important discriminator between large classes of gravity theories at cosmological and astrophysical scales. In this work we use a combination of simulated information of galaxy cluster mass profiles, inferred by Strong+Weak lensing analyses and by the study of the dynamics of the cluster member galaxies, to reconstruct the gravitational slip parameter $\eta$ and predict the accuracy with which it can be constrained with current and future galaxy cluster surveys. Performing a full-likelihood statistical analysis, we show that galaxy cluster observations can constrain $\eta$ down to the percent level already with a few tens of clusters. We discuss the significance of possible systematics, and show that the cluster masses and numbers of galaxy members used to reconstruct the dynamics mass profile have a mild effect on the predicted constraints.

A well-known problem of the $\Lambda$CDM model is the tension between the relatively high level of clustering, as quantified by the parameter $\sigma_8$, found in cosmic microwave background experiments and the smaller one obtained from large-scale observations in the late Universe. In this paper we show that coupled quintessence, i.e. a single dark energy scalar field conformally coupled to dark matter through a constant coupling, can solve this problem if the background is taken to be identical to the $\Lambda$CDM one. We show that two competing effects arise. On one hand, the additional scalar force is attractive, and is therefore expected to increase the clustering. On the other, in order to obtain the same background as $\Lambda$CDM, coupled quintessence must have a smaller amount of dark matter near the present epoch. We show that the second effect is dominating today and leads to an overall slower growth. Comparing to redshift distortion data, we find that coupled quintessence with $\Lambda$CDM background solves the tension between early and late clustering. We find for the coupling $\beta$ and for $\sigma_8$ the best fit values $|\beta| = 0.079^{+ 0.059}_{- 0.067}$ and $\sigma_8 = 0.818^{+0.115}_{-0.088}$. These values also fit the lensing data from the KiDS-450 survey. We also estimate that the future missions SKA and Euclid will constrain $\beta$ with an error of $\pm\, 1.5\times10^{-3}$ and for $\sigma_8$ of $\pm \,1.8\times10^{-3}$ at $1\sigma$ level.

Line intensity mapping opens up a new and exciting window for probing cosmology and fundamental physics during the Epoch of Reionisation, extending to redshifts previously untested by galaxy surveys. The power spectra of these line fluctuations are a promising tool to test gravity over a large range of scales and redshifts. We simulate cosmological volumes of 21cm fluctuations in general parametrisations of modified gravity, in order to calculate the corresponding power spectra, where additional parameters are the initial condition of matter perturbations $\alpha$ and the scale-dependent modified gravity parameter $Y$ (also known as $G_\mathrm{eff}$) that measures deviations from GR in the Poisson equation. We show the impact of these model-independent modifications of gravity, to either delay or expedite reionisation. For the 21cm intensity mapping survey to be performed by the SKA mission, we forecast the ability of line intensity mapping to constrain the parameters $Y$ and $\alpha$ at redshifts $z=6-11$, where $Y$ is assumed constant during this epoch (but without requiring constancy at all times). In our most conservative scenario, the $Y$ parameter can be constrained at the tens of percent level, while for improved modelling of foregrounds as well as of the (mildly) non-linear regime, up to sub-percent level constraints are attainable. We show the impact of jointly estimating reionisation model parameters and corresponding parameter correlations, as well as of foreground removal. We note, that tomography is crucial to break degeneracies and for constraints not to degrade significantly when adding reionisation model parameters, with most constraining power coming from the redshift bins $z=7-10$ where the shape of the 21cm power spectrum is evolving fastest.

The recent discovery of a $\gamma$-ray counterpart to a gravitational wave event has put extremely stringent constraints on the speed of gravitational waves at the present epoch. In turn, these constraints place strong theoretical pressure on potential modifications of gravity, essentially allowing only a conformally-coupled scalar to be active in the present Universe. In this paper, we show that direct detection of gravitational waves from optically identified sources can also measure or constrain the strength of the conformal coupling in scalar--tensor models through the time variation of the Planck mass. As a first rough estimate, we find that the LISA satellite can measure the dimensionless time variation of the Planck mass (the so-called parameter $\alpha_M$) at redshift around 1.5 with an error of about 0.03 to 0.13, depending on the assumptions concerning future observations. Stronger constraints can be achieved once reliable distance indicators at $z>2$ are developed, or with GW detectors that extend the capabilities of LISA, like the proposed Big Bang Observer. We emphasize that, just like the constraints on the gravitational speed, the bound on $\alpha_M$ is independent of the cosmological model.

The goal of this short report is to summarise some key results based on our previous works on model independent tests of gravity at large scales in the Universe, their connection with the properties of gravitational waves, and the implications of the recent measurement of the speed of tensors for the phenomenology of general families of gravity models for dark energy.

In this work, we use recent data on the Hubble expansion rate $H(z)$, the quantity $f\sigma_8(z)$ from redshift space distortions and the statistic $E_g$ from clustering and lensing observables to constrain in a model-independent way the linear anisotropic stress parameter $\eta$. This estimate is free of assumptions about initial conditions, bias, the abundance of dark matter and the background expansion. We denote this observable estimator as $\eta_{{\rm obs}}$. If $\eta_{{\rm obs}}$ turns out to be different from unity, it would imply either a modification of gravity or a non-perfect fluid form of dark energy clustering at sub-horizon scales. Using three different methods to reconstruct the underlying model from data, we report the value of $\eta_{{\rm obs}}$ at three redshift values, $z=0.29, 0.58, 0.86$. Using the method of polynomial regression, we find $\eta_{{\rm obs}}=0.57\pm1.05$, $\eta_{{\rm obs}}=0.48\pm0.96$, and $\eta_{{\rm obs}}=-0.11\pm3.21$, respectively. Assuming a constant $\eta_{{\rm obs}}$ in this range, we find $\eta_{{\rm obs}}=0.49\pm0.69$. We consider this method as our fiducial result, for reasons clarified in the text. The other two methods give for a constant anisotropic stress $\eta_{{\rm obs}}=0.15\pm0.27$ (binning) and $\eta_{{\rm obs}}=0.53 \pm 0.19$ (Gaussian Process). We find that all three estimates are compatible with each other within their $1\sigma$ error bars. While the polynomial regression method is compatible with standard gravity, the other two methods are in tension with it.

The current paper is addressing the possibility of the Dark Energy scalar field potential reconstruction from the SNe Ia data and the problems arising during the process. We describe the method and test its limits, stability of the reconstruction with respect to the $\Omega_m$ and $H_0$ parameters, as well as several issues connected to the errors propagation, with use of synthetic data. After that, we test the method with real Union2.1 and JLA, as well as recent PANTHEON SNe Ia datasets. It worths mentioning that in our approach we assume no {\it Ans\"atzen} on the dynamical variables (e.g., $H(z)$), which makes our method free of any degeneracies and biases which arise if one assume one or another way to parameterize $H(z)$, or the equation of state, or some other variable. On the other hand, the price we pay for this freedom is immense -- although the scheme demonstrates perfect reconstruction in the case of synthetic data, the reconstruction from the real data seems almost impossible, at least on the current data precision level. The errors of the resulting potential and of the kinetic term are huge; for the Hubble parameter they are smaller but still severalfold larger than those obtained with parameterization. Finally, we test the method with non-SNe $H(z)$ datasets, and the results are also disappointing -- similar to the SNe case, the reconstruction is not possible as the reconstructed kinetic term enters the negative values. This could be treated as a manifestation of the insufficient data precision, or as the overestimation of $(H_0, \Omega_m)$ values -- we discuss both possibilities.

Intensity mapping enables us to test cosmology and fundamental physics during the Epoch of Reionisation, with power spectra of line fluctuations probing a large range of scales and redshifts. Cosmological volumes of 21cm line fluctuations in general modified gravity scenarios are presented, where additional parameters are the initial condition parameter $\alpha$ of matter perturbations and the scale-dependent modified gravity parameter $Y$ that measures deviations from General Relativity in the Poisson equation. For upcoming surveys like the SKA we forecast, using the power spectra derived from our simulations, the ability of intensity mapping to constrain modifications of gravity during reionisation, as well as investigate correlations with astrophysical parameters.

Primordial black holes can be produced by a long range attractive fifth force stronger than gravity, mediated by a light scalar field interacting with nonrelativistic "heavy" particles. As soon as the energy fraction of heavy particles reaches a threshold, the fluctuations rapidly become nonlinear. The overdensities collapse into black holes or similar screened objects, without the need for any particular feature in the spectrum of primordial density fluctuations generated during inflation. We discuss whether such primordial black holes can constitute the total dark matter component in the Universe.

In this paper we present new constraints on the Hubble parameter $H_0$ using: (i) the available data on $H(z)$ obtained from cosmic chronometers (CCH); (ii) the Hubble rate data points extracted from the supernovae of Type Ia (SnIa) of the Pantheon compilation and the Hubble Space Telescope (HST) CANDELS and CLASH Multy-Cycle Treasury (MCT) programs; and (iii) the local HST measurement of $H_0$ provided by Riess et al. (2018), $H_0^{{\rm HST}}=(73.45\pm1.66)$ km/s/Mpc. Various determinations of $H_0$ using the Gaussian processes (GPs) method and the most updated list of CCH data have been recently provided by Yu, Ratra and Wang (2018). Using the Gaussian kernel they find $H_0=(67.42\pm 4.75)$ km/s/Mpc. Here we extend their analysis to also include the most released and complete set of SnIa data, which allows us to reduce the uncertainty by a factor $\sim 3$ with respect to the result found by only considering the CCH information. We obtain $H_0=(67.06\pm 1.68)$ km/s/Mpc, which favors again the lower range of values for $H_0$ and is in tension with $H_0^{{\rm HST}}$. The tension reaches the $2.71\sigma$ level. We round off the GPs determination too by taking also into account the error propagation of the kernel hyperparameters when the CCH with and without $H_0^{{\rm HST}}$ are used in the analysis. In addition, we present a novel method to reconstruct functions from data, which consists in a weighted sum of polynomial regressions (WPR). We apply it from a cosmographic perspective to reconstruct $H(z)$ and estimate $H_0$ from CCH and SnIa measurements. The result obtained with this method, $H_0=(68.90\pm 1.96)$ km/s/Mpc, is fully compatible with the GPs ones. Finally, a more conservative GPs+WPR value is also provided, $H_0=(68.45\pm 2.00)$ km/s/Mpc, which is still almost $2\sigma$ away from $H_0^{{\rm HST}}$.

We propose a simple, nonlocal modification to general relativity (GR) on large scales, which provides a model of late-time cosmic acceleration in the absence of the cosmological constant and with the same number of free parameters as in standard cosmology. The model is motivated by adding to the gravity sector an extra spin-2 field interacting nonlocally with the physical metric coupled to matter. The form of the nonlocal interaction is inspired by the simplest form of the Deser-Woodard (DW) model, $\alpha R\frac{1}{\Box}R$, with one of the Ricci scalars being replaced by a constant $m^{2}$, and gravity is therefore modified in the infrared by adding a simple term of the form $m^2\frac{1}{\Box}R$ to the Einstein-Hilbert term. We study cosmic expansion histories, and demonstrate that the new model can provide background expansions consistent with observations if $m$ is of the order of the Hubble expansion rate today, in contrast to the simple DW model with no viable cosmology. The model is best fit by $w_0\sim-1.075$ and $w_a\sim0.045$. We also compare the cosmology of the model to that of Maggiore and Mancarella (MM), $m^2R\frac{1}{\Box^2}R$, and demonstrate that the viable cosmic histories follow the standard-model evolution more closely compared to the MM model. We further demonstrate that the proposed model possesses the same number of physical degrees of freedom as in GR. Finally, we discuss the appearance of ghosts in the local formulation of the model, and argue that they are unphysical and harmless to the theory, keeping the physical degrees of freedom healthy.

In this paper we perform a reconstruction of the scalar field potential responsible for cosmic acceleration using SNe Ia data. After describing the method, we test it with real SNe Ia data-Union2.1 and JLA SNe datasets. We demonstrate that with the current data precision level, the full reconstruction is not possible. We discuss the problems which arise during the reconstruction process and the ways to overcome them.

The recent report that the gravitational wave speed equals the light speed puts strong constraints on the anisotropic stress parameter of many modified gravity models, a quantity that is directly observable through large-scale structure. We show here that models without a mass scale completely escape these constraints. We discuss a few relevant cases in detail: Brans-Dicke theory, nonlocal models, and Galileon Lagrangian.

The coincident detection of gravitational waves (GW) and a gamma-ray burst from a merger of neutron stars has placed an extremely stringent bound on the speed of GW. We showed previously that the presence of gravitational slip ($\eta$) in cosmology is intimately tied to modifications of GW propagation. This new constraint implies that the only remaining viable source of gravitational slip is a conformal coupling to gravity in scalar-tensor theories, while viable vector-tensor theories cannot now generate gravitational slip at all. We discuss structure formation in the remaining viable models, demonstrating that (i) the dark-matter growth rate must now be at least as fast as in GR, with the possible exception of the beyond Horndeski model. (ii) If there is any scale-dependence at all in the slip parameter, it is such that it takes the GR value at large scales. We show a consistency relation which must be violated if gravity is modified.

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at $z=0.44$ and RXC J2248.7-4431 at $z=0.35$ to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou\'e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range $\lambda$ for models in which the dependence of this parameter on the enviroment is negligible at the scale considered (i.e. $\lambda=const$) and fixing the fifth force strength to the value predicted in $f(R)$ gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range $\lambda\le1.61\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$ on the marginalized distribution. For RXJ 2248 instead a possible tension with the $\Lambda$CDM model appears when adding lensing information, with a lower limit $\lambda\ge0.14\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained.[...]

What are the fundamental limitations of reconstructing the properties of dark energy, given cosmological observations in the quasilinear regime in a range of redshifts, to be as precise as required? The aim of this paper is to address this question by constructing model-independent observables, whilst completely ignoring practical problems of real-world observations. Non-Gaussianities already present in the initial conditions are not directly accessible from observations, because of a perfect degeneracy with the non-Gaussianities arising from the nonlinear matter evolution in generalized dark energy models. By imposing a specific set of evolution equations that should cover a range of dark energy cosmologies, we however find a constraint equation for the linear structure growth rate $f_1$ expressed in terms of model-independent observables. Entire classes of dark energy models which do not satisfy this constraint equation could be ruled out, and for models satisfying it we could reconstruct e.g. the nonlocal bias parameters $b_1$ and $b_2$.

The Fisher matrix is a widely used tool to forecast the performance of future experiments and approximate the likelihood of large data sets. Most of the forecasts for cosmological parameters in galaxy clustering studies rely on the Fisher matrix approach for large-scale experiments like DES, Euclid, or SKA. Here we improve upon the standard method by taking into account three effects: the finite window function, the correlation between redshift bins, and the uncertainty on the bin redshift. The first two effects are negligible only in the limit of infinite surveys. The third effect, on the contrary, is negligible for infinitely small bins. Here we show how to take into account these effects and what the impact on forecasts of a Euclid-type experiment will be. The main result of this article is that the windowing and the bin cross-correlation induce a considerable change in the forecasted errors, of the order of 10-30% for most cosmological parameters, while the redshift bin uncertainty can be neglected for bins smaller than $\Delta z = 0.1$ roughly.

In many generalized models of gravity, perfect fluids in cosmology give rise to gravitational slip. Simultaneously, in very broad classes of such models, the propagation of gravitational waves is altered. We investigate the extent to which there is a one-to-one relationship between these two properties in three classes of models with one extra degree of freedom: scalar (Horndeski and beyond), vector (Einstein-Aether) and tensor (bimetric). We prove that in bimetric gravity and Einstein-Aether, it is impossible to dynamically hide the gravitational slip on all scales whenever the propagation of gravitational waves is modified. Horndeski models are much more flexible, but it is nonetheless only possible to hide gravitational slip dynamically when the action for perturbations is tuned to evolve in time toward a divergent kinetic term. These results provide an explicit, theoretical argument for the interpretation of future observations if they disfavoured the presence of gravitational slip.

In this work, we extend previous analyses of the structure formation in the $f(\Box^{-1}R)$ model of nonlocal gravity proposed by Deser and Woodard (DW), which reproduces the background expansion of $\Lambda$CDM with no need of a cosmological constant nor of any dimensional constant beside Newton's one. A previous analysis based on redshift-space distortions (RSD) data concluded that the model was ruled out. In this work we revisit the issue and find that, when recast in a localized model, the DW model is not ruled out and actually gives a better fit to RSD data than $\Lambda$CDM. At the same time, the model predicts a slightly lower value of $\sigma_{8}$ than $\Lambda$CDM, in agreement with recent estimates based on lensing. We also produce analytical approximations of the two modified gravity functions and of $f\sigma_{8}(z)$ as a function of redshift. Finally, we also show how much the fit depends on initial conditions when these are generalized with respect to a standard matter-dominated era.

We discuss the cosmological implications of nonlocal modifications of general relativity containing tensorial structures. Assuming the presence of standard radiation- and matter-dominated eras, we show that, except in very particular cases, the nonlocal terms contribute a rapidly growing energy density. These models therefore generically do not have a stable cosmological evolution.

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.

The aim of this paper is to constrain modified gravity with redshift space distortion observations and supernovae measurements. Compared with a standard LCDM analysis, we include three additional free parameters, namely the initial conditions of the matter perturbations,the overall perturbation normalization, and a scale-dependent Y parameter modifying the Poisson equation, in an attempt to perform a more model-independent analysis. First, we constrain the Poisson parameter Y (also called Geff ) by using currently available fsigma_8 data and the recent SN catalog JLA. We find that the inclusion of the additional free parameters makes the constraints significantly weaker than when fixing them to the standard cosmological value. Second, we constrain Y by using forecast growth-rate data for Euclid and SKA missions. Here again we point out the weakening of the constraints when the additional parameters are included. Finally, we adopt as modified gravity Poisson parameter the specific Horndeski form, and use scale-dependent forecasts to build an exclusion plot for the Yukawa potential akin to the ones realized in laboratory experiments, both for the Euclid and the SKA surveys.

Many theories of modified gravity with higher order derivatives are usually ignored because of serious problems that appear due to an additional ghost degree of freedom. Most dangerously, it causes an immediate decay of the vacuum. However, breaking Lorentz invariance can cure such abominable behavior. By analyzing a model that describes a massive graviton together with a remaining Boulware-Deser ghost mode we show that even ghostly theories of modified gravity can yield models that are viable at both classical and quantum levels and, therefore, they should not generally be ruled out. Furthermore, we identify the most dangerous quantum scattering process that has the main impact on the decay time and find differences to simple theories that only describe an ordinary scalar field and a ghost. Additionally, constraints on the parameters of the theory including some upper bounds on the Lorentz-breaking cutoff scale are presented. In particular, for a simple theory of massive gravity we find that a breaking of Lorentz invariance is allowed to happen even at scales above the Planck mass. Finally, we discuss the relevance to other theories of modified gravity.

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range $z=0-1.6$ and wave modes below $k=10 \text{h/Mpc}$. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, $\beta^{2}$, with precision smaller than 4\% and all other cosmological parameters better than 1\%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications.

Theories of massive gravity with one or two dynamical metrics generically lack stable and observationally-viable cosmological solutions that are distinguishable from $\Lambda$CDM. We consider an extension to trimetric gravity, with three interacting spin-2 fields which are not plagued by the Boulware-Deser ghost. We systematically explore every combination with two free parameters in search of background cosmologies that are competitive with $\Lambda$CDM. For each case we determine whether the expansion history satisfies viability criteria, and whether or not it contains beyond-$\Lambda$CDM phenomenology. Among the many models we consider, there are only three cases that seem to be both viable and distinguishable from standard cosmology. One of the models has only one free parameter and displays a crossing from above to below the phantom divide. The other two provide scaling behavior, although they contain future singularities that need to be studied in more detail. These models possess interesting features that make them compelling targets for a full comparison to observations of both cosmological expansion history and structure formation.

We discuss the cosmological implications of the $R~\Box^{-2}R$ nonlocal modification to standard gravity. We relax the assumption of special initial conditions in the local formulation of the theory, perform a full phase-space analysis of the system, and show that the late-time cosmology of the model exhibits two distinct evolution paths, on which a large range of values for the present equation of state can be reached. We then compare the general solutions to supernovae data and place constraints on the parameters of the model. In particular, we find that the mass parameter of the theory should be smaller than 1.2 in Hubble units.

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.

In the current state of cosmology, where cosmological parameters are being measured to percent accuracy, it is essential to understand all sources of error to high precision. In this paper we present the results of a study of the local variations in the Hubble constant measured at the distance scale of the Coma Cluster, and test the validity of correcting for the peculiar velocities predicted by gravitational instability theory. The study is based on N-body simulations, and includes models featuring a coupling between dark energy and dark matter, as well as two $\Lambda$CDM simulations with different values of $\sigma_8$. It is found that the variance in the local flows is significantly larger in the coupled models, which increases the uncertainty in the local measurements of the Hubble constant in these scenarios. By comparing the results from the different simulations, it is found that most of the effect is caused by the higher value of $\sigma_8$ in the coupled cosmologies, though this cannot account for all of the additional variance. Given the discrepancy between different estimates of the Hubble constant in the universe today, cosmological models causing a greater cosmic variance is something that we should be aware of.

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 use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

In a Bayesian context, theoretical parameters are correlated random variables. Then, the constraints on one parameter can be improved by either measuring this parameter more precisely - or by measuring the other parameters more precisely. Especially in the case of many parameters, a lengthy process of guesswork is then needed to determine the most efficient way to improve one parameter's constraints. In this short article, we highlight an extremely simple analytical expression that replaces the guesswork and that facilitates a deeper understanding of optimization with interdependent parameters.

Massive bigravity models are interesting alternatives to standard cosmology. In most cases, however, these models have been studied for a simplified scenario in which both metrics take homogeneous and isotropic forms [Friedmann-Lema\^itre-Robertson-Walker (FLRW)] with the same spatial curvatures. The interest to consider more general geometries arises, in particular, in view of the difficulty so far encountered in building stable cosmological solutions with homogeneous and isotropic metrics. Here we consider a number of cases in which the two metrics take more general forms, namely FLRW with different spatial curvatures---Lema\^itre, Lema\^itre-Tolman-Bondi (LTB), and Bianchi I---as well as cases where only one metric is linearly perturbed. We discuss possible consistent combinations and find that only some special cases of FLRW--Lema\^itre, LTB--LTB, and FLRW--Bianchi I combinations give consistent, nontrivial solutions.

We consider the general scalar field Horndeski Lagrangian coupled to matter. Within this class of models, we present two results that are independent of the particular form of the model. First, we show that in a Friedmann-Robertson-Walker metric the Horndeski Lagrangian coincides with the pressure of the scalar field. Second, we employ the previous result to identify the most general form of the Lagrangian that allows for cosmological scaling solutions, i.e. solutions where the ratio of matter to field density and the equation of state remain constant. Scaling solutions of this kind may help solving the coincidence problem since in this case the presently observed ratio of matter to dark energy does not depend on initial conditions, but rather on the theoretical parameters.

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.

From the galaxy power spectrum in redshift space, we derive semi-analytical results on the generic degeneracy of galaxy clustering measurements. Defining the observables $\bar{A}= Gb\sigma_8$ and $\bar{R} = Gf\sigma_8$, (being $G$ the growth function, $b$ the bias, $f$ the growth rate, and $\sigma_8$ the amplitude of the power spectrum), we perform a Fisher matrix formalism to forecast the expected precision of these quantities for a Euclid-like survey. Among the results we found that galaxy surveys have generically a slightly negative correlation between $\bar{A}$ and $\bar{R}$, and they can always measure $\bar{R}$ about 3.7 to 4.7 times better than $\bar{A}$.

Dark energy is often assumed to be composed by a single scalar field. The background cosmic expansion is not sufficient to determine whether this is true or not. We study multi-field scalar-tensor models with a general dark matter source and write the observable modified gravity parameters (effective gravitational constant and anisotropic stress) in the form of a ratio of polynomials in the Fourier wavenumber k of order 2N, where N is the number of scalar fields. By comparing these observables to real data it is in principle possible to determine the number of dark energy scalar fields coupled to gravity. We also show that there are no realistic non-trivial cases in which the order of the polynomials is reduced.

Modified gravity theories predict in general a non standard equation for the propagation of gravitational waves. Here we discuss the impact of modified friction and speed of tensor modes on cosmic microwave polarization B modes. We show that the non standard friction term, parametrized by $\alpha_{M}$, is degenerate with the tensor-to-scalar ratio $r$, so that small values of $r$ can be compensated by negative constant values of $\alpha_M$. We quantify this degeneracy and its dependence on the epoch at which $\alpha_{M}$ is different from the standard, zero, value and on the speed of gravitational waves $c_{T}$. In the particular case of scalar-tensor theories, $\alpha_{M}$ is constant and strongly constrained by background and scalar perturbations, $0\le \alpha_{M}< 0.01$ and the degeneracy with $r$ is removed. In more general cases however such tight bounds are weakened and the B modes can provide useful constraints on early-time modified gravity.

The theory of bigravity offers one of the simplest possibilities to describe a massive graviton while having self-accelerating cosmological solutions without a cosmological constant. However, it has been shown recently that bigravity is affected by early-time fast growing modes on the tensor sector. Here we argue that we can only trust the linear analysis up to when perturbations are in the linear regime and use a cut-off to stop the growing of the metric perturbations. This analysis, although more consistent, still leads to growing tensor modes that are unacceptably large for the theory to be compatible with measurements of the cosmic microwave background (CMB), both in temperature and polarization spectra. In order to suppress the growing modes and make the model compatible with CMB spectra, we find it necessary to either fine-tune the initial conditions, modify the theory or set the cut-off for the tensor perturbations of the second metric much lower than unity. Initial conditions such that the growing mode is sufficiently suppresed can be achieved in scenarios in which inflation ends at the GeV scale.

A recently proposed technique allows one to constrain both the background and perturbation cosmological parameters through the distribution function of supernova Ia apparent magnitudes. Here we extend this technique to alternative cosmological scenarios, in which the growth of structure does not follow the $\Lambda$CDM prescription. We apply the method first to the supernova data provided by the JLA catalog combined with all the current independent redshift distortion data and with low-redshift cluster data from Chandra and show that although the supernovae alone are not very constraining, they help in reducing the confidence regions. Then we apply our method to future data from LSST and from a survey that approximates the Euclid satellite mission. In this case we show that the combined data are nicely complementary and can constrain the normalization $\sigma_8$ and the growth rate index $\gamma$ to within $0.6\%$ and $7\%$, respectively. In particular, the LSST supernova catalog is forecast to give the constraint $\gamma (\sigma_8/0.83)^{6.7} = 0.55 \pm 0.1$. We also report on constraints relative to a step-wise parametrization of the growth rate of structures. These results show that supernova lensing serves as a good cross-check on the measurement of perturbation parameters from more standard techniques.

Euclid is the next ESA mission devoted to cosmology. It aims at observing most of the extragalactic sky, studying both gravitational lensing and clustering over $\sim$15,000 square degrees. The mission is expected to be launched in year 2020 and to last six years. The sheer amount of data of different kinds, the variety of (un)known systematic effects and the complexity of measures require efforts both in sophisticated simulations and techniques of data analysis. We review the mission main characteristics, some aspects of the the survey and highlight some of the areas of interest to this meeting

Nonlinear, ghost-free massive gravity has two tensor fields; when both are dynamical, the mass of the graviton can lead to cosmic acceleration that agrees with background data, even in the absence of a cosmological constant. Here the question of the stability of linear perturbations in this bimetric theory is examined. Instabilities are presented for several classes of models, and simple criteria for the cosmological stability of massive bigravity are derived. In this way, we identify a particular self-accelerating bigravity model, infinite-branch bigravity (IBB), which exhibits both viable background evolution and stable linear perturbations. We discuss the modified gravity parameters for IBB, which do not reduce to the standard $\Lambda$CDM result at early times, and compute the combined likelihood from measured growth data and type Ia supernovae. IBB predicts a present matter density $\Omega_{m0}=0.18$ and an equation of state $w(z)=-0.79+0.21z/(1+z)$. The growth rate of structure is well-approximated at late times by $f(z)\approx\Omega_{m}^{0.47}[1+0.21z/(1+z)]$. The implications of the linear instability for other bigravity models are discussed: the instability does not necessarily rule these models out, but rather presents interesting questions about how to extract observables from them when linear perturbation theory does not hold.

The Time Renormalization Group (TRG) is an effective method for accurate calculations of the matter power spectrum at the scale of the first baryonic acoustic oscillations. By using a particular variable transformation in the TRG formalism, we can reduce the 2D integral in the source term of the equations of motion for the power spectrum into a series of 1D integrals. The shape of the integrand allows us to pre-compute only thirteen antiderivatives numerically, which can then be reused when evaluating the outer integral. While this introduces a few challenges to keep numerical noise under control, we find that the computation time for nonlinear corrections to the matter power spectrum decreases by a factor of 50. This opens up the possibility to use of TRG for mass production as in Markov Chain Monte Carlo methods. A Fortran code demonstrating this new algorithm has been made publicly available.

We make precise the heretofore ambiguous statement that anisotropic stress is a sign of a modification of gravity. We show that in cosmological solutions of very general classes of models extending gravity --- all scalar-tensor theories (Horndeski), Einstein-Aether models and bimetric massive gravity --- a direct correspondence exists between perfect fluids apparently carrying anisotropic stress and a modification in the propagation of gravitational waves. Since the anisotropic stress can be measured in a model-independent manner, a comparison of the behavior of gravitational waves from cosmological sources with large-scale-structure formation could in principle lead to new constraints on the theory of gravity.

Most cosmological constraints on modified gravity are obtained assuming that the cosmic evolution was standard $\Lambda$CDM in the past and that the present matter density and power spectrum normalization are the same as in a $\Lambda$CDM model. Here we examine how the constraints change when these assumptions are lifted. We focus in particular on the parameter $Y$ (also called $G_{\mathrm{eff}}$) that quantifies the deviation from the Poisson equation. This parameter can be estimated by comparing with the model-independent growth rate quantity $f\sigma_{8}(z)$ obtained through redshift distortions. We reduce the model dependency in evaluating $Y$ by marginalizing over $\sigma_{8}$ and over the initial conditions, and by absorbing the degenerate parameter $\Omega_{m,0}$ into $Y$. We use all currently available values of $f\sigma_{8}(z)$. We find that the combination $\hat{Y}=Y\Omega_{m,0}$, assumed constant in the observed redshift range, can be constrained only very weakly by current data, $\hat{Y}=0.28_{-0.23}^{+0.35}$ at 68\% c.l. We also forecast the precision of a future estimation of $\hat{Y}$ in a Euclid-like redshift survey. We find that the future constraints will reduce substantially the uncertainty, $\hat{Y}=0.30_{-0.09}^{+0.08}$ , at 68\% c.l., but the relative error on $\hat{Y}$ around the fiducial remains quite high, of the order of 30\%. The main reason for these weak constraints is that $\hat{Y}$ is strongly degenerate with the initial conditions, so that large or small values of $\hat{Y}$ are compensated by choosing non-standard initial values of the derivative of the matter density contrast. Finally, we produce a forecast of a cosmological exclusion plot on the Yukawa strength and range parameters, which complements similar plots on laboratory scales but explores scales and epochs reachable only with large-scale galaxy surveys. (abridged)

We explore the impact of modified gravity on B-modes, identifying two main separate effects: lensing and propagation of tensor modes. The location of the inflationary peak of the BB spectrum depends on the speed of gravitational waves; the amplitude of the lensing contribution depends on the anisotropic stress. We single out these effects using the quasi-static regime and considering models for which the background and the growth of matter perturbations are standard. Using available data we obtain that the gravitational wave speed is compatible with the speed of light and constrained to within about 10%.

We discuss in detail a particularly simple example of a bimetric massive gravity model which seems to offer an alternative to the standard cosmological model at background level. For small redshifts, its equation of state is $w(z)\approx-1.22_{-0.02}^{+0.02}-0.64_{-0.04}^{+0.05}z/(1+z)$. Just like $\Lambda$CDM, it depends on a single parameter, has an analytical background expansion law and fits the expansion cosmological data well. However, confirming previous results, we find that the model is unstable at early times at small scales and speculate over possible ways to cure the instability. In the regime in which the model is stable, we find that it fits the linear perturbation observations well and has a growth index of approximately $\gamma=0.47$.

A range of Bayesian tools has become widely used in cosmological data treatment and parameter inference (see Kunz, Bassett & Hlozek (2007), Trotta (2008), Amendola, Marra & Quartin (2013)). With increasingly big datasets and higher precision, tools that enable us to further enhance the accuracy of our measurements gain importance. Here we present an approach based on internal robustness, introduced in Amendola, Marra & Quartin (2013) and adopted in Heneka, Marra & Amendola (2014), to identify biased subsets of data and hidden correlation in a model independent way.

We study models of quintessence consisting of a number of scalar fields coupled to several dark matter components. In the case of exponential potentials the scaling solutions can be described in terms of a single field. The corresponding effective logarithmic slope and effective coupling can be written in a simple form in terms of the individual slopes and couplings of the original fields. We also investigate solutions where the scalar potential is negligible, in particular those leading to transient matter dominated solutions. Finally, we compute the evolution equations for the linear perturbations which will allow these models to be tested against current and future observational data.

We find the general conditions for viable cosmological solution at the background level in bigravity models. Furthermore, we constrain the parameters by comparing to the Union 2.1 supernovae catalog and identify, in some cases analytically, the best fit parameter or the degeneracy curve among pairs of parameters. We point out that a bimetric model with a single free parameter predicts a simple relation between the equation of state and the density parameter, fits well the supernovae data and is a valid and testable alternative to $\Lambda$CDM. Additionally, we identify the conditions for a phantom behavior and show that viable bimetric cosmologies cannot cross the phantom divide.

We study a general scalar field Lagrangian coupled with matter and linear in $\Box\phi$ (also called KGB model). Within this class of models, we find the most general form of the Lagrangian that allows for cosmological scaling solutions, i.e. solutions where the ratio of matter to field density and the equation of state remain constant. Scaling solutions of this kind may help solving the coincidence problem since in this case the presently observed ratio of matter to dark energy does not depend on initial conditions, but rather on the theoretical parameters. Extending previous results we find that it is impossible to join in a single solution a matter era and the scaling attractor. This is an additional step towards finding the most general scaling Lagrangian within the Horndeski class, i.e. general scalar-tensor models with second order equations of motion.

The Multi-coupled Dark Energy (McDE) scenario has been recently proposed as a specific example of a cosmological model characterized by a non-standard physics of the dark sector of the universe that nevertheless gives an expansion history which does not significantly differ from the one of the standard $\Lambda $CDM model. In this work, we present the first constraints on the McDE scenario obtained by comparing the predicted evolution of linear density perturbations with a large compilation of recent data sets for the growth rate $f\sigma_{8}$, including 6dFGS, LRG, BOSS, WiggleZ and VIPERS. Confirming qualitative expectations, growth rate data provide much tighter bounds on the model parameters as compared to the extremely loose bounds that can be obtained when only the background expansion history is considered. In particular, the $95\%$ confidence level on the coupling strength $|\beta |$ is reduced from $|\beta |\leq 83$ (background constraints only) to $|\beta |\leq 0.88$ (background and linear perturbation constraints). We also investigate how these constraints further improve when using data from future wide-field surveys such as supernova data from LSST and growth rate data from Euclid-type missions. In this case the $95\%$ confidence level on the coupling further reduce to $|\beta |\leq 0.85$. Such constraints are in any case still consistent with a scalar fifth-force of gravitational strength, and we foresee that tighter bounds might be possibly obtained from the investigation of nonlinear structure formation in McDE cosmologies.[Abridged]

Soon the number of type Ia supernova (SN) measurements should exceed 100,000. Understanding the effect of weak lensing by matter structures on the supernova brightness will then be more important than ever. Although SN lensing is usually seen as a source of systematic noise, we will show that it can be in fact turned into signal. More precisely, the non-Gaussianity introduced by lensing in the SN Hubble diagram dispersion depends rather sensitively on the amplitude sigma8 of the matter power spectrum. By exploiting this relation, we are able to predict constraints on sigma8 of 7% (3%) for a catalog of 100,000 (500,000) SNe of average magnitude error 0.12 without having to assume that such intrinsic dispersion is known a priori. The intrinsic dispersion has been assumed to be Gaussian; possible intrinsic non-Gaussianities in the dataset (due to the SN themselves and/or to other transients) could be potentially dealt with by means of additional nuisance parameters describing higher moments of the intrinsic dispersion distribution function. This method is independent of and complementary to the standard methods based on CMB, cosmic shear or cluster abundance observables.

The use of advanced statistical analysis tools is crucial in order to improve cosmological parameter estimates via removal of systematic errors and identification of previously unaccounted for cosmological signals. Here we demonstrate the application of a new fully-Bayesian method, the internal robustness formalism, to scan for systematics and new signals in the recent supernova Ia Union compilations. Our analysis is tailored to maximize chances of detecting the anomalous subsets by means of a variety of sorting algorithms. We analyse supernova Ia distance moduli for effects depending on angular separation, redshift, surveys and hemispherical directions. The data have proven to be robust within 2 sigma, giving an independent confirmation of successful removal of systematics-contaminated supernovae. Hints of new cosmology, as for example the anisotropies reported by Planck, do not seem to be reflected in the supernova Ia data.

We present a new analysis on how to distinguish between isotropic and anisotropic cosmological models based on tracking the angular displacements of a large number of distant quasars over an extended period of time, and then performing a multipole-vector decomposition of the resulting displacement maps. We find that while the GAIA mission operating at its nominal specifications does not have sufficient angular resolution to resolve anisotropic universes from isotropic ones using this method within a reasonable timespan of ten years, a next-generation GAIA-like survey with a resolution ten times better should be equal to the task. Distinguishing between different anisotropic models is however more demanding. Keeping the observational timespan to ten years, we find that the angular resolution of the survey will need to be of order 0.1 micro-arcsecond in order for certain rotating anisotropic models to produce a detectable signature that is also unique to models of this class. However, should such a detection become possible, it would immediately allow us to rule out large local void models.

The effective anisotropic stress or gravitational slip $\eta=-\Phi/\Psi$ is a key variable in the characterisation of the physical origin of the dark energy, as it allows to test for a non-minimal coupling of the dark sector to gravity in the Jordan frame. It is however important to use a fully model-independent approach when measuring $\eta$ to avoid introducing a theoretical bias into the results. In this paper we forecast the precision with which future large surveys can determine $\eta$ in a way that only relies on directly observable quantities. In particular, we do not assume anything concerning the initial spectrum of perturbations, nor on its evolution outside the observed redshift range, nor on the galaxy bias. We first leave $\eta$ free to vary in space and time and then we model it as suggested in Horndeski models of dark energy. Among our results, we find that a future large scale lensing and clustering survey can constrain $\eta$ to within 10% if $k$-independent, and to within 60% or better at $k=0.1 h/$Mpc if it is restricted to follow the Horndeski model.

We consider the consequences of having no prior knowledge of the true dark energy model for the interpretation of cosmological observations. The magnitude of redshift-space distortions and weak-lensing shear is determined by the metric on the geodesics of which galaxies and light propagate. We show that, given precise enough observations, we can use these data to completely reconstruct the metric on our past lightcone and therefore to measure the scale- and time-dependence of the anisotropic stress and the evolution of the gravitational potentials in a model-independent manner. Since both dark matter and dark energy affect the visible sector only through the gravitational field they produce, they are inseparable without a model for dark energy: galaxy bias cannot be measured and therefore the distribution of dark matter determined; the peculiar velocity of dark matter can be identified with that of the galaxies only when the equivalence principle holds. Given these limitations, we show how one can nonetheless build tests for classes of dark energy models which depend on making measurements at multiple scales at a particular redshift. They are null tests on the model-independent observables, do not require modeling evolution in time and do not require any parametrization of the free functions of these models, such as the sound speed. We show how one can rule out or constrain the whole class of the most-general scalar-tensor theories even without assuming the quasi-static limit.

With the availability of thousands of type Ia supernovae in the near future the magnitude scatter induced by lensing will become a major issue as it affects parameter estimation. Current N-body simulations are too time consuming to be integrated in the likelihood analyses used for estimating the cosmological parameters. In this paper we show that in the weak lensing regime a statistical numerical approximation produces accurate results orders of magnitude faster. We write down simple fits to the second, third and fourth central moments of the lensing magnification probability distribution as a function of redshift, of the power spectrum normalization and of the present-day matter density. We also improve upon existing models of lensing variance and show that a shifted lognormal distribution fits well the numerical one. These fits can be easily employed in cosmological likelihood analyses. Moreover, our theoretical predictions make it possible to invert the problem and begin using supernovae lensing to constrain the cosmological parameters.

The persisting consistency of ever more accurate observational data with the predictions of the standard LCDM cosmological model puts severe constraints on possible alternative scenarios, but still does not shed any light on the fundamental nature of the cosmic dark sector.As large deviations from a LCDM cosmology are ruled out by data, the path to detect possible features of alternative models goes necessarily through the definition of cosmological scenarios that leave almost unaffected the background and -- to a lesser extent -- the linear perturbations evolution of the universe. In this context,the Multi-coupled DE (McDE) model was proposed by Baldi 2012 as a particular realization of an interacting Dark Energy field characterized by an effective screening mechanism capable of suppressing the effects of the coupling at the background and linear perturbation level. In the present paper, for the first time, we challenge the McDE scenario through a direct comparison with real data, in particular with the luminosity distance of Type Ia supernovae. By studying the existence and stability conditions of the critical points of the associated background dynamical system, we select only the cosmologically consistent solutions, and confront their background expansion history with data. Confirming previous qualitative results, the McDE scenario appears to be fully consistent with the adopted sample of Type Ia supernovae, even for coupling values corresponding to an associated scalar fifth-force about four orders of magnitude stronger than standard gravity. Our analysis demonstrates the effectiveness of the McDE background screening, and shows some new non-trivial asymptotic solutions for the future evolution of the universe. Our results show how the background expansion history might be highly insensitive to the fundamental nature and to the internal complexity of the dark sector. [Abridged]

There is an approximately 9% discrepancy, corresponding to 2.4sigma, between two independent constraints on the expansion rate of the universe: one indirectly arising from the cosmic microwave background and baryon acoustic oscillations, and one more directly obtained from local measurements of the relation between redshifts and distances to sources. We argue that by taking into account the local gravitational potential at the position of the observer this tension - strengthened by the recent Planck results - is partially relieved and the concordance of the standard model of cosmology increased. We estimate that measurements of the local Hubble constant are subject to a cosmic variance of about 2.4% (limiting the local sample to redshifts z>0.010) or 1.3% (limiting it to z>0.023), a more significant correction than that taken into account already. Nonetheless, we show that one would need a very rare fluctuation to fully explain the offset in the Hubble rates. If this tension is further strengthened, a cosmology beyond the standard model may prove necessary.

We investigate constraints on early dark energy (EDE) from the Cosmic Microwave Background (CMB) anisotropy, taking into account data from WMAP9 combined with latest small scale measurements from the South Pole Telescope (SPT). For a constant EDE fraction we propose a new parametrization with one less parameter but still enough to provide similar results to the ones previously studied in literature. The main emphasis of our analysis, however, compares a new set of different EDE parametrizations that reveal how CMB constraints depend on the redshift epoch at which Dark Energy was non negligible. We find that bounds on EDE get substantially weaker if dark energy starts to be non-negligible later, with early dark energy fraction Omega_e free to go up to about 5% at 2 sigma if the onset of EDE happens at z < 100. Tight bounds around 1-2% are obtained whenever dark energy is present at last scattering, even if its effects switch off afterwards. We show that the CMB mainly constrains the presence of Dark Energy at the time of its emission, while EDE-modifications of the subsequent growth of structure are less important.

A great deal of effort is currently being devoted to understanding, estimating and removing systematic errors in cosmological data. In the particular case of type Ia supernovae, systematics are starting to dominate the error budget. Here we propose a Bayesian tool for carrying out a systematic search for systematic contamination. This serves as an extension to the standard goodness-of-fit tests and allows not only to cross-check raw or processed data for the presence of systematics but also to pin-point the data that are most likely contaminated. We successfully test our tool with mock catalogues and conclude that the Union2.1 data do not possess a significant amount of systematics. Finally, we show that if one includes in Union2.1 the supernovae that originally failed the quality cuts, our tool signals the presence of systematics at over 3.8-sigma confidence level.

We present a new prescription for analysing cosmological perturbations in a more-general class of scalar-field dark-energy models where the energy-momentum tensor has an imperfect-fluid form. This class includes Brans-Dicke models, f(R) gravity, theories with kinetic gravity braiding and generalised galileons. We employ the intuitive language of fluids, allowing us to explicitly maintain a dependence on physical and potentially measurable properties. We demonstrate that hydrodynamics is not always a valid description for describing cosmological perturbations in general scalar-field theories and present a consistent alternative that nonetheless utilises the fluid language. We apply this approach explicitly to a worked example: k-essence non-minimally coupled to gravity. This is the simplest case which captures the essential new features of these imperfect-fluid models. We demonstrate the generic existence of a new scale separating regimes where the fluid is perfect and imperfect. We obtain the equations for the evolution of dark-energy density perturbations in both these regimes. The model also features two other known scales: the Compton scale related to the breaking of shift symmetry and the Jeans scale which we show is determined by the speed of propagation of small scalar-field perturbations, i.e. causality, as opposed to the frequently used definition of the ratio of the pressure and energy-density perturbations.

The aim of this paper is to answer the following two questions: (1) Given cosmological observations of the expansion history and linear perturbations in a range of redshifts and scales as precise as is required, which of the properties of dark energy could actually be reconstructed without imposing any parameterization? (2) Are these observables sufficient to rule out not just a particular dark energy model, but the entire general class of viable models comprising a single scalar field? This paper bears both good and bad news. On one hand, we find that the goal of reconstructing dark energy models is fundamentally limited by the unobservability of the present values of the matter density Omega_m0, the perturbation normalization sigma_8 as well as the present matter power spectrum. On the other, we find that, under certain conditions, cosmological observations can nonetheless rule out the entire class of the most general single scalar-field models, i.e. those based on the Horndeski Lagrangian.

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 discuss methods based on Principal Component Analysis to constrain the dark energy equation of state using a combination of Type Ia supernovae at low redshift and spectroscopic measurements of varying fundamental couplings at higher redshifts. We discuss the performance of this method when future better-quality datasets are available, focusing on two forthcoming ESO spectrographs - ESPRESSO for the VLT and CODEX for the E-ELT - which include these measurements as a key part of their science cases. These can realize the prospect of a detailed characterization of dark energy properties almost all the way up to redshift 4.

In recent years the possibility of measuring the temporal change of radial and transverse position of sources in the sky in real time have become conceivable thanks to the thoroughly improved technique applied to new astrometric and spectroscopic experiments, leading to the research domain we call Real-time cosmology. We review for the first time great part of the work done in this field, analysing both the theoretical framework and some endeavor to foresee the observational strategies and their capability to constrain models. We firstly focus on real time measurements of the overall redshift drift and angular separation shift in distant source, able to trace background cosmic expansion and large scale anisotropy, respectively. We then examine the possibility of employing the same kind of observations to probe peculiar and proper acceleration in clustered systems and therefore the gravitational potential. The last two sections are devoted to the short time future change of the cosmic microwave background, as well as to the temporal shift of the temperature anisotropy power spectrum and maps. We conclude revisiting in this context the effort made to forecast the power of upcoming experiments like CODEX, GAIA and PLANCK in providing these new observational tools.

We consider the case of a coupling in the dark cosmological sector, where a dark energy scalar field modifies the gravitational attraction between dark matter particles. We find that the strength of the coupling {\beta} is constrained using current Cosmic Microwave Background (CMB) data, including WMAP7 and SPT, to be less than 0.063 (0.11) at 68% (95%) confidence level. Further, we consider the additional effect of the CMB-lensing amplitude, curvature, effective number of relativistic species and massive neutrinos and show that the bound from current data on {\beta} is already strong enough to be rather stable with respect to any of these variables. The strongest effect is obtained when we allow for massive neutrinos, in which case the bound becomes slightly weaker, {\beta} < 0.084(0.14). A larger value of the effective number of relativistic degrees of freedom favors larger couplings between dark matter and dark energy as well as values of the spectral index closer to 1. Adding the present constraints on the Hubble constant, as well as from baryon acoustic oscillations and supernovae Ia, we find {\beta} < 0.050(0.074). In this case we also find an interesting likelihood peak for {\beta} = 0.041 (still compatible with 0 at 1{\sigma}). This peak comes mostly from a slight difference between the Hubble parameter HST result and the WMAP7+SPT best fit. Finally, we show that forecasts of Planck+SPT mock data can pin down the coupling to a precision of better than 1% and detect whether the marginal peak we find at small non zero coupling is a real effect.

Stable massive neutral particles emitted by astrophysical sources undergo deflection under the gravitational potential of our own galaxy. The deflection angle depends on the particle velocity and therefore non-relativistic particles will be deflected more than relativistic ones. If these particles can be detected through neutrino telescopes, cosmic ray detectors or directional dark matter detectors, their arrival directions would appear aligned on the sky along the source-lens direction. On top of this deflection, the arrival direction of non-relativistic particles is displaced with respect to the relativistic counterpart also due to the relative motion of the source with respect to the observer; this induces an alignment of detections along the sky projection of the source trajectory. The final alignment will be given by a combination of the directions induced by lensing and source proper motion. We derive the deflection-velocity relation for the Milky Way halo and suggest that searching for alignments on detection maps of particle telescopes could be a way to find new particles or new astrophysical phenomena.

We compute the number density of massive Black Holes (BHs) at the centre of galaxies at z=6 in different Dynamical Dark Energy (DDE) cosmologies, and compare it with existing observational lower limits, to derive constraints on the evolution of the Dark Energy equation of state parameter w. Our approach only assumes the canonical scenario for structure formation from the collapse of overdense regions of the Dark Matter dominated primordial density field on progressively larger scales; the Black Hole accretion and merging rate have been maximized in the computation so as to obtain robust constraints on w and on its look-back time derivative w_a. Our results provide independent constraints complementary to those obtained by combining Supernovae, Cosmic Microwave Background and Baryonic Acoustic Oscillations; while the latter concern combinations of w_0 and w_a leaving the time evolution of the state parameter w_a highly unconstrained, the BH abundance mainly provide upper limits on w_a, only weakly depending on w_0. Combined with the existing constraints, our results significantly restrict the allowed region in DDE parameter space, ruling out DDE models not providing cosmic time and fast growth factor large enough to allow for the building up of the observed abundance of BHs; in particular, models with -1.2 \leq w_0 \leq -1 and positive redshift evolution w_a > 0.8 - completely consistent with previous constraints - are strongly disfavoured by our independent constraints from BH abundance. Such range of parameters corresponds to "Quintom" DDE models, with w crossing -1 starting from larger values.

Coupling dark energy to dark matter provides one of the simplest way to effectively modify gravity at large scales without strong constraints from local (i.e. solar system) observations. Models of coupled dark energy have been studied several times in the past and are already significantly constrained by cosmic microwave background experiments. In this paper we estimate the constraints that future large-scale observations will be able to put on the coupling and in general on all the parameters of the model. We combine cosmic microwave background, tomographic weak lensing, redshift distortions and power spectrum probes. We show that next-generation observations can improve the current constraint on the coupling to dark matter by two orders of magnitude; this constraint is complementary to the current solar-system bounds on a coupling to baryons.

Growing Neutrino quintessence describes a form of dynamical dark energy that could explain why dark energy dominates the universe only in recent cosmological times. This scenario predicts the formation of large scale neutrino lumps which could allow for observational tests. We perform for the first time N-body simulations of the nonlinear growth of structures for cold dark matter and neutrino fluids in the context of Growing Neutrino cosmologies. Our analysis shows a pulsation - increase and subsequent decrease - of the neutrino density contrast. This could lead to interesting observational signatures, as an enhanced bulk flow in a situation where the dark matter density contrast only differs very mildly from the standard LCDM scenario. We also determine for the first time the statistical distribution of neutrino lumps as a function of mass at different redshifts. Such determination provides an essential ingredient for a realistic estimate of the observational signatures of Growing Neutrino cosmologies. Due to a breakdown of the non-relativistic Newtonian approximation our results are limited to redshifts z > 1.

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.

Our peculiar velocity with respect to the CMB rest frame is known to induce a large dipole in the CMB. However, the motion of an observer has also the effect of distorting the anisotropies at all scales, as shown by Challinor and Van Leeuwen (2002), due to aberration and Doppler effects. We propose to measure independently our local motion by using off-diagonal two-point correlation functions for high multipoles. We study the observability of the signal for temperature and polarization anisotropies. We point out that Planck can measure the velocity $\beta$ with an error of about 30% and the direction with an error of about 20 degrees. This method constitutes a cross-check, which can be useful to verify that our CMB dipole is due mainly to our velocity or to disentangle the velocity from other possible intrinsic sources. Although in this paper we focus on our peculiar velocity, a similar effect would result also from other intrinsic vectorial distortion of the CMB which would induce a dipolar lensing. Measuring the off-diagonal correlation terms is therefore a test for a preferred direction on the CMB sky.

We investigate the impact of non-linear corrections on dark energy parameter estimation from weak lensing probes. We find that using halofit expressions, suited to LCDM models, implies substantial discrepancies with respect to results directly obtained from N-body simulations, when w(z)\neq-1. Discrepancies appear strong when using models with w'(z=0)>0, as fiducial models; they are however significant even in the neighborhood of LCDM, where neglecting the degrees of freedom associated with the DE state equation can lead to a misestimate of the matter density parameter \Omega_m.

We extend the Figure of Merit formalism usually adopted to quantify the statistical performance of future dark energy probes to assess the robustness of a future mission to plausible systematic bias. We introduce a new robustness Figure of Merit which can be computed in the Fisher Matrix formalism given arbitrary systematic biases in the observable quantities. We argue that robustness to systematics is an important new quantity that should be taken into account when optimizing future surveys. We illustrate our formalism with toy examples, and apply it to future type Ia supernova (SNIa) and baryonic acoustic oscillation (BAO) surveys. For the simplified systematic biases that we consider, we find that SNIa are a somewhat more robust probe of dark energy parameters than the BAO. We trace this back to a geometrical alignement of systematic bias direction with statistical degeneracy directions in the dark energy parameter space.

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 consider a novel class of $f(\R)$ gravity theories where the connection is related to the conformally scaled metric $\hat g_{\mu\nu}=C(\R)g_{\mu\nu}$ with a scaling that depends on the scalar curvature $\R$ only. We call them C-theories and show that the Einstein and Palatini gravities can be obtained as special limits. In addition, C-theories include completely new physically distinct gravity theories even when $f(\R)=\R$. With nonlinear $f(\R)$, C-theories interpolate and extrapolate the Einstein and Palatini cases and may avoid some of their conceptual and observational problems. We further show that C-theories have a scalar-tensor formulation, which in some special cases reduces to simple Brans-Dicke-type gravity. If matter fields couple to the connection, the conservation laws in C-theories are modified. The stability of perturbations about flat space is determined by a simple condition on the lagrangian.

The interaction between the cosmon and neutrinos may solve the "why now problem" for dark energy cosmologies. Within growing neutrino quintessence it leads to the formation of nonlinear neutrino lumps. For a test of such models by the integrated Sachs-Wolfe effect for the cosmic microwave background (CMB) we estimate the size and time evolution of the gravitational potential induced by these lumps. A population of lumps with size of 100 Mpc or more could lead to observable effects on the CMB anisotropies for low angular momenta. The linear approximation is found to be invalid for the relevant length scales. Quantitative estimates depend strongly on the details of the transition between the linear and nonlinear regimes. In particular, important backreaction effects arise from the nonlinearities of the cosmon interactions. At the present stage the uncertainties of the estimate make it difficult to constrain the parameter space of growing neutrino models. We explicitly discuss scenarios and models that are compatible with the CMB observations.

We reanalyze the supernovae data from the Union Compilation including the weak lensing effects caused by inhomogeneities. We compute the lensing probability distribution function for each background solution described by the parameters Omega_M, Omega_L and w in the presence of inhomogeneities, approximately modeled with a single-mass population of halos. We then perform a likelihood analysis in the space of FLRW-parameters and compare our results with the standard approach. We find that the inclusion of lensing can move the best-fit model significantly towards the cosmic concordance of the flat LCDM model, improving the agreement with the constraints coming from the cosmic microwave background and baryon acoustic oscillations.

The characterization of dark energy is a central task of cosmology. To go beyond a cosmological constant, we need to introduce at least an equation of state and a sound speed and consider observational tests that involve perturbations. If dark energy is not completely homogeneous on observable scales then the Poisson equation is modified and dark matter clustering is directly affected. One can then search for observational effects of dark energy clustering using dark matter as a probe. In this paper we exploit an analytical approximate solution of the perturbation equations in a general dark energy cosmology to analyze the performance of next-decade large scale surveys in constraining equation of state and sound speed. We find that tomographic weak lensing and galaxy redshift surveys can constrain the sound speed of the dark energy only if the latter is small, of the order of $c_{s}\lesssim0.01$ (in units of $c$). For larger sound speeds the error grows to 100% and more. We conclude that large scale structure observations contain very little information about the perturbations in canonical scalar field models with a sound speed of unity. Nevertheless, they are able to detect the presence of "cold" dark energy, i.e. a dark energy with non-relativistic speed of sound.

Two recently proposed techniques, involving the measurement of the cosmic parallax and redshift drift, provide novel ways of directing probing (over a time-span of several years) the background metric of the universe and therefore shed light on the dark energy conundrum. The former makes use of upcoming high-precision astrometry measurements to either observe or put tight constraints on cosmological anisotropy for off-center observers, while the latter employs high-precision spectroscopy to give an independent test of the present acceleration of the universe. In this paper, we show that both methods can break the degeneracy between LTB void models and more traditional dark energy theories. Using the near-future observational missions Gaia and CODEX we show that this distinction might be made with high confidence levels in the course of a decade.

Cosmic parallax is the change of angular separation between pair of sources at cosmological distances induced by an anisotropic expansion. An accurate astrometric experiment like Gaia could observe or put constraints on cosmic parallax. Examples of anisotropic cosmological models are Lemaitre-Tolman-Bondi void models for off-center observers (introduced to explain the observed acceleration without the need for dark energy) and Bianchi metrics. If dark energy has an anisotropic equation of state, as suggested recently, then a substantial anisotropy could arise at $z \lesssim 1$ and escape the stringent constraints from the cosmic microwave background. In this paper we show that such models could be constrained by the Gaia satellite or by an upgraded future mission.

In this paper we place new constraints on a f(R) modified gravity model recently proposed by Hu and Sawicki. After checking that the Hu and Sawicki model produces a viable cosmology, i.e. a matter dominated epoch followed by a late-time acceleration, we constrain some of its parameters by using recent observations from the UNION compilation of luminosity distances of Supernovae type Ia, including complementary information from Baryonic Acoustic Oscillations, Hubble expansion, and age data. We found that the data considered is unable to place significant constraints on the model parameters and we discuss the impact of a different assumption of the background model in cosmic parameters inference.

Refined astrometry measurements allow us to detect large-scale deviations from isotropy through real-time observations of changes in the angular separation between sources at cosmic distances. This "cosmic parallax" effect is a powerful consistency test of FRW metric and may set independent constraints on cosmic anisotropy. We apply this novel general test to LTB cosmologies with off-center observers and show that future satellite missions such as Gaia might achieve accuracies that would put limits on the off-center distance which are competitive with CMB dipole constraints.

These are the findings of the Joint Dark Energy Mission (JDEM) Figure of Merit (FoM) Science Working Group (SWG), the FoMSWG. JDEM is a space mission planned by NASA and the DOE for launch in the 2016 time frame. The primary mission is to explore the nature of dark energy. In planning such a mission, it is necessary to have some idea of knowledge of dark energy in 2016, and a way to quantify the performance of the mission. In this paper we discuss these issues.

In order for a modified gravity model to be a candidate for cosmological dark energy it has to pass stringent local gravity experiments. We find that a Brans-Dicke (BD) theory with well-defined second order corrections that include the Gauss-Bonnet term possess this feature. We construct the generic second order theory that gives, to linear order, a BD metric solution for a point-like mass source. We find that these theories interpolate between general relativity (GR) and BD gravity. In particular it is found that the relevant Eddington parameter, that is commonly heavily constrained by time delay experiments, can be arbitrarily close to the GR value of 1, with an arbitrary BD parameter. We find the region where the solution is stable to small timelike perturbations.

The Baryon Acoustic Oscillations (BAOs) or baryon wiggles which are present in the galaxy power spectrum at scales 100-150Mpc/h are powerful features with which to constrain cosmology. The potential of these probes is such that these are now included as primary science goals in the planning of several future galaxy surveys. However, there is not a uniquely defined BAO Method in the literature but a range of implementations. We study the assumptions and cosmological performances of three different BAO methods: the full Fourier space power spectrum [P(k)], the `wiggles only' in Fourier space and the spherical harmonics power spectrum [C(l)]. We contrast the power of each method to constrain cosmology for two fiducial surveys taken from the Dark Energy Task Force (DETF) report and equivalent to future ground and space based spectroscopic surveys. We find that, depending on the assumptions used, the dark energy Figure of Merit (FoM) can change by up to a factor of 35 for a given fiducial model and survey. We compare our results with the DETF implementation and, discuss the robustness of each probe, by quantifying the dependence of the FoM with the wavenumber range. The more information used by a method, the higher its statistical performance, but the higher its sensitivity to systematics and implementations details.

It has been suggested recently that the change in cosmological redshift (the Sandage test of expansion) could be observed in the next generation of large telescopes and ultra-stable spectrographs. In a recent paper we estimated the change of peculiar velocity, i.e. the peculiar acceleration, in nearby galaxies and clusters and shown it to be of the same order of magnitude as the typical cosmological signal. Mapping the acceleration field allows for a reconstruction of the galactic gravitational potential without assuming virialization. In this paper we focus on the peculiar acceleration in our own Galaxy, modeled as a Kuzmin disc and a dark matter spherical halo. We estimate the peculiar acceleration for all known Galactic globular clusters and find some cases with an expected velocity shift in excess of 20 cm/sec for observations fifteen years apart, well above the typical cosmological acceleration. We then compare the predicted signal for a MOND (modified Newtonian dynamics) model in which the spherical dark matter halo is absent. We find that the signal pattern is qualitatively different, showing that the peculiar acceleration field could be employed to test competing theories of gravity. However the difference seems too small to be detectable in the near future.

We discuss the phenomenology of the dark energy in first order perturbation theory, demonstrating that the dark energy cannot be fully constrained unless the dark matter is found, and that there are two functions that characterise the observational properties of the dark sector for cosmological probes. We argue that measuring these two functions should be an important goal for observational cosmology in the next decades.

In a recent paper [1] (PRL 98,131302,2007) we have shown that f(R)=R + mu R^{n} modified gravity dark energy models are not cosmologically viable because during the matter era that precedes the accelerated stage the cosmic expansion is given by a sim t^{1/2} rather than a sim t^{2/3}, where a is a scale factor and t is the cosmic time. A recent work [2] (PLB 639,135-143,2006) by Capozziello et al. criticised our results presenting some apparent counter-examples to our claim in f(R)= mu R^{n} models. We show here that those particular R^{n} models can produce an expansion as a sim t^{2/3} but this does not connect to a late-time acceleration. Hence, though acceptable f(R) dark energy models may exist, the R^{n} models presented in Capozziello et al. are not cosmologically viable, confirming our previous results in Ref. [1].

We identify the class of f(R) dark energy models which have a viable cosmology, i.e. a matter dominated epoch followed by a late-time acceleration. The deviation from a LambdaCDM model (f=R-Lambda) is quantified by the function m=Rf_{,RR}/f_{,R}. The matter epoch corresponds to m(r=-1) simeq +0 (where r=-Rf_{,R}/f) while the accelerated attractor exists in the region 0<m<1. We find that the equation of state w_DE of all such ``viable'' f(R) models exhibits two features: w_DE diverges at some redshift z_{c} and crosses the cosmological constant boundary (``phantom crossing'') at a redshift z_{b} smaller than z_{c}. Using the observational data of Supernova Ia and Cosmic Microwave Background, we obtain the constraint m<O(0.1) and we find that the phantom crossing could occur at z_{b}>1, i.e. within reach of observations. If we add local gravity constraints, the bound on m becomes very stringent, with m several orders of magnitude smaller than unity in the region whose density is much larger than the present cosmological density. The representative models that satisfy both cosmological and local gravity constraints take the asymptotic form m(r)=C(-r-1)^p with p>1 as r approaches -1.

A model of non-local dilaton interactions, motivated by string duality symmetries, is applied to a scenario of "coupled quintessence" in which the dilaton dark energy is non-locally coupled to the dark-matter sources. It is shown that the non-local effects tend to generate a backreaction which -- for strong enough coupling -- can automatically compensate the acceleration due to the negative pressure of the dilaton potential, thus asymptotically restoring the standard (dust-dominated) decelerated regime. This result is illustrated by analytical computations and numerical examples.

Although the Gauss-Bonnet term is a topological invariant for general relativity, it couples naturally to a quintessence scalar field, modifying gravity at solar system scales. We determine the solar system constraints due to this term by evaluating the post-Newtonian metric for a distributional source. We find a mass dependent, 1/r^7 correction to the Newtonian potential, and also deviations from the Einstein gravity prediction for light-bending. We constrain the parameters of the theory using planetary orbits, the Cassini spacecraft data, and a laboratory test of Newton's law, always finding extremely tight bounds on the energy associated to the Gauss-Bonnet term. We discuss the relevance of these constraints to late-time cosmological acceleration.

The growth factor of linear fluctuations is probably one of the least known quantity in observational cosmology. Here we discuss the constraints that baryon oscillations in galaxy power spectra from future surveys can put on a conveniently parametrized growth factor. We find that spectroscopic surveys of $5000 deg^2$ extending to $z \approx 3$ could estimate the growth index $\gamma$ within 0.06; a similar photometric survey would give $\Delta\gamma\approx 0.15$. This test provides an important consistency check for standard cosmological model and could constrain modified gravity models. We discuss the errors and the figure of merit for various combinations of redshift errors and survey size.

It has been proposed recently to observe the change in cosmological redshift of distant galaxies or quasars with the next generation of large telescope and ultra-stable spectrographs (the so-called Sandage-Loeb test). Here we investigate the possibility of observing the change in peculiar velocity in nearby clusters and galaxies. This ``peculiar acceleration'' could help reconstructing the gravitational potential without assuming virialization. We show that the expected effect is of the same order of magnitude of the cosmological velocity shift. Finally, we discuss how to convert the theoretical predictions into quantities directly related to observations.

Several experiments in the near future will test dark energy through its effects on the linear growth of matter perturbations. It is therefore important to find simple and at the same time general parametrizations of the linear growth rate. We show that a simple fitting formula that generalizes previous expressions reproduces the growth function in models that allow for a growth faster than standard, as for instance in scalar-tensor models. We use data from galaxy and Lyman-$\alpha$ power spectra to constrain the linear growth rate. We find $\gamma=0.6_{-0.3}^{+0.4}$ for the growth rate index and $\eta=0.0_{-0.2}^{+0.3}$ for the additional growth parameter we introduce.

We investigate quintessence cosmologies with a matter component consisting of particles with an increasing mass. While negligible in early cosmology, the appearance of a growing matter component has stopped the evolution of the cosmon field at a redshift around six. In turn, this has triggered the accelerated expansion of the Universe. We propose to associate growing matter with neutrinos. Then the presently observed dark energy density and its equation of state are determined by the neutrino mass.

We introduce a convenient parametrization of dark energy models that is general enough to include several modified gravity models and generalized forms of dark energy. In particular we take into account the linear perturbation growth factor, the anisotropic stress and the modified Poisson equation. We discuss the sensitivity of large scale weak lensing surveys like the proposed DUNE satellite to these parameters. We find that a large-scale weak-lensing tomographic survey is able to easily distinguish the Dvali-Gabadadze-Porrati model from LCDM and to determine the perturbation growth index to an absolute error of 0.02-0.03.

We derive the conditions under which dark energy models whose Lagrangian densities f are written in terms of the Ricci scalar R are cosmologically viable. We show that the cosmological behavior of f(R) models can be understood by a geometrical approach consisting in studying the m(r) curve on the (r, m) plane, where m=Rf_{,RR}/f_{,R} and r=-Rf_{,R}/f with f_{,R}=df/dR. This allows us to classify the f(R) models into four general classes, depending on the existence of a standard matter epoch and on the final accelerated stage. The existence of a viable matter dominated epoch prior to a late-time acceleration requires that the variable m satisfies the conditions m(r) approx+0 and dm/dr>-1 at r approx-1. For the existence of a viable late-time acceleration we require instead either (i) m=-r-1, (sqrt{3}-1)/2<m<1 and dm/dr<-1 or (ii) 0<m<1 at r=-2. These conditions identify two regions in the (r, m) space, one for the matter era and the other for the acceleration. Only models with a m(r) curve that connects these regions and satisfy the requirements above lead to an acceptable cosmology. The models of the type f(R)=alpha R^{-n} and f=R+alpha R^{-n} do not satisfy these conditions for any n>0 and n<-1 and are thus cosmologically unacceptable. Similar conclusions can be reached for many other examples discussed in the text. In most cases the standard matter era is replaced by a cosmic expansion with scale factor a propto t^{1/2}. We also find that f(R) models can have a strongly phantom attractor but in this case there is no acceptable matter era.

All $f(R)$ modified gravity theories are conformally identical to models of quintessence in which matter is coupled to dark energy with a strong coupling. This coupling induces a cosmological evolution radically different from standard cosmology. We find that in all $f(R)$ theories that behave as a power of $R$ at large or small $R$ (which include most of those proposed so far in the literature) the scale factor during the matter phase grows as $t^{1/2}$ instead of the standard law $t^{2/3}$. This behaviour is grossly inconsistent with cosmological observations (e.g. WMAP), thereby ruling out these models even if they pass the supernovae test and can escape the local gravity constraints.

Models where the dark matter component of the universe interacts with the dark energy field have been proposed as a solution to the cosmic coincidence problem, since in the attractor regime both dark energy and dark matter scale in the same way. In these models the mass of the cold dark matter particles is a function of the dark energy field responsible for the present acceleration of the universe, and different scenarios can be parameterized by how the mass of the cold dark matter particles evolves with time. In this article we study the impact of a constant coupling $\delta$ between dark energy and dark matter on the determination of a redshift dependent dark energy equation of state $w_{DE}(z)$ and on the dark matter density today from SNIa data. We derive an analytical expression for the luminosity distance in this case. In particular, we show that the presence of such a coupling increases the tension between the cosmic microwave background data from the analysis of the shift parameter in models with constant $w_{DE}$ and SNIa data for realistic values of the present dark matter density fraction. Thus, an independent measurement of the present dark matter density can place constraints on models with interacting dark energy.

Models with a scalar field coupled to the Gauss-Bonnet Lagrangian appear naturally from Kaluza-Klein compactifications of pure higher-dimensional gravity. We study linear, cosmological perturbations in the limits of weak coupling and slow-roll, and derive simple expressions for the main observable sub-horizon quantities: the anisotropic stress factor, the time-dependent gravitational constant, and the matter perturbation growth factor. Using present observational data, and assuming slow-roll for the dark energy field, we find that the fraction of energy density associated with the coupled Gauss-Bonnet term cannot exceed 15%. The bound should be treated with caution, as there are significant uncertainies in the data used to obtain it. Even so, it indicates that the future prospects for constraining the coupled Gauss-Bonnet term with cosmological observations are encouraging.

We point out that the Type Ia supernovae in the SNLS dataset are consistent with an early beginning of the cosmic acceleration if dark energy interacts strongly with dark matter. We find that the acceleration could have started at redshift z =3 and higher.

The equation of state (EOS) in quartessence models interpolates between two stages: $p\simeq 0$ at high energy densities and $p\approx -\rho$ at small ones. In the quartessence models analyzed up to now, the EOS is convex, implying increasing adiabatic sound speed ($c_{s}^{2}$) as the energy density decreases in an expanding Universe. A non-negligible $c_{s}^{2}$ at recent times is the source of the matter power spectrum problem that plagued all convex (non-silent) quartessence models. Viability for these cosmologies is only possible in the limit of almost perfect mimicry to $\Lambda$CDM. In this work we investigate if similarity to $\Lambda$CDM is also required in the class of quartessence models whose EOS changes concavity as the Universe evolves. We focus our analysis in the simple case in which the EOS has a step-like shape, such that at very early times $p\simeq0$, and at late times $p\simeq const<0$. For this class of models a non-negligible $c_{s}^{2}$ is a transient phenomenon, and could be relevant only at a more early epoch. We show that agreement with a large set of cosmological data requires that the transition between these two asymptotic states would have occurred at high redshift ($z_t\gtrsim38$). This leads us to conjecture that the cosmic expansion history of any successful non-silent quartessence is (practically) identical to the $\Lambda$CDM one.

A cosmological model that aims at solving the coincidence problem should show that dark energy and dark matter follow the same scaling solution from some time onward. At the same time, the model should contain a sufficiently long matter-dominated epoch that takes place before acceleration in order to guarantee a decelerated epoch and structure formation. So a successful cosmological model requires the occurrence of a sequence of epochs, namely a radiation era, a matter-dominated era and a final accelerated scaling attractor with $\Omega_{\phi} \simeq 0.7$. In this paper we derive the generic form of a scalar-field Lagrangian that possesses scaling solutions in the case where the coupling $Q$ between dark energy and dark matter is a free function of the field $\phi$. We then show, rather surprisingly, that the aforementioned sequence of epochs cannot occur for a vast class of generalized coupled scalar field Lagrangians that includes, to our knowledge, all scaling models in the current literature.

Dark Matter (DM) and Dark Energy (DE) can be both described in terms of ultra-light Pseudo-Goldstone-Bosons (PGB) with masses m_{DM} ~ 10^{-23}eV and m_{DE} <= 10^{-33}eV respectively. Following Barbieri et al, we entertain the possibility that a PGB exists with mass m_I intermediate between these two limits, giving a partial contribution to DM. We evaluate the related effects on the power spectrum of the matter density perturbations and on the cosmic microwave background and we derive the bounds on the density fraction, f_I, of this intermediate field from current data, with room for a better sensitivity on f_I in the near future. We also give a simple and unified analytic description of the free streaming effects both for an ultra-light scalar and for a massive neutrino.

We derive new constraints set by SNIa experiments (`gold' data sample of Riess et al.), X-ray galaxy cluster data (Allen et al. Chandra measurements of the X-ray gas mass fraction in 26 clusters), large scale structure (Sloan Digital Sky Survey spectrum) and cosmic microwave background (WMAP) on the quartessence Chaplygin model. We consider both adiabatic perturbations and intrinsic non-adiabatic perturbations such that the effective sound speed vanishes (Silent Chaplygin). We show that for the adiabatic case, only models with equation of state parameter $ |\alpha |\lesssim 10^{-2}$ are allowed: this means that the allowed models are very close to \LambdaCDM. In the Silent case, however, the results are consistent with observations in a much broader range, -0.3<\alpha<0.7.

Cosmological scaling solutions are particularly important in solving the coincidence problem of dark energy. We derive the equations of sub-Hubble linear matter perturbations for a general scalar-field Lagrangian--including quintessence, tachyon, dilatonic ghost condensate and k-essence--and solve them analytically for scaling solutions. We find that matter perturbations are always damped if a phantom field is coupled to dark matter and identify the cases in which the gravitational potential is constant. This provides an interesting possibility to place stringent observational constraints on scaling dark energy models.

We present the initial conditions for a decaying cosmological perturbation and study its signatures in the CMB anisotropies and matter power spectra. An adiabatic decaying mode in presence of components which are not described as perfect fluids (such as collisionless matter) decays slower than in a perfect-fluid dominated universe and displays super-Hubble oscillations. By including a correlated decaying mode with a red or a scale invariant spectrum, the anisotropy pattern shows super-imposed oscillations before the first Dopplear peak, while with a blue spectrum the amplitude of the secondary peaks relative to the first one and the matter power spectrum can be altered. WMAP first year data constrain the decaying to growing ratio of scale invariant adiabatic fluctuations at the matter-radiation equality to less than 10%.

We estimate the constraints that the recent high-redshift sample of supernovae type Ia put on a phenomenological interaction between dark energy and dark matter. The interaction can be interpreted as arising from the time variation of the mass of dark matter particles. We find that the coupling correlates with the equation of state: roughly speaking, a negative coupling (in our sign convention) implies phantom energy ($w_{\phi}<-1$) while a positive coupling implies ``ordinary'' dark energy. The constraints from the current supernovae Ia Hubble diagram favour a negative coupling and an equation of state $w_{\phi}<-1$. A zero or positive coupling is in fact unlikely at 99% c.l. (assuming constant equation of state); at the same time non-phantom values ($w_{\phi}>-1$) are unlikely at 95%. We show also that the usual bounds on the energy density weaken considerably when the coupling is introduced: values as large as $\Omega_{m0}=0.7$ become acceptable for as concerns SNIa. We find that the rate of change of the mass $\dot{m}/m$ of the dark matter particles is constrained to be $\delta_{0}$ in a Hubble time, with $-10<\delta_{0}<-1$ to 95% c.l.. We show that a large positive coupling might in principle avoid the future singularity known as ``big rip'' (occurring for $w_{\phi}<-1$) but the parameter region for this to occur is almost excluded by the data. We also forecast the constraints that can be obtained from future experiments, focusing on supernovae and baryon oscillations in the power spectra of deep redshift surveys. We show that the method of baryon oscillations holds the best potential to contrain the coupling.

Scalar field models with non-standard kinetic terms have been proposed in the context of k-inflation, of Born-Infeld lagrangians, of phantom energy and, more in general, of low-energy string theory. In general, scalar fields are expected to couple to matter inducing a new interaction. In this paper I derive the cosmological perturbation equations and the Yukawa correction to gravity for such general models. I find three interesting results: first, when the field behaves as phantom energy (equation of state less than -1) then the coupling strength is negative, inducing a long-range repulsive force; second, the dark energy field might cluster on astrophysical scales; third, applying the formalism to a Brans-Dicke theory with general kinetic term it is shown that its Newtonian effects depend on a single parameter that generalizes the Brans-Dicke constant.

We discuss the possible consistency of the recently discovered Type Ia supernovae at z>1 with models in which dark energy is strongly coupled to a significant fraction of dark matter, and in which an (asymptotic) accelerated phase exists where dark matter and dark energy scale in the same way. Such a coupling has been suggested for a possible solution of the coincidence problem, and is also motivated by string cosmology models of "late time" dilaton interactions. Our analysis shows that, for coupled dark energy models, the recent data are still consistent with acceleration starting as early as at $z=3$ (to within 90% c.l.), although at the price of a large "non-universality" of the dark energy coupling to different matter fields. Also, as opposed to uncoupled models which seem to prefer a ``phantom'' dark energy, we find that a large amount of coupled dark matter is compatible with present data only if the dark energy field has a conventional equation of state w>-1.

With the WMAP data we can now begin to test realistic models of inflation involving multiple scalar fields. These naturally lead to correlated adiabatic and isocurvature (entropy) perturbations with a running spectral index. We present the first full (9 parameter) likelihood analysis of double inflation with WMAP data and find that despite the extra freedom, supersymmetric hybrid potentials are strongly constrained with less than 7% correlated isocurvature component allowed when standard priors are imposed on the cosomological parameters. As a result we also find that Akaike & Bayesian model selection criteria rather strongly prefer single-field inflation, just as equivalent analysis prefers a cosmological constant over dynamical dark energy in the late universe. It appears that simplicity is the best guide to our universe.

We discuss the constraints that future photometric and spectroscopic redshift surveys can put on dark energy through the baryon oscillations of the power spectrum. We model the dark energy either with a perfect fluid or a scalar field and take into account the information contained in the linear growth function. We show that the growth function helps to break the degeneracy in the dark energy parameters and reduce the errors on $w_0,w_1$ roughly by 30% making more appealing multicolor surveys based on photometric redshifts. We find that a 200 square degrees spectroscopic survey reaching $z = 3$ can constrain $w_0,w_1$ to within $\Delta w_0=0.21,\Delta w_1=0.26$ and to $\Delta w_0=0.39,\Delta w_1=0.54$ using photometric redshifts with absolute uncertainty of 0.02. In the scalar field case we show that the slope $n$ of the inverse power-law potential for dark energy can be constrained to $\Delta n=0.26$ (spectroscopic redshifts) or $\Delta n=0.40$ (photometric redshifts), i.e. better than with future ground-based supernovae surveys or CMB data.

The skewness of the large scale distribution of matter has long been known to be a probe of gravitational clustering. Here we show that the skewness is also a probe of violation of the equivalence principle between dark matter and baryons. The predicted level of violation can be tested with the forecoming data from the Sloan Digital Sky Survey.

I review the linear and second-order perturbation theory in dark energy models with explicit interaction to matter in view of applications to N-body simulations and non-linear phenomena. Several new or generalized results are obtained: the general equations for the linear perturbation growth; an analytical expression for the bias induced by a species-dependent interaction; the Yukawa correction to the gravitational potential due to dark energy interaction; the second-order perturbation equations in coupled dark energy and their Newtonian limit. I also show that a density-dependent effective dark energy mass arises if the dark energy coupling is varying.

We present the results of a series of N-body simulations in cosmologies where dark matter (DM) is coupled to dark energy (DE), so easing the cosmic coincidence problem. The dark-dark coupling introduces two novel effects in N-body dynamics: (i) DM particle masses vary with time; (ii) gravity between DM particles is ruled by a constant $G^{*}$, greater than Newton's constant $G$, holding in other 2-body interactions. As a consequence, baryons and DM particle distributions develop a large scale bias. Here we investigate DE models with Ratra-Peebles (RP) potentials; the dark-dark coupling is set in a parametric range compatible with observations, for as concern background and linear perturbation properties. We study the halo mass function, the halo density profile and the behavior of the non-linear bias. We find that non-linear dynamics puts additional constraints to the coupling parameter. They mostly arise from density profiles, that we find to yield higher concentrations, in coupled RP models, with respect to (uncoupled) dynamical DE cosmologies. Such enhancement, although being a strong effect in some coupling parameter range, is just a minor change for smaller but significant values of the coupling parameter. With these further restrictions, coupled DE models with RP potential are consistent with non-linear observables.

We compare the WMAP temperature power spectrum and SNIa data to models with a generalized Chaplygin gas as dark energy. The generalized Chaplygin gas is a component with an exotic equation of state, p_X=-A/\rho^\alpha_X (a polytropic gas with negative constant and exponent). Our main result is that, restricting to a flat universe and to adiabatic pressure perturbations for the generalized Chaplygin gas, the constraints at 95% CL to the present equation of state w_X = p_X / \rho_X and to the parameter \alpha are -1\leq w_X < -0.8, 0 \leq \alpha <0.2, respectively. Moreover, we show that a Chaplygin gas (\alpha =1) as a candidate for dark energy is ruled out by our analysis at more than the 99.99% CL. A generalized Chaplygin gas as a unified dark matter candidate (\Omega_{CDM}=0) appears much less likely than as a dark energy model, although its \chi^2 is only two sigma away from the expected value.

The satellite experiment WMAP has produced for the first time a high-coverage, high resolution survey of the microwave sky, releasing publicly available data that are likely to remain unrivalled for years to come. Here we compare the WMAP temperature power spectrum, along with an exhautive compilation of previous experiments, to models of dark energy that allow for a tracking epoch at the present, deriving updated bounds on the dark energy equation of state and the other cosmological parameters. Moreover, we complement the analysis by including a coupling of the dark energy to dark matter. The main results are: a) the WMAP data alone constrains the equation of state of tracking dark energy to be w_\phi<-0.67(-0.49) to 68%(95%) (confining the analysis to w_\phi>-1), which implies for an inverse power law potential an exponent \alpha<0.99(2.08); b) the dimensionless coupling to dark matter is |\beta|<0.07(0.13). Including the results from the supernovae Ia further constrains the dark energy equation of state.

Most models of dark energy predict the beginning of the accelerated epoch at z <= 1. However, there are no observational or theoretical evidences in favor of such a recent start of the cosmic acceleration. In fact, a model of dark energy coupled to dark matter is explicitely constructed that a) is accelerated even at high z; b) allows structure formation during acceleration; and c) is consistent with the type Ia supernovae Hubble diagram, including the farthest known supernova SN1997ff at z~1.7 . It is shown that the accelerated epoch in this model could have started as early as z~5.

We estimate the growth of matter perturbations in a class of recently proposed dark-energy models based on the (loop-corrected) gravi-dilaton string effective action, and characterized by a global attractor epoch in which dark-matter and dark-energy density scale with the same effective equation of state. Unlike most dark-energy models, we find that the accelerated phase might start even at redshifts as high as z~5 (thus relaxing the coincidence problem), while still producing at present an acceptable level of matter fluctuations. We also show that such an early acceleration is not in conflict with the recently discovered supernova SN1997ff at z~1.7 . The comparison of the predicted value of sigma_{8} with the observational data provides interesting constraints on the fundamental parameters of the given model of dilaton-dark matter interactions.

In most models of dark energy the structure formation stops when the accelerated expansion begins. In contrast, we show that the coupling of dark energy to dark matter may induce the growth of perturbations even in the accelerated regime. In particular, we show that this occurs in the models proposed to solve the cosmic coincidence problem, in which the ratio of dark energy to dark matter is constant. Moreover, if the dark energy couples only to dark matter and not to baryons, as requested by the constraints imposed by local gravity measurements, the baryon fluctuations develop a constant, scale-independent, large-scale bias which is in principle directly observable.

We compare the latest cosmic microwave background data with theoretical predictions including correlated adiabatic and CDM isocurvature perturbations with a simple power-law dependence. We find that there is a degeneracy between the amplitude of correlated isocurvature perturbations and the spectral tilt. A negative (red) tilt is found to be compatible with a larger isocurvature contribution. Estimates of the baryon and CDM densities are found to be almost independent of the isocurvature amplitude. The main result is that current microwave background data do not exclude a dominant contribution from CDM isocurvature fluctuations on large scales.

It is well-known that even high quality cosmic microwave background (CMB) observations are not sufficient on their own to determine the equation of state of the dark energy, due to the effect of the so-called geometric degeneracy at large multipoles and the cosmic variance at small ones. In contrast, we find that CMB data can put tight constraints on another fundamental property of the dark energy, namely its coupling to dark matter. We compare the current high-resolution CMB data to models of dark energy characterized by an inverse power law or exponential potential and by the coupling to dark matter. We determine the curve of degeneracy between the dark energy equation of state and the dimensionless Hubble parameter h and show that even an independent perfect determination of h may be insufficient to distinguish dark energy from a pure cosmological constant with the current dataset. On the other hand, we find that the interaction with dark matter is firmly bounded, regardless of the potential. In terms of the dimensionless ratio \beta of the dark energy interaction to gravity, we find \beta <0.16 (95% c.l.). This implies that the effective equation of state between equivalence and tracking has been close to the pure matter equation of state within 1% and that scalar gravity is at least 40 times weaker than tensor gravity. Further, we show that an experiment limited by cosmic variance only, like the Planck mission, can put an upper bound \beta < 0.05 (95% c.l.). This shows that CMB observations have a strong potentiality not only as a test of cosmic kinematics but also as a gravitational probe.

We show that all cosmological models with an accelerated stationary global attractor reduce asymptotically to a dark energy field with an exponential potential coupled linearly to a perfect fluid dark matter. In such models the abundance of the dark components reaches a stationary value and therefore the problem of their present coincidence is solved. The requirement of a vanishing coupling of the baryons in order to pass local gravity experiments induces the existence of an intermediate baryon-dominated era. We discuss in detail the properties of these models and show that to accomodate standard nucleosynthesis they cannot produce a microwave background consistent with observations. We conclude that, among stationary models, only a time-dependent coupling or equation of state might provide a realistic cosmology.

In most models of dark energy the structure formation stops when the accelerated expansion begins. In contrast, we show that the coupling of dark energy to dark matter may induce the growth of perturbations even in the accelerated regime. In particular, we show that this occurs in the models proposed to solve the cosmic coincidence problem, in which the ratio of dark energy to dark matter is constant. Depending on the parameters, the growth may be much faster than in a standard matter-dominated era. Moreover, if the dark energy couples only to dark matter and not to baryons, as requested by the constraints imposed by local gravity measurements, the baryon fluctuations develop a constant, scale-independent, large-scale bias which is in principle directly observable. We find that a lower limit to the baryon bias b>0.5 requires the total effective parameter of state w_e=1+p/rho to be larger than 0.6 while a limit b>0.73 would rule out the model.

The estimation of cosmological parameters from cosmic microwave experiments has almost always been performed assuming gaussian data. In this paper the sensitivity of the parameter estimation to different assumptions on the probability distribution of the fluctuations is tested. Specifically, adopting the Edgeworth expansion, I show how the cosmological parameters depend on the skewness of the C_l spectrum. In the particular case of skewness independent of the multipole number I find that the primordial slope, the baryon density and the cosmological constant increase with the skewness.

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.

We present the results from the correlation analysis of two galactic redshift surveys, the extended CfA2+SSRS2 and the slice centered on $\delta=-12^{\circ}$ of Las Campanas. Furthermore, we evaluate the likelihood confidence regions for the CDM model and for the fractal model parameters. Our results indicate that, although the CDM model is in good agreement with the data, the fractal description cannot be ruled out, because of its intrinsic high variance.

We propose a cosmological model that makes a significant step toward solving the coincidence problem of the near similarity at the present of the dark energy and dark matter components. Our cosmology has the following properties: a) among flat and homogeneous spaces, the present universe is a global attractor: all the possible initial conditions lead to the observed proportion of dark energy and dark matter; once reached, it remains fixed forever; b) the expansion is accelerated at the present, as requested by the supernovae observations; c) the model is consistent with the large-scale structure and microwave background data; d) the dark energy and the dark matter densities always scale similarly after equivalence and are close to within two orders of magnitude. The model makes use of a non-linear coupling of the dark energy to the dark matter that switches on after structure formation.

We calculate analytically and numerically the Dyer-Roeder distance in perfect fluid quintessence models and give an accurate fit to the numerical solutions for all the values of the density parameter and the quintessence equation of state. Then we apply our solutions to the estimation of $H_{0}$ from multiple image time delays and find that the inclusion of quintessence modifies sensibly the likelihood distribution of $H_{0}$, generally reducing the best estimate with respect to a pure cosmological constant. Marginalizing over the other parameters ($\Omega_{m}$ and the quintessence equation of state), we obtain $H_{0}=71\pm 6$ km/sec/Mpc for an empty beam and $H_{0}=64\pm 4$ km/sec/Mpc for a filled beam. We also discuss the future prospects for distinguishing quintessence from a cosmological constant with time delays.

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.

The recent high-quality Boomerang data allow to test many competing cosmological models. Here I present a seven-parameter likelihood analysis of dark energy models with exponential potential and explicit coupling to dark matter. Such a model is conformally equivalent to a scalar field with non-minimal coupling to gravity. So far, the constraints on a dark energy - dark matter coupling were extremely weak. The Boomerang data constrain the dimensionless coupling $\beta $ to be smaller than 0.1, an order of magnitude better than previous limits. In terms of the constant $\xi $ of non-minimally coupled theories, this amounts to $\xi <0.01.$ On the other hand, Boomerang has not enough sensitivity to put constraints on the potential slope.

A new component of the cosmic medium, a light scalar field or ''quintessence '', has been proposed recently to explain cosmic acceleration with a dynamical cosmological constant. Such a field is expected to be coupled explicitely to ordinary matter, unless some unknown symmetry prevents it. I investigate the cosmological consequences of such a coupled quintessence (CQ) model, assuming an exponential potential and a linear coupling. This model is conformally equivalent to Brans-Dicke Lagrangians with power-law potential. I evaluate the density perturbations on the cosmic microwave background and on the galaxy distribution at the present and derive bounds on the coupling constant from the comparison with observational data. A novel feature of CQ is that during the matter dominated era the scalar field has a finite and almost constant energy density. This epoch, denoted as $\phi $MDE, is responsible of several differences with respect to uncoupled quintessence: the multipole spectrum of the microwave background is tilted at large angles, the acoustic peaks are shifted, their amplitude is changed, and the present 8Mpc$/h$ density variance is diminished. The present data constrain the dimensionless coupling constant to $|\beta |\leq 0.1$.

We investigate to which extent a time variation of the gravitational constant or other fundamental constants affects the best fit of the Hubble diagram of type Ia supernovae. In particular, we show that a slow increase of $G$ in the past, below experimental constraints, can reconcile the SNe Ia observations with an open zero-$\Lambda$ universe.

I analyze the density perturbations in a cosmological model with a scalar field coupled to ordinary matter, such as one obtains in string theory and in conformally transformed scalar-tensor theories. The spectrum of multipoles on the last scattering surface and the power spectrum at the present are compared with observations to derive bounds on the coupling constant and on the exponential potential slope. It is found that the acoustic peaks and the power spectrum are strongly sensitive to the model parameters. The models that best fit the galaxy spectrum and satisfy the cluster abundance test have energy density $\Omega_{\phi}\simeq 0.1$ and a scale factor expansion law $a\sim t^{p}, p\simeq 0.68$.

We investigate the clustering properties of 13 QSO lines of sight in flat space, with average redshifts from z~2 to 4. We estimate the 1-D power spectrum and the integral density of neighbours, and discuss their variation with respect to redshift and column density. We compare the results with standard CDM models, and estimate the power spectrum of Lyman-alpha clustering as a function both of redshift and column density. We find that a) there is no significant periodicity or characteristic scale; b) the clustering depends both on column density and redshift; c) the clustering increases linearly only if at the same time the HI column density decreases strongly with redshift. The results remain qualitatively the same assuming an open cosmological model.

A class of generalized non-minimal coupling theories is investigated, in search of scaling attractors able to provide an accelerated expansion at the present time. Solutions are found in the strong coupling regime and when the coupling function and the potential verify a simple relation. In such cases, which include power law and exponential functions, the dynamics is independent of the exact form of the coupling and the potential. The constraint from the time variability of $G$, however, limits the fraction of energy in the scalar field to less than 4% of the total energy density, and excludes accelerated solutions at the present.

We analyse the Las Campanas Redshift Survey using the integrated conditional density (or density of neighbors) in volume-limited subsamples up to unprecedented scales (200 Mpc/$h$) in order to determine without ambiguity the behavior of the density field. We find that the survey is well described by a fractal up to 20-30 Mpc/$h$, but flattens toward homogeneity at larger scales. Although the data are still insufficient to establish with high significance the expected homogeneous behavior, and therefore to rule out a fractal trend to larger scales, a fit with a CDM-like spectrum with high normalization well represents the data.

We consider the prospects for detecting weak gravitational lensing by underdensities (voids) in the large-scale matter distribution. We derive the basic expressions for magnification and distortion by spherical voids. Clustering of the background sources and cosmic variance are the main factors which limit in principle the detection of lensing by voids. We conclude that only voids with radii larger than $\sim 100$ \hm have lensing signal to noise larger than unity.

I present here a generalization of the maximum likelihood method and the $\chi^2$ method to the cases in which the data are {\it not} assumed to be Gaussian distributed. The method, based on the multivariate Edgeworth expansion, can find several astrophysical applications. I mention only two of them. First, in the microwave background analysis, where it cannot be excluded that the initial perturbations are non-Gaussian. Second, in the large scale structure statistics, as we already know that the galaxy distribution deviates from Gaussianity on the scales at which non-linearity is important. As a first interesting result I show here how the confidence regions are modified when non-Gaussianity is taken into account.

We present the statistical properties of large scale galaxy distribution in the LEDA redshift database. This catalog contains more than 40,000 redshifts over all the sky. We find that LEDA, although seriously affected by incompleteness, shows quite stable statistical properties. In particular, we have considered the behaviour of the two points correlation function and of the power spectrum of the density fluctuations, and we have done several tests to check whether the incompleteness of the catalog affect these statistical quantities. Our conclusion is that galaxy distribution in this catalog has fractal properties up to the distance 200 Mpc with D = 2.1 \pm 0.2, with no sign towards homogeneity: this result is statistically stable. We analyze also the properties of the angular distributions, finding a complete agreement with the results obtained in redshift space. Finally, we compare these results with those obtained in other redshift surveys, finding that this sample well reproduces the properties of galaxy distribution found in the different catalogs.

We estimate the 1--D power spectrum and the integral density of neighbors in 14 QSO lines of sight. We discuss their variation with respect to redshift and column density and compare the results with standard CDM models.

We generalize in several ways the results existing in the literature: a) we make use of an exact general relativistic solution for a spherical, nearly empty cavity in the matter dominated era to evaluate the null geodesics and the Sachs-Wolfe effect; b) we evaluate the magnitude of the adiabatic fluctuations of the photon-baryon plasma; c) we study the influence of the shell profile; and d) we take into account the finite thickness of the last scattering surface (LSS) and the influence of its position with respect to the void center. We find empirically an analytic approximation to the Sachs-Wolfe effect for all crossing geometries and we derive an upper limit of $\approx$ 25 $h^{-1}$ Mpc for the comoving radii of voids sitting on the LSS in order to achieve compatibility with COBE's data. As a nearly empty void has an overcomoving expansion of a factor of $\approx$ 4 between decoupling and the present, the maximum allowed size at present is $\approx$ 100 $h^{-1}$ Mpc. On the other hand, the smallness of the comoving size relative to the sound horizon reduces strongly the adiabatic effect by Silk damping and makes it negligible. Most of the signature of primordial voids comes therefore from metric effects and consists of subdegree spots blue or red depending on whether the center lies beyond or within the LSS. In conclusion we refine and confirm earlier constraints on a power law void spectrum originated in an inflationary phase transition and capable of generating the observed large scale structure.

We calculate analytically the bubble nucleation rate in a model of first order inflation which is able to produce large scale structure. The computation includes the first-order departure from the thin-wall limit, the explicit derivation of the pre-exponential factor, and the gravitational correction. The resulting bubble spectrum is then compared with constraints from the large scale structure and the microwave background. We show that there are models which pass all the constraints and produce bubble-like perturbations of interesting size. Furthermore, we show that it is in principle possible to reconstruct completely the inflationary two-field potential from observations.

We study the behaviour of the power spectrum (PS) in the case of fractal structures. We show that in this case the main observational features of the PS, the large scale flattening and the scaling of the amplitude with sample depth, are related to finite size effects, due to the fractal nature of galaxy distribution in the sample. Comparing with the recent results of the PS for the CfA2 we conclude that this catalog is consistent with the fractal description.

A crucial issue in cosmology is the determination of the fluctuation power spectrum.The standard picture of the matter clustering, the Cold Dark Matter model (and its variant), assumes that,on scales smaller than a certain ``flattening scale'' $\lambda_f$, the power spectrum increases with the scale, while on much larger scales it decreases so to match the tiny fluctuations observed in the microwave background. The standard picture also assumes that, once a statistically homogeneous sample is reached, the power spectrum amplitude is fixed, and any major variation should be attributed to luminosity segregation. However, the determination of $\lambda_f$ and of the absolute amplitude, if any, is still matter of debate. In particular, there is no consensus on whether the turnaround has been detected or not, and on the actual importance of the luminosity segregation effect. We show that, due to the finiteness of the sample the standard analysis of self-similar (fractal) distributions yields a turnaround for scales close to the survey scale, and a systematic amplitude shift with the survey scale. We point out that both features, bending and scaling, are in agreement with recent determination of the power spectrum, in particular with the CfA2 spectrum. We remark that the standard power spectrum is not a well defined statistical tool to characterize the galaxy distribution when a homogeneity scale has not been reached. In order to perform an analysis that does not imply any a priori assumption one should study the PS of the density, rather than the PS of the density contrast.

It is already understood that the increasing observational evidence for an open Universe may be reconciled with inflation if our horizon is contained inside one single huge bubble nucleated during the inflationary phase transition. In the scenario we present here, the Universe consists of infinitely many superhorizon bubbles, like our own, the distribution of which can be made to peak at $\Omega_0\approx 0.2$. Therefore, unlike the existing literature, we do not have to rely upon the anthropic principle nor upon special initial conditions.

We generalize the maximum likelihood method to non-Gaussian distribution functions by means of the multivariate Edgeworth expansion. We stress the potential interest of this technique in all those cosmological problems in which the determination of a non-Gaussian signature is relevant, e.g. in the analysis of large scale structure and cosmic microwave background. A first important result is that the asymptotic confidence limits on the parameters are systematically widened when the higher order correlation functions are positive, with respect to the Gaussian case.

Following a suggestion by Gasperini and Veneziano, that String Cosmology can be reconciled with Inflation and, hence, with the Standard Big Bang, we display an analytical solution which possesses four interesting properties: (1) it is non-singular; (2) it distinguishes the dynamics of the external scale factor, $a(t)$, from that of the internal one, $b(t)$; (3) it exhibits a non-monotonic behavior of $a(t)$; and (4) it stabilizes both Newton's constant and $b(t)$ (the latter to a finite, non-vanishing value). The interest of the non-monotonic evolution of $a(t)$ consists in the fact that it contains three phases of accelerated expansion, contraction and expansion before the final decelerated expansion which eventually becomes the Standard Big Bang. The total number of e-folds of the three accelerated eras can be calculated and tuned to fit the requests of observational astronomy.

To test the primordial power spectra predicted by a double inflationary model with a break of amplitude $\Delta=3$ at a scale of $2\pi/k\approx 10 \hm$ and CDM as dominant matter content, we perform PM simulations with 128$^3$ particles on a 256$^3$ grid. The broken scale invariance of the power spectra explains the extra power observed in the large-scale matter distribution. COBE-normalized spectra and a linear biasing with $b\approx 2$ are shown to reproduce the reconstructed power spectra from the CfA catalog. Identifying galactic halos with overdensity of approximately two times the cell variance, we can fit the angular correlation function using both the Limber equation and creating a APM-like angular projection with the observed luminosity function. Finally, the higher order moments of the galaxy distribution are shown to fit reasonably well the observed values.

We study the effect of the non-Gaussian clustering of galaxies on the statistics of pencil beam surveys. We find that the higher order moments of the galaxy distribution play an important role in the probability distribution for the power spectrum peaks. Taking into account the observed values for the kurtosis of galaxy distribution we derive the general probability distribution for the power spectrum modes in non-Gaussian models and show that the probability to obtain the 128$\hm$ periodicity found in pencil beam surveys is raised by roughly one order of magnitude. The non-Gaussianity of the galaxy distribution is however still insufficient to explain the reported peak-to-noise ratio of the periodicity, so that extra power on large scales seems required.

We analyze the statistical properties of bubble models for the large-scale distribution of galaxies. To this aim, we realize static simulations, in which galaxies are mostly randomly arranged in the regions surrounding bubbles. As a first test, we realize simulations of the Lick map, by suitably projecting the three-dimensional simulations. In this way, we are able to safely compare the angular correlation function implied by a bubbly geometry to that of the APM sample. We find that several bubble models provide an adequate amount of large-scale correlation, which nicely fits that of APM galaxies. Further, we apply the statistics of the count-in-cell moments to the three-dimensional distribution and compare them with available observational data on variance, skewness and kurtosis. Based on our purely geometrical constructions, we find that a well defined hierarchical scaling of higher order moments up to scales $\sim 70\hm$. The overall emerging picture is that the bubbly geometry is well suited to reproduce several aspects of large-scale clustering.

A generic outcome of theories with scalar-tensor coupling is the existence of inflationary attractors, either power-law or de Sitter. The fluctuations arising during this phase are Gaussian and their spectrum depends on the wavenumber $k$ according to the power-law $k^{1/(1-p)}$, where $p$ is the inflationary power-law exponent. We investigate to which extent these properties depend on the coupling function and on the potential. We find the class of models in which viable attractors exist. Within this class, we find that the cosmic expansion and the scaling of the fluctuation spectrum are independent of the coupling function. Further, the analytical solution of the Fokker-Planck equation shows that the deviations from Gaussianity are negligible.