We study a metric cubic gravity theory considering odd-parity modes of linear inhomogeneous perturbations on a spatially homogeneous Bianchi type I manifold close to the isotropic de Sitter spacetime. We show that in the regime of small anisotropy, the theory possesses new degrees of freedom compared to General Relativity, whose kinetic energy vanishes in the limit of exact isotropy. From the mass dispersion relation we show that such theory always possesses at least one ghost mode as well as a very short-time-scale (compared to the Hubble time) classical tachyonic (or ghost-tachyonic) instability. In order to confirm our analytic analysis, we also solve the equations of motion numerically and we find that this instability is developed well before a single e-fold of the scale factor. This shows that this gravity theory, as it is, cannot be used to construct viable cosmological models.

We study the phenomenology of a class of minimally modified gravity theories called $f(\mathcal{H})$ theories, in which the usual general relativistic Hamiltonian constraint is replaced by a free function of it. After reviewing the construction of the theory and a consistent matter coupling, we analyze the dynamics of cosmology at the levels of both background and perturbations, and present a concrete example of the theory with a $3$-parameter family of the function $f$. Finally, we compare this example model to Planck data as well as some later-time probes, showing that such a realization of $f(\mathcal{H})$ theories fits the data significantly better than the standard $\Lambda$CDM model, in particular by modifying gravity at intermediate redshifts, $z\simeq743$.

We investigate the cosmological viability of the generalized proca theory. We first implement the background and linear perturbation equations of motion in the Boltzmann code and then study the constraints on the parameters of the generalized proca theory after running MCMC against the cosmological data set. With Planck + HST data, we obtain the constraint $h=0.7334_{-0.0269}^{+0.0246}$, which indicates that the tension between early universe and late time universe within this theory is removed. By adding other late-time data sets (BAO, RSD, etc.) we show that the tension is reduced, as the 2$\sigma$ allowed region for $h$ in Proca, $h=0.7041_{-0.0087}^{+0.0094}$, overlaps with the 2$\sigma$ region of the HST data.

We provide two derivations of the baryonic equations that can be straightforwardly implemented in existing Einstein--Boltzmann solvers. One of the derivations begins with an action principle, while the other exploits the conservation of the stress-energy tensor. While our result is manifestly covariant and satisfies the Bianchi identities, we point out that this is not the case for the implementation of the seminal work by Ma and Bertschinger and in the existing Boltzmann codes. We also study the tight coupling approximation up to the second order without choosing any gauge using the covariant full baryon equations. We implement the improved baryon equations in a Boltzmann code and investigate the change in the estimate of cosmological parameters by performing an MCMC analysis. With the covariantly correct baryon equations of motion, we find 1% deviation for the best fit values of the cosmological parameters that should be taken into account. While in this paper, we study the Lambda-CDM model only, our baryon equations can be easily implemented in other models and various modified gravity theories.