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.

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.

Theories of gravity in the beyond Horndeski class recover the predictions of general relativity in the solar system whilst admitting novel cosmologies, including late-time de Sitter solutions in the absence of a cosmological constant. Deviations from Newton's law are predicted inside astrophysical bodies, which allow for falsifiable, smoking-gun tests of the theory. In this work we study the pulsations of stars by deriving and solving the wave equation governing linear adiabatic oscillations to find the modified period of pulsation. Using both semi-analytic and numerical models, we perform a preliminary survey of the stellar zoo in an attempt to identify the best candidate objects for testing the theory. Brown dwarfs and Cepheid stars are found to be particularly sensitive objects and we discuss the possibility of using both to test the theory.