The minimal theory of quasidilaton massive gravity with or without a Horndeski-type kinetic term for the quasidilaton field propagates only three physical modes: the two massive tensor polarizations and one scalar mode. This reduced number of degrees of freedom is realized by a Lorentz symmetry violation at cosmological scales and the presence of appropriate constraints that remove unwanted modes. Vacuum cosmological solutions have been considered in a previous work, and it has been shown that the late-time de Sitter attractor is stable under inhomogeneous perturbations. In this work, we explore the stability of cosmological solutions in the presence of matter fields. Assuming for simplicity that the quasidilaton scalar is on an attractor at the level of the background, we derive stability conditions in the subhorizon limit, and find the scalar sound speeds, as well as the modification with respect to general relativity to the gravitational potential in the quasistatic approximation. We also find that the speed limit of gravitational waves coincides with the speed of light for any homogeneous and isotropic cosmological background, on or away from the attractor.
I. INTRODUCTIONOne of the most important questions in modern cosmology remains the elucidation of the origin of the current accelerated expansion of the Universe [1]. One of the studied avenues to address this so-called dark energy puzzle lies in large-distance, or infra-red (IR), modifications to general relativity (GR). For example, gravity could be weaker at distances larger than the separation between galaxy clusters so that a large cosmological constant would not gravitate as strongly as in GR and would lead to a modest acceleration as observed. Alternatively, in the absence of a cosmological constant an IR modified behavior of gravity could dynamically source the acceleration of the expansion.Among the IR modifications to Einstein gravity, adding a mass to the gravitational field is often presented as one of the simplest options, at least conceptually. In practice, it has however been difficult to construct a theory of massive gravity without compromising basic assumptions. In particular, the de Rham Gabadadze Tolley (dRGT) massive gravity [2], unique healthy non-linear theory of pure massive gravity in four dimensions, does not accommodate Friedmann Lemaître Robertson Walker (FLRW) solutions, meaning that either isotropy or homogeneity at large scales need to be considered as approximate [3]. Although this does not compromise the validity of the theory as long as the breaking of the approximate symmetry occurs at large enough distances [4], it effectively renders the analysis more tenuous. In order to evade the practical issues, several extensions of dRGT massive gravity have thus since then been proposed, for instance: massive bigravity [5], mass-varying massive gravity [6], modified matter couplings [7], scalar field extensions [8][9][10][11][12], and Lorentz symmetry violating extensions [13][14][15][16]. This has proven to be a fruitful path: se...