In the context of Finsler-Randers theory we consider, for a first time, the cosmological scenario of the varying vacuum. In particular, we assume the existence of a cosmological fluid source described by an ideal fluid and the varying vacuum terms. We determine the cosmological history of this model by performing a detailed study on the dynamics of the field equations. We determine the limit of General Relativity, while we find new eras in the cosmological history provided by the geometrodynamical terms provided by the Finsler-Randers theory.

We perform for the first time a dynamical system analysis of both the background and perturbation equations, of $\Lambda$CDM cosmology and quintessence scenario with an exponential potential. In the former case the perturbations do not change the stability of the late-time attractor of the background equations, and the system still results in the dark-energy dominated, de Sitter solution, having passed from the correct dark-matter era with $\gamma\approx6/11$. However, in the case of quintessence the incorporation of perturbations changes the stability and properties of the background evolution, and the only conditionally stable points present either an exponentially increasing matter clustering not favored by observations, or Laplacian instabilities, and thus not physically interesting. This result is a severe disadvantage of quintessence cosmology comparing to $\Lambda$CDM paradigm.

We study the cosmological evolution of the field equations in the context of Einstein-Aether cosmology by including a scalar field in a spatially flat Friedmann--Lema\^{\i}tre--Robertson--Walker spacetime. Our analysis is separated into two separate where a pressureless fluid source is included or absent. In particular, we determine the critical points of the field equations and we study the stability of the specific solutions. The limit of general relativity is fully recovered, while the dynamical system admits de Sitter solutions which can describe the past inflationary era and the future late-time attractor. Results for generic scalar field potentials are presented while some numerical behaviours are given for specific potential forms.

In the context of $f\left( R\right) $-gravity with a spatially flat FLRW metric containing an ideal fluid, we use the method of invariant transformations to specify families of models which are integrable. We find three families of $f(R)$ theories for which new analytical solutions are given and closed-form solutions are provided.

We study for the first time the dynamical properties and the growth index of linear matter perturbations of the Finsler-Randers (FR) cosmological model, for which we consider that the cosmic fluid contains matter, radiation and a scalar field. Initially, for various FR scenarios we implement a critical point analysis and we find solutions which provide cosmic acceleration and under certain circumstances we can have de-Sitter points as stable late-time attractors. Then we derive the growth index of matter fluctuations in various Finsler-Randers cosmologies. Considering cold dark matter and neglecting the scalar field component from the perturbation analysis we find that the asymptotic value of the growth index is $\gamma_{\infty}^{(FR)}\approx\frac {9}{16}$, which is close to that of the concordance $\Lambda$ cosmology, $\gamma^{(\Lambda)} \approx\frac{6}{11}$. In this context, we show that the current FR model provides the same Hubble expansion with that of Dvali, Gabadadze and Porrati (DGP) gravity model. However, the two models can be distinguished at the perturbation level since the growth index of FR model is $\sim18.2\%$ lower than that of the DPG gravity $\gamma^{(DGP)} \approx \frac{11}{16}$. If we allow pressure in the matter fluid then we obtain $\gamma_{\infty}^{(FR)}\approx\frac{9(1+w_{m})(1+2w_{m})}{2[8+3w_{m}% (5+3w_{m})]}$, where $w_{m}$ is the matter equation of state parameter. Finally, we extend the growth index analysis by using the scalar field and we find that the evolution of the growth index in FR cosmologies is affected by the presence of scalar field.

In the context of generalised Brans-Dicke cosmology we use the Killing tensors of the minisuperspace in order to determine the unspecified potential of a scalar-tensor gravity theory. Specifically, based on the existence of contact symmetries of the field equations, we find four types of potentials which provide exactly integrable dynamical systems. We investigate the dynamical properties of these potentials by using a critical point analysis and we find solutions which lead to cosmic acceleration and under specific conditions we can have de-Sitter points as stable late-time attractors.