We study perturbation theory for large-scale structure in the most general scalar-tensor theories propagating a single scalar degree of freedom, which include Horndeski theories and beyond. We model the parameter space using the effective field theory of dark energy. For Horndeski theories, the gravitational field and fluid equations are invariant under a combination of time-dependent transformations of the coordinates and fields. This symmetry allows one to construct a physical adiabatic mode which fixes the perturbation-theory kernels in the squeezed limit and ensures that the well-known consistency relations for large-scale structure, originally derived in general relativity, hold in modified gravity as well. For theories beyond Horndeski, instead, one generally cannot construct such an adiabatic mode. Because of this, the perturbation-theory kernels are modified in the squeezed limit and the consistency relations for large-scale structure do not hold. We show, however, that the modification of the squeezed limit depends only on the linear theory. We investigate the observational consequences of this violation by computing the matter bispectrum. In the squeezed limit, the largest effect is expected when considering the cross-correlation between different tracers. Moreover, the individual contributions to the 1-loop matter power spectrum do not cancel in the infrared limit of the momentum integral, modifying the power spectrum on non-linear scales.

We introduce $k$-evolution, a relativistic $N$-body code based on $\textit{gevolution}$, which includes clustering dark energy among its cosmological components. To describe dark energy, we use the effective field theory approach. In particular, we focus on $k$-essence with a speed of sound much smaller than unity but we lay down the basis to extend the code to other dark energy and modified gravity models. We develop the formalism including dark energy non-linearities but, as a first step, we implement the equations in the code after dropping non-linear self-coupling in the $k$-essence field. In this simplified setup, we compare $k$-evolution simulations with those of $\texttt{CLASS}$ and $\textit{gevolution}$ 1.2, showing the effect of dark matter and gravitational non-linearities on the power spectrum of dark matter, of dark energy and of the gravitational potential. Moreover, we compare $k$-evolution to Newtonian $N$-body simulations with back-scaled initial conditions and study how dark energy clustering affects massive halos.

We point out that dark-energy perturbations may become unstable in the presence of a gravitational wave of sufficiently large amplitude. We study this effect for the cubic Horndeski operator (braiding), proportional to $\alpha_{\rm B}$. The scalar that describes dark-energy fluctuations features ghost and/or gradient instabilities for gravitational-wave amplitudes that are produced by typical binary systems. Taking into account the populations of binary systems, we conclude that the instability is triggered in the whole Universe for $|\alpha_{\rm B} |\gtrsim 10^{-2}$, i.e. when the modification of gravity is sizeable. The instability is triggered by massive black-hole binaries down to frequencies corresponding to $10^{10}$ km: the instability is thus robust, unless new physics enters on even longer wavelengths. The fate of the instability and the subsequent time-evolution of the system depend on the UV completion, so that the theory may end up in a state very different from the original one. The same kind of instability is present in beyond-Horndeski theories for $|\alpha_{\rm H}| \gtrsim 10^{-20}$. In conclusion, the only dark-energy theories with sizeable cosmological effects that avoid these problems are $k$-essence models, with a possible conformal coupling with matter.

We develop an approach to parametrize cosmological perturbations beyond linear order for general dark energy and modified gravity models characterized by a single scalar degree of freedom. We derive the full nonlinear action, focusing on Horndeski theories. In the quasi-static, non-relativistic limit, there are a total of six independent relevant operators, three of which start at nonlinear order. The new nonlinear couplings modify, beyond linear order, the generalized Poisson equation relating the Newtonian potential to the matter density contrast. We derive this equation up to cubic order in perturbations and, in a companion article, we apply it to compute the one-loop matter power spectrum. Within this approach, we also discuss the Vainshtein regime around spherical sources and the relation between the Vainshtein scale and the nonlinear scale for structure formation.

We develop an approach to compute observables beyond the linear regime of dark matter perturbations for general dark energy and modified gravity models. We do so by combining the Effective Field Theory of Dark Energy and Effective Field Theory of Large-Scale Structure approaches. In particular, we parametrize the linear and nonlinear effects of dark energy on dark matter clustering in terms of the Lagrangian terms introduced in a companion paper, focusing on Horndeski theories and assuming the quasi-static approximation. The Euler equation for dark matter is sourced, via the Newtonian potential, by new nonlinear vertices due to modified gravity and, as in the pure dark matter case, by the effects of short-scale physics in the form of the divergence of an effective stress tensor. The effective fluid introduces a counterterm in the solution to the matter continuity and Euler equations, which allows a controlled expansion of clustering statistics on mildly nonlinear scales. We use this setup to compute the one-loop dark-matter power spectrum.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, taking into account the large occupation numbers. We describe dark energy using the effective field theory approach, in the context of generalized scalar-tensor theories. When the $m_3^3$ (cubic Horndeski) and $\tilde m_4^2$ (beyond Horndeski) operators are present, the gravitational wave acts as a classical background for $\pi$ and modifies its dynamics. In particular, $\pi$ fluctuations are described by a Mathieu equation and feature instability bands that grow exponentially. Focusing on the regime of small gravitational-wave amplitude, corresponding to narrow resonance, we calculate analytically the produced $\pi$, its energy and the change of the gravitational-wave signal. The resonance is affected by $\pi$ self-interactions in a way that we cannot describe analytically. This effect is very relevant for the operator $m_3^3$ and it limits the instability. In the case of the $\tilde m_4^2$ operator self-interactions can be neglected, at least in some regimes. The modification of the gravitational-wave signal is observable for $3 \times 10^{-20} \lesssim \alpha_{\rm H} \lesssim 10^{-17}$ with a LIGO/Virgo-like interferometer and for $10^{-16} \lesssim \alpha_{\rm H} \lesssim 10^{-10}$ with a LISA-like one.

We study the screening mechanism in the most general scalar-tensor theories that leave gravitational waves unaffected and are thus compatible with recent LIGO/Virgo observations. Using the effective field theory of dark energy approach, we consider the general action for perturbations beyond linear order, focussing on the quasi-static limit. When restricting to the subclass of theories that satisfy the gravitational wave constraints, the fully nonlinear effective Lagrangian contains only three independent parameters. One of these, $\beta_1$, is uniquely present in degenerate higher-order theories. We compute the two gravitational potentials for a spherically symmetric matter source and we find that for $\beta_1 \ge 0$ they decrease as the inverse of the distance, as in standard gravity, while the case $\beta_1 < 0$ is ruled out. For $\beta_1 > 0$, the two potentials differ and their gravitational constants are not the same on the inside and outside of the body. Generically, the bound on anomalous light bending in the Solar System constrains $\beta_1 \lesssim 10^{-5}$. Standard gravity can be recovered outside the body by tuning the parameters of the model, in which case $\beta_1 \lesssim 10^{-2}$ from the Hulse-Taylor pulsar.

A light scalar degree of freedom, as the one possibly responsible for the accelerated expansion of the Universe, could leave observable traces in the inspiral gravitational wave signal of binary systems. In order to study these effects, we extend the effective field theory formalism of Goldberger and Rothstein to minimal scalar-tensor theories of gravity. This class of models is still very broad, because the couplings of the scalar to matter are far less constrained than those a massless spin-2 field. In most of the paper we focus on conformal couplings. Using the effective field theory approach, we discuss the emergence of violations of the strong equivalence principle even in models that exhibit universality of free fall at the microscopic level. Our results on the conservative dynamics of the binary and its power emitted in gravitational and scalar radiation agree with those obtained with the standard post-Newtonian formalism. We also compare them to more recent work. Finally, we discuss the implications of considering a disformal coupling to matter.

We show that theories of mimetic gravity can be viewed as degenerate higher-order scalar-tensor (DHOST) theories that admit an extra local (gauge) symmetry in addition to the usual diffeomorphism invariance. We reformulate and classify mimetic theories in this perspective. Using the effective theory of dark energy, recently extended to include DHOST theories, we then investigate the linear perturbations about a homogeneous and isotropic background for all mimetic theories. We also include matter, in the form of a $k$-essence scalar field, and we derive the quadratic action for linear perturbations in this case.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, through the processes $\gamma \to \pi\pi$ and $\gamma \to \gamma \pi$, made possible by the spontaneous breaking of Lorentz invariance. Within the EFT of Dark Energy (or Horndeski/beyond Horndeski theories) the first process is large for the operator $\frac12 \tilde m_4^2(t) \, \delta g^{00}\, \left( {}^{(3)}\! R + \delta K_\mu^\nu \delta K^\mu_\nu -\delta K^2 \right)$, so that the recent observations force $\tilde m_4 =0$ (or equivalently $\alpha_{\rm H}=0$). This constraint, together with the requirement that gravitational waves travel at the speed of light, rules out all quartic and quintic GLPV theories. Additionally, we study how the same couplings affect the propagation of gravitons at loop order. The operator proportional to $\tilde m_4^2$ generates a calculable, non-Lorentz invariant higher-derivative correction to the graviton propagation. The modification of the dispersion relation provides a bound on $\tilde m_4^2$ comparable to the one of the decay. Conversely, operators up to cubic Horndeski do not generate sizeable higher-derivative corrections.

The observation of GW170817 and its electromagnetic counterpart implies that gravitational waves travel at the speed of light, with deviations smaller than a few $\times 10^{-15}$. We discuss the consequences of this experimental result for models of dark energy and modified gravity characterized by a single scalar degree of freedom. To avoid tuning, the speed of gravitational waves must be unaffected not only for our particular cosmological solution, but also for nearby solutions obtained by slightly changing the matter abundance. For this to happen the coefficients of various operators must satisfy precise relations that we discuss both in the language of the Effective Field Theory of Dark Energy and in the covariant one, for Horndeski, beyond Horndeski and degenerate higher-order theories. The simplification is dramatic: of the three functions describing quartic and quintic beyond Horndeski theories, only one remains and reduces to a standard conformal coupling to the Ricci scalar for Horndeski theories. We show that the deduced relations among operators do not introduce further tuning of the models, since they are stable under quantum corrections.

We compare analytical computations with numerical simulations for dark-matter clustering, in general relativity and in the normal branch of DGP gravity (nDGP). Our analytical frameword is the Effective Field Theory of Large-Scale Structure (EFTofLSS), which we use to compute the one-loop dark-matter power spectrum, including the resummation of infrared bulk displacement effects. We compare this to a set of 20 COLA simulations at redshifts $z = 0$, $z=0.5$, and $z =1$, and fit the free parameter of the EFTofLSS, called the speed of sound, in both $\Lambda$CDM and nDGP at each redshift. At one-loop at $z = 0$, the reach of the EFTofLSS is $k_{\rm reach}\approx 0.14 \, h { \rm Mpc^{-1}}$ for both $\Lambda$CDM and nDGP. Along the way, we compare two different infrared resummation schemes and two different treatments of the time dependence of the perturbative expansion, concluding that they agree to approximately $1\%$ over the scales of interest. Finally, we use the ratio of the COLA power spectra to make a precision measurement of the difference between the speeds of sound in $\Lambda$CDM and nDGP, and verify that this is proportional to the modification of the linear coupling constant of the Poisson equation.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

We study generalized disformal transformations, including derivatives of the metric, in the context of the Effective Field Theory of Inflation. All these transformations do not change the late-time cosmological observables but change the coefficients of the operators in the action: some couplings are effectively redundant. At leading order in derivatives and up to cubic order in perturbations, one has 6 free functions that can be used to set to zero 6 of the 17 operators at this order. This is used to show that the tensor three-point function cannot be modified at leading order in derivatives, while the scalar-tensor-tensor correlator can only be modified by changing the scalar dynamics. At higher order in derivatives there are transformations that do not affect the Einstein-Hilbert action: one can find 6 additional transformations that can be used to simplify the inflaton action, at least when the dynamics is dominated by the lowest derivative terms. We also identify the leading higher-derivative corrections to the tensor power spectrum and bispectrum.

Most existing theories of dark energy and/or modified gravity, involving a scalar degree of freedom, can be conveniently described within the framework of the Effective Theory of Dark Energy, based on the unitary gauge where the scalar field is uniform. We extend this effective approach by allowing the Lagrangian in unitary gauge to depend on the time derivative of the lapse function. Although this dependence generically signals the presence of an extra scalar degree of freedom, theories that contain only one propagating scalar degree of freedom, in addition to the usual tensor modes, can be constructed by requiring the initial Lagrangian to be degenerate. Starting from a general quadratic action, we derive the dispersion relations for the linear perturbations around Minkowski and a cosmological background. Our analysis directly applies to the recently introduced Degenerate Higher-Order Scalar-Tensor (DHOST) theories. For these theories, we find that one cannot recover a Poisson-like equation in the static linear regime except for the subclass that includes the Horndeski and so-called "beyond Horndeski" theories. We also discuss Lorentz-breaking models inspired by Horava gravity.

We explore general scalar-tensor models in the presence of a kinetic mixing between matter and the scalar field, which we call Kinetic Matter Mixing. In the frame where gravity is de-mixed from the scalar this is due to disformal couplings of matter species to the gravitational sector, with disformal coefficients that depend on the gradient of the scalar field. In the frame where matter is minimally coupled, it originates from the so-called beyond Horndeski quadratic Lagrangian. We extend the Effective Theory of Interacting Dark Energy by allowing disformal coupling coefficients to depend on the gradient of the scalar field as well. In this very general approach, we derive the conditions to avoid ghost and gradient instabilities and we define Kinetic Matter Mixing independently of the frame metric used to described the action. We study its phenomenological consequences for a $\Lambda$CDM background evolution, first analytically on small scales. Then, we compute the matter power spectrum and the angular spectra of the CMB anisotropies and the CMB lensing potential, on all scales. We employ the public version of COOP, a numerical Einstein-Boltzmann solver that implements very general scalar-tensor modifications of gravity. Rather uniquely, Kinetic Matter Mixing weakens gravity on short scales, predicting a lower $\sigma_8$ with respect to the $\Lambda$CDM case. We propose this as a possible solution to the tension between the CMB best-fit model and low-redshift observables.

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

We explore the phenomenological consequences of general late-time modifications of gravity in the quasi-static approximation, in the case where cold dark matter is non-minimally coupled to the gravitational sector. Assuming spectroscopic and photometric surveys with configuration parameters similar to those of the Euclid mission, we derive constraints on our effective description from three observables: the galaxy power spectrum in redshift space, tomographic weak-lensing shear power spectrum and the correlation spectrum between the integrated Sachs-Wolfe effect and the galaxy distribution. In particular, with $\Lambda$CDM as fiducial model and a specific choice for the time dependence of our effective functions, we perform a Fisher matrix analysis and find that the unmarginalized $68\%$ CL errors on the parameters describing the modifications of gravity are of order $\sigma\sim10^{-2}$--$10^{-3}$. We also consider two other fiducial models. A nonminimal coupling of CDM enhances the effects of modified gravity and reduces the above statistical errors accordingly. In all cases, we find that the parameters are highly degenerate, which prevents the inversion of the Fisher matrices. Some of these degeneracies can be broken by combining all three observational probes.

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.

We consider the superpotential formalism to describe the evolution of scalar fields during inflation, generalizing it to include the case with non-canonical kinetic terms. We provide a characterization of the attractor behaviour of the background evolution in terms of first and second slow-roll parameters (which need not be small). We find that the superpotential is useful in justifying the separate universe approximation from the gradient expansion, and also in computing the spectra of primordial perturbations around attractor solutions in the $\delta N$ formalism. As an application, we consider a class of models where the background trajectories for the inflaton fields are derived from a product separable superpotential. In the perspective of the holographic inflation scenario, such models are dual to a deformed CFT boundary theory, with $D$ mutually uncorrelated deformation operators. We compute the bulk power spectra of primordial adiabatic and entropy cosmological perturbations, and show that the results agree with the ones obtained by using conformal perturbation theory in the dual picture.

We introduce a new class of scalar-tensor theories that extend Horndeski, or "generalized galileon", models. Despite possessing equations of motion of higher order in derivatives, we show that the true propagating degrees of freedom obey well-behaved second-order equations and are thus free from Ostrogradski instabilities, in contrast to the standard lore. Remarkably, the covariant versions of the original galileon Lagrangians-obtained by direct replacement of derivatives with covariant derivatives-belong to this class of theories. These extensions of Horndeski theories exhibit an uncommon, interesting phenomenology: the scalar degree of freedom affects the speed of sound of matter, even when the latter is minimally coupled to gravity.

Effective theories of a scalar $\phi$ invariant under the internal \textit{galileon symmetry} $\phi\to\phi+b_\mu x^\mu$ have been extensively studied due to their special theoretical and phenomenological properties. In this paper, we introduce the notion of \textit{weakly broken galileon invariance}, which characterizes the unique class of couplings of such theories to gravity that maximally retain their defining symmetry. The curved-space remnant of the galileon's quantum properties allows to construct (quasi) de Sitter backgrounds largely insensitive to loop corrections. We exploit this fact to build novel cosmological models with interesting phenomenology, relevant for both inflation and late-time acceleration of the universe.

We present a unifying treatment of dark energy and modified gravity that allows distinct conformal-disformal couplings of matter species to the gravitational sector. In this very general approach, we derive the conditions to avoid ghost and gradient instabilities. We compute the equations of motion for background quantities and linear perturbations. We illustrate our formalism with two simple scenarios, where either cold dark matter or a relativistic fluid is nonminimally coupled. This extends previous studies of coupled dark energy to a much broader spectrum of gravitational theories.

We have recently proposed a new class of gravitational scalar-tensor theories free from Ostrogradski instabilities, in arXiv:1404.6495. As they generalize Horndeski theories, or "generalized" galileons, we call them G$^3$. These theories possess a simple formulation when the time hypersurfaces are chosen to coincide with the uniform scalar field hypersurfaces. We confirm that they contain only three propagating degrees of freedom by presenting the details of the Hamiltonian formulation. We examine the coupling between these theories and matter. Moreover, we investigate how they transform under a disformal redefinition of the metric. Remarkably, these theories are preserved by disformal transformations that depend on the scalar field gradient, which also allow to map subfamilies of G$^3$ into Horndeski theories.

We review and extend a novel approach that we introduced recently, to describe general dark energy or scalar-tensor models. Our approach relies on an ADM formulation based on the hypersurfaces where the underlying scalar field is uniform. The advantage of this approach is that it can describe in the same language and in a minimal way a vast number of existing models, such as quintessence models, $F(R)$ theories, scalar tensor theories, their Horndeski extensions and beyond. It also naturally includes Horava-Lifshitz theories. As summarized in this review, our approach provides a unified treatment of the linear cosmological perturbations about a FLRW universe, obtained by a systematic expansion of our general action up to quadratic order. This shows that the behaviour of these linear perturbations is generically characterized by five time-dependent functions. We derive the full equations of motion in the Newtonian gauge, and obtain in particular the equation of state for dark energy perturbations, in the Horndeski case, in terms of these functions. Our unifying description thus provides the simplest and most systematic way to confront theoretical models with current and future cosmological observations.

We show that the prediction for the primordial tensor power spectrum cannot be modified at leading order in derivatives. Indeed, one can always set to unity the speed of propagation of gravitational waves during inflation by a suitable disformal transformation of the metric, while a conformal one can make the Planck mass time-independent. Therefore, the tensor amplitude unambiguously fixes the energy scale of inflation. Using the Effective Field Theory of Inflation, we check that predictions are independent of the choice of frame, as expected. The first corrections to the standard prediction come from two parity violating operators with three derivatives. Also the correlator $<\gamma\gamma\gamma>$ is standard and only receives higher derivative corrections. These results hold also in multifield models of inflation and in alternatives to inflation and make the connection between a (quasi) scale-invariant tensor spectrum and inflation completely robust.

The recently derived consistency relations for Large Scale Structure do not hold if the Equivalence Principle (EP) is violated. We show it explicitly in a toy model with two fluids, one of which is coupled to a fifth force. We explore the constraints that galaxy surveys can set on EP violation looking at the squeezed limit of the 3-point function involving two populations of objects. We find that one can explore EP violations of order 10^{-3} - 10^{-4} on cosmological scales. Chameleon models are already very constrained by the requirement of screening within the Solar System and only a very tiny region of the parameter space can be explored with this method. We show that no violation of the consistency relations is expected in Galileon models.

We generalize the recently derived single-field consistency relations of Large Scale Structure in two directions. First, we treat the effect of the long modes (with momentum q) on the short ones (with momentum k) non-perturbatively, by writing resummed consistency relations which do not require k/q \delta_q << 1. These relations do not make any assumptions on the short-scales physics and are extended to include (an arbitrary number of) multiple long modes, internal lines with soft momenta and soft loops. We do several checks of these relations in perturbation theory and we verify that the effect of soft modes always cancels out in equal-time correlators. Second, we write the relations directly in redshift space, without assuming the single-stream approximation: not only the long mode affects the short scales as a homogeneous gravitational field, but it also displaces them by its velocity along the line-of-sight. Redshift space consistency relations still vanish when short modes are taken at equal time: an observation of a signal in the squeezed limit would point towards multifield inflation or a violation of the equivalence principle.

We derive consistency relations for the late universe (CDM and \Lambda CDM): relations between an n-point function of the density contrast \delta and an (n+1)-point function in the limit in which one of the (n+1) momenta becomes much smaller than the others. These are based on the observation that a long mode, in single-field models of inflation, reduces to a diffeomorphism since its freezing during inflation all the way until the late universe, even when the long mode is inside the horizon (but out of the sound horizon). These results are derived in Newtonian gauge, at first and second order in the small momentum q of the long mode and they are valid non-perturbatively in the short-scale \delta. In the non-relativistic limit our results match with (Kehagias and Riotto '13) and (Peloso and Pietroni '13). These relations are a consequence of diffeomorphism invariance; they are not satisfied in the presence of extra degrees of freedom during inflation or violation of the Equivalence Principle (extra forces) in the late universe.

We compute the full cosmic microwave background temperature bispectrum generated by nonlinearities after single-field inflation. By integrating the photon temperature at second order along a perturbed geodesic in Newtonian gauge, we derive an expression for the observed temperature fluctuations that, for the first time, clarifies the separation of the gravitational lensing and time-delay effects from the purely second-order contributions. We then use the second-order Boltzmann code CosmoLib$2^{\rm nd}$ to calculate these contributions and their bispectrum. Including the perturbations in the photon path, the numerically computed bispectrum exactly matches the expected squeezed limit. Moreover, the analytic squeezed-limit formula reproduces well the signal-to-noise and shape of the full bispectrum, potentially facilitating the subtraction of the bias induced by second-order effects. For a cosmic-variance limited experiment with $l_{\rm max} = 2000$, the bias on a local signal is $f_{\rm NL}^{\rm loc} =0.73$ negligible for equilateral and orthogonal signals. The signal-to-noise ratio is unity at $l_{\rm max} \sim 3000$, suggesting that second-order effects may hopefully be measured in the future.

We use the so-called eikonal approximation, recently introduced in the context of cosmological perturbation theory, to compute power spectra for multi-component fluids. We demonstrate that, at any given order in standard perturbation theory, multipoint power spectra do not depend on the large-scale adiabatic modes. Moreover, we employ perturbation theories to decipher how nonadiabatic modes, such as a relative velocity between two different components, damp the small-scale matter power spectrum, a mechanism recently described in the literature. In particular, we do an explicit calculation at 1-loop order of this effect. While the 1-loop result eventually breaks down, we show how the damping effect can be fully captured by the help of the eikonal approximation. A relative velocity not only induces mode damping but also creates large-scale anisotropic modulations of the matter power spectrum amplitude. We illustrate this for the Local Group environment.

We introduce the eikonal approximation to study the effect of the large-scale motion of cosmic fluids on their small-scale evolution. This approach consists in collecting the impact of the long-wavelength displacement field into a single or finite number of random variables, whose statistical properties can be computed from the initial conditions. For a single dark matter fluid, we show that we can recover the nonlinear propagators of renormalized perturbation theory. These are obtained with no need to assume that the displacement field follows the linear theory. Then we extend the eikonal approximation to many fluids. In particular, we study the case of two non-relativistic components and we derive their resummed propagators in the presence of isodensity modes. Unlike the adiabatic case, where only the phase of small-scale modes is affected by the large-scale advection field, the isodensity modes change also the amplitude on small scales. We explicitly solve the case of cold dark matter-baryon mixing and find that the isodensity modes induce only very small corrections to the resummed propagators.

The effective field theory of cosmological perturbations stems from considering a cosmological background solution as a state displaying spontaneous breaking of time translations and (adiabatic) perturbations as the related Nambu-Goldstone modes. With this insight, one can systematically develop a theory for the cosmological perturbations during inflation and, with minor modifications, also describe in full generality the gravitational interactions of dark energy, which are relevant for late-time cosmology. The formalism displays a unique set of Lagrangian operators containing an increasing number of cosmological perturbations and derivatives. We give an introductory description of the unitary gauge formalism for theories with broken gauge symmetry---that allows to write down the most general Lagrangian---and of the Stueckelberg "trick"---that allows to recover gauge invariance and to make the scalar field explicit. We show how to apply this formalism to gravity and cosmology and we reproduce the detailed analysis of the action in the ADM variables. We also review some basic applications to inflation and dark energy.

We propose a minimal description of single field dark energy/modified gravity within the effective field theory formalism for cosmological perturbations, which encompasses most existing models. We start from a generic Lagrangian given as an arbitrary function of the lapse and of the extrinsic and intrinsic curvature tensors of the time hypersurfaces in unitary gauge, i.e. choosing as time slicing the uniform scalar field hypersurfaces. Focusing on linear perturbations, we identify seven Lagrangian operators that lead to equations of motion containing at most two (space or time) derivatives, the background evolution being determined by the time dependent coefficients of only three of these operators. We then establish a dictionary that translates any existing or future model whose Lagrangian can be written in the above form into our parametrized framework. As an illustration, we study Horndeski's-or generalized Galileon-theories and show that they can be described, up to linear order, by only six of the seven operators mentioned above. This implies, remarkably, that the dynamics of linear perturbations can be more general than that of Horndeski while remaining second order. Finally, in order to make the link with observations, we provide the entire set of linear perturbation equations in Newtonian gauge, the effective Newton constant in the quasi-static approximation and the ratio of the two gravitational potentials, in terms of the time-dependent coefficients of our Lagrangian.

We propose a universal description of dark energy and modified gravity that includes all single-field models. By extending a formalism previously applied to inflation, we consider the metric universally coupled to matter fields and we write in terms of it the most general unitary gauge action consistent with the residual unbroken symmetries of spatial diffeomorphisms. Our action is particularly suited for cosmological perturbation theory: the background evolution depends on only three operators. All other operators start at least at quadratic order in the perturbations and their effects can be studied independently and systematically. In particular, we focus on the properties of a few operators which appear in non-minimally coupled scalar-tensor gravity and galileon theories. In this context, we study the mixing between gravity and the scalar degree of freedom. We assess the quantum and classical stability, derive the speed of sound of fluctuations and the renormalization of the Newton constant. The scalar can always be de-mixed from gravity at quadratic order in the perturbations, but not necessarily through a conformal rescaling of the metric. We show how to express covariant field-operators in our formalism and give several explicit examples of dark energy and modified gravity models in our language. Finally, we discuss the relation with the covariant EFT methods recently appeared in the literature.

We compute the cosmic microwave background temperature bispectrum generated by nonlinearities at recombination on all scales. We use CosmoLib$2^{\rm nd}$, a numerical Boltzmann code at second-order to compute CMB bispectra on the full sky. We consistently include all effects except gravitational lensing, which can be added to our result using standard methods. The bispectrum is peaked on squeezed triangles and agrees with the analytic approximation in the squeezed limit at the few per cent level for all the scales where this is applicable. On smaller scales, we recover previous results on perturbed recombination. For cosmic-variance limited data to $l_{\rm max} =2000$, its signal-to-noise is $S/N=0.47$ and will bias a local signal by $f_{\rm NL}^{\rm loc}\simeq 0.82$.

We compute the matter bispectrum in the presence of primordial local non-Gaussianity over a wide range of scales, including the very small nonlinear ones. We use the Halo Model approach, considering non-Gaussian corrections to the halo profiles, the halo mass function and the bias functions. We compare our results in the linear and mildly nonlinear scales to a large ensemble of Gaussian and non-Gaussian numerical simulations. We consider both squeezed and equilateral configurations, at redshift z = 0 and z = 1. For z = 0, the deviations between the Halo Model and the simulations are smaller than 10% in the squeezed limit, both in the Gaussian and non-Gaussian cases. The Halo Model allows to make predictions on scales much smaller than those reached by numerical simulations. For local non-Gaussian initial conditions with a parameter fNL = 100, we find an enhancement of the bispectrum in the squeezed configuration k = k3 = k2 >> k1 \sim 0.01 h^{-1} Mpc, of \sim 15% and \sim 25% on scales k \sim 1 h^{-1} Mpc, at z = 0 and z = 1 respectively. This is mainly due to the non-Gaussian corrections in the linear bias. Finally we provide a very simple expression valid for any scenario, i.e. for any choice of the halo profile, mass and bias functions, which allow for a fast evaluation of the bispectrum on squeezed configurations.

With future galaxy surveys, a huge number of Fourier modes of the distribution of the large scale structures in the Universe will become available. These modes are complementary to those of the CMB and can be used to set constraints on models of the early universe, such as inflation. Using a MCMC analysis, we compare the power of the CMB with that of the combination of CMB and galaxy survey data, to constrain the power spectrum of primordial fluctuations generated during inflation. We base our analysis on the Planck satellite and a spectroscopic redshift survey with configuration parameters close to those of the Euclid mission as examples. We first consider models of slow-roll inflation, and show that the inclusion of large scale structure data improves the constraints by nearly halving the error bars on the scalar spectral index and its running. If we attempt to reconstruct the inflationary single-field potential, a similar conclusion can be reached on the parameters characterizing the potential. We then study models with features in the power spectrum. In particular, we consider ringing features produced by a break in the potential and oscillations such as in axion monodromy. Adding large scale structures improves the constraints on features by more than a factor of two. In axion monodromy we show that there are oscillations with small amplitude and frequency in momentum space that are undetected by CMB alone but can be measured by including galaxy surveys in the analysis.

Future lensing surveys will be nearly full-sky and reach an unprecedented depth, probing scales closer and closer to the Hubble radius. This motivates the study of the cosmic shear beyond the small-angle approximation and including general relativistic corrections that are usually suppressed on sub-Hubble scales. The complete expression of the reduced cosmic shear at second order including all relativistic effects was derived in [1]. In the present paper we compute the resulting cosmic shear bispectrum when all these effects are properly taken into account and we compare it to primordial non-Gaussianity of the local type. The new general relativistic effects are generically smaller than the standard non-linear couplings. However, their relative importance increases at small multipoles and for small redshifts of the sources. The dominant effect among these non standard corrections is due to the inhomogeneity of the source redshift. In the squeezed limit, its amplitude can become of the order of the standard couplings when the redshift of the sources is below 0.5. Moreover, while the standard non-linear couplings depend on the angle between the short and long mode, the relativistic corrections do not and overlap almost totally with local type non-Gaussianity. We find that they can contaminate the search for a primordial local signal by f_NL>10.

The CMB bispectrum generated by second-order effects at recombination can be calculated analytically when one of the three modes has a wavelength much longer than the other two and is outside the horizon at recombination. This was pointed out in \cite{Creminelli:2004pv} and here we correct their results. We derive a simple formula for the bispectrum, $f_{NL}^{loc} = - (1/6+ \cos 2 \theta) \cdot (1- 1/2 \cdot d \ln (l_S^2 C_{S})/d \ln l_S)$, where $C_S$ is the short scale spectrum and $\theta$ the relative orientation between the long and the short modes. This formula is exact and takes into account all effects at recombination, including recombination-lensing, but neglects all late-time effects such as ISW-lensing. The induced bispectrum in the squeezed limit is small and will negligibly contaminate the Planck search for a local primordial signal: this will be biased only by $f_{NL}^{loc}\approx-0.4$. The above analytic formula includes the primordial non-Gaussianity of any single-field model. It also represents a consistency check for second-order Boltzmann codes: we find substantial agreement with the CMBquick code.

We study large-scale structure formation in the presence of a quintessence component with zero speed of sound in the framework of Eulerian Perturbation Theory. Due to the absence of pressure gradients, quintessence and dark matter are comoving and can be studied as a unique fluid in terms of the total energy density contrast and the common velocity. In this description the clustering of quintessence enhances the linear term proportional to the velocity divergence in the continuity equation by a factor (1+w) Omega_Q / Omega_m. This is responsible for a rapid evolution of the growth rate at low redshifts, and modifies the standard relation between the velocity divergence and the growth factor. For the total fluid, the solutions for the linear growth function and growth rate can be written in integral forms and admit simple fitting formulae, as in the LambdaCDM case. At second order in perturbation theory, we derive an explicit expression for the kernels F_2 and G_2. They receive modifications of the order of the ratio between quintessence and total energy density perturbations, which affect the corresponding tree-level bispectra. We finally compute the cumulative signal-to-noise in the power spectrum, bispectrum and reduced bispectrum, expected for departures from a LambdaCDM cosmology both in the clustering and smooth quintessence scenarios. The reduced bispectrum, in particular, receives sensible modifications only in the clustering case and can potentially be used to detect or rule out the model.

We give a pedagogical review of a covariant and fully non-perturbative approach to study nonlinear perturbations in cosmology. In the first part, devoted to cosmological fluids, we define a nonlinear extension of the uniform-density curvature perturbation and derive its evolution equation. In the second part, we focus our attention on multiple scalar fields and present a nonlinear description in terms of adiabatic and entropy perturbations. In both cases, we show how the formalism presented here enables one to easily obtain equations up to second, third and higher orders.

We study the spherical collapse model in the presence of quintessence with negligible speed of sound. This case is particularly motivated for w<-1 as it is required by stability. As pressure gradients are negligible, quintessence follows dark matter during the collapse. The spherical overdensity behaves as a separate closed FLRW universe, so that its evolution can be studied exactly. We derive the critical overdensity for collapse and we use the extended Press-Schechter theory to study how the clustering of quintessence affects the dark matter mass function. The effect is dominated by the modification of the linear dark matter growth function. A larger effect occurs on the total mass function, which includes the quintessence overdensities. Indeed, here quintessence constitutes a third component of virialized objects, together with baryons and dark matter, and contributes to the total halo mass by a fraction ~ (1+w) Omega_Q / Omega_m. This gives a distinctive modification of the total mass function at low redshift.

We compute the reduced cosmic shear up to second order in the gravitational potential without relying on the small angle or thin-lens approximation. This is obtained by solving the Sachs equation which describes the deformation of the infinitesimal cross-section of light bundle in the optical limit, and maps galaxy intrinsic shapes into their angular images. The calculation is done in the Poisson gauge without a specific matter content, including vector and tensor perturbations generated at second order and taking account of the inhomogeneities of a fixed redshift source plane. Our final result is expressed in terms of spin-2 operators on the sphere and is valid on the full sky. Beside the well known lens-lens and Born corrections that dominate on small angular scales, we find new non-linear couplings. These are a purely general relativistic intrinsic contribution, a coupling between the gravitational potential at the source with the lens, couplings between the time delay with the lens, couplings between two photon deflections, as well as non-linear couplings due to the second-order vector and tensor components. The inhomogeneity in the redshift of the source induces a coupling between the photon redshift with the lens. All these corrections become important on large angular scales and should thus be included when computing higher-order observables such as the bispectrum, in full or partially full-sky surveys.

We calculate the Cosmic Microwave Background anisotropy bispectrum on large angular scales in the absence of primordial non-Gaussianities, assuming exact matter dominance and extending at second order the classic Sachs-Wolfe result \delta T/T=\Phi/3. The calculation is done in Poisson gauge. Besides intrinsic contributions calculated at last scattering, one must consider integrated effects. These are associated to lensing, and to the time dependence of the potentials (Rees-Sciama) and of the vector and tensor components of the metric generated at second order. The bispectrum is explicitly computed in the flat-sky approximation. It scales as l^(-4) in the scale invariant limit and the shape dependence of its various contributions is represented in 3d plots. Although all the contributions to the bispectrum are parametrically of the same order, the full bispectrum is dominated by lensing. In the squeezed limit it corresponds to f_NL^local = -1/6 - cos(2 \theta), where \theta is the angle between the short and the long modes; the angle dependent contribution comes from lensing. In the equilateral limit it corresponds to f_NL^equil ~ 3.13.

We compute the connected four-point correlation function of the primordial curvature perturbation generated during inflation with standard kinetic terms, where the correlation is established via exchange of a graviton between two pairs of scalar fluctuations. Any such correlation yields a contribution to the scalar trispectrum of the order of the tensor to scalar ratio r. This contribution is numerically one order of magnitude larger than the one previously calculated on the basis of scalar perturbations interacting at a point and satisfies a simple relation in the limit where the momentum of the graviton which is exchanged becomes much smaller than the external momenta. We conclude that the total non-linearity parameter generated by single-field models of slow-roll inflation is at maximum tauNL ~ r.

We compute the non-Gaussianity of the curvature perturbation generated by ekpyrotic collapse with multiple fields. The transition from the multi-field scaling solution to a single-field dominated regime converts initial isocurvature field perturbations to an almost scale-invariant comoving curvature perturbation. In the specific model of two fields, $\phi_1$ and $\phi_2$, with exponential potentials, $-V_i \exp (-c_i \phi_i)$, we calculate the bispectrum of the resulting curvature perturbation. We find that the non-Gaussianity is dominated by non-linear evolution on super-Hubble scales and hence is of the local form. The non-linear parameter of the curvature perturbation is given by $f_{NL} = 5 c_j^2 /12$, where $c_j$ is the exponent of the potential for the field which becomes sub-dominant at late times. Since $c_j^2$ must be large, in order to generate an almost scale invariant spectrum, the non-Gaussianity is inevitably large. By combining the present observational constraints on $f_{\rm NL}$ and the scalar spectral index, the specific model studied in this paper is thus ruled out by current observational data.

We study generic single-field dark energy models, by a parametrization of the most general theory of their perturbations around a given background, including higher derivative terms. In appropriate limits this approach reproduces standard quintessence, k-essence and ghost condensation. We find no general pathology associated to an equation of state w_Q < -1 or in crossing the phantom divide w_Q = -1. Stability requires that the w_Q < -1 side of dark energy behaves, on cosmological scales, as a k-essence fluid with a virtually zero speed of sound. This implies that one should set the speed of sound to zero when comparing with data models with w_Q < -1 or crossing the phantom divide. We summarize the theoretical and stability constraints on the quintessential plane (1+w_Q) vs. speed of sound squared.

We study non-linear primordial adiabatic and isocurvature perturbations and their non-Gaussianity. After giving a general formulation in the context of an extended delta N-formalism, we analyse in detail two illustrative examples. The first is a mixed curvaton-inflaton scenario in which fluctuations of both the inflaton and a curvaton (a light isocurvature field during inflation) contribute to the primordial density perturbation. The second example is that of double inflation involving two decoupled massive scalar fields during inflation. In the mixed curvaton-inflaton scenario we find that the bispectrum of primordial isocurvature perturbations may be large and comparable to the bispectrum of adiabatic curvature perturbations.

We present a covariant formalism for studying nonlinear perturbations of scalar fields. In particular, we consider the case of two scalar fields and introduce the notion of adiabatic and isocurvature covectors. We obtain differential equations governing the evolution of these two covectors, as well as the evolution equation for the covector associated with the curvature perturbation. The form of these equations is very close to the analogous equations obtained in the linear theory, but our equations are fully nonlinear and exact. As an application of our formalism, we expand these equations at second order in the perturbations. On large scales, we obtain a closed system of coupled scalar equations giving the evolution of the second-order adiabatic and entropy perturbations in terms of the first-order perturbations. These equations in general contain a nonlocal term which, however, rapidly decays in an expanding universe.

We study nonlinear cosmological perturbations during the post-inflationary evolution, using the equivalence between a perfect barotropic fluid and a derivatively coupled scalar field with Lagrangian [-(\partial \phi)^2]^[(1+w)/2w]. Since this Lagrangian is just a special case of k-inflation, this approach is analogous to the one employed in the study of non-Gaussianities from inflation. We use this method to derive the second order metric during matter dominance in the comoving gauge directly as function of the primordial inflationary perturbation \zeta. Going to Poisson gauge, we recover the metric previously derived in the literature.

We present an approach to cosmological perturbations based on a covariant perturbative expansion between two worldlines in the real inhomogeneous universe. As an application, at an arbitrary order we define an exact scalar quantity which describes the inhomogeneities in the number of e-folds on uniform density hypersurfaces and which is conserved on all scales for a barotropic ideal fluid. We derive a compact form for its conservation equation at all orders and assign it a simple physical interpretation. To make a comparison with the standard perturbation theory, we develop a method to construct gauge-invariant quantities in a coordinate system at arbitrary order, which we apply to derive the form of the n-th order perturbation in the number of e-folds on uniform density hypersurfaces and its exact evolution equation. On large scales, this provides the gauge-invariant expression for the curvature perturbation on uniform density hypersurfaces and its evolution equation at any order.

We study the bispectrum of the curvature perturbation on uniform energy density hypersurfaces in models of inflation with two scalar fields evolving simultaneously. In the case of a separable potential, it is possible to compute the curvature perturbation up to second order in the perturbations, generated on large scales due to the presence of non-adiabatic perturbations, by employing the $\delta N$-formalism, in the slow-roll approximation. In this case, we provide an analytic formula for the nonlinear parameter $f_{NL}$. We apply this formula to double inflation with two massive fields, showing that it does not generate significant non-Gaussianity; the nonlinear parameter at the end of inflation is slow-roll suppressed. Finally, we develop a numerical method for generic two-field models of inflation, which allows us to go beyond the slow-roll approximation and confirms our analytic results for double inflation.

We give a detailed and improved presentation of our recently proposed formalism for non-linear perturbations in cosmology, based on a covariant and fully non-perturbative approach. We work, in particular, with a covector combining the gradients of the energy density and of the local number of e-folds to obtain a non-linear generalization of the familiar linear uniform density perturbation. We show that this covector obeys a remarkably simple conservation equation which is exact, fully non-linear and valid at all scales. We relate explicitly our approach to the coordinate-based formalisms for linear perturbations and for second-order perturbations. We also consider other quantities, which are conserved on sufficiently large scales for adiabatic perturbations, and discuss the issue of gauge invariance.

We develop a covariant formalism to study nonlinear perturbations of dissipative and interacting relativistic fluids. We derive nonlinear evolution equations for various covectors defined as linear combinations of the spatial gradients of the local number of e-folds and of some scalar quantities characterizing the fluid, such as the energy density or the particle number density. For interacting fluids we decompose perturbations into adiabatic and entropy components and derive their coupled evolution equations, recovering and extending the results obtained in the context of the linear theory. For non-dissipative and noninteracting fluids, these evolution equations reduce to the conservation equations that we have obtained in recent works. We also illustrate geometrically the meaning of the covectors that we have introduced.

We define fully non-perturbative generalizations of the uniform density and comoving curvature perturbations, which are known, in the linear theory, to be conserved on sufficiently large scales for adiabatic perturbations. Our non-linear generalizations are defined geometrically, independently of any coordinate system. We give the equations governing their evolution on all scales. Also, in order to make contact with previous works on first and second order perturbations, we introduce a coordinate system and show that previous results can be recovered, on large scales, in a remarkably simple way, after restricting our definitions to first and second orders in a perturbative expansion.

We study the possibility that large scale cosmological perturbations have been generated during the domination and decay of a massive particle species whose mass depends on the expectation value of a light scalar field. We discuss the constraints that must be imposed on the field in order to remain light and on the annihilation cross section and decay rate of the massive particles in order for the mechanism to be efficient. We compute the resulting curvature perturbations after the mass domination, recovering the results of Dvali, Gruzinov, and Zaldarriaga in the limit of total domination. By comparing the amplitude of perturbations generated by the mass domination to those originally present from inflation, we conclude that this mechanism can be the primary source of perturbations only if inflation does not rely on slow-roll conditions.

We study the quantum fluctuations of scalar fields with a variable effective mass during an inflationary phase. We consider the situation where the effective mass depends on a background scalar field, which evolves during inflation from being frozen into a damped oscillatory phase when the Hubble parameter decreases below its mass. We find power spectra with suppressed amplitude on large scales, similar to the standard massless spectrum on small scales, and affected by modulations on intermediate scales. We stress the analogies and differences with the parametric resonance in the preheating scenario. We also discuss some potentially observable consequences when the scalar field behaves like a curvaton.

We discuss second-order cosmological perturbations on super-Hubble scales, in a scalar field dominated universe, such as during single field inflation. In this contest we show that the gauge-invariant curvature perturbations defined on uniform density and comoving hypersurfaces coincide and that perturbations are adiabatic in the large scale limit. Since it has been recently shown that the uniform density curvature perturbation is conserved on large scales if perturbations are adiabatic, we conclude that both the uniform density and comoving curvature perturbations at second-order, in a scalar field dominated universe, are conserved. Finally, in the light of this result, we comment on the variables recently used in the literature to compute non-Gaussianities.

A recent variant of the inflationary paradigm is that the ``primordial'' curvature perturbations come from quantum fluctuations of a scalar field, subdominant and effectively massless during inflation, called the ``curvaton'', instead of the fluctuations of the inflaton field. We consider the situation where the primordial curvature perturbations generated by the quantum fluctuations of an inflaton and of a curvaton field are of the same order of magnitude. We compute the curvature perturbation and its spectrum in this case and we discuss the observational consequences.

We compare the standard single scalar field inflationary predictions with those of an inflationary phase driven by a tachyon field. A slow-roll formalism is defined for tachyon inflation, and we derive the spectra of scalar and tensor perturbations as well as the consistency relations. At lowest order the predictions of standard and tachyon inflation are shown to be the same. Higher order deviations are present and their observational relevance is discussed. We then study some typical inflationary tachyon potentials, discuss their observational consequences and compare them with recent data. All the models predict a negative and very small running of the scalar spectral index, and they consistently lie within the 1\sigma contour of the data set. However, the regime of blue scalar spectral index and large gravity waves cannot be explored by these models. Finally, a new exact solution of the unperturbed and perturbed coupled gravity and tachyon equations is also presented.

We study the evolution of perturbations during the domination and decay of a massive particle species whose mass and decay rate are allowed to depend on the expectation value of a light scalar field. We specialize in the case where the light field is slow-rolling, showing that during a phase of inhomogeneous mass-domination and decay the isocurvature perturbation of the light field is converted into a curvature perturbation with an efficiency which is nine times larger than when the mass is fixed. We derive a condition on the annihilation cross section and on the decay rate for the domination of the massive particles and we show that standard model particles cannot dominate the universe before nucleosynthesis. We also compare this mechanism with the curvaton model. Finally, observational signatures are discussed. A cold dark matter isocurvature mode can be generated if the dark matter is produced out of equilibrium by both the inflaton and the massive particle species decay. Non-Gaussianities are present: they are chi-square deviations. However, they might be too small to be observable.

We compare the standard single scalar field inflationary predictions with those of an inflationary phase driven by a tachyon field. A slow-roll formalism is defined for tachyon inflation and we derive the spectra of scalar and tensor perturbations as well as the consistency relations. At lowest order the predictions of standard and tachyon inflation are the same. Higher order deviations are present and their observational relevance is discussed. We discuss the observational consequences of some typical inflationary tachyon potentials and compare them with recent data. All the models predict a negative and very small running of the scalar spectral index, and they consistently lie within the 1$\sigma$ contour of the data set. However, the regime of blue scalar spectral index and large gravity waves is not explored by these models.

At low energy, the four-dimensional effective action of the ekpyrotic model of the universe is equivalent to a slightly modified version of the pre big bang model. We discuss cosmological perturbations in these models. In particular we address the issue of matching the perturbations from a collapsing to an expanding phase in full generality. We show that, generically, one obtains $n=0$ for the spectrum of scalar perturbations in the original pre big model (with vanishing potential). When an exponential potential for the dilaton is included, a scale invariant spectrum ($n=1$) of adiabatic scalar perturbations is produced under very generic matching conditions, both in a modified pre big bang and ekpyrotic scenario. We also derive general results valid for power law scale factors matched to a radiation dominated era.

We derive the gauge invariant perturbation equations for a 5-dimensional bulk spacetime in the presence of a brane. The equations are derived in full generality, without specifying a particular energy content of the bulk or the brane. We do not assume Z_2 symmetry, and show that the degree of freedom associated with brane motion plays a crucial role. Our formalism may also be used in the Z_2 symmetric case where it simplifies considerably.

We reexamine Wald's no-hair theorem for global anisotropy in the brane world scenarios. We derive a set of sufficient conditions which must be satisfied by the brane matter and bulk metric so that a homogeneous and anisotropic brane asymptotically evolves to a de Sitter spacetime in the presence of a positive cosmological constant on the brane. We discuss the violations of these sufficient conditions and we show that a negative nonlocal energy density or the presence of strong anisotropic stress (i.e., a magnetic field) may lead the brane to collapse. We discuss the generality of these conditions.

We present an alternative scenario for cosmic structure formation where initial fluctuations are due to Kalb-Ramond axions produced during a pre-big bang phase of inflation. We investigate whether this scenario, where the fluctuations are induced by seeds and therefore are of isocurvature nature, can be brought in agreement with present observations by a suitable choice of cosmological parameters. We also discuss several observational signatures which can distinguish axion seeds from standard inflationary models. We finally discuss the gravitational wave background induced in this model and we show that it may be well within the range of future observations.

We present an alternative scenario for cosmic structure formation where initial fluctuations are due to Kalb-Ramond axions produced during a pre-big bang phase of inflation. We investigate whether this scenario, where the fluctuations are induced by seeds and therefore are of isocurvature nature, can be brought in agreement with present observations by a suitable choice of cosmological parameters. We also discuss several observational signatures which can distinguish axion seeds from standard inflationary models.

A recently proposed mechanism for large-scale structure in string cosmology --based on massless axionic seeds-- is further analyzed and extended to the acoustic-peak region. Existence, structure, and normalization of the peaks turn out to depend crucially on the overall evolution of extra dimensions during the pre-big bang phase: conversely, precise cosmic microwave background anisotropy data in the acoustic-peak region will provide, within the next decade, a window on string-theory's extra dimensions before their eventual compactification.