We use the Equation of State (EoS) approach to study the evolution of the dark sector in Horndeski models, the most general scalar-tensor theories with second order equations of motion. By including the effects of the dark sector into our code EoS\_class, we demonstrate the numerical stability of the formalism and excellent agreement with results from other publicly available codes for a range of parameters describing the evolution of the function characterising the perturbations for Horndeski models, $\alpha_{\rm x}$, with ${\rm x}=\{{\rm K}, {\rm B}, {\rm M}, {\rm T}\}$. After demonstrating that on sub-horizon scales ($k\gtrsim 10^{-3}~{\rm Mpc}^{-1}$ at $z=0$) velocity perturbations in both the matter and the dark sector are typically subdominant with respect to density perturbations in the equation of state for perturbations, we find an attractor solution for the dark sector gauge-invariant density perturbation $\Delta_{\rm ds}$ by neglecting its time derivatives in the equation describing its time evolution, as commonly done in the well-known quasi-static approximation. Using this result, we provide simplified expressions for the equation-of-state functions: the dark sector entropy perturbations $w_{\rm ds}\Gamma_{\rm ds}$ and anisotropic stress $w_{\rm ds}\Pi_{\rm ds}$. From this we derive a growth factor-like equation for both matter and dark sector and are able to capture the relevant physics for several observables with great accuracy. We finally present new analytical expressions for the well-known modified gravity phenomenological functions $\mu$, $\eta$ and $\Sigma$ for a generic Horndeski model as functions of $\alpha_{\rm x}$. We show that on small scales they reproduce expressions presented in previous works, but on large scales, we find differences with respect to other works.

We investigate generalized Einstein-Aether theories that are compatible with the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation, and lensing data. For a given dark energy equation of state, $w_{\rm de}$, we formulate a designer approach and we investigate their impact on the CMB temperature anisotropy and matter power spectra. We use the Equation of State approach to parametrize the perturbations and find that this approach is particularly useful in identifying the most suitable and numerically efficient parameters to explore in a Markov chain Monte Carlo (MCMC) analysis. We find the data constrains models with $w_{\rm de} = -1$ to be compatible with $\Lambda$CDM. For $w_{\rm de} \not= -1$ models, which avoid the gravitational waves constraint through the entropy perturbation, we constrain $w_{\rm de}$ to be $w_{\rm de } = -1.06^{+0.08}_{-0.03}$ (CMB) and $w_{\rm de } =-1.04^{+0.05}_{-0.02}$ (CMB+Lensing) at $68\%$C.L., and find that these models can be different from $\Lambda$CDM and still be compatible with the data. We also find that these models can ameliorate some anomalies in $\Lambda$CDM when confronted with data, such as the low-$\ell$ and high-$k$ power in the CMB temperature anisotropy and matter power spectra respectively, but not simultaneously. We also investigate the anomalous lensing amplitude, quantified by $A_{\rm lens}$, and find that that for $w_{\rm de} = -1$ models, $A_{\rm lens} = 1.15^{+0.07}_{-0.08}$ (CMB) and $A_{\rm lens} = 1.12\pm0.05$ (CMB+Lensing) at $68\%$C.L. $\sim$ 2$\sigma$ larger than expected, similar to previous analyses of $\Lambda$CDM.

We explore the constraints on dark sector models imposed by the recent observation of coincident gravitational waves and gamma rays from a binary neutron star merger, GW170817. Rather than focusing on specific models as has been considered by other authors, we explore this in the context of the equation of state approach of which the specific models are special cases. After confirming the strong constraints found by others for Horndeski, Einstein-Aether and massive gravity models, we discuss how it is possible to construct models which might evade the constraints from GW170817 but still leading to cosmologically interesting modifications to gravity. Possible examples are ``miracle cancellations" such as in $f(R)$ models, nonlocal models and higher-order derivatives. The latter two rely on the dimensionless ratio of the wave number of the observed gravitational waves to the Hubble expansion rate being very large ($\sim10^{19}$) which is used to suppress modifications to the speed of gravitational waves.

We investigate the effects of Population III stars on the sky-averaged 21-cm background radiation, which traces the collective emission from all sources of ultraviolet and X-ray photons before reionization is complete. While UV photons from PopIII stars can in principle shift the onset of radiative coupling of the 21-cm transition -- and potentially reionization -- to early times, we find that the remnants of PopIII stars are likely to have a more discernible impact on the 21-cm signal than PopIII stars themselves. The X-rays from such sources preferentially heat the IGM at early times, which elongates the epoch of reheating and results in a more gradual transition from an absorption signal to emission. This gradual heating gives rise to broad, asymmetric wings in the absorption signal, which stand in contrast to the relatively sharp, symmetric signals that arise in models treating PopII sources only. A stronger signature of PopIII, in which the position of the absorption minimum becomes inconsistent with PopII-only models, requires extreme star-forming events that may not be physically plausible, lending further credence to predictions of relatively high frequency absorption troughs, $\nu_{\min} \sim 100$ MHz. As a result, though the trough location alone may not be enough to indicate the presence of PopIII, the asymmetric wings should arise even if only a few PopIII stars form in each halo before the transition to PopII star formation occurs, provided that the PopIII IMF is sufficiently top-heavy and at least some PopIII stars form in binaries.

We investigate the evolution of cosmological perturbations in models of dark energy described by a time-like unit normalized vector field specified by a general function $\mathcal{F}(\mathcal{K})$, so-called Generalized Einstein-Aether models. First we study the background dynamics of such models via a designer approach in an attempt to model this theory as dark energy. We find that only one specific form of this designer approach matches $\Lambda$CDM at background order and we also obtain a differential equation which $\mathcal{F}(\mathcal{K})$ must satisfy for general $w$CDM cosmologies. We also present the equations of state for perturbations in Generalized Einstein-Aether models, which completely parametrize these models at the level of linear perturbations. A generic feature of modified gravity models is that they introduce new degrees of freedom. By fully eliminating these we are able to express the gauge invariant entropy perturbation and the scalar, vector, and tensor anisotropic stresses in terms of the perturbed fluid variables and metric perturbations only. These can then be used to study the evolution of perturbations in the scalar, vector, and tensor sectors and we use these to evolve the Newtonian gravitational potentials.