We provide further evidence that a massless cosmological scalar field with a non-minimal coupling to the Ricci curvature of the type $M^2_{\rm pl}(1+\xi \sigma^n/M_{\rm pl}^n) $ alleviates the existing tension between local measurements of the Hubble constant and its inference from CMB anisotropies and baryonic acoustic oscillations data in presence of a cosmological constant. In these models, the expansion history is modified compared to $\Lambda$CDM at early time, mimicking a change in the effective number of relativistic species, and gravity weakens after matter-radiation equality. Compared to $\Lambda$CDM, a quadratic ($n=2$) coupling increases the Hubble constant when {\em Planck} 2018 (alone or in combination with BAO and SH0ES) measurements data are used in the analysis. Negative values of the coupling, for which the scalar field decreases, seem favored and consistency with Solar System can be naturally achieved for a large portion of the parameter space without the need of any screening mechanism. We show that our results are robust to the choice of $n$, also presenting the analysis for $n=4$.

We propose new massive gravity theories with 5 dynamical degrees of freedom. We evade uniqueness theorems regarding the form of the kinetic and potential terms by adopting the "generalized massive gravity" framework, where a global translation invariance is broken. By exploiting the rotation symmetry in the field space, we determine two novel classes of theories. The first one is an extension of generalized massive gravity with a non-minimal coupling. On the other hand, the second theory produces a mass term that is different from de Rham, Gabadadze, Tolley construction and trivially has 5 degrees of freedom. Both theories allows for stable cosmological solutions without infinite strong coupling, which are free of ghost and gradient instabilities.

One of the most appealing approaches to ease the Hubble tension is the inclusion of an early dark energy (EDE) component that adds energy to the Universe in a narrow redshift window around the time of recombination and dilutes faster than radiation afterwards. In this paper, we analyze EDE in the framework of $\alpha$-attractor models. As well known, the success in alleviating the Hubble tension crucially depends on the shape of the energy injection. We show how different types of energy injections can be easily obtained, thanks to the freedom in choosing the functional form of the potential inspired by $\alpha$-attractor models. To confirm our intuition we perform an MCMC analysis for three representative cases and find indeed that $H_0$ is significantly larger than in $\Lambda$CDM like in other EDE models. Unlike axion-driven EDE models with super Planckian decay constant, the curvature of EDE models required by the data is natural in the context of recent theoretical developments in $\alpha$-attractors.

We present a thorough analysis of the cosmological perturbations in Generalised Massive Gravity. This is an extension of de Rham-Gabadadze-Tolley theory where the translation invariance in the St\"uckelberg field space is broken. This allows the mass parameters to be promoted to functions of the St\"uckelberg fields. We consider an exact cosmological background in this theory and study the stability of perturbations. We derive conditions to avoid ghost, gradient and tachyonic instability. The cosmology is an extension of the self-accelerating branch of the constant mass parameter theory, but now all five massive graviton polarisations propagate. For concreteness, we consider a minimal version of the theory where cosmology undergoes an accelerated expansion at late times and show that the perturbative stability is preserved for a range of parameters.

For field theories in curved spacetime, defining how matter gravitates is part of the theory building process. In this letter, we adopt Bekenstein's multiple geometries approach to allow part of the matter sector to follow the geodesics on a general pseudo-Riemannian geometry, constructed from a tensor and a $U(1)$ gauge field. This procedure allows us to generate a previously unknown corner of vector-tensor theories. In the Jordan frame, apparent high-derivative terms of the vector field are reduced by integrating out an auxiliary variable, at the cost of introducing new matter interactions. As a simple example, we consider a conformal relation between different geometries and demonstrate the presence of an auxiliary degree. We conclude with a discussion of applications, in particular for the early universe.

We revisit the question of viability of bigravity cosmology as a candidate for dark energy. In the context of the low energy limit model, where matter couples to a single metric, we study linear perturbations around homogeneous and isotropic backgrounds to derive the Poisson's equation for the Newtonian potential. Extending to second order perturbations, we identify the Vainshtein radius below which non-linear scalar self interactions conspire to reproduce GR on local scales. We combine all of these results to determine the parameter space that allows a late time de-Sitter attractor compatible with observations and a successful Vainsthein mechanism. We find that the requirement on having a successful Vainsthein mechanism is not compatible with the existence of cosmological solutions at early times.

The cuscuton is a scalar field with infinite speed of propagation. It was introduced in the context of cosmology but it has also been claimed to resemble Horava gravity in a certain limit. Here we revisit the cuscuton theory as a Lorentz-violating gravity theory. We clarify its relation with Horava gravity and Einstein-aether theory, analyze its causal structure, and consider its initial value formulation. Finally, we discuss to which extent the cuscuton theory can be used as a proxy for Horava gravity in the context of gravitational collapse and formation of universal horizons.

The recent detection of gravitational wave GW170817 has placed a severe bound on the deviation of the speed of gravitational waves from the speed of light. We explore the consequences of this detection for Horava gravity.

We present a novel cosmological solution in the framework of extended quasidilaton theory which underwent scrutiny recently. We only consider terms that do not generate the Boulware-Deser degree of freedom, hence the "ghost-free" quasidilaton theory, and show three new branches of cosmological evolution therein. One of the solutions passes the perturbative stability tests. This new solution exhibits a late time self-acceleration and all graviton polarizations acquire masses that converge to a constant in the asymptotic future. Moreover, all modes propagate at the speed of light. We propose that this solution can be used as a benchmark model for future phenomenological studies.

The extended quasidilaton theory is one of the simplest Lorentz-invariant massive gravity theories which can accommodate a stable self-accelerating vacuum solution. In this paper we revisit this theory and study the effect of matter fields. For a matter sector that couples minimally to the physical metric, we find hints of a Jeans type instability in the IR. In the analogue k-essence field set-up, this instability manifests itself as an IR ghost for the scalar field perturbation, but this can be interpreted as a classical instability that becomes relevant below some momentum scale in terms of matter density perturbations. We also consider the effect of the background evolution influenced by matter on the stability of the gravity sector perturbations. In particular, we address the previous claims of ghost instability in the IR around the late time attractor. We show that, although the matter-induced modification of the evolution potentially brings tension to the stability conditions, one goes beyond the regime of validity of the effective theory well before the solutions become unstable. We also draw attention to the fact that the IR stability conditions are also enforced by the existence requirements of consistent background solutions.

We study the cosmology of bimetric theory with a composite matter coupling. We find two possible branches of background evolution. We investigate the question of stability of cosmological perturbations. For the tensor and vector perturbations, we derive conditions on the absence of ghost and gradient instabilities. For the scalar modes, we obtain conditions for avoiding ghost degrees. In the first branch, we find that one of the scalar modes becomes a ghost at the late stages of the evolution. Conversely, this problem can be avoided in the second branch. However, we also find that the constraint for the second branch prevents the doubly coupled matter fields from being the standard ingredients of cosmology. We thus conclude that a realistic and stable cosmological model requires additional minimally coupled matter fields.

We study cosmological perturbations in bimetric theory with two fluids each of which is coupled to one of the two metrics. Focusing on a healthy branch of background solutions, we clarify the stability of the cosmological perturbations. For this purpose, we extend the condition for the absence of the so-called Higuchi ghost, and show that the condition is guaranteed to be satisfied on the healthy branch. We also calculate the squared propagation speeds of perturbations and derive the conditions for the absence of the gradient instability. To avoid the gradient instability, we find that the model parameters are weakly constrained.

Recently the first example of a unitary theory of Lorentz-invariant massive gravity allowing for stable self-accelerating de Sitter solutions was found, extending the quasidilaton theory. In this paper we further generalize this new action for the quasidilaton field by introducing general Lagrangian terms which are consistent with the quasidilaton symmetry while leading to second order equations of motion. We find that the structure of the theory, compared to the simplest stable example, does not change on introducing these new terms.

We present a review of cosmological solutions in non-linear massive gravity, focusing on the stability of perturbations. Although homogeneous and isotropic solutions have been found, these are now known to suffer from either Higuchi ghost or a new non-linear ghost instability. We discuss two approaches to alleviate this issue. By relaxing the symmetry of the background by e.g. breaking isotropy in the hidden sector, it is possible to accommodate a stable cosmological solution. Alternatively, extending the theory to allow for new dynamical degrees of freedom can also remove the conditions which lead to the instability. As examples for this case, we study the stability of self-accelerating solutions in the quasi-dilatonic extension and generic cosmological solutions in the varying mass extension. While the quasi-dilaton case turns out to be unstable, the varying mass case allows stable regimes of parameters. Viable self-accelerating solutions in the varying mass theory yet remain to be found.

We investigate nonlinear stability of two classes of cosmological solutions in massive gravity: isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions and anisotropic FLRW solutions. For this purpose we construct the linear cosmological perturbation theory around axisymmetric Bianchi type--I backgrounds. We then expand the background around the two classes of solutions, which are fixed points of the background evolution equation, and analyze linear perturbations on top of it. This provides a consistent truncation of nonlinear perturbations around these fixed point solutions and allows us to analyze nonlinear stability in a simple way. In particular, it is shown that isotropic FLRW solutions exhibit nonlinear ghost instability. On the other hand, anisotropic FLRW solutions are shown to be ghost-free for a range of parameters and initial conditions.

We study cosmological perturbations around self-accelerating solutions to two extensions of nonlinear massive gravity: the quasi-dilaton theory and the mass-varying theory. We examine stability of the cosmological solutions, and the extent to which the vanishing of the kinetic terms for scalar and vector perturbations of self-accelerating solutions in massive gravity is generic when the theory is extended. We find that these kinetic terms are in general non-vanishing in both extensions, though there are constraints on the parameters and background evolution from demanding that they have the correct sign. In particular, the self-accelerating solutions of the quasi-dilaton theory are always unstable to scalar perturbations with wavelength shorter than the Hubble length.

We discuss the detectability of gravitational waves with a time dependent mass contribution, by means of the stochastic gravitational wave observations. Such a mass term typically arises in the cosmological solutions of massive gravity theories. We conduct the analysis based on a general quadratic action, and thus the results apply universally to any massive gravity theories in which modification of general relativity appears primarily in the tensor modes. The primary manifestation of the modification in the gravitational wave spectrum is a sharp peak. The position and height of the peak carry information on the present value of the mass term, as well as the duration of the inflationary stage. We also discuss the detectability of such a gravitational wave signal using the future-planned gravitational wave observatories.

We argue that all homogeneous and isotropic solutions in nonlinear massive gravity are unstable. For this purpose, we study the propagating modes on a Bianchi type--I manifold. We analyze their kinetic terms and dispersion relations as the background manifold approaches the homogeneous and isotropic limit. We show that in this limit, at least one ghost always exists and that its frequency tends to vanish for large scales, meaning that it cannot be integrated out from the low energy effective theory. This ghost mode is interpreted as a leading nonlinear perturbation around a homogeneous and isotropic background.

In the scope of the nonlinear massive gravity, we study fixed points of evolution equations for a Bianchi type--I universe. We find a new attractor solution with non-vanishing anisotropy, on which the physical metric is isotropic but the Stuckelberg configuration is anisotropic. As a result, at the background level, the solution describes a homogeneous and isotropic universe, while a statistical anisotropy is expected from perturbations, suppressed by smallness of the graviton mass.

We present a fully nonlinear study of long wavelength cosmological perturbations within the framework of the projectable Horava-Lifshitz gravity, coupled to a single scalar field. Adopting the gradient expansion technique, we explicitly integrate the dynamical equations up to any order of the expansion, then restrict the integration constants by imposing the momentum constraint. While the gradient expansion relies on the long wavelength approximation, amplitudes of perturbations do not have to be small. When the $\lambda\to 1$ limit is taken, the obtained nonlinear solutions exhibit a continuous behavior at any order of the gradient expansion, recovering general relativity in the presence of a scalar field and the "dark matter as an integration constant". This is in sharp contrast to the results in the literature based on the "standard" (and naive) perturbative approach where in the same limit, the perturbative expansion of the action breaks down and the scalar graviton mode appears to be strongly coupled. We carry out a detailed analysis on the source of these apparent pathologies and determine that they originate from an improper application of the perturbative approximation in the momentum constraint. We also show that there is a new branch of solutions, valid in the regime where $|\lambda-1|$ is smaller than the order of perturbations. In the limit $\lambda\to1$, this new branch allows the theory to be continuously connected to general relativity (plus "dark matter").

We study cosmological perturbations of self-accelerating universe solutions in the recently proposed nonlinear theory of massive gravity, with general matter content. While the broken diffeomorphism invariance implies that there generically are 2 tensor, 2 vector and 2 scalar degrees of freedom in the gravity sector, we find that the scalar and vector degrees have vanishing kinetic terms and nonzero mass terms. Depending on their nonlinear behavior, this indicates either nondynamical nature of these degrees or strong couplings. Assuming the former, we integrate out the 2 vector and 2 scalar degrees of freedom. We then find that in the scalar and vector sectors, gauge-invariant variables constructed from metric and matter perturbations have exactly the same quadratic action as in general relativity. The difference from general relativity arises only in the tensor sector, where the graviton mass modifies the dispersion relation of gravitational waves, with a time-dependent effective mass. This may lead to modification of stochastic gravitational wave spectrum.

In the context of a recently proposed nonlinear massive gravity with Lorentz-invariant mass terms, we investigate open Friedmann-Robertson-Walker (FRW) universes driven by arbitrary matter source. While the flat FRW solutions were recently shown to be absent, the proof does not extend to the open universes. We find three independent branches of solutions to the equations of motion for the St\"uckelberg scalars. One of the branches does not allow any nontrivial FRW cosmologies, as in the previous no-go result. On the other hand, both of the other two branches allow general open FRW universes governed by the Friedmann equation with the matter source, the standard curvature term and an effective cosmological constant $\Lambda_{\pm}=c_{\pm}m_g^2$. Here, $m_g$ is the graviton mass, +and - represent the two branches, and $c_{\pm}$ are constants determined by the two dimensionless parameters of the theory. Since an open FRW universe with a sufficiently small curvature constant can approximate a flat FRW universe but there is no exactly flat FRW solution, the theory exhibits a discontinuity at the flat FRW limit.

In the framework of the power-counting renormalizable theory of gravitation, recently proposed by Ho\v{r}ava, we study the limit $\lambda\to\infty$, which is arguably the most natural candidate for the ultraviolet fixed point of the renormalization group flow. In the projectable version with dynamical critical exponent $z=3$, we analyze the Friedmann-Robertson-Walker background driven by the so-called "dark matter as integration constant" and perturbations around it. We show that amplitudes of quantum fluctuations for both scalar and tensor gravitons remain finite in the limit and that the theory is weakly coupled under a certain condition.

We compute the phenomenological signatures of a model (Watanabe et al' 09) of anisotropic inflation driven by a scalar and a vector field. The action for the vector is U(1) invariant, and the model is free of ghost instabilities. A suitable coupling of the scalar to the kinetic term of the vector allows for a slow roll evolution of the vector vev, and hence for a prolonged anisotropic expansion; this provides a counter example to the cosmic no hair conjecture. We compute the nonvanishing two point correlation functions between physical modes of the system, and express them in terms of power spectra with angular dependence. The anisotropy parameter g_* for the scalar-scalar spectrum (defined as in the Ackerman et al '07 parametrization) turns out to be negative in the simplest realization of the model, which, therefore, cannot account for the angular dependence emerged in some analyses of the WMAP data. A g_* of order -0.1 is achieved when the energy of the vector is about 6-7 orders of magnitude smaller than that of the scalar during inflation. For such values of the parameters, the scalar-tensor correlation (which is in principle a distinctive signature of anisotropic spaces) is smaller than the tensor-tensor correlation.

The Minimal Supersymmetric Standard Model has several flat directions, which can naturally be excited during inflation. If they have a slow (perturbative) decay, they may affect the thermalization of the inflaton decay products. In the present paper, we consider the system of udd and QLd flat directions, which breaks the U(1)xSU(2)xSU(3) symmetry completely. In the unitary gauge and assuming a general soft breaking mass configuration, we show that for a range of parameters, the background condensate of flat directions can undergo a fast non-perturbative decay, due to non-adiabatic evolution of the eigenstates. We find that both the background evolution and part of the decay can be described accurately by previously studied gauged toy models of flat direction decay.

We show that expanding or contracting Kasner universes are unstable due to the amplification of gravitational waves (GW). As an application of this general relativity effect, we consider a pre-inflationary anisotropic geometry characterized by a Kasner-like expansion, which is driven dynamically towards inflation by a scalar field. We investigate the evolution of linear metric fluctuations around this background, and calculate the amplification of the long-wavelength GW of a certain polarization during the anisotropic expansion (this effect is absent for another GW polarization, and for scalar fluctuations). These GW are superimposed to the usual tensor modes of quantum origin from inflation, and are potentially observable if the total number of inflationary e-folds exceeds the minimum required to homogenize the observable universe only by a small margin. Their contribution to the temperature anisotropy angular power spectrum decreases with the multipole l as l^(-p), where p depends on the slope of the initial GW power-spectrum. Constraints on the long-wavelength GW can be translated into limits on the total duration of inflation and the initial GW amplitude. The instability of classical GW (and zero-vacuum fluctuations of gravitons) during Kasner-like expansion (or contraction) may have other interesting applications. In particular, if GW become non-linear, they can significantly alter the geometry before the onset of inflation.

We compute the nonperturbative decay of supersymmetric flat directions due to their D-term potential. Flat directions can develop large vacuum expectation values (vevs) during inflation, and, if they are long-lived, this can strongly affect the reheating and thermalization stages after the inflation. We study a generic system of two U(1) or SU(2) flat directions which are cosmologically evolving after inflation. After proper gauge fixing, we show that the excitations of the fields around this background can undergo exponential amplification, at the expense of the energy density of the flat directions. We compute this effect for several values of the masses and the initial vevs of the two flat directions, through a combination of analytical methods and extensive numerical simulations. For a wide range of parameters the flat directions decay within their first few rotations.

We extend the standard theory of cosmological perturbations to homogeneous but anisotropic universes. We present an exhaustive computation for the case of a Bianchi I model, with a residual isotropy between two spatial dimensions, which is undergoing complete isotropization at the onset of inflation; we also show how the computation can be further extended to more general backgrounds. In presence of a single inflaton field, there are three physical perturbations (precisely as in the isotropic case), which are obtained (i) by removing gauge and nondynamical degrees of freedom, and (ii) by finding the combinations of the remaining modes in terms of which the quadratic action of the perturbations is canonical. The three perturbations, which later in the isotropic regime become a scalar mode and two tensor polarizations (gravitational wave), are coupled to each other already at the linearized level during the anisotropic phase. This generates nonvanishing correlations between different modes of the CMB anisotropies, which can be particularly relevant at large scales (and, potentially, be related to the large scale anomalies in the WMAP data). As an example, we compute the spectrum of the perturbations in this Bianchi I geometry, assuming that the inflaton is in a slow roll regime also in the anisotropic phase. For this simple set-up, fixing the initial conditions for the perturbations appears more difficult than in the standard case, and additional assumptions seem to be needed to provide predictions for the CMB anisotropies.

Most of the analysis of the Cosmic Microwave Background relies on the assumption of statistical isotropy. However, given some recent evidence pointing against isotropy, as for instance the observed alignment of different multipoles on large scales, it is worth testing this assumption against the increasing amount of available data. As a pivot model, we assume that the spectrum of the primordial perturbations depends also on their directionality (rather than just on the magnitude of their momentum, as in the standard case). We explicitly compute the correlation matrix for the temperature anisotropies in the simpler case in which there is a residual isotropy between two spatial directions. As a concrete example, we consider a different initial expansion rate along one direction, and the following isotropization which takes place during inflation. Depending on the amount of inflation, this can lead to broken statistical isotropy on the largest observable scales.