The Stochastic Gravitational-Wave Background (SGWB) is expected to be a key observable for Gravitational-Wave (GW) interferometry. Its detection will open a new window on early Universe cosmology, on the astrophysics of compact objects and, as shown in this Letter, on the particle physics content of the Universe. In this Letter we show that, besides their effects on the Cosmic Microwave Background (CMB) and on Large Scale Structure (LSS), relativistic particles in the early Universe leave a clear imprint on the anisotropies of the SGWB. In particular we show that a change in the number of decoupled relativistic particles shifts the angular power spectrum of the SGWB, as both the Sachs-Wolfe (SW) and the Integrated Sachs-Wolfe (ISW) terms are affected. Being very large-angle effects, these lead to new testable predictions for future GW interferometers.

The detection of gravitational waves (GWs) propagating through cosmic structures can provide invaluable information on the geometry and content of our Universe, as well as on the fundamental theory of gravity. In order to test possible departures from General Relativity, it is essential to analyse, in a modified gravity setting, how GWs propagate through a perturbed cosmological space-time. Working within the framework of geometrical optics, we develop tools to address this topic for a broad class of scalar-tensor theories, including scenarios with non-minimal, derivative couplings between scalar and tensor modes. We determine the corresponding evolution equations for the GW amplitude and polarization tensor. The former satisfies a generalised evolution equation that includes possible effects due to a variation of the effective Planck scale; the latter can fail to be parallely transported along GW geodesics unless certain conditions are satisfied. We apply our general formulas to specific scalar-tensor theories with unit tensor speed, and then focus on GW propagation on a perturbed space-time. We determine corrections to standard formulas for the GW luminosity distance and for the evolution of the polarization tensor, which depend both on modified gravity and on the effects of cosmological perturbations. Our results can constitute a starting point to disentangle among degeneracies from different sectors that can influence GW propagation through cosmological space-times.

The detection and characterization of the Stochastic Gravitational Wave Background (SGWB) is one of the main goals of Gravitational Wave (GW) experiments. The observed SGWB will be the combination of GWs from cosmological (as predicted by many models describing the physics of the early Universe) and astrophysical origins, which will arise from the superposition of GWs from unresolved sources whose signal is too faint to be detected. Therefore, it is important to have a proper modeling of the astrophysical SGWB (ASGWB) in order to disentangle the two signals; moreover, this will provide additional information on astrophysical properties of compact objects. Applying the Cosmic Rulers formalism, we compute the observed ASGWB angular power spectrum, hence using gauge invariant quantities, accounting for all effects intervening between the source and the observer. These are the so-called projection effects, which include Kaiser, Doppler and gravitational potentials effect. Our results show that these projection effects are the most important at the largest scales, and they contribute to up to tens of percent of the angular power spectrum amplitude, with the Kaiser term being the largest at all scales. While the exact impact of these results will depend on instrumental and astrophysical details, a precise theoretical modeling of the ASGWB will necessarily need to include all these projection effects.

The Einstein Telescope (ET), a proposed European ground-based gravitational-wave detector of third-generation, is an evolution of second-generation detectors such as Advanced LIGO, Advanced Virgo, and KAGRA which could be operating in the mid 2030s. ET will explore the universe with gravitational waves up to cosmological distances. We discuss its main scientific objectives and its potential for discoveries in astrophysics, cosmology and fundamental physics.

An observable stochastic background of gravitational waves is generated whenever primordial black holes are created in the early universe thanks to a small-scale enhancement of the curvature perturbation. We calculate the anisotropies and non-Gaussianity of such stochastic gravitational waves background which receive two contributions, the first at formation time and the second due to propagation effects. The former contribution can be generated if the distribution of the curvature perturbation is characterized by a local and scale-invariant shape of non-Gaussianity. Under such an assumption, we conclude that a sizeable magnitude of anisotropy and non-Gaussianity in the gravitational waves would suggest that primordial black holes may not comply the totality of the dark matter.

A future detection of the Stochastic Gravitational Wave Background (SGWB) with GW experiments is expected to open a new window on early universe cosmology and on the astrophysics of compact objects. In this paper we study SGWB anisotropies, that can offer new tools to discriminate between different sources of GWs. In particular, the cosmological SGWB inherits its anisotropies both (i) at its production and (ii) during its propagation through our perturbed universe. Concerning (i), we show that it typically leads to anisotropies with order one dependence on frequency. We then compute the effect of (ii) through a Boltzmann approach, including contributions of both large-scale scalar and tensor linearized perturbations. We also compute for the first time the three-point function of the SGWB energy density, which can allow one to extract information on GW non-Gaussianity with interferometers. Finally, we include non-linear effects associated with long wavelength scalar fluctuations, and compute the squeezed limit of the 3-point function for the SGWB density contrast. Such limit satisfies a consistency relation, conceptually similar to what found in the literature for the case of CMB perturbations.

We study wide-angle correlations in the galaxy power spectrum in redshift space, including all general relativistic effects and the Kaiser Rocket effect in general relativity. We find that the Kaiser Rocket effect becomes important on large scales and at high redshifts, and leads to new contributions in wide-angle correlations. We believe this effect might be very important for future large volume surveys.

Foreground components in the Cosmic Microwave Background (CMB) are sparse in a needlet representation, due to their specific morphological features (anisotropy, non-Gaussianity). This leads to the possibility of applying needlet thresholding procedures as a component separation tool. In this work, we develop algorithms based on different needlet-thresholding schemes and use them as extensions of existing, well-known component separation techniques, namely ILC and template-fitting. We test soft- and hard-thresholding schemes, using different procedures to set the optimal threshold level. We find that thresholding can be useful as a denoising tool for internal templates in experiments with few frequency channels, in conditions of low signal-to-noise. We also compare our method with other denoising techniques, showing that thresholding achieves the best performance in terms of reconstruction accuracy and data compression while preserving the map resolution. The best results in our tests are in particular obtained when considering template-fitting in an LSPE like experiment, especially for B-mode spectra.

In the model where Primordial Black Holes (PBHs) form from large primordial curvature (C) perturbations, i.e., CPBHs, constraints on PBH abundance provide in principle constraints on the primordial curvature power spectrum. This connection however depends necessarily on the details of PBH formation mechanism. In this paper we provide, for the first time, constraints on the primordial curvature power spectrum from the latest limits on PBH abundance, taking into account all the steps from gravitational collapse in real space to PBH formation. In particular, we use results from numerical relativity simulations and peak theory to study the conditions for PBH formation for a range of perturbation shapes, including non-linearities, perturbation profile and a careful treatment of smoothing and filtering scales. We then obtain updated PBH formation conditions and translate that into primordial spectrum constraints for a wide range of shapes and abundances. These updated constraints cover a range of scales not probed by other cosmological observables. Our results show that the correct and accurate modelling of non-linearities, filtering and typical perturbation profile, is crucial for deriving meaningful cosmological implications.

Waveforms of gravitational waves provide information about a variety of parameters for the binary system merging. However, standard calculations have been performed assuming a FLRW universe with no perturbations. In reality this assumption should be dropped: we show that the inclusion of cosmological perturbations translates into corrections to the estimate of astrophysical parameters derived for the merging binary systems. We compute corrections to the estimate of the luminosity distance due to velocity, volume, lensing and gravitational potential effects. Our results show that the amplitude of the corrections will be negligible for current instruments, mildly important for experiments like the planned DECIGO, and very important for future ones such as the Big Bang Observer.

The Stochastic Gravitational Wave Background (SGWB) is expected to be a key observable for Gravitational Wave (GW) interferometry. Its detection will open a new window on early universe cosmology and on the astrophysics of compact objects. Using a Boltzmann approach, we study the angular anisotropies of the GW energy density, which is an important tool to disentangle the different cosmological and astrophysical contributions to the SGWB. Anisotropies in the cosmological background are imprinted both at its production, and by GW propagation through the large-scale scalar and tensor perturbations of the universe. The first contribution is not present in the Cosmic Microwave Background (CMB) radiation (as the universe is not transparent to photons before recombination), causing an order one dependence of the anisotropies on frequency. Moreover, we provide a new method to characterize the cosmological SGWB through its possible deviation from a Gaussian statistics. In particular, the SGWB will become a new probe of the primordial non-Gaussianity of the large-scale cosmological perturbations.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

This Astro-2020 White Paper deals with what we might learn from future gravitational wave observations about the early universe phase transitions and their energy scales, primordial black holes, Hubble parameter, dark matter and dark energy, modified theories of gravity and extra dimensions.

Tensions between cosmic microwave background observations and the growth of the large-scale structure inferred from late-time probes pose a serious challenge to the concordance $\Lambda$CDM cosmological model. State-of-the-art data from the Planck satellite predicts a higher rate of structure growth than what preferred by low-redshift observables. Such tension has hitherto eluded conclusive explanations in terms of straightforward modifications to $\Lambda$CDM, e.g. the inclusion of massive neutrinos or a dynamical dark energy component. Here, we investigate models of 'quartessence' -- a single dark component mimicking both dark matter and dark energy -- whose non-vanishing sound speed inhibits structure growth at late times on scales smaller than its corresponding Jeans' length. In principle, this could reconcile high- and low-redshift observations. We put this hypothesis to test against temperature and polarisation spectra from the latest Planck release, SDSS DR12 measurements of baryon acoustic oscillations and redshift-space distortions, and cosmic shear correlation functions from KiDS. This the first time that any specific model of quartessence is applied to actual data. We show that, if we naively apply $\Lambda$CDM nonlinear prescription to quartessence, the combined data sets allow for tight constraints on the model parameters. Apparently, quartessence alleviates the tension between the total matter fraction and late-time structure clustering, although in fact the tension is transferred from the latter to the quartessence sound speed parameter. However, we found that this strongly depends upon information from nonlinear scales. Indeed, if we relax this assumption, quartessence models appear still viable. For this reason, we argue that the nonlinear behaviour of quartessence deserves further investigation and may lead to a deeper understanding of the physics of the dark Universe.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We perform for the first time N-body simulations of Interacting Dark Energy assuming non-Gaussian initial conditions, with the aim of investigating possible degeneracies of these two theoretically independent phenomena in different observational probes. We focus on the large-scale matter distribution, as well as on the statistical and structural properties of collapsed halos and cosmic voids. On very large scales, we show that it is possible to choose the Interaction and non-Gaussian parameters such that their effects on the halo power spectrum cancel, and the power spectrum is indistinguishable from a $\mathrm{\Lambda CDM}$ model. On small scales, measurements of the non-linear matter power spectrum, halo-matter bias, halo and subhalo mass function and cosmic void number function validate the degeneracy determined on large scales. However, the internal structural properties of halos and cosmic voids, namely halo concentration-mass relation and void density profile, are very different from those measured in the $\mathrm{\Lambda CDM}$ model, thereby breaking the degeneracy. In practice, the values of $\mathrm{f_{NL}}$ required to cancel the effect of interaction are already ruled by observations. Our results show in principle that the combination of large- and small-scale probes is needed to constrain Interacting Dark Energy and Primordial non-Gaussianity separately.

Using the consistency relation in Fourier space, we derive the observed galaxy bispectrum from single-field inflation in the squeezed limit, in which one of the three modes has a wavelength much longer than the other two. This provides a non-trivial check of the full computation of the bispectrum based on second-order cosmological perturbation theory in this limit. We show that gauge modes need to be carefully removed in the second-order cosmological perturbations in order to calculate the observed galaxy bispectrum in the squeezed limit. We then give an estimate of the effective non-Gaussianity due to general relativistic lightcone effects that could mimic a primordial non-Gaussian signal.

Given the important role that the galaxy bispectrum has recently acquired in cosmology and the scale and precision of forthcoming galaxy clustering observations, it is timely to derive the full expression of the large-scale bispectrum going beyond approximated treatments which neglect integrated terms or higher-order bias terms or use the Limber approximation. On cosmological scales, relativistic effects that arise from observing on the past light-cone alter the observed galaxy number counts, therefore leaving their imprints on N-point correlators at all orders. In this paper we compute for the first time the bispectrum including all general relativistic, local and integrated, effects at second order, the tracers' bias at second order, geometric effects as well as the primordial non-Gaussianity contribution. This is timely considering that future surveys will probe scales comparable to the horizon where approximations widely used currently may not hold; neglecting these effects may introduce biases in estimation of cosmological parameters as well as primordial non-Gaussianity.

We study up to second order the galaxy number over-density that depends on magnification in redshift space on cosmological scales for a concordance model. The result contains all general relativistic effects up to second order which arise from observing on the past light cone, including all redshift and lensing distortions, contributions from velocities, Sachs-Wolfe, integrated SW and time-delay terms. We find several new terms and contributions that could be potentially important for an accurate calculation of the bias on estimates of non-Gaussianity and on precision parameter estimates.

We study the detectability of large-scale velocity effects on galaxy clustering, by simulating galaxy surveys and combining the clustering of different types of tracers of large-scale structure. We employ a set of lognormal mocks that simulate a $20.000$ deg$^2$ near-complete survey up to $z=0.8$, in which each galaxy mock traces the spatial distribution of dark matter of that mock with a realistic bias prescription. We find that the ratios of the monopoles of the power spectra of different types of tracers carry most of the information that can be extracted from a multi-tracer analysis. In particular, we show that with a multi-tracer technique it will be possible to detect velocity effects with $\gtrsim 3 \sigma$. Finally, we investigate the degeneracy of these effects with the (local) non-Gaussianity parameter $f_{\rm NL}$, and how large-scale velocity contributions could be mistaken for the signatures of primordial non-Gaussianity.

We introduce a method to create mock galaxy catalogues in redshift space including general relativistic effects to linear order in the cosmological perturbations. We dub our method LIGER, short for `light cones with general relativity'. LIGER takes a (N-body or hydrodynamic) Newtonian simulation as an input and outputs the distribution of galaxies in comoving redshift space. This result is achieved making use of a coordinate transformation and simultaneously accounting for lensing magnification. The calculation includes both local corrections and terms that have been integrated along the line of sight. Our fast implementation allows the production of many realizations that can be used to forecast the performance of forthcoming wide-angle surveys and to estimate the covariance matrix of the observables. To facilitate this use, we also present a variant of LIGER designed for large-volume simulations with low mass resolution. In this case, the galaxy distribution on large scales is obtained by biasing the matter-density field. Finally, we present two sample applications of LIGER. First, we discuss the impact of weak gravitational lensing onto the angular clustering of galaxies in a Euclid-like survey. In agreement with previous analytical studies, we find that magnification bias can be measured with high confidence. Second, we focus on two generally neglected Doppler-induced effects: magnification and the change of number counts with redshift. We show that the corresponding redshift-space distortions can be detected at 5.5$\sigma$ significance with the completed Square Kilometre Array.

We study the nonlinear gravitational dynamics of a universe filled with a pressureless fluid and a cosmological constant $\Lambda$ in the context of Newtonian gravity, and in the relativistic post-Friedmann approach proposed in paper I [I. Milillo et al., Phys. Rev. D 92, 023519 (2015).]. The post-Friedmann approximation scheme is based on the $1/c$ expansion of the space-time metric and the energy-momentum tensor, and includes nonlinear Newtonian cosmology. Here we establish the nonlinear post-Friedmann framework in the Lagrangian-coordinates approach for structure formation. For this we first identify a Lagrangian gauge which is suitable for incorporating nonzero vorticity. We analyze our results in two limits: at the leading order we recover the fully nonlinear Newtonian cosmological equations in the Lagrangian formulation, and we provide a space-time metric consistent from the perspective of general relativity. We then linearize our expressions and recover the relativistic results at first order in cosmological perturbation theory. Therefore, the introduced approximation scheme provides a unified treatment for the two leading-order regimes, from the small scales described by Newtonian gravity to the large linear scale, where first-order relativistic cosmological perturbation theory gives a very good description of structure formation.

We present the galaxy number overdensity up to second order in redshift space on cosmological scales for a concordance model. The result contains all general relativistic effects up to second order that arise from observing on the past light cone, including all redshift effects, lensing distortions from convergence and shear, and contributions from velocities, Sachs-Wolfe, integrated SW and time-delay terms. This result will be important for accurate calculation of the bias on estimates of non-Gaussianity and on precision parameter estimates, introduced by nonlinear projection effects.

We present a detailed derivation of the observed galaxy number over-density on cosmological scales up to second order in perturbation theory. We include all relativistic effects that arise from observing on the past lightcone. The derivation is in a general gauge, and applies to all dark energy models (including interacting dark energy) and many modified gravity models. The result will be important for accurate cosmological parameter estimation, including non-Gaussianity, since all projection effects need to be taken into account. It also offers the potential for new probes of General Relativity, dark energy and modified gravity. This paper accompanies Paper I which presents the key results for the concordance model in Poisson gauge.

Galaxy clustering on very large scales can be probed via the 2-point correlation function in the general case of wide and deep separations, including all the lightcone and relativistic effects. Using our recently developed formalism, we analyze the behavior of the local and integrated contributions and how these depend on redshift range, linear and angular separations and luminosity function. Relativistic corrections to the local part of the correlation can be non-negligible but they remain generally sub-dominant. On the other hand, the additional correlations arising from lensing convergence and time-delay effects can become very important and even dominate the observed total correlation function. We investigate different configurations formed by the observer and the pair of galaxies, and we find that the case of near-radial large-scale separations is where these effects will be the most important.

In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.

We study the parity-odd part (that we shall call Doppler term) of the linear galaxy two-point correlation function that arises from wide-angle, velocity, Doppler lensing and cosmic acceleration effects. As it is important at low redshift and at large angular separations, the Doppler term is usually neglected in the current generation of galaxy surveys. For future wide-angle galaxy surveys such as Euclid, SPHEREx and SKA, however, we show that the Doppler term must be included. The effect of these terms is dominated by the magnification due to relativistic aberration effects and the slope of the galaxy redshift distribution and it generally mimics the effect of the local type primordial non-Gaussianity with the effective nonlinearity parameter $f_{\rm NL}^{\rm eff}$ of a few, we show that this would affect forecasts on measurements of $f_{\rm NL}$ at low-redshift. Our results show that a survey at low redshift with large number density over a wide area of the sky could detect the Doppler term with a signal-to-noise ratio of $\sim 1-20$, depending on survey specifications.

We discuss the question of gauge choice when analysing relativistic density perturbations at second order. We compare Newtonian and General Relativistic approaches. Some misconceptions in the recent literature are addressed. We show that the comoving-synchronous gauge is the unique gauge in General Relativity that corresponds to the Lagrangian frame and is entirely appropriate to describe the matter overdensity at second order. The comoving-synchronous gauge is the simplest gauge in which to describe Lagrangian bias at second order.

We extend previous analyses of wide-angle correlations in the galaxy power spectrum in redshift space to include all general relativistic effects. These general relativistic corrections to the standard approach become important on large scales and at high redshifts, and they lead to new terms in the wide-angle correlations. We show that in principle the new terms can produce corrections of nearly 10 % on Gpc scales over the usual Newtonian approximation. General relativistic corrections will be important for future large-volume surveys such as SKA and Euclid, although the problem of cosmic variance will present a challenge in observing this.

We investigate how well future large-scale radio surveys could measure different shapes of primordial non-Gaussianity; in particular we focus on angle-dependent non-Gaussianity arising from primordial anisotropic sources, whose bispectrum has an angle dependence between the three wavevectors that is characterized by Legendre polynomials $\mathcal{P}_L$ and expansion coefficients $c_L$. We provide forecasts for measurements of galaxy power spectrum, finding that Large-Scale Structure (LSS) data could allow measurements of primordial non-Gaussianity competitive or improving upon current constraints set by CMB experiments, for all the shapes considered. We argue that the best constraints will come from the possibility to assign redshift information to radio galaxy surveys, and investigate a few possible scenarios for the EMU and SKA surveys. A realistic (futuristic) modeling could provide constraints of $f_{\rm NL}^{\rm loc} \approx 1 (0.5)$ for the local shape, $f_{\rm NL}$ of $\mathcal{O}(10) (\mathcal{O}(1))$ for the orthogonal, equilateral and folded shapes, and $c_{L=1} \approx 80 (2)$, $c_{L=2} \approx 400 (10)$ for angle-dependent non-Gaussianity. The more futuristic forecasts show the potential of LSS analyses to considerably improve current constraints on non-Gaussianity, and so on models of the primordial Universe. Finally, we find the minimum requirements that would be needed to reach $\sigma(c_{L=1})=10$, which can be considered as a typical (lower) value predicted by some (inflationary) models.

We investigate the cosmological dependence and the constraining power of large-scale galaxy correlations, including all redshift-distortions, wide-angle, lensing and gravitational potential effects on linear scales. We analyze the cosmological information present in the lensing convergence and in the gravitational potential terms describing the so-called "relativistic effects," and we find that, while smaller than the information contained in intrinsic galaxy clustering, it is not negligible. We investigate how neglecting them does bias cosmological measurements performed by future spectroscopic and photometric large-scale surveys such as SKA and Euclid. We perform a Fisher analysis using the CLASS code, modified to include scale-dependent galaxy bias and redshift-dependent magnification and evolution bias. Our results show that neglecting relativistic terms introduces an error in the forecasted precision in measuring cosmological parameters of the order of a few tens of percent, in particular when measuring the matter content of the Universe and primordial non-Gaussianity parameters. Therefore, we argue that radial correlations and integrated relativistic terms need to be taken into account when forecasting the constraining power of future large-scale number counts of galaxy surveys.

We present a nonlinear post-Friedmann framework for structure formation, generalizing to cosmology the weak-field (post-Minkowskian) approximation, unifying the treatment of small and large scales. We consider a universe filled with a pressureless fluid and a cosmological constant $\Lambda$, the theory of gravity is Einstein's general relativity and the background is the standard flat $\Lambda$CDM cosmological model. We expand the metric and the energy-momentum tensor in powers of $1/c$, keeping the matter density and peculiar velocity as exact fundamental variables. We assume the Poisson gauge, including scalar and tensor modes up to $1/c^4$ order and vector modes up to $1/c^5$ terms. Through a redefinition of the scalar potentials as a resummation of the metric contributions at different orders, we obtain a complete set of nonlinear equations, providing a unified framework to study structure formation from small to superhorizon scales, from the nonlinear Newtonian to the linear relativistic regime. We explicitly show the validity of our scheme in the two limits: at leading order we recover the fully nonlinear equations of Newtonian cosmology; when linearized, our equations become those for scalar and vector modes of first-order relativistic perturbation theory in the Poisson gauge. Tensor modes are nondynamical at the $1/c^4$ order we consider: they are purely nonlinear and describe a distortion of the spatial slices determined at this order by a constraint, quadratic in the scalar and vector variables. The main results of our analysis are as follows: (a) at leading order a purely Newtonian nonlinear energy current sources a frame-dragging gravitomagnetic vector potential, and (b) in the leading-order Newtonian regime and in the linear relativistic regime the two scalar metric potentials are the same, while the nonlinearity of general relativity makes them different.

The observed galaxy power spectrum acquires relativistic corrections from lightcone effects, and these corrections grow on very large scales. Future galaxy surveys in optical, infrared and radio bands will probe increasingly large wavelength modes and reach higher redshifts. In order to exploit the new data on large scales, an accurate analysis requires inclusion of the relativistic effects. This is especially the case for primordial non-Gaussianity and for extending tests of dark energy models to horizon scales. Here we investigate the latter, focusing on models where the dark energy interacts non-gravitationally with dark matter. Interaction in the dark sector can also lead to large-scale deviations in the power spectrum. If the relativistic effects are ignored, the imprint of interacting dark energy will be incorrectly identified and thus lead to a bias in constraints on interacting dark energy on very large scales.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

If dark energy and dark matter interact via exchange of energy and momentum, then this may affect the galaxy power spectrum on large scales. When this happens, it may be degenerate with the signal from primordial non-Gaussianity via scale-dependent bias. We consider a class of interacting dark energy models and show that the matter overdensity is scale-dependent on large scales. We estimate the effective non-Gaussianity arising from the large-scale effects of interaction in the dark sector. The signal of dark sector interaction can be disentangled from a primordial non-Gaussian signal by measuring the power at two redshifts.

We investigate the properties of the 2-point galaxy correlation function at very large scales, including all geometric and local relativistic effects -- wide-angle effects, redshift space distortions, Doppler terms and Sachs-Wolfe type terms in the gravitational potentials. The general three-dimensional correlation function has a nonzero dipole and octupole, in addition to the even multipoles of the flat-sky limit. We study how corrections due to primordial non-Gaussianity and General Relativity affect the multipolar expansion, and we show that they are of similar magnitude (when f_NL is small), so that a relativistic approach is needed. Furthermore, we look at how large-scale corrections depend on the model for the growth rate in the context of modified gravity, and we discuss how a modified growth can affect the non-Gaussian signal in the multipoles.

Quintessence can cluster only on horizon scales. What is the effect on the observed matter distribution? To answer this, we need a relativistic approach that goes beyond the standard Newtonian calculation and deals properly with large scales. Such an approach has recently been developed for the case when dark energy is vacuum energy, which does not cluster at all. We extend this relativistic analysis to deal with dynamical dark energy. Using three quintessence potentials as examples, we compute the angular power spectrum for the case of an HI intensity map survey. Compared to the concordance model with the same small-scale power at z=0, quintessence boosts the angular power by up to ~15% at high redshifts, while power in the two models converges at low redshifts. The difference is mainly due to the background evolution, driven mostly by the normalization of the power spectrum today. The dark energy perturbations make only a small contribution on the largest scales, and a negligible contribution on smaller scales. Ironically, the dark energy perturbations remove the false boost of large-scale power that arises if we impose the (unphysical) assumption that the dark energy is smooth.

In this paper we obtain the bispectrum of dark matter density perturbations in the frame of covariant cubic Galileon theories. This result is obtained by means of a semi-analytic approach to second-order perturbations in Galileon cosmologies, assuming Gaussian initial conditions. In particular, we show that, even in the presence of large deviations of the linear growth-rate w.r.t. the {\Lambda}CDM one, at the bispectrum level such deviations are reduced to a few percent.

We propose an analytic procedure that allows to determine quantitatively the deviation in the behavior of cosmological perturbations between a given f(R) modified gravity model and a LCDM reference model. Our method allows to study structure formation in these models from the largest scales, of the order of the Hubble horizon, down to scales deeply inside the Hubble radius, without employing the so-called "quasi-static" approximation. Although we restrict our analysis here to linear perturbations, our technique is completely general and can be extended to any perturbative order.

We use Luminous Red Galaxies from the Sloan Digital Sky Survey II to test the cosmological structure growth in two alternatives to the standard LCDM+GR cosmological model. We compare observed three-dimensional clustering in SDSS DR7 with theoretical predictions for the standard vanilla LCDM+GR model, Unified Dark Matter cosmologies and the normal branch DGP. In computing the expected correlations in UDM cosmologies, we derive a parameterized formula for the growth factor in these models. For our analysis we apply the methodology tested in Raccanelli et al. 2010 and use the measurements of Samushia et al. 2011, that account for survey geometry, non-linear and wide-angle effects and the distribution of pair orientation. We show that the estimate of the growth rate is potentially degenerate with wide-angle effects, meaning that extremely accurate measurements of the growth rate on large scales will need to take such effects into account. We use measurements of the zeroth and second order moments of the correlation function from SDSS DR7 data and the Large Suite of Dark Matter Simulations, and perform a likelihood analysis to constrain the parameters of the models. Using information on the clustering up to r_max = 120 Mpc/h, and after marginalizing over the bias, we find, for UDM models, a speed of sound < 6.1e-4, and, for the nDGP model, a cross-over scale r_c > 340 Mpc, at 95% confidence level.

We present parameter estimation forecasts for future 3D cosmic shear surveys for a class of Unified Dark Matter (UDM) models, where a single scalar field mimics both Dark Matter (DM) and Dark Energy (DE). These models have the advantage that they can describe the dynamics of the Universe with a single matter component providing an explanation for structure formation and cosmic acceleration. A crucial feature of the class of UDM models we use in this work is characterised by a parameter, c_inf (in units of the speed of light c=1), that is the value of the sound speed at late times, and on which structure formation depends. We demonstrate that the properties of the DM-like behaviour of the scalar field can be estimated with very high precision with large-scale, fully 3D weak lensing surveys. We found that 3D weak lensing significantly constrains c_inf, and we find minimal errors 0.00003, for the fiducial value c_inf=0.001, and 0.000026, for c_inf=0.012. Moreover, we compute the Bayesian evidence for UDM models over the LCDM model as a function of c_inf. For this purpose, we can consider the LCDM model as a UDM model with c_inf=0. We find that the expected evidence clearly shows that the survey data would unquestionably favour UDM models over the LCDM model, for the values c_inf>0.001.

We present an analysis of the cross-correlation between the CMB and the large-scale structure (LSS) of the Universe in Unified Dark Matter (UDM) scalar field cosmologies. We work out the predicted cross-correlation function in UDM models, which depends on the speed of sound of the unified component, and compare it with observations from six galaxy catalogues (NVSS, HEAO, 2MASS, and SDSS main galaxies, luminous red galaxies, and quasars). We sample the value of the speed of sound and perform a likelihood analysis, finding that the UDM model is as likely as the LambdaCDM, and is compatible with observations for a range of values of c_\infinity (the value of the sound speed at late times) on which structure formation depends. In particular, we obtain an upper bound of c_\infinity^2 \leq 0.009 at 95% confidence level, meaning that the LambdaCDM model, for which c_\infinity^2 = 0, is a good fit to the data, while the posterior probability distribution peaks at the value c_\infinity^2=10^(-4) . Finally, we study the time dependence of the deviation from LambdaCDM via a tomographic analysis using a mock redshift distribution and we find that the largest deviation is for low-redshift sources, suggesting that future low-z surveys will be best suited to constrain UDM models.

In this work we analyze and review cosmological models in which the dynamics of a single scalar field accounts for a unified description of the Dark Matter and Dark Energy sectors, dubbed Unified Dark Matter (UDM) models. In this framework, we consider the general Lagrangian of k-essence, which allows to find solutions around which the scalar field describes the desired mixture of Dark Matter and Dark Energy. We also discuss static and spherically symmetric solutions of Einstein's equations for a scalar field with non-canonical kinetic term, in connection with galactic halo rotation curves.

We investigate the time evolution of the gravitational potential Phi for a special class of non-adiabatic Unified Dark matter Models described by scalar field lagrangians. These models predict the same background evolution as in the LambdaCDM and possess a non-vanishing speed of sound. We provide a very accurate approximation of Phi, valid after the recombination epoch, in the form of a Bessel function of the first kind. This approximation may be useful for a future deeper analysis of Unified Dark Matter scalar field models.

We investigate the general properties of Unified Dark Matter (UDM) scalar field models with Lagrangians with a non-canonical kinetic term, looking specifically for models that can produce a fast transition between an early Einstein-de Sitter CDM-like era and a later Dark Energy like phase, similarly to the barotropic fluid UDM models in JCAP1001(2010)014. However, while the background evolution can be very similar in the two cases, the perturbations are naturally adiabatic in fluid models, while in the scalar field case they are necessarily non-adiabatic. The new approach to building UDM Lagrangians proposed here allows to escape the common problem of the fine-tuning of the parameters which plague many UDM models. We analyse the properties of perturbations in our model, focusing on the the evolution of the effective speed of sound and that of the Jeans length. With this insight, we can set theoretical constraints on the parameters of the model, predicting sufficient conditions for the model to be viable. An interesting feature of our models is that what can be interpreted as w_{DE} can be <-1 without violating the null energy conditions.

We investigate the general properties of Unified Dark Matter (UDM) fluid models where the pressure and the energy density are linked by a barotropic equation of state (EoS) $p = p(\rho)$ and the perturbations are adiabatic. The EoS is assumed to admit a future attractor that acts as an effective cosmological constant, while asymptotically in the past the pressure is negligible. UDM models of the dark sector are appealing because they evade the so-called "coincidence problem" and "predict" what can be interpreted as $w_{\rm DE} \approx -1$, but in general suffer the effects of a non-negligible Jeans scale that wreak havoc in the evolution of perturbations, causing a large Integrated Sachs-Wolfe effect and/or changing structure formation at small scales. Typically, observational constraints are violated, unless the parameters of the UDM model are tuned to make it indistinguishable from $\Lambda$CDM. Here we show how this problem can be avoided, studying in detail the functional form of the Jeans scale in adiabatic UDM perturbations and introducing a class of models with a fast transition between an early Einstein-de Sitter CDM-like era and a later $\Lambda$CDM-like phase. If the transition is fast enough, these models may exhibit satisfactory structure formation and CMB fluctuations. To consider a concrete case, we introduce a toy UDM model and show that it can predict CMB and matter power spectra that are in agreement with observations for a wide range of parameter values.

Weak gravitational lensing is a powerful tool for studying both the geometry and the dynamics of the Universe. Its power spectrum contains information on the sources emitting photons and on the large--scale structures that these hotons cross. We calculate the weak lensing cosmic convergence and shear power spectra, in linear theory and Limber's approximation, for two different classes of cosmological models: the standard \LambdaCDM and Unified Dark Matter (UDM) models. The latter models attempt to solve the problems of the dark matter in the dynamics of galaxies and galaxy clusters and of the late-time acceleration of the Universe expansion by introducing a scalar field that mimics both dark matter and dark energy. A crucial feature of the UDM models is the speed of sound c_\infty, that is the value of the sound speed at late times, on which structure formation depends. In this paper, we provide the predictions of the UDM models for the weak lensing signal,with various values of c_\infty. We consider both the Cosmic Microwave Background and background galaxies at different redshifts as sources for the weak lensing power spectra. We find that the power spectra in UDM models are more sensitive to the variations of c_\infty for sources located at low redshifts. Moreover, we find that for c_\infty>10^{-3} (in units of the speed of light), the UDM weak lensing convergence power spectrum C^{\kappa\kappa}(l) for background galaxies is strongly suppressed with respect to the \LambdaCDM spectrum, particularly at small angular scales l\gtrsim100.

We use scalar-field Lagrangians with a non-canonical kinetic term to obtain unified dark matter models where both the dark matter and the dark energy, the latter mimicking a cosmological constant, are described by the scalar field itself. In this framework, we propose a technique to reconstruct models where the effective speed of sound is small enough that the scalar field can cluster. These models avoid the strong time evolution of the gravitational potential and the large Integrated Sachs-Wolfe effect which have been a serious drawback of previously considered models. Moreover, these unified dark matter scalar field models can be easily generalized to behave as dark matter plus a dark energy component behaving like any type of quintessence fluid.

We perform an analytical study of the Integrated Sachs-Wolfe (ISW) effect within the framework of Unified Dark Matter models based on a scalar field which aim at a unified description of dark energy and dark matter. Computing the temperature power spectrum of the Cosmic Microwave Background anisotropies we are able to isolate those contributions that can potentially lead to strong deviations from the usual ISW effect occurring in a $\Lambda$CDM universe. This helps to highlight the crucial role played by the sound speed in the Unified Dark Matter models. Our treatment is completely general in that all the results depend only on the speed of sound of the dark component and thus it can be applied to a variety of unified models, including those which are not described by a scalar field but relies on a single dark fluid.

We investigate the static and spherically symmetric solutions of Einstein's equations for a scalar field with non-canonical kinetic term, assumed to provide both the dark matter and dark energy components of the Universe. In particular, we give a prescription to obtain solutions (dark halos) whose rotation curve v_c(r) is in good agreement with observational data. We show that there exist suitable scalar field Lagrangians that allow to describe the cosmological background evolution and the static solutions with a single dark fluid.

Considering the general Lagrangian of k-essence models, we study and classify them through variables connected to the fluid equation of state parameter w_\kappa. This allows to find solutions around which the scalar field describes a mixture of dark matter and cosmological constant-like dark energy, an example being the purely kinetic model proposed by Scherrer. Making the stronger assumption that the scalar field Lagrangian is exactly constant along solutions of the equation of motion, we find a general class of k-essence models whose classical trajectories directly describe a unified dark matter/dark energy (cosmological constant) fluid. While the simplest case of a scalar field with canonical kinetic term unavoidably leads to an effective sound speed c_s=1, thereby inhibiting the growth of matter inhomogeneities, more general non-canonical k-essence models allow for the possibility that c_s << 1 whenever matter dominates.