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.

In the context of Einstein-aether scalar field cosmology we solve the field equations and determine exact and analytic solutions. In particular, we consider a model proposed by Kanno and Soda where the aether and the scalar fields interact through the aether coefficient parameters, which are promoted to be functions of the scalar field. For this model, we write the field equations by using the minisuperspace approach and we determine the scalar field potentials which leads to Liouville--integrable systems. We solve the field equations for five families of scalar field potentials and, whether it is feasible, we write down the analytic solutions by using closed-form functions.

In the present work we perform a systematic analysis of a new dark energy parametrization and its various corrections at first and higher orders around the presence epoch $z=0$, where the first order correction of this dark energy parametrization recovers the known Chevallier-Polarski-Linder model. We have considered up to the third order corrections of this parametrization and investigate the models at the level of background and perturbations. The models have been constrained using the latest astronomical datasets from a series of potential astronomical data, such as the cosmic microwave background observations, baryon acoustic oscillations measurements, recent Pantheon sample of the supernova type Ia and the Hubble parameter measurements. From the analyses we found that all parametrization favor the quintessential character of the dark energy equation of state where the phantom crossing is marginally allowed (within 68\% CL). Finally, we perform the Bayesian analysis using \texttt{MCEvidence} to quantify the statistical deviations of the parametrizations compared to the standard $\Lambda$CDM cosmology. The Bayesian analysis reports that $\Lambda $CDM is favored over all the DE parametrizations.

We perform a detailed analysis for the dynamics of Chiral cosmology in a spatially flat Friedmann-Lema\^{\i}tre-Robertson-Walker universe with a mixed potential term. The stationary points are categorized in four families. Previous results in the literature are recovered while new phases in the cosmological evolution are found. From our analysis we find nine different cosmological solutions, the eight describe scaling solutions, where the one is that of a pressureless fluid, while only one de Sitter solution is recovered.

Interacting dark energy models are widely renowned for giving an explanation to the cosmic coincidence problem as well as several observational issues. According to the recent observational data, and so far we are concerned with the literature, the choice of the interaction function between dark matter and dark energy is always questionable since there is no such underlying theory that could derive it. Thus, in this work we have raised this issue by proposing two new nonlinear interaction functions and constrain them using cosmic microwave background (CMB) from Planck 2018, baryon acoustic oscillations (BAO), dark energy survey and a measurement of the Hubble constant $H_0$ from Hubble Space Telescope (HST) 2019. The dark energy equation of state is considered to be constant throughout the work and the geometry of the universe is assumed to be homogeneous and isotropic with zero spatial curvature. Our analyses report that a non-zero interaction is always allowed by the observational data and the dark energy equation of state is bent towards the phantom regime. In particular, when $H_0$ from HST is added to Planck 2018+BAO, we find an evidence for a non-zero coupling at more than $2\sigma$ confidence level. Our analyses also report that for both the models, $H_0$ is close to its local measurements and thus alleviating the $H_0$ tension. In particular, one of the interacting models perfectly solves the $H_0$ tension.

We investigate the recently introduced metastable dark energy (DE) models after the final Planck 2018 legacy release. The essence of the present work is to analyze their evolution at the level of perturbations. Our analyses show that both the metastable dark energy models considered in this article, are excellent candidates to alleviate the $H_0$ tension. In particular, for the present models, Planck 2018 alone can alleviate the $H_0$ tension within 68\% CL. Along with the final cosmic microwave background data from the Planck 2018 legacy release, we also include external cosmological datasets in order to asses the robustness of our findings.

In the present article we study the cosmological evolution of a two-scalar field gravitational theory defined in the Jordan frame. Specifically, we assume one of the scalar fields to be minimally coupled to gravity, while the second field which is the Brans-Dicke scalar field is nonminimally coupled to gravity and also coupled to the other scalar field. In the Einstein frame this theory reduces to a two-scalar field theory where the two fields can interact only in the potential term, which means that the quintom theory is recovered. The cosmological evolution is studied by analyzing the equilibrium points of the field equations in the Jordan frame. We find that the theory can describe the cosmological evolution in large scales, while inflationary solutions are also provided.

In a spatially flat Friedmann-Lema\^{i}tre-Robertson-Walker universe, we investigate a unified cosmic fluid scenario endowed with bulk viscosity in which the coefficient of the bulk viscosity has a power law evolution. The power law in the bulk viscous coefficient is a general case in this study which naturally includes several choices as special cases. Considering such a general bulk viscous scenario, in the present work we have extracted the observational constraints using the latest cosmological datasets and examine their behaviour at the level of both background and perturbations. From the observational analyses, we find that a non-zero bulk viscous coefficient is always favored and some of the models in this series are able to weaken the current tension on $H_0$ for some dataset. However, from the Bayesian evidence analysis, $\Lambda$CDM is favored over the bulk viscous model.

We consider a class of interacting dark energy models in a flat and nonflat FLRW universe where the interaction is characterized by the modified evolution of the pressureless dark matter as $a^{-3+\delta (a)}$, $a$ being the FLRW scale factor and $\delta (a)$ quantifies the interaction rate. By assuming the most natural and nonsingular parametrization for $\delta (a)$ as $\delta \left( a\right) =\sum_{i} \delta_i (1-a)^i $, where $\delta _{i}$'s ($i=0,1,2,3,..$) are constants, we reconstruct the expansion history of the universe for three particular choices of the DE sector using different cosmological datasets. Our analyses show that the non-interacting scenario is consistent with the observations while the interaction is not strictly ruled out. We reconstruct in the following way. We start with the first two terms of $\delta (a)$ above and constrain $\delta_0$, $\delta_1$. Then we consider up to the second order terms in $\delta (a)$ but fix $\delta_0$, $\delta_1$ to their constrained values and constrain $\delta_2$; similarly we constrain $\delta_3$, and finally we constrain $(\delta_0, \delta_1, \delta_1, \delta_3)$ by keeping all of them to be free as a generalized case. Our reconstruction technique shows that the constraints on $\delta_2$ (fixing $\delta_0$ and $\delta_1$) and $\delta_3$ (fixing $\delta_0$, $\delta_1$ and $\delta_2$) are almost zero for any interaction model and thus the effective scenario is well described by the linear parametrization $\delta (a)\simeq \delta_{0}+\delta _{1}(1-a)$. Additionally, a strong negative correlation between $\delta_0$, $\delta_1$ is observed independently of the dark energy fluid and the curvature of our universe.

Unified cosmological models have received a lot of attention in astrophysics community for explaining both the dark matter and dark energy evolution. The Chaplygin cosmologies, a well known name in this group have been investigated matched with observations from different sources. Obviously, Chaplygin cosmologies have to obey restrictions in order to be consistent with the observational data. As a consequence, alternative unified models, differing from Chaplygin model, are of special interest. In the present work we consider a specific example of such a unified cosmological model, that is quantified by only a single parameter $\mu$, that can be considered as a minimal extension of the $\Lambda$-cold dark matter cosmology. We investigate its observational boundaries together with an analysis of the universe at large scale. Our study shows that at early time the model behaves like a dust, and as time evolves, it mimics a dark energy fluid depicting a clear transition from the early decelerating phase to the late cosmic accelerating phase. Finally, the model approaches the cosmological constant boundary in an asymptotic manner. We remark that for the present unified model, the estimations of $H_0$ are slightly higher than its local estimation and thus alleviating the $H_0$ tension.

Lie symmetries are applied to classify the source of the magnetic field for the Pulsar equation near to the surface of the neutron star. We find that there are six possible different admitted Lie algebras. We apply the corresponding Lie invariants to reduce the Pulsar equation close to the surface to an ordinary differential equation. This equation is solved either with the use of Lie symmetries or the application of the ARS algorithm for singularity analysis to write the analytic solution as a Laurent expansion. These solutions are called inner solutions.

We examine the phase space of Ho\v{r}ava-Lifshitz cosmology for a wide range of self-interacting potentials for the scalar field under the detailed-balance condition and without imposing it, by means of the powerful method of $f$-devisers. A compactification approach is performed for the exponential potential and for potentials beyond the exponential one, extending the previous findings in the literature. By using this approach it is possible to describe the finite region of the phase space and the region where the phase-space variables becomes infinity. Furthermore, we present several results concerning the stability of the \emph{de Sitter} solution in Ho\v{r}ava-Lifshitz cosmology using Center Manifold theory. The advantages of this procedure are unveiled immediately when it is compared with the Normal Forms Calculations presented before in the literature.

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.

We investigate the cosmological dynamics of interacting dark energy models in which the interaction function is a nonlinear in terms of the energy densities. Considering explicitly the interaction between a pressureless dark matter and a scalar field, minimally coupled to Einstein gravity, we explore the dynamics of the spatially flat FLRW universe for the exponential potential of the scalar field.. We perform the stability analysis for the three nonlinear interaction models of our consideration through the analysis of critical points and we investigate the cosmological parameters and we discuss the physical behaviour at the critical points. From the analysis of the critical points we find a number of possibilities that include the stable late time accelerated solution, $w$CDM-like solution, radiation-like solution and moreover the unstable inflationary solution as well.

Cosmological evolution driven incorporating continuous particle creation by the time-varying gravitational field is investigated. We consider a spatially flat, homogeneous and isotropic universe with two matter fluids in the context of general relativity. One fluid is endowed with gravitationally induced `adiabatic' particle creation, while the second fluid simply satisfies the conservation of energy. We show that the dynamics of the two fluids is entirely controlled by a single nonlinear differential equation involving the particle creation rate, $\Gamma (t)$. We consider a very general particle creation rate, $\Gamma (t)$, that reduces to several special cases of cosmological interest, including $\Gamma =$ constant, $% \Gamma \propto 1/H^{n}$ ($n\in \mathbb{N}$), $\Gamma \propto \exp (1/H)$. Finally, we present singular algebraic solutions of the gravitational field equations for the two-fluid particle creation models and discuss their stability.

Cosmological models where dark matter (DM) and dark energy (DE) interact with each other are the general scenarios in compared to the non-interacting models. The interaction is usually motivated from the phenomenological ground and thus there is no such rule to prefer a particular interaction between DM and DE. Being motivated, in this work, allowing an exponential interaction between DM and DE in a spatially flat homogeneous and isotropic universe, we explore the dynamics of the universe through the constraints of the free parameters where the strength of the interaction is characterized by the dimensionless coupling parameter $\xi$ and the equation of state (EoS) for DE, $w_x$, is supposed to be a constant. The interaction scenario is fitted using the latest available observational data. Our analyses report that the observational data permit a non-zero value of $\xi$ but it is very small and consistent with $\xi =0$. From the constraints on $w_x$, we find that both phantom ($w_x< -1$) and quintessence ($w_x> -1$) regimes are equally allowed but $w_x$ is very close to `$-1$'. The overall results indicate that at the background level, the interaction model cannot be distinguished from the base $\Lambda $-cold dark matter model while from the perturbative analyses, the interaction model mildly deviates from the base model. We highlight that, even if we allow DM and DE to interact in an exponential manner, but according to the observational data, the evidence for a non-zero coupling is very small.

We consider the existence of an "inflaton" described by an homogeneous scalar field in the Szekeres cosmological metric. The gravitational field equations are reduced to two families of solutions which describe the homogeneous Kantowski-Sachs spacetime and an inhomogeneous FLRW(-like) spacetime with spatial curvature a constant. The main differences with the original Szekeres spacetimes containing only pressure-free matter are discussed. We investigate the stability of the two families of solution by studying the critical points of the field equations. We find that there exist stable solutions which describe accelerating spatially-flat FLRW geometries.

The cosmological viability of varying $G\left( t\right) $ and $\Lambda \left( t\right) $ cosmology is discussed by determining the cosmological eras provided by the theory. Such a study is performed with the determination of the critical points while stability analysis is performed. The application of Renormalization group in the ADM formalism of General Relativity provides a modified second-order theory of gravity where varying $G\left( t\right) $ plays the role of a minimally coupled field, different from that of Scalar-tensor theories, while $\Lambda\left( t\right) =\Lambda\left( G\left( t\right) \right) $ is a potential term. We find that the theory provides two de Sitter phases and a tracking solution. In the presence of matter source, two new critical points are introduced, where the matter source contributes to the universe. One of those points describes the $\Lambda$CDM cosmology and in order for the solution at the point to be cosmologically viable, it has to be unstable. Moreover, the second point, where matter exists, describes a universe where the dark energy parameter for the equation of state has a different value from that of the cosmological constant.

We use a mathematical approach based on the constraints systems in order to reconstruct the equation of state and the inflationary potential for the inflaton field from the observed spectral indices for the density perturbations $n_{s}$ and the tensor to scalar ratio $r$. From the astronomical data, we can observe that the measured values of these two indices lie on a two-dimensional surface. We express these indices in terms of the Hubble slow-roll parameters and we assume that $n_{s}-1=h\left( r\right) $. For the function $h\left( r\right) $, we consider three cases, where $h\left( r\right) $ is constant, linear and quadratic, respectively. From this, we derive second-order equations whose solutions provide us with the explicit forms for the expansion scale-factor, the scalar-field potential, and the effective equation of state for the scalar field. Finally, we show that for there exist mappings which transform one cosmological solution to another and allow new solutions to be generated from existing ones.

In a higher-order modified teleparallel theory cosmological we present analytical cosmological solutions. In particular we determine forms of the unknown potential which drives the scalar field such that the field equations form a Liouville integrable system. For the determination of the conservation laws we apply the Cartan symmetries. Furthermore, inspired from our solutions, a toy model is studied and it is shown that it can describe the Supernova data, while at the same time introduces dark matter components in the Hubble function. When the extra matter source is a stiff fluid then we show how analytical solutions for Bianchi I universes can be constructed from our analysis. Finally, we perform a global dynamical analysis of the field equations by using variables different from that of the Hubble-normalization.

We consider nonlinear redshift-dependent equation of state parameters as dark energy models in a spatially flat Friedmann-Lema\^itre-Robertson-Walker universe. To depict the expansion history of the universe in such cosmological scenarios, we take into account the large scale behaviour of such parametric models and fit them using a set of latest observational data with distinct origin that includes cosmic microwave background radiation, Supernove Type Ia, baryon acoustic oscillations, redshift space distortion, weak gravitational lensing, Hubble parameter measurements from cosmic chronometers and finally the local Hubble constant from Hubble space telescope. The fitting technique avails the publicly available code Cosmological Monte Carlo (CosmoMC), to extract the cosmological information out of these parametric dark energy models. From our analysis it follows that those models could describe the late time accelerating phase of the universe, while they are distinguished from the $\Lambda-$cosmology.

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 show that in scalar-field cosmology, a dust fluid follows as quantum corrections from solutions of the Wheeler-DeWitt equation generated by Lie symmetries. The energy density of the dust fluid is related with the frequency of the wavefunction.

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.

The existence and the stability conditions for some exact relativistic solutions of special interest are studied in a higher-order modified teleparallel gravitational theory. The theory with the use of a Lagrange multiplier is equivalent with that of General Relativity with a minimally coupled noncanonical field. The conditions for the existence of de Sitter solutions and ideal gas solutions in the case of vacuum are studied as also the stability criteria. Furthermore, in the presence of matter the behaviour of scaling solutions is given. Finally, we discuss the degrees of freedom of the field equations and we reduce the field equations in an algebraic equation, where in order to demonstrate our result we show how this noncanonical scalar field can reproduce the Hubble function of $\Lambda$-cosmology.

In a homogeneous and isotropic universe with non-zero spatial curvature we consider the effects of gravitational particle production in the dynamics of the universe. We show that the dynamics of the universe in such a background is characterized by a single nonlinear differential equation which is significantly dependent on the rate of particle creation and whose solutions can be dominated by the curvature effects at early times. For different particle creation rates we apply the singularity test in order to find the analytic solutions of the background dynamics. We describe the behavior of the cosmological solutions for both open and closed universes. We also show how the effects of curvature can be produced by the presence of a second perfect fluid with an appropriate equation of state. By combining that results with the analysis of the critical points we find that our consideration can be related with the pre-inflationary era. Specifically we find that for negative spatial curvature small changes of the Milne spacetime leads to a de Sitter universe.

A fourth-order theory of gravity is considered which in terms of dynamics has the same degrees of freedom and number of constraints as those of scalar-tensor theories. In addition it admits a canonical point-like Lagrangian description. We study the critical points of the theory and we show that it can describe the matter epoch of the universe and that two accelerated phases can be recovered one of which describes a de Sitter universe. Finally for some models exact solutions are presented.

A class of generalized Galileon cosmological models, which can be described by a point-like Lagrangian, is considered in order to utilize Noether's Theorem to determine conservation laws for the field equations. In the Friedmann-Lema\^itre-Robertson-Walker universe, the existence of a nontrivial conservation law indicates the integrability of the field equations. Due to the complexity of the latter, we apply the differential invariants approach in order to construct special power-law solutions and study their stability.

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.

We investigate Kantowski-Sachs models in Einstein-{\ae}ther theory with a perfect fluid source using the singularity analysis to prove the integrability of the field equations and dynamical system tools to study the evolution. We find an inflationary source at early times, and an inflationary sink at late times, for a wide region in the parameter space. The results by A. A. Coley, G. Leon, P. Sandin and J. Latta (JCAP 12, 010, 2015), are then re-obtained as particular cases. Additionally, we select other values for the non-GR parameters which are consistent with current constraints, getting a very rich phenomenology. In particular, we find solutions with infinite shear, zero curvature, and infinite matter energy density in comparison with the Hubble scalar. We also have stiff-like future attractors, anisotropic late-time attractors, or both, in some special cases. Such results are developed analytically, and then verified by numerics. Finally, the physical interpretation of the new critical points is discussed.

An algorithm is used to generate new solutions of the scalar field equations in homogeneous and isotropic universes. Solutions can be found for pure scalar fields with various potentials in the absence and presence of spatial curvature and other perfect fluids. A series of generalisations of the Chaplygin gas and bulk viscous cosmological solutions for inflationary universes are found. Furthermore other closed-form solutions which provide inflationary universes are presented. We also show how the Hubble slow-roll parameters can be calculated using the solution algorithm and we compare these inflationary solutions with the observational data provided by the Planck 2015 collaboration in order to constraint and rule out some of these models.

In the cosmological scenario in $f\left( T\right) $ gravity, we find analytical solutions for an isotropic and homogeneous universe containing a dust fluid and radiation and for an empty anisotropic Bianchi I universe. The method that we apply is that of movable singularities of differential equations. For the isotropic universe, the solutions are expressed in terms of a Laurent expansion, while for the anisotropic universe we find a family of exact Kasner-like solutions in vacuum. Finally, we discuss when a nonlinear $f\left( T\right) $-gravity theory provides solutions for the teleparallel equivalence of general relativity and derive conditions for exact solutions of general relativity to solve the field equations of an $f(T)$ theory.

In a flat Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) geometry, we consider the expansion of the universe powered by the gravitationally induced `adiabatic' matter creation. To demonstrate how matter creation works well with the expanding universe, we have considered a general creation rate and analyzed this rate in the framework of dynamical analysis. The dynamical analysis hints the presence of a non-singular universe (without the big bang singularity) with two successive accelerated phases, one at the very early phase of the universe (i.e. inflation), and the other one describes the current accelerating universe, where this early, late accelerated phases are associated with an unstable fixed point (i.e. repeller) and a stable fixed (attractor) points, respectively. We have described this phenomena by analytic solutions of the Hubble function and the scale factor of the FLRW universe. Using Jacobi Last multiplier method, we have found a Lagrangian for this matter creation rate describing this scenario of the universe. To match with our early physics results, we introduce an equivalent dynamics driven by a single scalar field and discussed the associated observable parameters compared them with the latest PLANCK data sets. Finally, introducing the teleparallel modified gravity, we have established an equivalent gravitational theory in the framework of matter creation.

We show that the recent results of \ [Int. J. Mod. Phys. D 25 (2016) 1650051] on the application of Lie/Noether symmetries in scalar field cosmology are well-known in the literature while the problem could have been solved easily under a coordinate transformation. That follows from the property, that the admitted group of invariant transformations of dynamical system is independent on the coordinate system.

We consider the application of group invariant transformations in order to constrain a flat isotropic and homogeneous cosmological model, containing of a Brans-Dicke scalar field and a perfect fluid with a constant equation of state parameter $w$, where the latter is not interacting with the scalar field in the gravitational action integral. The requirement that the Wheeler-DeWitt equation be invariant under one-parameter point transformations provides us with two families of power-law potentials for the Brans-Dicke field, in which the powers are functions of the Brans-Dicke parameter $\omega_{BD}$ and the parameter $w$. The existence of the Lie symmetry in the Wheeler-DeWitt equation is equivalent to the existence of a conserved quantity in field equations and with oscillatory terms in the wavefunction of the universe. This enables us to solve the field equations. For a specific value of the conserved quantity, we find a closed-form solution for the Hubble factor, which is equivalent to a cosmological model in general relativity containing two perfect fluids. This provides us with different models for specific values of the parameters $\omega_{BD},$ and $w$. Finally, the results hold for the specific case where the Brans-Dicke parameter $\omega_{BD}$ is zero, that is, for the O'Hanlon massive dilaton theory, and consequently for $f\left( R\right) $ gravity in the metric formalism.

We discuss the application of the Hojman's Symmetry Approach for the determination of conservation laws in Cosmology, which has been recently applied by various authors in different cosmological models. We show that Hojman's method for regular Hamiltonian systems, where the Hamiltonian function is one of the involved equations of the system, is equivalent to the application of Noether's Theorem for generalized transformations. That means that for minimally-coupled scalar field cosmology or other modified theories which are conformally related with scalar-field cosmology, like $f(R)$ gravity, the application of Hojman's method provide us with the same results with that of Noether's theorem. Moreover we study the special Ansatz. $\phi\left( t\right) =\phi\left( a\left( t\right) \right) $, which has been introduced for a minimally-coupled scalar field, and we study the Lie and Noether point symmetries for the reduced equation. We show that under this Ansatz, the unknown function of the model cannot be constrained by the requirement of the existence of a conservation law and that the Hojman conservation quantity which arises for the reduced equation is nothing more than the functional form of Noetherian conservation laws for the free particle. On the other hand, for $f(T)$ teleparallel gravity, it is not the existence of Hojman's conservation laws which provide us with the special function form of $f(T)$ functions, but the requirement that the reduced second-order differential equation admits a Jacobi Last multiplier, while the new conservation law is nothing else that the Hamiltonian function of the reduced equation.

Lie symmetries are discussed for the Wheeler-De Witt equation in Bianchi Class A cosmologies. In particular, we consider General Relativity, minimally coupled scalar field gravity and Hybrid Gravity as paradigmatic examples of the approach. Several invariant solutions are determined and classified according to the form of the scalar field potential. The approach gives rise to a suitable method to select classical solutions and it is based on the first principle of the existence of symmetries.

We apply as selection rule to determine the unknown functions of a cosmological model the existence of Lie point symmetries for the Wheeler-DeWitt equation of quantum gravity. Our cosmological setting consists of a flat Friedmann-Robertson-Walker metric having the scale factor $a(t)$, a scalar field with potential function $V(\phi)$ minimally coupled to gravity and a vector field of its kinetic energy is coupled with the scalar field by a coupling function $f(\phi)$. Then, the Lie symmetries of this dynamical system are investigated by utilizing the behavior of the corresponding minisuperspace under the infinitesimal generator of the desired symmetries. It is shown that by applying the Lie symmetry condition the form of the coupling function and also the scalar field potential function may be explicitly determined so that we are able to solve the Wheeler-DeWitt equation. Finally, we show how we can use the Lie symmetries in order to construct conservation laws and exact solutions for the field equations.

We consider quintessence scalar field cosmology in which the Lagrangian of the scalar field is modified by the Generalized Uncertainty Principle. We show that the perturbation terms which arise from the deformed algebra are equivalent with the existence of a second scalar field, where the two fields interact in the kinetic part. Moreover, we consider a spatially flat Friedmann-Lema\^{\i}tre-Robertson-Walker spacetime (FLRW), and we derive the gravitational field equations. We show that the modified equation of state parameter $w_{GUP}$ can cross the phantom divide line; that is $w_{GUP}<-1$. Furthermore, we derive the field equations in the dimensionless parameters, the dynamical system which arises is a singular perturbation system in which we study the existence of the fixed points in the slow manifold. Finally, we perform numerical simulations for some well known models and we show that for these models with the specific initial conditions, the parameter $w_{GUP}$ crosses the phantom barrier.

We give a general method to find exact cosmological solutions for scalar-field dark energy in the presence of perfect fluids. We use the existence of invariant transformations for the Wheeler De Witt (WdW) equation. We show that the existence of a point transformation under which the WdW equation is invariant is equivalent to the existence of conservation laws for the field equations, which indicates the existence of analytical solutions. We extend previous work by providing exact solutions for the Hubble parameter and the effective dark-energy equation of state parameter for cosmologies containing a combination of perfect fluid and a scalar field whose self-interaction potential is a power of hyperbolic functions. We find solutions explicitly when the perfect fluid is radiation or cold dark matter and determine the effects of non-zero spatial curvature. Using the Planck 2015 data, we determine the evolution of the effective equation of state of the dark energy. Finally, we study the global dynamics using dimensionless variables. We find that if the current cosmological model is Liouville integrable (admits conservation laws) then there is a unique stable point which describes the de-Sitter phase of the universe.

We comment upon the application of Hojman's method for the determination of conservation laws in Cosmology, which has been introduced by Capozziello \& Roshan (Phys. Lett. B 726 (2013) 471 (arXiv:1308.3910)), and has been applied recently in the cosmological scenario of a nonminimally coupled scalar field by Paolella \& Capozziello (Phys. Lett. A (2015), in press (arXiv:1503.00098)). We apply the Ansatz, $\phi\left( t\right) =\phi\left( a\left( t\right) \right) $, which was introduced by the cited authors for a minimally-coupled scalar field, and we study the Lie and Noether point symmetries for the reduced equation. We show that under this Ansatz the unknown function of the model cannot be constrained by the requirement of the existence of a conservation law and that the Hojman conservation quantity which arises for the reduced equation is nothing more than the functional form of the Noether conservation law of momentum for the free particle. Finally we show that Hojman's method for Hamiltonian systems, in which the Hamiltonian function is one of the involved equations of the system, is equivalent with the application of Noether's Theorem for generalized transformations.

In this thesis, we study the one parameter point transformations which leave invariant the differential equations. In particular we study the Lie and the Noether point symmetries of second order differential equations. We establish a new geometric method which relates the point symmetries of the differential equations with the collineations of the underlying manifold where the motion occurs. This geometric method is applied in order the two and three dimensional Newtonian dynamical systems to be classified in relation to the point symmetries; to generalize the Newtonian Kepler-Ermakov system in Riemannian spaces; to study the symmetries between classical and quantum systems and to investigate the geometric origin of the Type II hidden symmetries for the homogeneous heat equation and for the Laplace equation in Riemannian spaces. At last but not least, we apply this geometric approach in order to determine the dark energy models by use the Noether symmetries as a geometric criterion in modified theories of gravity.

We propose to use dynamical symmetries of the field equations, in order to classify the dark energy models in the context of scalar field (quintessence or phantom) FLRW cosmologies. Practically, symmetries provide a useful mathematical tool in physical problems since they can be used to simplify a given system of differential equations as well as to determine the integrability of the physical system. The requirement that the field equations admit dynamical symmetries results in two potentials one of which is the well known Unified Dark Matter (UDM) potential and another new potential. For each hyperbolic potential we obtain the corresponding analytic solution of the field equations. The proposed analysis suggests that the requirement of the contact symmetry appears to be very competitive to other independent tests used to probe the functional form of a given potential and thus the associated nature of dark energy. Finally, in order to test the viability of the above scalar field models we perform a joint likelihood analysis using some of the latest cosmological data.

We consider the two scalar field cosmology in a FRW spatially flat spacetime where the scalar fields interact both in the kinetic part and the potential. We apply the Noether point symmetries in order to define the interaction of the scalar fields. We use the point symmetries in order to write the field equations in the normal coordinates and we find that the Lagrangian of the field equations which admits at least three Noether point symmetries describes linear Newtonian systems. Furthermore, by using the corresponding conservation laws we find exact solutions of the field equations. Finally, we generalize our results to the case of N scalar fields interacting both in their potential and their kinematic part in a flat FRW background.

Symmetries play a crucial role in physics and, in particular, the Noether symmetries are a useful tool both to select models motivated at a fundamental level, and to find exact solutions for specific Lagrangians. In this work, we consider the application of point symmetries in the recently proposed metric-Palatini Hybrid Gravity in order to select the $f({\cal R})$ functional form and to find analytical solutions for the field equations and for the related Wheeler-DeWitt (WDW) equation. We show that, in order to find out integrable $f({\cal R})$ models, conformal transformations in the Lagrangians are extremely useful. In this context, we explore two conformal transformations of the forms $d\tau=N(a) dt$ and $d\tau=N(\phi) dt$. For the former conformal transformation, we found two cases of $f({\cal R})$ functions where the field equations admit Noether symmetries. In the second case, the Lagrangian reduces to a Brans-Dicke-like theory with a general coupling function. For each case, it is possible to transform the field equations by using normal coordinates to simplify the dynamical system and to obtain exact solutions. Furthermore, we perform quantization and derive the WDW equation for the minisuperspace model. The Lie point symmetries for the WDW equation are determined and used to find invariant solutions.

Spherically symmetric solutions for f(T) gravity models are derived by the so called Noether Symmetry Approach. First, we present a full set of Noether symmetries for some minisuperspace models. Then, we compute analytical solutions and find that spherically symmetric solutions in f(T) gravity can be recast in terms of Schwarzschild-like solutions modified by a distortion function depending on a characteristic radius. The obtained solutions are more general than those obtained by the usual solution methods.

In this paper, we present a complete Noether Symmetry analysis in the framework of scalar-tensor cosmology. Specifically, we consider a non-minimally coupled scalar field action embedded in the FLRW spacetime and provide a full set of Noether symmetries for related minisuperspaces. The presence of symmetries implies that the dynamical system becomes integrable and then we can compute cosmological analytical solutions for specific functional forms of coupling and potential functions selected by the Noether Approach.

We investigate the main features of the flat Friedmann-Lema{\i}tre-Robertson-Walker cosmological models in the f(T) teleparallel gravity. In particular, a general approach to find out exact cosmological solutions in f (T) gravity is discussed. Instead of taking into account phenomenological models, we consider as a selection criterion, the existence of Noether symmetries in the cosmological f(T) point-like Lagrangian. We find that only power-law models admit extra Noether symmetries. A complete analysis of such cosmological models is developed.

Conformally related metrics and Lagrangians are considered in the context of scalar-tensor gravity cosmology. After the discussion of the problem, we pose a lemma in which we show that the field equations of two conformally related Lagrangians are also conformally related if and only if the corresponding Hamiltonian vanishes. Then we prove that to every non-minimally coupled scalar field, we may associate a unique minimally coupled scalar field in a conformally related space with an appropriate potential. The latter result implies that the field equations of a non-minimally coupled scalar field are the same at the conformal level with the field equations of the minimally coupled scalar field. This fact is relevant in order to select physical variables among conformally equivalent systems. Finally, we find that the above propositions can be extended to a general Riemannian space of n-dimensions.

We perform a detailed study of the modified gravity $f(R)$ models in the light of the basic geometrical symmetries, namely Lie and Noether point symmetries, which serve to illustrate the phenomenological viability of the modified gravity paradigm as a serious alternative to the traditional scalar field approaches. In particular, we utilize a model-independent selection rule based on first integrals, due to Noether symmetries of the equations of motion, in order to identify the viability of $f(R)$ models in the context of flat FLRW cosmologies. The Lie/Noether point symmetries are computed for six modified gravity models that include also a cold dark matter component. As it is expected, we confirm that all the proposed modified gravity models admit the trivial first integral namely energy conservation. We find that only the $f(R)=(R^{b}-2\Lambda)^{c}$ model, which generalizes the concordance $\Lambda$ cosmology, accommodates extra Lie/Noether point symmetries. For this $f(R)$ model the existence of non-trivial Noether (first) integrals can be used to determine the integrability of the model. Indeed within this context we solve the problem analytically and thus we provide for the first time the evolution of the main cosmological functions such as the scale factor of the universe and the Hubble expansion rate.

We employ a three fluid model in order to construct a cosmological model in the Friedmann Robertson Walker flat spacetime, which contains three types of matter dark energy, dark matter and a perfect fluid with a linear equation of state. Dark matter is described by dust and dark energy with a scalar field with potential V({\phi}). In order to fix the scalar field potential we demand Lie symmetry invariance of the field equations, which is a model-independent assumption. The requirement of an extra Lie symmetry selects the exponential scalar field potential. The further requirement that the analytic solution is invariant under the point transformation generated by the Lie symmetry eliminates dark matter and leads to a quintessence and a phantom cosmological model containing a perfect fluid and a scalar field. Next we assume that the Lagrangian of the system admits an extra Noether symmetry. This new assumption selects the scalar field potential to be exponential and forces the perfect fluid to be stiff. Furthermore the existence of the Noether integral allows for the integration of the dynamical equations. We find new analytic solutions to quintessence and phantom cosmologies which contain all three fluids. Using these solutions one is able to compute analytically all main cosmological functions, such as the scale factor, the scalar field, the Hubble expansion rate, the deceleration parameter etc.

We propose to use a model-independent criterion based on first integrals of motion, due to Noether symmetries of the equations of motion, in order to classify the dark energy models in the context of scalar field (quintessence or phantom) FLRW cosmologies. In general, the Noether symmetries play an important role in physics because they can be used to simplify a given system of differential equations as well as to determine the integrability of the system. The Noether symmetries are computed for nine distinct accelerating cosmological scenarios that contain a homogeneous scalar field associated with different types of potentials. We verify that all the scalar field potentials, presented here, admit the trivial first integral namely energy conservation, as they should. We also find that the exponential potential inspired from scalar field cosmology, as well as some types of hyperbolic potentials, include extra Noether symmetries. This feature suggests that these potentials should be preferred along the hierarchy of scalar field potentials. Finally, using the latter potentials, in the framework of either quintessence or phantom scalar field cosmologies that contain also a non-relativistic matter(dark matter) component, we find that the main cosmological functions, such as the scale factor of the universe, the scalar field, the Hubble expansion rate and the metric of the FRLW space-time, are computed analytically. Interestingly, under specific circumstances the predictions of the exponential and hyperbolic scalar field models are equivalent to those of the $\Lambda$CDM model, as far as the global dynamics and the evolution of the scalar field are concerned. The present analysis suggests that our technique appears to be very competitive to other independent tests used to probe the functional form of a given potential and thus the associated nature of dark energy.